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.

A comparison of anatomic and cellular transcriptome structures across 40 human brain diseases.

Genes associated with risk for brain disease exhibit characteristic expression patterns that reflect both anatomical and cell type relationships. Brain-wide transcriptomic patterns of disease risk genes provide a molecular-based signature, based on differential co-expression, that is often unique to that disease. Brain diseases can be compared and aggregated based on the similarity of their signatures which often associates diseases from diverse phenotypic classes. Analysis of 40 common human brain diseases identifies 5 major transcriptional patterns, representing tumor-related, neurodegenerative, psychiatric and substance abuse, and 2 mixed groups of diseases affecting basal ganglia and hypothalamus. Further, for diseases with enriched expression in cortex, single-nucleus data in the middle temporal gyrus (MTG) exhibits a cell type expression gradient separating neurodegenerative, psychiatric, and substance abuse diseases, with unique excitatory cell type expression differentiating psychiatric diseases. Through mapping of homologous cell types between mouse and human, most disease risk genes are found to act in common cell types, while having species-specific expression in those types and preserving similar phenotypic classification within species. These results describe structural and cellular transcriptomic relationships of disease risk genes in the adult brain and provide a molecular-based strategy for classifying and comparing diseases, potentially identifying novel disease relationships.

Atypical co-development of the thalamus and cortex in autism: Evidence from age-related white-gray contrast change.

Recent studies have shown that white-gray contrast (WGC) of either cortical or subcortical gray matter provides for accurate predictions of age in typically developing (TD) children, and that, at least for the cortex, it changes differently with age in subjects with autism spectrum disorder (ASD) compared to their TD peers. Our previous study showed different patterns of contrast change between ASD and TD in sensorimotor and association cortices. While that study was confined to the cortex, we hypothesized that subcortical structures, particularly the thalamus, were involved in the observed cortical dichotomy between lower and higher processing. The current paper investigates that hypothesis using the WGC measures from the thalamus in addition to those from the cortex. We compared age-related WGC changes in the thalamus to those in the cortex. To capture the simultaneity of this change across the two structures, we devised a metric capturing the co-development of the thalamus and cortex (CoDevTC), proportional to the magnitude of cortical and thalamic age-related WGC change. We calculated this metric for each of the subjects in a large homogeneous sample taken from the Autism Brain Imaging Data Exchange (ABIDE) (N = 434). We used structural MRI data from the largest high-quality cross-sectional sample (NYU) as well as two other large high-quality sites, GU and OHSU, all three using Siemens 3T scanners. We observed that the co-development features in ASD and TD exhibit contrasting patterns; specifically, some higher-order thalamic nuclei, such as the lateral dorsal nucleus, exhibited reduction in codevelopment with most of the cortex in ASD compared to TD. Moreover, this difference in the CoDevTC pattern correlates with a number of behavioral measures across multiple cognitive and physiological domains. The results support previous notions of altered connectivity in autism, but add more specific evidence about the heterogeneity in thalamocortical development that elucidates the mechanisms underlying the clinical features of ASD.

A critical role of brain network architecture in a continuum model of autism spectrum disorders spanning from healthy individuals with genetic liability to individuals with ASD.

Studies have shown cortical alterations in individuals with autism spectrum disorders (ASD) as well as in individuals with high polygenic risk for ASD. An important addition to the study of altered cortical anatomy is the investigation of the underlying brain network architecture that may reveal brain-wide mechanisms in ASD and in polygenic risk for ASD. Such an approach has been proven useful in other psychiatric disorders by revealing that brain network architecture shapes (to an extent) the disorder-related cortical alterations. This study uses data from a clinical dataset-560 male subjects (266 individuals with ASD and 294 healthy individuals, CTL, mean age at 17.2 years) from the Autism Brain Imaging Data Exchange database, and data of 391 healthy individuals (207 males, mean age at 12.1 years) from the Pediatric Imaging, Neurocognition and Genetics database. ASD-related cortical alterations (group difference, ASD-CTL, in cortical thickness) and cortical correlates of polygenic risk for ASD were assessed, and then statistically compared with structural connectome-based network measures (such as hubs) using spin permutation tests. Next, we investigated whether polygenic risk for ASD could be predicted by network architecture by building machine-learning based prediction models, and whether the top predictors of the model were identified as disease epicenters of ASD. We observed that ASD-related cortical alterations as well as cortical correlates of polygenic risk for ASD implicated cortical hubs more strongly than non-hub regions. We also observed that age progression of ASD-related cortical alterations and cortical correlates of polygenic risk for ASD implicated cortical hubs more strongly than non-hub regions. Further investigation revealed that structural connectomes predicted polygenic risk for ASD (r = 0.30, p < 0.0001), and two brain regions (the left inferior parietal and left suparmarginal) with top predictive connections were identified as disease epicenters of ASD. Our study highlights a critical role of network architecture in a continuum model of ASD spanning from healthy individuals with genetic risk to individuals with ASD. Our study also highlights the strength of investigating polygenic risk scores in addition to multi-modal neuroimaging measures to better understand the interplay between genetic risk and brain alterations associated with ASD.

