Summary Human functional MRI (fMRI) research primarily focuses on analyzing data averaged across groups, which limits the detail, specificity, and clinical utility of fMRI resting-state functional connectivity (RSFC) and task activation maps. To push our understanding of functional brain organization to the level of individual humans, we assembled a novel MRI dataset containing five hours of RSFC data, six hours of task fMRI, multiple structural MRIs, and neuropsychological tests from each of ten adults. Using these data, we generated ten high fidelity, individual-specific functional connectomes. This individual connectome approach revealed several new types of spatial and organizational variability in brain networks, including unique network features and topologies that corresponded with structural and task-derived brain features. We are releasing this highly-sampled, individual-focused dataset as a resource for neuroscientists, and we propose precision individual connectomics as a model for future work examining the organization of healthy and diseased individual human brains.

Magnetic resonance imaging (MRI) has transformed our understanding of the human brain through well-replicated mapping of abilities to specific structures (for example, lesion studies) and functions1–3 (for example, task functional MRI (fMRI)). Mental health research and care have yet to realize similar advances from MRI. A primary challenge has been replicating associations between inter-individual differences in brain structure or function and complex cognitive or mental health phenotypes (brain-wide association studies (BWAS)). Such BWAS have typically relied on sample sizes appropriate for classical brain mapping4 (the median neuroimaging study sample size is about 25), but potentially too small for capturing reproducible brain–behavioural phenotype associations5,6. Here we used three of the largest neuroimaging datasets currently available—with a total sample size of around 50,000 individuals—to quantify BWAS effect sizes and reproducibility as a function of sample size. BWAS associations were smaller than previously thought, resulting in statistically underpowered studies, inflated effect sizes and replication failures at typical sample sizes. As sample sizes grew into the thousands, replication rates began to improve and effect size inflation decreased. More robust BWAS effects were detected for functional MRI (versus structural), cognitive tests (versus mental health questionnaires) and multivariate methods (versus univariate). Smaller than expected brain–phenotype associations and variability across population subsamples can explain widespread BWAS replication failures. In contrast to non-BWAS approaches with larger effects (for example, lesions, interventions and within-person), BWAS reproducibility requires samples with thousands of individuals.

SUMMARY The cerebellum contains the majority of neurons in the human brain and is unique for its uniform cytoarchitecture, absence of aerobic glycolysis, and role in adaptive plasticity. Despite anatomical and physiological differences between the cerebellum and cerebral cortex, group-average functional connectivity studies have identified networks related to specific functions in both structures. Recently, precision functional mapping of individuals revealed that functional networks in the cerebral cortex exhibit measurable individual specificity. Using the highly-sampled Midnight Scan Club (MSC) dataset, we found the cerebellum contains reliable, individual-specific network organization that is significantly more variable than the cerebral cortex. The frontoparietal network, thought to support adaptive control, was the only network overrepresented in the cerebellum compared to the cerebral cortex (2.3-fold). Temporally, all cerebellar resting state signals lagged behind the cerebral cortex (125-380ms), supporting the hypothesis that the cerebellum engages in a domain-general function in the adaptive control of all cortical processes.

Magnetic resonance imaging (MRI) continues to drive many important neuroscientific advances. However, progress in uncovering reproducible associations between individual differences in brain structure/function and behavioral phenotypes (e.g., cognition, mental health) may have been undermined by typical neuroimaging sample sizes (median N=25)1,2. Leveraging the Adolescent Brain Cognitive Development (ABCD) Study3 (N=11,878), we estimated the effect sizes and reproducibility of these brain wide associations studies (BWAS) as a function of sample size. The very largest, replicable brain wide associations for univariate and multivariate methods were r=0.14 and r=0.34, respectively. In smaller samples, typical for brain wide association studies, irreproducible, inflated effect sizes were ubiquitous, no matter the method (univariate, multivariate). Until sample sizes started to approach consortium levels, BWAS were underpowered and statistical errors assured. Multiple factors contribute to replication failures4,5,6; here, we show that the pairing of small brain behavioral phenotype effect sizes with sampling variability is a key element in widespread BWAS replication failure. Brain behavioral phenotype associations stabilize and become more reproducible with sample sizes of N>2,000. While investigator initiated brain behavior research continues to generate hypotheses and propel innovation, large consortia are needed to usher in a new era of reproducible human brain wide association studies.

Highlights d Precision functional mapping identifies subcortical integration zones in individuals d Systematic connectivity reveals motor, cognitive, and visual attention integration zones d Subcortical functional organization has individually variable and conserved features d Integration zones map onto variably effective DBS sites, suggesting clinical utility