Neuroimaging has been facing a data deluge characterized by the exponential growth of both raw and processed data. As a result, mining the massive quantities of digital data collected in these studies offers unprecedented opportunities and has become paramount for today's research. As the neuroimaging community enters the world of "Big Data", there has been a concerted push for enhanced sharing initiatives, whether within a multisite study, across studies, or federated and shared publicly. This article will focus on the database and processing ecosystem developed at the Montreal Neurological Institute (MNI) to support multicenter data acquisition both nationally and internationally, create database repositories, facilitate data-sharing initiatives, and leverage existing software toolkits for large-scale data processing.

The neural correlates of human personality have been of longstanding interest; however, most studies in the field have relied on modest sample sizes and few replicable results have been reported to date. We investigated relationships between personality and brain gray matter in a sample of generally healthy, older (mean age 73 years) adults from Scotland drawn from the Lothian Birth Cohort 1936. Participants (N = 578) completed a brain MRI scan and self-reported Big Five personality trait measures. Conscientiousness trait scores were positively related to brain cortical thickness in a range of regions, including bilateral parahippocampal gyrus, bilateral fusiform gyrus, left cingulate gyrus, right medial orbitofrontal cortex, and left dorsomedial prefrontal cortex. These associations - most notably in frontal regions - were modestly-to-moderately attenuated by the inclusion of biomarker variables assessing allostatic load and smoking status. None of the other personality traits showed robust associations with brain cortical thickness, nor did we observe any personality trait associations with cortical surface area and gray matter volume. These findings indicate that brain cortical thickness is associated with conscientiousness, perhaps partly accounted for by allostatic load and smoking status.

Generative models focused on multifactorial causal mechanisms in brain disorders are scarce and generally based on limited data. Despite the biological importance of the multiple interacting processes, their effects remain poorly characterized from an integrative analytic perspective. Here, we propose a spatiotemporal multifactorial causal model (MCM) of brain (dis)organization and therapeutic intervention that accounts for local causal interactions, effects propagation via physical brain networks, cognitive alterations, and identification of optimum therapeutic interventions. In this article, we focus on describing the model and applying it at the population-based level for studying late onset Alzheimer's disease (LOAD). By interrelating six different neuroimaging modalities and cognitive measurements, this model accurately predicts spatiotemporal alterations in brain amyloid-β (Aβ) burden, glucose metabolism, vascular flow, resting state functional activity, structural properties, and cognitive integrity. The results suggest that a vascular dysregulation may be the most-likely initial pathologic event leading to LOAD. Nevertheless, they also suggest that LOAD it is not caused by a unique dominant biological factor (e.g. vascular or Aβ) but by the complex interplay among multiple relevant direct interactions. Furthermore, using theoretical control analysis of the identified population-based multifactorial causal network, we show the crucial advantage of using combinatorial over single-target treatments, explain why one-target Aβ based therapies might fail to improve clinical outcomes, and propose an efficiency ranking of possible LOAD interventions. Although still requiring further validation at the individual level, this work presents the first analytic framework for dynamic multifactorial brain (dis)organization that may explain both the pathologic evolution of progressive neurological disorders and operationalize the influence of multiple interventional strategies.

Verbal and non-verbal intelligence in children is highly correlated, and thus, it has been difficult to differentiate their neural substrates. Nevertheless, recent studies have shown that verbal and non-verbal intelligence can be dissociated and focal cortical regions corresponding to each have been demonstrated. However, the pattern of structural covariance corresponding to verbal and non-verbal intelligence remains unexplored. In this study, we used 586 longitudinal anatomical MRI scans of subjects aged 6-18 years, who had concurrent intelligence quotient (IQ) testing on the Wechsler Abbreviated Scale of Intelligence. Structural covariance networks (SCNs) were constructed using interregional correlations in cortical thickness for low-IQ (Performance IQ=100±8, Verbal IQ=100±7) and high-IQ (PIQ=121±8, VIQ=120±9) groups. From low- to high-VIQ group, we observed constrained patterns of anatomical coupling among cortical regions, complemented by observations of higher global efficiency and modularity, and lower local efficiency in high-VIQ group, suggesting a shift towards a more optimal topological organization. Analysis of nodal topological properties (regional efficiency and participation coefficient) revealed greater involvement of left-hemispheric language related regions including inferior frontal and superior temporal gyri for high-VIQ group. From low- to high-PIQ group, we did not observe significant differences in anatomical coupling patterns, global and nodal topological properties. Our findings indicate that people with higher verbal intelligence have structural brain differences from people with lower verbal intelligence - not only in localized cortical regions, but also in the patterns of anatomical coupling among widely distributed cortical regions, possibly resulting to a system-level reorganization that might lead to a more efficient organization in high-VIQ group.

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.