The Canadian Open Neuroscience Platform-An open science framework for the neuroscience community.

The Canadian Open Neuroscience Platform (CONP) takes a multifaceted approach to enabling open neuroscience, aiming to make research, data, and tools accessible to everyone, with the ultimate objective of accelerating discovery. Its core infrastructure is the CONP Portal, a repository with a decentralized design, where datasets and analysis tools across disparate platforms can be browsed, searched, accessed, and shared in accordance with FAIR principles. Another key piece of CONP infrastructure is NeuroLibre, a preprint server capable of creating and hosting executable and fully reproducible scientific publications that embed text, figures, and code. As part of its holistic approach, the CONP has also constructed frameworks and guidance for ethics and data governance, provided support and developed resources to help train the next generation of neuroscientists, and has fostered and grown an engaged community through outreach and communications. In this manuscript, we provide a high-level overview of this multipronged platform and its vision of lowering the barriers to the practice of open neuroscience and yielding the associated benefits for both individual researchers and the wider community.

A 3D atlas of functional human brain energetic connectome based on neuropil distribution.

The human brain is energetically expensive, yet the key factors governing its heterogeneous energy distributions across cortical regions to support its diversity of functions remain unexplored. Here, we built up a 3D digital cortical energy atlas based on the energetic costs of all neuropil activities into a high-resolution stereological map of the human cortex with cellular and synaptic densities derived, respectively, from ex vivo histological staining and in vivo PET imaging. The atlas was validated with PET-measured glucose oxidation at the voxel level. A 3D cortical activity map was calculated to predict the heterogeneous activity rates across all cortical regions, which revealed that resting brain is indeed active with heterogeneous neuronal activity rates averaging around 1.2 Hz, comprising around 70% of the glucose oxidation of the cortex. Additionally, synaptic density dominates spatial patterns of energetics, suggesting that the cortical energetics rely heavily on the distribution of synaptic connections. Recent evidence from functional imaging studies suggests that some cortical areas act as hubs (i.e., interconnecting distinct and functionally active regions). An inverse allometric relationship was observed between hub metabolic rates versus hub volumes. Hubs with smaller volumes have higher synapse density, metabolic rate, and activity rates compared to nonhubs. The open-source BrainEnergyAtlas provides a granular framework for exploring revealing design principles in energy-constrained human cortical circuits across multiple spatial scales.

Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain.

BACKGROUND: Higher-order cognition is hypothesized to be implemented via distributed cortical networks that are linked via long-range connections. However, it is unknown how computational advantages of long-range connections reflect cortical microstructure and microcircuitry. METHODS: We investigated this question by (i) profiling long-range cortical connectivity using resting-state functional magnetic resonance imaging (MRI) and cortico-cortical geodesic distance mapping, (ii) assessing how long-range connections reflect local brain microarchitecture, and (iii) examining the microarchitectural similarity of regions connected through long-range connections. RESULTS: Analysis of 2 independent datasets indicated that sensory/motor areas had more clustered short-range connections, while transmodal association systems hosted distributed, long-range connections. Meta-analytical decoding suggested that this topographical difference mirrored shifts in cognitive function, from perception/action towards emotional/social processing. Analysis of myelin-sensitive in vivo MRI as well as postmortem histology and transcriptomics datasets established that gradients in functional connectivity distance are paralleled by those present in cortical microarchitecture. Notably, long-range connections were found to link spatially remote regions of association cortex with an unexpectedly similar microarchitecture. CONCLUSIONS: By mapping covarying topographies of long-range functional connections and cortical microcircuits, the current work provides insights into structure-function relations in human neocortex.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations.

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Accurate and Efficient Simulation of Very High-Dimensional Neural Mass Models with Distributed-Delay Connectome Tensors.

