In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth.

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth.

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

The Brain's Topographical Organization Shapes Dynamic Interaction Patterns That Support Flexible Behavior Based on Rules and Long-Term Knowledge.

Adaptive behavior relies both on specific rules that vary across situations and stable long-term knowledge gained from experience. The frontoparietal control network (FPCN) is implicated in the brain's ability to balance these different influences on action. Here, we investigate how the topographical organization of the cortex supports behavioral flexibility within the FPCN. Functional properties of this network might reflect its juxtaposition between the dorsal attention network (DAN) and the default mode network (DMN), two large-scale systems implicated in top-down attention and memory-guided cognition, respectively. Our study tests whether subnetworks of FPCN are topographically proximal to the DAN and the DMN, respectively, and how these topographical differences relate to functional differences: the proximity of each subnetwork is anticipated to play a pivotal role in generating distinct cognitive modes relevant to working memory and long-term memory. We show that FPCN subsystems share multiple anatomical and functional similarities with their neighboring systems (DAN and DMN) and that this topographical architecture supports distinct interaction patterns that give rise to different patterns of functional behavior. The FPCN acts as a unified system when long-term knowledge supports behavior but becomes segregated into discrete subsystems with different patterns of interaction when long-term memory is less relevant. In this way, our study suggests that the topographical organization of the FPCN and the connections it forms with distant regions of cortex are important influences on how this system supports flexible behavior.

MRI economics: Balancing sample size and scan duration in brain wide association studies.

A pervasive dilemma in neuroimaging is whether to prioritize sample size or scan duration given fixed resources. Here, we systematically investigate this trade-off in the context of brain-wide association studies (BWAS) using resting-state functional magnetic resonance imaging (fMRI). We find that total scan duration (sample size × scan duration per participant) robustly explains individual-level phenotypic prediction accuracy via a logarithmic model, suggesting that sample size and scan duration are broadly interchangeable. The returns of scan duration eventually diminish relative to sample size, which we explain with principled theoretical derivations. When accounting for fixed costs associated with each participant (e.g., recruitment, non-imaging measures), we find that prediction accuracy in small-scale BWAS might benefit from much longer scan durations (>50 min) than typically assumed. Most existing large-scale studies might also have benefited from smaller sample sizes with longer scan durations. Both logarithmic and theoretical models of the relationships among sample size, scan duration and prediction accuracy explain well-predicted phenotypes better than poorly-predicted phenotypes. The logarithmic and theoretical models are also undermined by individual differences in brain states. These results replicate across phenotypic domains (e.g., cognition and mental health) from two large-scale datasets with different algorithms and metrics. Overall, our study emphasizes the importance of scan time, which is ignored in standard power calculations. Standard power calculations inevitably maximize sample size at the expense of scan duration. The resulting prediction accuracies are likely lower than would be produced with alternate designs, thus impeding scientific discovery. Our empirically informed reference is available for future study design: WEB_APPLICATION_LINK.

Macroscale intrinsic dynamics are associated with microcircuit function in focal and generalized epilepsies.

Epilepsies are a group of neurological disorders characterized by abnormal spontaneous brain activity, involving multiscale changes in brain functional organizations. However, it is not clear to what extent the epilepsy-related perturbations of spontaneous brain activity affect macroscale intrinsic dynamics and microcircuit organizations, that supports their pathological relevance. We collect a sample of patients with temporal lobe epilepsy (TLE) and genetic generalized epilepsy with tonic-clonic seizure (GTCS), as well as healthy controls. We extract massive temporal features of fMRI BOLD time-series to characterize macroscale intrinsic dynamics, and simulate microcircuit neuronal dynamics used a large-scale biological model. Here we show whether macroscale intrinsic dynamics and microcircuit dysfunction are differed in epilepsies, and how these changes are linked. Differences in macroscale gradient of time-series features are prominent in the primary network and default mode network in TLE and GTCS. Biophysical simulations indicate reduced recurrent connection within somatomotor microcircuits in both subtypes, and even more reduced in GTCS. We further demonstrate strong spatial correlations between differences in the gradient of macroscale intrinsic dynamics and microcircuit dysfunction in epilepsies. These results emphasize the impact of abnormal neuronal activity on primary network and high-order networks, suggesting a systematic abnormality of brain hierarchical organization.