Insular cortex is a core hub involved in multiple cognitive and socio-affective processes. Yet, the anatomical mechanisms that explain how it is involved in such a diverse array of functions remain incompletely understood. Here, we tested the hypothesis that changes in myeloarchitecture across the insular cortex explain how it can be involved in many different facets of cognitive function. Detailed intracortical profiling, performed across hundreds of insular locations on the basis of myelin-sensitive magnetic resonance imaging (MRI), was compressed into a lower-dimensional space uncovering principal axes of myeloarchitectonic variation. Leveraging two datasets with different high-resolution MRI contrasts, we obtained robust support for two principal dimensions of insular myeloarchitectonic differentiation in vivo, one running from ventral anterior to posterior banks and one radiating from dorsal anterior towards both ventral anterior and posterior subregions. Analyses of post mortem 3D histological data showed that the antero-posterior axis was mirrored in cytoarchitectural markers, even when controlling for sulco-gyral folding. Resting-state functional connectomics in the same individuals and ad hoc meta-analyses showed that myelin gradients in the insula relate to diverse affiliation to macroscale intrinsic functional systems, showing differential shifts in functional network embedding across each myelin-derived gradient. Collectively, our findings offer a novel approach to capture structure-function interactions of a key node of the limbic system, and suggest a multidimensional structural basis underlying the diverse functional roles of the insula.

The maturational schedule of typical brain development is tightly constrained; deviations from it are associated with cognitive atypicalities, and are potentially predictive of developmental disorders. Previously, we have shown that the white/gray contrast at the inner border of the cortex is a good predictor of chronological age, and is sensitive to aspects of brain development that reflect cognitive performance. Here we extend that work to include the white/gray contrast at the border of subcortical structures. We show that cortical and subcortical contrast together yield better age-predictions than any non-kernel-based method based on a single image-type, and that the residuals of the improved predictions provide new insight into unevenness in cognitive performance. We demonstrate the improvement in age predictions in two large datasets: the NIH Pediatric Data, with 831 scans of typically developing individuals between 4 and 22 years of age; and the Pediatric Imaging, Neurocognition, and Genetics data, with 909 scans of individuals in a similar age-range. Assessment of the relation of the residuals of these age predictions to verbal and performance IQ revealed correlations in opposing directions, and a principal component analysis of the residuals of the model that best fit the contrast data produced components related to either performance IQ or verbal IQ. Performance IQ was associated with the first principle component, reflecting increased cortical contrast, broadly, with almost no subcortical presence; verbal IQ was associated with the second principle component, reflecting reduced contrast in the basal ganglia and increased contrast in the bilateral arcuate fasciculi.

Several reports have described cortical thickness (CTh) developmental trajectories, with conflicting results. Some studies have reported inverted-U shape curves with peaks of CTh in late childhood to adolescence, while others suggested predominant monotonic decline after age 6. In this study, we reviewed CTh developmental trajectories in the NIH MRI Study of Normal Brain Development, and in a second step, evaluated the impact of post-processing quality control (QC) procedures on identified trajectories. The quality-controlled sample included 384 individual subjects with repeated scanning (1-3 per subject, total scans n=753) from 4.9 to 22.3years of age. The best-fit model (cubic, quadratic, or first-order linear) was identified at each vertex using mixed-effects models. The majority of brain regions showed linear monotonic decline of CTh. There were few areas of cubic trajectories, mostly in bilateral temporo-parietal areas and the right prefrontal cortex, in which CTh peaks were at, or prior to, age 8. When controlling for total brain volume, CTh trajectories were even more uniformly linear. The only sex difference was faster thinning of occipital areas in boys compared to girls. The best-fit model for whole brain mean thickness was a monotonic decline of 0.027mm per year. QC procedures had a significant impact on identified trajectories, with a clear shift toward more complex trajectories (i.e., quadratic or cubic) when including all scans without QC (n=954). Trajectories were almost exclusively linear when using only scans that passed the most stringent QC (n=598). The impact of QC probably relates to decreasing the inclusion of scans with CTh underestimation secondary to movement artifacts, which are more common in younger subjects. In summary, our results suggest that CTh follows a simple linear decline in most cortical areas by age 5, and all areas by age 8. This study further supports the crucial importance of implementing post-processing QC in CTh studies of development, aging, and neuropsychiatric disorders.

