A central goal in systems neuroscience is the parcellation of the cerebral cortex into discrete neurobiological "atoms". Resting-state functional magnetic resonance imaging (rs-fMRI) offers the possibility of in-vivo human cortical parcellation. Almost all previous parcellations relied on one of two approaches. The local gradient approach detects abrupt transitions in functional connectivity patterns. These transitions potentially reflect cortical areal boundaries defined by histology or visuotopic fMRI. By contrast, the global similarity approach clusters similar functional connectivity patterns regardless of spatial proximity, resulting in parcels with homogeneous (similar) rs-fMRI signals. Here we propose a gradient-weighted Markov Random Field (gwMRF) model integrating local gradient and global similarity approaches.Using task-fMRI and rs-fMRI across diverse acquisition protocols, we found gwMRF parcellations to be more homogeneous than four previously published parcellations.Furthermore, gwMRF parcellations agreed with the boundaries of certain cortical areas defined using histology and visuotopic fMRI. Some parcels captured sub-areal (somatotopic and visuotopic) features that likely reflect distinct computational units within known cortical areas. These results suggest that gwMRF parcellations reveal neurobiologically meaningful features of brain organization and are potentially useful for future applications requiring dimensionality reduction of voxel-wise fMRI data. Multi-resolution parcellations generated from 1489 participants are available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Scha efer2018_LocalGlobal)

Resting-state functional magnetic resonance imaging (rs-fMRI) offers the opportunity to delineate individual-specific brain networks. A major question is whether individual-specific network topography (i.e., location and spatial arrangement) is behaviorally relevant. Here, we propose a multi-session hierarchical Bayesian model (MS-HBM) for estimating individualspecific cortical networks and investigate whether individual-specific network topography can predict human behavior. The multiple layers of the MS-HBM explicitly differentiate intra-subject (within-subject) from inter-subject (between-subject) network variability. By ignoring intra-subject variability, previous network mappings might confuse intra-subject variability for inter-subject differences. Compared with other approaches, MS-HBM parcellations generalized better to new rs-fMRI and task-fMRI data from the same subjects.More specifically, MS-HBM parcellations estimated from a single rs-fMRI session (10 minutes) showed comparable generalizability as parcellations estimated by two state-of-theart methods using five sessions (50 minutes). We also showed that behavioral phenotypes across cognition, personality and emotion could be predicted by individual-specific network topography with modest accuracy, comparable to previous reports predicting phenotypes based on connectivity strength. Network topography estimated by MS-HBM was more effective for behavioral prediction than network size, as well as network topography estimated by other parcellation approaches. Thus, similar to connectivity strength, individualspecific network topography might also serve as a fingerprint of human behavior. substantially across participants (Harrison et al., 2015;Laumann et al., 2015;Wang et al., 2015;Glasser et al., 2016;Gordon et al., 2017aGordon et al., , 2017bGordon et al., , 2017cBraga and Buckner, 2017).Yet, the possible behavioral relevance of individual differences in network size and network topography (location and spatial arrangement) remains unknown.We proposed a multi-session hierarchical Bayesian model (MS-HBM) for estimating individual-specific network parcellations of the cerebral cortex and investigated whether individual-specific network topography and size are associated with human behavior. The multiple layers of the MS-HBM allowed explicit separation of inter-subject (betweensubject) and intra-subject (within-session) functional connectivity variability. Previous individual-specific network mappings only accounted for inter-subject variability, but not intra-subject variability. However, inter-subject and intra-subject RSFC variability can be markedly different across regions (Mueller et al., 2013;Chen et al., 2015;Laumann et al., 2015). For example, the motor cortex exhibits high intra-subject functional connectivity variability, but low inter-subject functional connectivity variability .Therefore, observed RSFC variability in the motor cortex might be incorrectly attributed to inter-subject spatial variability of brain networks, rather than just intra-...

There is significant interest in the development and application of deep neural networks (DNNs) to neuroimaging data. A growing literature suggests that DNNs outperform their classical counterparts in a variety of neuroimaging applications, yet there are few direct comparisons of relative utility. Here, we compared the performance of three DNN architectures and a classical machine learning algorithm (kernel regression) in predicting individual phenotypes from whole-brain resting-state functional connectivity (RSFC) patterns. One of the DNNs was a generic fully-connected feedforward neural network, while the other two DNNs were recently published approaches specifically designed to exploit the structure of connectome data. By using a combined sample of almost 10,000 participants from the Human Connectome Project (HCP) and UK Biobank, we showed that the three DNNs and kernel regression achieved similar performance across a wide range of behavioral and demographic measures. Furthermore, the generic feedforward neural network exhibited similar performance to the two state-of-the-art connectome-specific DNNs. When predicting fluid intelligence in the UK Biobank, performance of all algorithms dramatically improved when sample size increased from 100 to 1000 subjects. Improvement was smaller, but still significant, when sample size increased from 1000 to 5000 subjects. Importantly, kernel regression was competitive across all sample sizes. Overall, our study suggests that kernel regression is as effective as DNNs for RSFC-based behavioral prediction, while incurring significantly lower computational costs. Therefore, kernel regression might serve as a useful baseline algorithm for future studies.