This paper introduces methods and a novel toolbox that efficiently integrates high-dimensional Neural Mass Models (NMMs) specified by two essential components. The first is the set of nonlinear Random Differential Equations (RDEs) of the dynamics of each neural mass. The second is the highly sparse three-dimensional Connectome Tensor (CT) that encodes the strength of the connections and the delays of information transfer along the axons of each connection. To date, simplistic assumptions prevail about delays in the CT, often assumed to be Dirac-delta functions. In reality, delays are distributed due to heterogeneous conduction velocities of the axons connecting neural masses. These distributed-delay CTs are challenging to model. Our approach implements these models by leveraging several innovations. Semi-analytical integration of RDEs is done with the Local Linearization (LL) scheme for each neural mass, ensuring dynamical fidelity to the original continuous-time nonlinear dynamic. This semi-analytic LL integration is highly computationally-efficient. In addition, a tensor representation of the CT facilitates parallel computation. It also seamlessly allows modeling distributed delays CT with any level of complexity or realism. This ease of implementation includes models with distributed-delay CTs. Consequently, our algorithm scales linearly with the number of neural masses and the number of equations they are represented with, contrasting with more traditional methods that scale quadratically at best. To illustrate the toolbox's usefulness, we simulate a single Zetterberg-Jansen and Rit (ZJR) cortical column, a single thalmo-cortical unit, and a toy example comprising 1000 interconnected ZJR columns. These simulations demonstrate the consequences of modifying the CT, especially by introducing distributed delays. The examples illustrate the complexity of explaining EEG oscillations, e.g., split alpha peaks, since they only appear for distinct neural masses. We provide an open-source Script for the toolbox.

BigBrain 3D atlas of cortical layers: Cortical and laminar thickness gradients diverge in sensory and motor cortices.

Histological atlases of the cerebral cortex, such as those made famous by Brodmann and von Economo, are invaluable for understanding human brain microstructure and its relationship with functional organization in the brain. However, these existing atlases are limited to small numbers of manually annotated samples from a single cerebral hemisphere, measured from 2D histological sections. We present the first whole-brain quantitative 3D laminar atlas of the human cerebral cortex. It was derived from a 3D histological atlas of the human brain at 20-micrometer isotropic resolution (BigBrain), using a convolutional neural network to segment, automatically, the cortical layers in both hemispheres. Our approach overcomes many of the historical challenges with measurement of histological thickness in 2D, and the resultant laminar atlas provides an unprecedented level of precision and detail. We utilized this BigBrain cortical atlas to test whether previously reported thickness gradients, as measured by MRI in sensory and motor processing cortices, were present in a histological atlas of cortical thickness and which cortical layers were contributing to these gradients. Cortical thickness increased across sensory processing hierarchies, primarily driven by layers III, V, and VI. In contrast, motor-frontal cortices showed the opposite pattern, with decreases in total and pyramidal layer thickness from motor to frontal association cortices. These findings illustrate how this laminar atlas will provide a link between single-neuron morphology, mesoscale cortical layering, macroscopic cortical thickness, and, ultimately, functional neuroanatomy.

Therapeutic effect of tempo in Mozart's "Sonata for two pianos" (K. 448) in patients with epilepsy: An electroencephalographic study.

BACKGROUND: Mozart's "Sonata for two pianos" (Köchel listing 448) has proven effective as music therapy for patients with epilepsy, but little is understood about the mechanism of which feature in it impacted therapeutic effect. This study explored whether tempo in that piece is important for its therapeutic effect. METHODS: We measured the effects of tempo in Mozart's sonata on clinical and electroencephalographic parameters of 147 patients with epilepsy who listened to the music at slow, original, or accelerated speed. As a control, patients listened to Haydn's Symphony no. 94 at original speed. RESULTS: Listening to Mozart's piece at original speed significantly reduced the number of interictal epileptic discharges. It decreased beta power in the frontal, parietal, and occipital regions, suggesting increased auditory attention and reduced visual attention. It also decreased functional connectivity among frontal, parietal, temporal, and occipital brain regions, also suggesting increased auditory attention and reduced visual attention. No such effects were observed after patients listened to the slow or fast version of Mozart's piece, or to Haydn's symphony at normal speed. CONCLUSIONS: These results suggest that Mozart's "Sonata for two pianos" may exert therapeutic effects by regulating attention when played at its original tempo, but not slower or faster. These findings may help guide the design and optimization of music therapy against epilepsy.

Early brain development in infants at high risk for autism spectrum disorder.