Idiopathic Parkinson's disease (PD) is a neurodegenerative disorder diagnosed on the basis of motor symptoms, but that also includes cognitive and visuo-spatial deficits. Though PD is known to initially affect subcortical regions, the cortex also exhibits neuronal loss in the course of the disease as post mortem studies have shown. So far, PD-related pattern of cortical damage remains unclear, because of disease-caused heterogeneity, and also in part because of methodological issues such as the limitations of Voxel Based Morphometry. Here corticometry was used, a technique that decouples local surface from thickness, to obtain a better picture of PD corticomorphometric patterns. We acquired MRI volumes for 33 healthy controls (HC) and 49 PD patients, extracted local cortical thickness and surface area and modeled both of them as a function of group and age for each participant. Cortical thickness averaged on the whole cortex did not differ between the two groups while mean surface area was significantly larger in the PD group. The bilateral parietal lobule, the right superior frontal gyrus, the left cingulate cortex and the left insular cortex exhibited larger local surface area in the PD group. The right precuneus exhibited cortical thinning associated with age in the PD group and not in the HC group. Furthermore, cortical thinning was observed in the PD group compared with the control group in the left medial supplementary motor area (SMA) and in the right dorsal pre-SMA. Finally, we found the left temporal pole thickness to correlate with disease duration, as well as the bilateral occipital cortex and Broca's area. These results suggest that PD etiology is associated with specific cortical alterations, which could account for cognitive deficits that arise as the disease evolves. Finally, our results observed in the occipital cortex as a function of disease duration may indicate the increase in PD-related visuo-spatial deficits, which can sometimes result in hallucinations later on in the disease. In the future, MRI-generated corticometry, combined with additional behavioral markers, may prove to be a useful diagnosis tool to characterize the evolution of motor and cognitive deficits in PD.

Knowing the maturational schedule of typical brain development is critical to our ability to identify deviations from it; such deviations have been related to cognitive performance and even developmental disorders. Chronological age can be predicted from brain images with considerable accuracy, but with limited spatial specificity, particularly in the case of the cerebral cortex. Methods using multi-modal data have shown the greatest accuracy, but have made limited use of cortical measures. Methods using complex measures derived from voxels throughout the brain have also shown great accuracy, but are difficult to interpret in terms of cortical development. Measures based on cortical surfaces have yielded less accurate predictions, suggesting that perhaps cortical maturation is less strongly related to chronological age than is maturation of deep white matter or subcortical structures. We question this suggestion. We show that a simple metric based on the white/gray contrast at the inner border of the cortex is a good predictor of chronological age. We demonstrate this in two large datasets: the NIH Pediatric Data, with 832 scans of typically developing children, adolescents, and young adults; and the Pediatric Imaging, Neurocognition, and Genetics data, with 760 scans of individuals in a similar age-range. Further, our usage of an elastic net penalized linear regression model reveals the brain regions which contribute most to age-prediction. Moreover, we show that the residuals of age-prediction based on this white/gray contrast metric are not merely random errors, but are strongly related to IQ, suggesting that this metric is sensitive to aspects of brain development that reflect cognitive performance.