Translating phenotypic prediction models from big to small anatomical MRI data using meta-matching.

Individualized phenotypic prediction based on structural MRI is an important goal in neuroscience. Prediction performance increases with larger samples, but small-scale datasets with fewer than 200 participants are often unavoidable. We have previously proposed a "meta-matching" framework to translate models trained from large datasets to improve the prediction of new unseen phenotypes in small collection efforts. Meta-matching exploits correlations between phenotypes, yielding large improvement over classical machine learning when applied to prediction models using resting-state functional connectivity as input features. Here, we adapt the two best performing meta-matching variants ("meta-matching finetune" and "meta-matching stacking") from our previous study to work with T1-weighted MRI data by changing the base neural network architecture to a 3D convolution neural network. We compare the two meta-matching variants with elastic net and classical transfer learning using the UK Biobank (N = 36,461), Human Connectome Project Young Adults (HCP-YA) dataset (N = 1,017) and HCP-Aging dataset (N = 656). We find that meta-matching outperforms elastic net and classical transfer learning by a large margin, both when translating models within the same dataset, as well as translating models across datasets with different MRI scanners, acquisition protocols and demographics. For example, when translating a UK Biobank model to 100 HCP-YA participants, meta-matching finetune yielded a 136% improvement in variance explained over transfer learning, with an average absolute gain of 2.6% (minimum = -0.9%, maximum = 17.6%) across 35 phenotypes. Overall, our results highlight the versatility of the meta-matching framework.

Multilayer meta-matching: translating phenotypic prediction models from multiple datasets to small data.

Resting-state functional connectivity (RSFC) is widely used to predict phenotypic traits in individuals. Large sample sizes can significantly improve prediction accuracies. However, for studies of certain clinical populations or focused neuroscience inquiries, small-scale datasets often remain a necessity. We have previously proposed a "meta-matching" approach to translate prediction models from large datasets to predict new phenotypes in small datasets. We demonstrated large improvement of meta-matching over classical kernel ridge regression (KRR) when translating models from a single source dataset (UK Biobank) to the Human Connectome Project Young Adults (HCP-YA) dataset. In the current study, we propose two meta-matching variants ("meta-matching with dataset stacking" and "multilayer meta-matching") to translate models from multiple source datasets across disparate sample sizes to predict new phenotypes in small target datasets. We evaluate both approaches by translating models trained from five source datasets (with sample sizes ranging from 862 participants to 36,834 participants) to predict phenotypes in the HCP-YA and HCP-Aging datasets. We find that multilayer meta-matching modestly outperforms meta-matching with dataset stacking. Both meta-matching variants perform better than the original "meta-matching with stacking" approach trained only on the UK Biobank. All meta-matching variants outperform classical KRR and transfer learning by a large margin. In fact, KRR is better than classical transfer learning when less than 50 participants are available for finetuning, suggesting the difficulty of classical transfer learning in the very small sample regime. The multilayer meta-matching model is publicly available at GITHUB_LINK.

Comparison between gradients and parcellations for functional connectivity prediction of behavior.

Resting-state functional connectivity (RSFC) is widely used to predict behavioral measures. To predict behavioral measures, representing RSFC with parcellations and gradients are the two most popular approaches. Here, we compare parcellation and gradient approaches for RSFC-based prediction of a broad range of behavioral measures in the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. Among the parcellation approaches, we consider group-average "hard" parcellations (Schaefer et al., 2018), individual-specific "hard" parcellations (Kong et al., 2021a), and an individual-specific "soft" parcellation (spatial independent component analysis with dual regression; Beckmann et al., 2009). For gradient approaches, we consider the well-known principal gradients (Margulies et al., 2016) and the local gradient approach that detects local RSFC changes (Laumann et al., 2015). Across two regression algorithms, individual-specific hard-parcellation performs the best in the HCP dataset, while the principal gradients, spatial independent component analysis and group-average "hard" parcellations exhibit similar performance. On the other hand, principal gradients and all parcellation approaches perform similarly in the ABCD dataset. Across both datasets, local gradients perform the worst. Finally, we find that the principal gradient approach requires at least 40 to 60 gradients to perform as well as parcellation approaches. While most principal gradient studies utilize a single gradient, our results suggest that incorporating higher order gradients can provide significant behaviorally relevant information. Future work will consider the inclusion of additional parcellation and gradient approaches for comparison.