How individual differences in brain network organization track behavioral variability is a fundamental question in systems neuroscience. Recent work suggests that resting-state and task-state functional connectivity can predict specific traits at the individual level. However, most studies focus on single behavioral traits, thus not capturing broader relationships across behaviors. In a large sample of 1858 typically developing children from the Adolescent Brain Cognitive Development (ABCD) study, we show that predictive network features are distinct across the domains of cognitive performance, personality scores and mental health assessments. On the other hand, traits within each behavioral domain are predicted by similar network features. Predictive network features and models generalize to other behavioral measures within the same behavioral domain. Although tasks are known to modulate the functional connectome, predictive network features are similar between resting and task states. Overall, our findings reveal shared brain network features that account for individual variation within broad domains of behavior in childhood.

Individuals vary widely in their tendency to seek stimulation and act impulsively, early developing traits with genetic origins. Failures to regulate these behaviors increase risk for maladaptive outcomes including substance abuse. Here, we explored the neuroanatomical correlates of sensation seeking and impulsivity in healthy young adults. Our analyses revealed links between sensation seeking and reduced cortical thickness that were preferentially localized to regions implicated in cognitive control, including anterior cingulate and middle frontal gyrus (n ϭ 1015). These associations generalized to self-reported motor impulsivity, replicated in an independent group (n ϭ 219), and correlated with heightened alcohol, tobacco, and caffeine use. Critically, the relations between sensation seeking and brain structure were evident in participants without a historyofalcoholortobaccouse,suggestingthatobservedassociationswithanatomyarenotsolelyaconsequenceofsubstanceuse.Theseresults demonstrate that individual differences in the tendency to seek stimulation, act on impulse, and engage in substance use are correlated with the anatomical structure of cognitive control circuitry. Our findings suggest that, in healthy populations, covariation across these complex multidimensional behaviors may in part originate from a common underlying biology.

The human brain is comprised of a complex web of functional networks that link anatomically distinct regions. However, the biological mechanisms supporting network organization remain elusive, particularly across cortical and subcortical territories with vastly divergent cellular and molecular properties. Here, using human and primate brain transcriptional atlases, we demonstrate that spatial patterns of gene expression show strong correspondence with limbic and somato/motor cortico-striatal functional networks. Network-associated expression is consistent across independent human datasets and evolutionarily conserved in non-human primates. Genes preferentially expressed within the limbic network (encompassing nucleus accumbens, orbital/ventromedial prefrontal cortex, and temporal pole) relate to risk for psychiatric illness, chloride channel complexes, and markers of somatostatin neurons. Somato/motor associated genes are enriched for oligodendrocytes and markers of parvalbumin neurons. These analyses indicate that parallel cortico-striatal processing channels possess dissociable genetic signatures that recapitulate distributed functional networks, and nominate molecular mechanisms supporting cortico-striatal circuitry in health and disease.

Major depressive disorder emerges from the complex interactions of biological systems that span genes and molecules through cells, networks, and behavior. Establishing how neurobiological processes coalesce to contribute to depression requires a multiscale approach, encompassing measures of brain structure and function as well as genetic and cell-specific transcriptional data. Here, we examine anatomical (cortical thickness) and functional (functional variability, global brain connectivity) correlates of depression and negative affect across three population-imaging datasets: UK Biobank, Brain Genomics Superstruct Project, and Enhancing NeuroImaging through Meta Analysis (ENIGMA; combined n ≥ 23,723). Integrative analyses incorporate measures of cortical gene expression, postmortem patient transcriptional data, depression genome-wide association study (GWAS), and single-cell gene transcription. Neuroimaging correlates of depression and negative affect were consistent across three independent datasets. Linking ex vivo gene down-regulation with in vivo neuroimaging, we find that transcriptional correlates of depression imaging phenotypes track gene down-regulation in postmortem cortical samples of patients with depression. Integrated analysis of single-cell and Allen Human Brain Atlas expression data reveal somatostatin interneurons and astrocytes to be consistent cell associates of depression, through both in vivo imaging and ex vivo cortical gene dysregulation. Providing converging evidence for these observations, GWAS-derived polygenic risk for depression was enriched for genes expressed in interneurons, but not glia. Underscoring the translational potential of multiscale approaches, the transcriptional correlates of depression-linked brain function and structure were enriched for disorder-relevant molecular pathways. These findings bridge levels to connect specific genes, cell classes, and biological pathways to in vivo imaging correlates of depression.

Converging evidence indicates that groups of patients with nominally distinct psychiatric diagnoses are not separated by sharp or discontinuous neurobiological boundaries. In healthy populations, individual differences in behavior are reflected in variability across the collective set of functional brain connections (functional connectome). These data suggest that the spectra of transdiagnostic symptom profiles observed in psychiatric patients may map onto detectable patterns of network function. To examine the manner through which neurobiological variation might underlie clinical presentation, we obtained fMRI data from over 1,000 individuals, including 210 diagnosed with a primary psychotic disorder or affective psychosis (bipolar disorder with psychosis and schizophrenia or schizoaffective disorder), 192 presenting with a primary affective disorder without psychosis (unipolar depression, bipolar disorder without psychosis), and 608 demographically matched healthy comparison participants recruited through a large-scale study of brain imaging and genetics. Here, we examine variation in functional connectomes across psychiatric diagnoses, finding striking evidence for disease connectomic “fingerprints” that are commonly disrupted across distinct forms of pathology and appear to scale as a function of illness severity. The presence of affective and psychotic illnesses was associated with graded disruptions in frontoparietal network connectivity (encompassing aspects of dorsolateral prefrontal, dorsomedial prefrontal, lateral parietal, and posterior temporal cortices). Conversely, other properties of network connectivity, including default network integrity, were preferentially disrupted in patients with psychotic illness, but not patients without psychotic symptoms. This work allows us to establish key biological and clinical features of the functional connectomes of severe mental disease.