Brain enlargement has been observed in children with autism spectrum disorder (ASD), but the timing of this phenomenon, and the relationship between ASD and the appearance of behavioural symptoms, are unknown. Retrospective head circumference and longitudinal brain volume studies of two-year olds followed up at four years of age have provided evidence that increased brain volume may emerge early in development. Studies of infants at high familial risk of autism can provide insight into the early development of autism and have shown that characteristic social deficits in ASD emerge during the latter part of the first and in the second year of life. These observations suggest that prospective brain-imaging studies of infants at high familial risk of ASD might identify early postnatal changes in brain volume that occur before an ASD diagnosis. In this prospective neuroimaging study of 106 infants at high familial risk of ASD and 42 low-risk infants, we show that hyperexpansion of the cortical surface area between 6 and 12 months of age precedes brain volume overgrowth observed between 12 and 24 months in 15 high-risk infants who were diagnosed with autism at 24 months. Brain volume overgrowth was linked to the emergence and severity of autistic social deficits. A deep-learning algorithm that primarily uses surface area information from magnetic resonance imaging of the brain of 6-12-month-old individuals predicted the diagnosis of autism in individual high-risk children at 24 months (with a positive predictive value of 81% and a sensitivity of 88%). These findings demonstrate that early brain changes occur during the period in which autistic behaviours are first emerging.

Association of cumulative prenatal adversity with infant subcortical structure volumes and child problem behavior and its moderation by a coexpression polygenic risk score of the serotonin system.

Prenatal adversity has been linked to later psychopathology. Yet, research on cumulative prenatal adversity, as well as its interaction with offspring genotype, on brain and behavioral development is scarce. With this study, we aimed to address this gap. In Finnish mother-infant dyads, we investigated the association of a cumulative prenatal adversity sum score (PRE-AS) with (a) child emotional and behavioral problems assessed with the Strengths and Difficulties Questionnaire at 4 and 5 years (N = 1568, 45.3% female), (b) infant amygdalar and hippocampal volumes (subsample N = 122), and (c) its moderation by a hippocampal-specific coexpression polygenic risk score based on the serotonin transporter (SLC6A4) gene. We found that higher PRE-AS was linked to greater child emotional and behavioral problems at both time points, with partly stronger associations in boys than in girls. Higher PRE-AS was associated with larger bilateral infant amygdalar volumes in girls compared to boys, while no associations were found for hippocampal volumes. Further, hyperactivity/inattention in 4-year-old girls was related to both genotype and PRE-AS, the latter partially mediated by right amygdalar volumes as preliminary evidence suggests. Our study is the first to demonstrate a dose-dependent sexually dimorphic relationship between cumulative prenatal adversity and infant amygdalar volumes.

EEG effective connectivity during the first year of life mirrors brain synaptogenesis, myelination, and early right hemisphere predominance.

INTRODUCTION: The maturation of electroencephalogram (EEG) effective connectivity in healthy infants during the first year of life is described. METHODS: Participants: A cross-sectional sample of 125 healthy at-term infants, from 0 to 12 months of age, underwent EEG in a state of quiet sleep. PROCEDURES: The EEG primary currents at the source were described with the sLoreta method. An unmixing algorithm was applied to reduce the leakage, and the isolated effective coherence, a direct and directed measurement of information flow, was calculated. RESULTS AND DISCUSSION: Initially, the highest indices of connectivity are at the subcortical nuclei, continuing to the parietal lobe, predominantly the right hemisphere, then expanding to temporal, occipital, and finally the frontal areas, which is consistent with the myelination process. Age-related connectivity changes were mostly long-range and bilateral. Connections increased with age, mainly in the right hemisphere, while they mainly decreased in the left hemisphere. Increased connectivity from 20 to 30 Hz, mostly at the right hemisphere. These findings were consistent with right hemisphere predominance during the first three years of life. Theta and alpha connections showed the greatest changes with age. Strong connectivity was found between the parietal, temporal, and occipital regions to the frontal lobes, responsible for executive functions and consistent with behavioral development during the first year. The thalamus exchanges information bidirectionally with all cortical regions and frequency bands. CONCLUSIONS: The maturation of EEG connectivity during the first year in healthy infants is very consistent with synaptogenesis, reductions in synaptogenesis, myelination, and functional and behavioral development.

Early protein energy malnutrition impacts life-long developmental trajectories of the sources of EEG rhythmic activity.