Homotopic local-global parcellation of the human cerebral cortex from resting-state functional connectivity.

Resting-state fMRI is commonly used to derive brain parcellations, which are widely used for dimensionality reduction and interpreting human neuroscience studies. We previously developed a model that integrates local and global approaches for estimating areal-level cortical parcellations. The resulting local-global parcellations are often referred to as the Schaefer parcellations. However, the lack of homotopic correspondence between left and right Schaefer parcels has limited their use for brain lateralization studies. Here, we extend our previous model to derive homotopic areal-level parcellations. Using resting-fMRI and task-fMRI across diverse scanners, acquisition protocols, preprocessing and demographics, we show that the resulting homotopic parcellations are as homogeneous as the Schaefer parcellations, while being more homogeneous than five publicly available parcellations. Furthermore, weaker correlations between homotopic parcels are associated with greater lateralization in resting network organization, as well as lateralization in language and motor task activation. Finally, the homotopic parcellations agree with the boundaries of a number of cortical areas estimated from histology and visuotopic fMRI, while capturing sub-areal (e.g., somatotopic and visuotopic) features. Overall, these results suggest that the homotopic local-global parcellations represent neurobiologically meaningful subdivisions of the human cerebral cortex and will be a useful resource for future studies. Multi-resolution parcellations estimated from 1479 participants are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Yan2023_homotopic).

Comparison of individualized behavioral predictions across anatomical, diffusion and functional connectivity MRI.

A fundamental goal across the neurosciences is the characterization of relationships linking brain anatomy, functioning, and behavior. Although various MRI modalities have been developed to probe these relationships, direct comparisons of their ability to predict behavior have been lacking. Here, we compared the ability of anatomical T1, diffusion and functional MRI (fMRI) to predict behavior at an individual level. Cortical thickness, area and volume were extracted from anatomical T1 images. Diffusion Tensor Imaging (DTI) and approximate Neurite Orientation Dispersion and Density Imaging (NODDI) models were fitted to the diffusion images. The resulting metrics were projected to the Tract-Based Spatial Statistics (TBSS) skeleton. We also ran probabilistic tractography for the diffusion images, from which we extracted the stream count, average stream length, and the average of each DTI and NODDI metric across tracts connecting each pair of brain regions. Functional connectivity (FC) was extracted from both task and resting-state fMRI. Individualized prediction of a wide range of behavioral measures were performed using kernel ridge regression, linear ridge regression and elastic net regression. Consistency of the results were investigated with the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. In both datasets, FC-based models gave the best prediction performance, regardless of regression model or behavioral measure. This was especially true for the cognitive component. Furthermore, all modalities were able to predict cognition better than other behavioral components. Combining all modalities improved prediction of cognition, but not other behavioral components. Finally, across all behaviors, combining resting and task FC yielded prediction performance similar to combining all modalities. Overall, our study suggests that in the case of healthy children and young adults, behaviorally-relevant information in T1 and diffusion features might reflect a subset of the variance captured by FC.

Proportional intracranial volume correction differentially biases behavioral predictions across neuroanatomical features, sexes, and development.