The current evolution of 'cloud neuroscience' leads to more efforts with the large-scale EEG applications, by using EEG pipelines to handle the rapidly accumulating EEG data. However, there are a few specific cloud platforms that seek to address the cloud computational challenges of EEG big data analysis to benefit the EEG community. In response to the challenges, a WeBrain cloud platform (https://webrain.uestc.edu.cn/) is designed as a web-based brainformatics platform and computational ecosystem to enable large-scale EEG data storage, exploration and analysis using cloud high-performance computing (HPC) facilities. WeBrain connects researchers from different fields to EEG and multimodal tools that have become the norm in the field and the cloud processing power required to handle those large EEG datasets. This platform provides an easy-to-use system for novice users (even no computer programming skills) and provides satisfactory maintainability, sustainability and flexibility for IT administrators and tool developers. A range of resources are also available on https://webrain.uestc.edu.cn/, including documents, manuals, example datasets related to WeBrain, and collected links to open EEG datasets and tools. It is not necessary for users or administrators to install any software or system, and all that is needed is a modern web browser, which reduces the technical expertise required to use or manage WeBrain. The WeBrain platform is sponsored and driven by the China-Canada-Cuba international brain cooperation project (CCC-Axis, http://ccc-axis.org/), and we hope that WeBrain will be a promising cloud brainformatics platform for exploring brain information in large-scale EEG applications in the EEG community.

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.

Normal aging is accompanied by various cognitive functional declines. Recent studies have revealed disruptions in the coordination of large-scale functional brain networks such as the default mode network in advanced aging. However, organizational alterations of the structural brain network at the system level in aging are still poorly understood. Here, using cortical thickness, we investigated the modular organization of the cortical structural networks in 102 young and 97 normal aging adults. Brain networks for both cohorts displayed a modular organization overlapping with functional domains such as executive and auditory/language processing. However, compared with the modular organization of young adults, the aging group demonstrated a significantly reduced modularity that might be indicative of reduced functional segregation in the aging brain. More importantly, the aging brain network exhibited reduced intra-/inter-module connectivity in modules corresponding to the executive function and the default mode network of young adults, which might be associated with the decline of cognitive functions in aging. Finally, we observed age-associated alterations in the regional characterization in terms of their intra/inter-module connectivity. Our results indicate that aging is associated with an altered modular organization in the structural brain networks and provide new evidence for disrupted integrity in the large-scale brain networks that underlie cognition.

Prevailing psychometric theories of intelligence posit that individual differences in cognitive performance are attributable to three main sources of variance: the general factor of intelligence (g), cognitive ability domains, and specific test requirements and idiosyncrasies. Cortical thickness has been previously associated with g. In the present study, we systematically analyzed associations between cortical thickness and cognitive performance with and without adjusting for the effects of g in a representative sample of children and adolescents (N=207, Mean age=11.8; SD=3.5; Range=6 to 18.3 years). Seven cognitive tests were included in a measurement model that identified three first-order factors (representing cognitive ability domains) and one second-order factor representing g. Residuals of the cognitive ability domain scores were computed to represent g-independent variance for the three domains and seven tests. Cognitive domain and individual test scores as well as residualized scores were regressed against cortical thickness, adjusting for age, gender and a proxy measure of brain volume. g and cognitive domain scores were positively correlated with cortical thickness in very similar areas across the brain. Adjusting for the effects of g eliminated associations of domain and test scores with cortical thickness. Within a psychometric framework, cortical thickness correlates of cognitive performance on complex tasks are well captured by g in this demographically representative sample.

Previous studies of cortical asymmetry have relied mainly on voxel-based morphometry (VBM), or manual segmentation of regions of interest. This study uses fully automated, surface-based techniques to analyse position and surface area asymmetry for the mid-surfaces of 112 right-handed subjects' cortical hemispheres from a cohort of young adults. Native space measurements of local surface area asymmetry and vertex position asymmetry were calculated from surfaces registered to a previously validated hemisphere-unbiased surface-based template. Our analysis confirms previously identified hemispheric asymmetries (Yakovlevian torque, frontal and occipital petalia) in enhanced detail. It does not support previous findings of gender/asymmetry interactions or rightward planum parietale areal increase. It reveals several new findings, including a striking leftward increase in surface area of the supramarginal gyrus (peak effect 18%), compared with a smaller areal increase in the left Heschl's gyrus and planum temporale region (peak effect 8%). A second finding was rightward increase in surface area (peak effect 10%) in a band around the medial junction between the occipital lobe, and parietal and temporal lobes. By clearly separating out the effects of structural translocation and surface area change from those of thickness and curvature, this study resolves the confound of these variables inherent in VBM studies.

A variety of methods have been developed to identify brain networks with spontaneous, coherent activity in resting-state functional magnetic resonance imaging (fMRI). We propose here a generic statistical framework to quantify the stability of such resting-state networks (RSNs), which was implemented with k-means clustering. The core of the method consists in bootstrapping the available datasets to replicate the clustering process a large number of times and quantify the stable features across all replications. This bootstrap analysis of stable clusters (BASC) has several benefits: (1) it can be implemented in a multi-level fashion to investigate stable RSNs at the level of individual subjects and at the level of a group; (2) it provides a principled measure of RSN stability; and (3) the maximization of the stability measure can be used as a natural criterion to select the number of RSNs. A simulation study validated the good performance of the multi-level BASC on purely synthetic data. Stable networks were also derived from a real resting-state study for 43 subjects. At the group level, seven RSNs were identified which exhibited a good agreement with the previous findings from the literature. The comparison between the individual and group-level stability maps demonstrated the capacity of BASC to establish successful correspondences between these two levels of analysis and at the same time retain some interesting subject-specific characteristics, e.g. the specific involvement of subcortical regions in the visual and fronto-parietal networks for some subjects.

Asymmetries of multiple components of the rodent cerebrum have been described at various levels of organization. Yet, despite its ubiquitous nature, many confusing and sometimes contradictory reports regarding structural asymmetries in the rodent brain have been published. There is a need, therefore, for a whole-brain imaging analysis technique for asymmetry studies that is both accurate, reproducible and robust. To this end, a comprehensive three-dimensional examination of differences in brain structure in an inbred mouse strain was undertaken. The goal of this study was thus to use high-resolution magnetic resonance imaging to assess structural asymmetries in the adult C57Bl/6J mouse brain. Fixed brain T2-weighted images of 20 male C57Bl/6J mice were acquired on a 7T scanner at 32 microm isotropic resolution. We used voxel-based analyses to examine structural asymmetries throughout the whole mouse brain. The striatum, medial-posterior regions of the thalamus, and motor, sensorimotor, and visual cortex were found to be asymmetrical. The most significant asymmetry was found in the hippocampus and, specifically, the dentate gyrus. In each case, the left region was larger than the right. No other regions of the mouse brain showed structural asymmetry. The results in the dentate gyrus were confirmed using stereology, revealing a correlation of r=0.61 between magnetic resonance and stereological measures. Hippocampal, along with cortical asymmetry, has been discussed repeatedly in the literature, yet a clear pattern of directionality, until this point, has not been described. The findings of asymmetry in the striatum and absence of asymmetry in the rest of the brain are novel and show the advantage of using the whole-brain three-dimensional techniques developed herein for assessing asymmetry.

Spatial normalization, registration, and segmentation techniques for Magnetic Resonance Imaging (MRI) often use a target or template volume to facilitate processing, take advantage of prior information, and define a common coordinate system for analysis. In the neuroimaging literature, the MNI305 Talairach-like coordinate system is often used as a standard template. However, when studying pediatric populations, variation from the adult brain makes the MNI305 suboptimal for processing brain images of children. Morphological changes occurring during development render the use of age-appropriate templates desirable to reduce potential errors and minimize bias during processing of pediatric data. This paper presents the methods used to create unbiased, age-appropriate MRI atlas templates for pediatric studies that represent the average anatomy for the age range of 4.5-18.5 years, while maintaining a high level of anatomical detail and contrast. The creation of anatomical T1-weighted, T2-weighted, and proton density-weighted templates for specific developmentally important age-ranges, used data derived from the largest epidemiological, representative (healthy and normal) sample of the U.S. population, where each subject was carefully screened for medical and psychiatric factors and characterized using established neuropsychological and behavioral assessments. Use of these age-specific templates was evaluated by computing average tissue maps for gray matter, white matter, and cerebrospinal fluid for each specific age range, and by conducting an exemplar voxel-wise deformation-based morphometry study using 66 young (4.5-6.9 years) participants to demonstrate the benefits of using the age-appropriate templates. The public availability of these atlases/templates will facilitate analysis of pediatric MRI data and enable comparison of results between studies in a common standardized space specific to pediatric research.

Epileptic activity underlying idiopathic generalized epilepsy (IGE) is related to abnormal thalamo-cortical interactions. Our purpose was to map in vivo the organization of the thalamo-cortical network in IGE. We measured cortical thickness and thalamic volumes on MRI in 23 IGE patients with generalized tonic-clonic seizures only and 46 healthy controls. Significant correlations between thalamic volumes and cortical thickness were interpreted as thalamo-cortical network connections. In controls, thickness of frontal, limbic, and occipital regions was positively correlated with the thalamic volumes, corresponding to known anatomical connections from sacrificial tracer studies in primates and human in vivo DTI data. In patients, thalamo-cortical network correlations increased in fronto-central and parietal regions, but decreased in limbic areas. Group analysis revealed that, compared to controls, IGE patients had bilateral thalamic atrophy and widespread cortical thinning that was most prominent in fronto-central areas, with a prevalence of up to 40%. Duration of epilepsy affected negatively thalamic volumes and thickness of fronto-central and limbic cortices. These effects were significantly different from aging in controls. Patients with poorly controlled seizures showed an even faster progression in these neocortical regions. Fronto-centro-parietal atrophy in IGE is likely the effect of generalized seizure activity inducing thalamo-cortical network remodeling. On the other hand, limbic abnormalities may take place secondary to thalamic disconnection.

A recent study found differences in localised regions of the cortex between the YAC128 mouse model of Huntington's Disease (HD) and wild-type mice. There are, however, few tools to automatically examine shape differences in the cortices of mice. This paper describes an algorithm for automatically measuring cortical thickness across the entire cortex from MRI of fixed mouse brain specimens. An analysis of the variance of the method showed that, on average, a 50 microm (0.05 mm) localised difference in cortical thickness can be measured using MR scans. Applying these methods to 8-month-old YAC128 mouse model mice representing an early stage of HD, we found an increase in cortical thickness in the sensorimotor cortex, and also revealed regions wherein decreasing striatal volume correlated with increasing cortical thickness, indicating a potential compensatory response.