Protein Energy Malnutrition (PEM) has lifelong consequences on brain development and cognitive function. We studied the lifelong developmental trajectories of resting-state EEG source activity in 66 individuals with histories of Protein Energy Malnutrition (PEM) limited to the first year of life and in 83 matched classmate controls (CON) who are all participants of the 49 years longitudinal Barbados Nutrition Study (BNS). qEEGt source z-spectra measured deviation from normative values of EEG rhythmic activity sources at 5-11 years of age and 40 years later at 45-51 years of age. The PEM group showed qEEGt abnormalities in childhood, including a developmental delay in alpha rhythm maturation and an insufficient decrease in beta activity. These profiles may be correlated with accelerated cognitive decline.

Harmonized-Multinational qEEG norms (HarMNqEEG).

This paper extends frequency domain quantitative electroencephalography (qEEG) methods pursuing higher sensitivity to detect Brain Developmental Disorders. Prior qEEG work lacked integration of cross-spectral information omitting important functional connectivity descriptors. Lack of geographical diversity precluded accounting for site-specific variance, increasing qEEG nuisance variance. We ameliorate these weaknesses. (i) Create lifespan Riemannian multinational qEEG norms for cross-spectral tensors. These norms result from the HarMNqEEG project fostered by the Global Brain Consortium. We calculate the norms with data from 9 countries, 12 devices, and 14 studies, including 1564 subjects. Instead of raw data, only anonymized metadata and EEG cross-spectral tensors were shared. After visual and automatic quality control, developmental equations for the mean and standard deviation of qEEG traditional and Riemannian DPs were calculated using additive mixed-effects models. We demonstrate qEEG "batch effects" and provide methods to calculate harmonized z-scores. (ii) We also show that harmonized Riemannian norms produce z-scores with increased diagnostic accuracy predicting brain dysfunction produced by malnutrition in the first year of life and detecting COVID induced brain dysfunction. (iii) We offer open code and data to calculate different individual z-scores from the HarMNqEEG dataset. These results contribute to developing bias-free, low-cost neuroimaging technologies applicable in various health settings.

A sub+cortical fMRI-based surface parcellation.

Both cortical and subcortical structures are organized into a large number of distinct areas reflecting functional and cytoarchitectonic differences. Mapping these areas is of fundamental importance to neuroscience. A central obstacle to this task is the inaccuracy associated with bringing results from individuals into a common space. The vast individual differences in morphology pose a serious problem for volumetric registration. Surface-based approaches fare substantially better, but have thus far been used only for cortical parcellation, leaving subcortical parcellation in volumetric space. We extend the surface-based approach to include also the subcortical deep gray-matter structures, thus achieving a uniform representation across both cortex and subcortex, suitable for use with surface-based metrics that span these structures, for example, white/gray contrast. Using data from the Enhanced Nathan Klein Institute-Rockland Sample, limited to individuals between 19 and 69 years of age, we generate a functional parcellation of both the cortical and subcortical surfaces. To assess this extended parcellation, we show that (a) our parcellation provides greater homogeneity of functional connectivity patterns than do arbitrary parcellations matching in the number and size of parcels; (b) our parcels align with known cortical and subcortical architecture; and (c) our extended functional parcellation provides an improved fit to the complexity of life-span (6-85 years) changes in white/gray contrast data compared to arbitrary parcellations matching in the number and size of parcels, supporting its use with surface-based measures. We provide our extended functional parcellation for the use of the neuroimaging community.

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.

Convolutional neural networks for cytoarchitectonic brain mapping at large scale.

Human brain atlases provide spatial reference systems for data characterizing brain organization at different levels, coming from different brains. Cytoarchitecture is a basic principle of the microstructural organization of the brain, as regional differences in the arrangement and composition of neuronal cells are indicators of changes in connectivity and function. Automated scanning procedures and observer-independent methods are prerequisites to reliably identify cytoarchitectonic areas, and to achieve reproducible models of brain segregation. Time becomes a key factor when moving from the analysis of single regions of interest towards high-throughput scanning of large series of whole-brain sections. Here we present a new workflow for mapping cytoarchitectonic areas in large series of cell-body stained histological sections of human postmortem brains. It is based on a Deep Convolutional Neural Network (CNN), which is trained on a pair of section images with annotations, with a large number of un-annotated sections in between. The model learns to create all missing annotations in between with high accuracy, and faster than our previous workflow based on observer-independent mapping. The new workflow does not require preceding 3D-reconstruction of sections, and is robust against histological artefacts. It processes large data sets with sizes in the order of multiple Terabytes efficiently. The workflow was integrated into a web interface, to allow access without expertise in deep learning and batch computing. Applying deep neural networks for cytoarchitectonic mapping opens new perspectives to enable high-resolution models of brain areas, introducing CNNs to identify borders of brain areas.

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.