Individual differences in brain anatomy can be used to predict variations in cognitive ability. Most studies to date have focused on broad population-level trends, but the extent to which the observed predictive features are shared across sexes and age groups remains to be established. While it is standard practice to account for intracranial volume (ICV) using proportion correction in both regional and whole-brain morphometric analyses, in the context of brain-behavior predictions the possible differential impact of ICV correction on anatomical features and subgroups within the population has yet to be systematically investigated. In this work, we evaluate the effect of proportional ICV correction on sex-independent and sex-specific predictive models of individual cognitive abilities across multiple anatomical properties (surface area, gray matter volume, and cortical thickness) in healthy young adults (Human Connectome Project; n = 1013, 548 females) and typically developing children (Adolescent Brain Cognitive Development study; n = 1823, 979 females). We demonstrate that ICV correction generally reduces predictive accuracies derived from surface area and gray matter volume, while increasing predictive accuracies based on cortical thickness in both adults and children. Furthermore, the extent to which predictive models generalize across sexes and age groups depends on ICV correction: models based on surface area and gray matter volume are more generalizable without ICV correction, while models based on cortical thickness are more generalizable with ICV correction. Finally, the observed neuroanatomical features predictive of cognitive abilities are unique across age groups regardless of ICV correction, but whether they are shared or unique across sexes (within age groups) depends on ICV correction. These findings highlight the importance of considering individual differences in ICV, and show that proportional ICV correction does not remove the effects of cranial volume from anatomical measurements and can introduce ICV bias where previously there was none. ICV correction choices affect not just the strength of the relationships captured, but also the conclusions drawn regarding the neuroanatomical features that underlie those relationships.

Sensory-motor cortices shape functional connectivity dynamics in the human brain.

Large-scale biophysical circuit models provide mechanistic insights into the micro-scale and macro-scale properties of brain organization that shape complex patterns of spontaneous brain activity. We developed a spatially heterogeneous large-scale dynamical circuit model that allowed for variation in local synaptic properties across the human cortex. Here we show that parameterizing local circuit properties with both anatomical and functional gradients generates more realistic static and dynamic resting-state functional connectivity (FC). Furthermore, empirical and simulated FC dynamics demonstrates remarkably similar sharp transitions in FC patterns, suggesting the existence of multiple attractors. Time-varying regional fMRI amplitude may track multi-stability in FC dynamics. Causal manipulation of the large-scale circuit model suggests that sensory-motor regions are a driver of FC dynamics. Finally, the spatial distribution of sensory-motor drivers matches the principal gradient of gene expression that encompasses certain interneuron classes, suggesting that heterogeneity in excitation-inhibition balance might shape multi-stability in FC dynamics.

Individual-Specific Areal-Level Parcellations Improve Functional Connectivity Prediction of Behavior.

Resting-state functional magnetic resonance imaging (rs-fMRI) allows estimation of individual-specific cortical parcellations. We have previously developed a multi-session hierarchical Bayesian model (MS-HBM) for estimating high-quality individual-specific network-level parcellations. Here, we extend the model to estimate individual-specific areal-level parcellations. While network-level parcellations comprise spatially distributed networks spanning the cortex, the consensus is that areal-level parcels should be spatially localized, that is, should not span multiple lobes. There is disagreement about whether areal-level parcels should be strictly contiguous or comprise multiple noncontiguous components; therefore, we considered three areal-level MS-HBM variants spanning these range of possibilities. Individual-specific MS-HBM parcellations estimated using 10 min of data generalized better than other approaches using 150 min of data to out-of-sample rs-fMRI and task-fMRI from the same individuals. Resting-state functional connectivity derived from MS-HBM parcellations also achieved the best behavioral prediction performance. Among the three MS-HBM variants, the strictly contiguous MS-HBM exhibited the best resting-state homogeneity and most uniform within-parcel task activation. In terms of behavioral prediction, the gradient-infused MS-HBM was numerically the best, but differences among MS-HBM variants were not statistically significant. Overall, these results suggest that areal-level MS-HBMs can capture behaviorally meaningful individual-specific parcellation features beyond group-level parcellations. Multi-resolution trained models and parcellations are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Kong2022_ArealMSHBM).

Correction: Time of day is associated with paradoxical reductions in global signal fluctuation and functional connectivity.

[This corrects the article DOI: 10.1371/journal.pbio.3000602.].

Heritability of individualized cortical network topography.

Human cortex is patterned by a complex and interdigitated web of large-scale functional networks. Recent methodological breakthroughs reveal variation in the size, shape, and spatial topography of cortical networks across individuals. While spatial network organization emerges across development, is stable over time, and is predictive of behavior, it is not yet clear to what extent genetic factors underlie interindividual differences in network topography. Here, leveraging a nonlinear multidimensional estimation of heritability, we provide evidence that individual variability in the size and topographic organization of cortical networks are under genetic control. Using twin and family data from the Human Connectome Project (n = 1,023), we find increased variability and reduced heritability in the size of heteromodal association networks (h2 : M = 0.34, SD = 0.070), relative to unimodal sensory/motor cortex (h2 : M = 0.40, SD = 0.097). We then demonstrate that the spatial layout of cortical networks is influenced by genetics, using our multidimensional estimation of heritability (h2-multi; M = 0.14, SD = 0.015). However, topographic heritability did not differ between heteromodal and unimodal networks. Genetic factors had a regionally variable influence on brain organization, such that the heritability of network topography was greatest in prefrontal, precuneus, and posterior parietal cortex. Taken together, these data are consistent with relaxed genetic control of association cortices relative to primary sensory/motor regions and have implications for understanding population-level variability in brain functioning, guiding both individualized prediction and the interpretation of analyses that integrate genetics and neuroimaging.

The detailed organization of the human cerebellum estimated by intrinsic functional connectivity within the individual.

Distinct regions of the cerebellum connect to separate regions of the cerebral cortex forming a complex topography. Although cerebellar organization has been examined in group-averaged data, study of individuals provides an opportunity to discover features that emerge at a higher spatial resolution. Here, functional connectivity MRI was used to examine the cerebellum of two intensively sampled individuals (each scanned 31 times). Connectivity to somatomotor cortex showed the expected crossed laterality and topography of the body maps. A surprising discovery was connectivity to the primary visual cortex along the vermis with evidence for representation of the central field. Within the hemispheres, each individual displayed a hierarchical progression from the inverted anterior lobe somatomotor map through to higher-order association zones. The hierarchy ended at Crus I/II and then progressed in reverse order through to the upright somatomotor map in the posterior lobe. Evidence for a third set of networks was found in the most posterior extent of the cerebellum. Detailed analysis of the higher-order association networks revealed robust representations of two distinct networks linked to the default network, multiple networks linked to cognitive control, as well as a separate representation of a language network. Although idiosyncratic spatial details emerged between subjects, each network could be detected in both individuals, and seed regions placed within the cerebellum recapitulated the full extent of the spatially specific cerebral networks. The observation of multiple networks in juxtaposed regions at the Crus I/II apex confirms the importance of this zone to higher-order cognitive function and reveals new organizational details.NEW & NOTEWORTHY Stable, within-individual maps of cerebellar organization reveal orderly macroscale representations of the cerebral cortex with local juxtaposed zones representing distinct networks. In addition, individuals reveal idiosyncratic organizational features.

Convergent molecular, cellular, and cortical neuroimaging signatures of major depressive disorder.

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.

Time of day is associated with paradoxical reductions in global signal fluctuation and functional connectivity.

The brain exhibits substantial diurnal variation in physiology and function, but neuroscience studies rarely report or consider the effects of time of day. Here, we examined variation in resting-state functional MRI (fMRI) in around 900 individuals scanned between 8 AM and 10 PM on two different days. Multiple studies across animals and humans have demonstrated that the brain's global signal (GS) amplitude (henceforth referred to as "fluctuation") increases with decreased arousal. Thus, in accord with known circadian variation in arousal, we hypothesised that GS fluctuation would be lowest in the morning, increase in the midafternoon, and dip in the early evening. Instead, we observed a cumulative decrease in GS fluctuation as the day progressed. Although respiratory variation also decreased with time of day, control analyses suggested that this did not account for the reduction in GS fluctuation. Finally, time of day was associated with marked decreases in resting-state functional connectivity across the whole brain. The magnitude of decrease was significantly stronger than associations between functional connectivity and behaviour (e.g., fluid intelligence). These findings reveal time of day effects on global brain activity that are not easily explained by expected arousal state or physiological artefacts. We conclude by discussing potential mechanisms for the observed diurnal variation in resting brain activity and the importance of accounting for time of day in future studies.

Deep neural networks and kernel regression achieve comparable accuracies for functional connectivity prediction of behavior and demographics.

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.