brainlife.io: a decentralized and open-source cloud platform to support neuroscience research.

Neuroscience is advancing standardization and tool development to support rigor and transparency. Consequently, data pipeline complexity has increased, hindering FAIR (findable, accessible, interoperable and reusable) access. brainlife.io was developed to democratize neuroimaging research. The platform provides data standardization, management, visualization and processing and automatically tracks the provenance history of thousands of data objects. Here, brainlife.io is described and evaluated for validity, reliability, reproducibility, replicability and scientific utility using four data modalities and 3,200 participants.

Relation of connectome topology to brain volume across 103 mammalian species.

The brain connectome is an embedded network of anatomically interconnected brain regions, and the study of its topological organization in mammals has become of paramount importance due to its role in scaffolding brain function and behavior. Unlike many other observable networks, brain connections incur material and energetic cost, and their length and density are volumetrically constrained by the skull. Thus, an open question is how differences in brain volume impact connectome topology. We address this issue using the MaMI database, a diverse set of mammalian connectomes reconstructed from 201 animals, covering 103 species and 12 taxonomy orders, whose brain size varies over more than 4 orders of magnitude. Our analyses focus on relationships between volume and modular organization. After having identified modules through a multiresolution approach, we observed how connectivity features relate to the modular structure and how these relations vary across brain volume. We found that as the brain volume increases, modules become more spatially compact and dense, comprising more costly connections. Furthermore, we investigated how spatial embedding shapes network communication, finding that as brain volume increases, nodes' distance progressively impacts communication efficiency. We identified modes of variation in network communication policies, as smaller and bigger brains show higher efficiency in routing- and diffusion-based signaling, respectively. Finally, bridging network modularity and communication, we found that in larger brains, modular structure imposes stronger constraints on network signaling. Altogether, our results show that brain volume is systematically related to mammalian connectome topology and that spatial embedding imposes tighter restrictions on larger brains.

Partial entropy decomposition reveals higher-order information structures in human brain activity.

The standard approach to modeling the human brain as a complex system is with a network, where the basic unit of interaction is a pairwise link between two brain regions. While powerful, this approach is limited by the inability to assess higher-order interactions involving three or more elements directly. In this work, we explore a method for capturing higher-order dependencies in multivariate data: the partial entropy decomposition (PED). Our approach decomposes the joint entropy of the whole system into a set of nonnegative atoms that describe the redundant, unique, and synergistic interactions that compose the system's structure. PED gives insight into the mathematics of functional connectivity and its limitation. When applied to resting-state fMRI data, we find robust evidence of higher-order synergies that are largely invisible to standard functional connectivity analyses. Our approach can also be localized in time, allowing a frame-by-frame analysis of how the distributions of redundancies and synergies change over the course of a recording. We find that different ensembles of regions can transiently change from being redundancy-dominated to synergy-dominated and that the temporal pattern is structured in time. These results provide strong evidence that there exists a large space of unexplored structures in human brain data that have been largely missed by a focus on bivariate network connectivity models. This synergistic structure is dynamic in time and likely will illuminate interesting links between brain and behavior. Beyond brain-specific application, the PED provides a very general approach for understanding higher-order structures in a variety of complex systems.

Functional hierarchies in brain dynamics characterized by signal reversibility in ferret cortex.

Brain signal irreversibility has been shown to be a promising approach to study neural dynamics. Nevertheless, the relation with cortical hierarchy and the influence of different electrophysiological features is not completely understood. In this study, we recorded local field potentials (LFPs) during spontaneous behavior, including awake and sleep periods, using custom micro-electrocorticographic (µECoG) arrays implanted in ferrets. In contrast to humans, ferrets remain less time in each state across the sleep-wake cycle. We deployed a diverse set of metrics in order to measure the levels of complexity of the different behavioral states. In particular, brain irreversibility, which is a signature of non-equilibrium dynamics, captured by the arrow of time of the signal, revealed the hierarchical organization of the ferret's cortex. We found different signatures of irreversibility and functional hierarchy of large-scale dynamics in three different brain states (active awake, quiet awake, and deep sleep), showing a lower level of irreversibility in the deep sleep stage, compared to the other. Irreversibility also allowed us to disentangle the influence of different cortical areas and frequency bands in this process, showing a predominance of the parietal cortex and the theta band. Furthermore, when inspecting the embedded dynamic through a Hidden Markov Model, the deep sleep stage was revealed to have a lower switching rate and lower entropy production. These results suggest functional hierarchies in organization that can be revealed through thermodynamic features and information theory metrics.

Structural-functional brain network coupling predicts human cognitive ability.

Individual differences in general cognitive ability (GCA) have a biological basis within the structure and function of the human brain. Network neuroscience investigations revealed neural correlates of GCA in structural as well as in functional brain networks. However, whether the relationship between structural and functional networks, the structural-functional brain network coupling (SC-FC coupling), is related to individual differences in GCA remains an open question. We used data from 1030 adults of the Human Connectome Project, derived structural connectivity from diffusion weighted imaging, functional connectivity from resting-state fMRI, and assessed GCA as a latent g-factor from 12 cognitive tasks. Two similarity measures and six communication measures were used to model possible functional interactions arising from structural brain networks. SC-FC coupling was estimated as the degree to which these measures align with the actual functional connectivity, providing insights into different neural communication strategies. At the whole-brain level, higher GCA was associated with higher SC-FC coupling, but only when considering path transitivity as neural communication strategy. Taking region-specific variations in the SC-FC coupling strategy into account and differentiating between positive and negative associations with GCA, allows for prediction of individual cognitive ability scores in a cross-validated prediction framework (correlation between predicted and observed scores: r = 0.25, p < .001). The same model also predicts GCA scores in a completely independent sample (N = 567, r = 0.19, p < .001). Our results propose structural-functional brain network coupling as a neurobiological correlate of GCA and suggest brain region-specific coupling strategies as neural basis of efficient information processing predictive of cognitive ability.

Reliability of variability and complexity measures for task and task-free BOLD fMRI.

Brain activity continuously fluctuates over time, even if the brain is in controlled (e.g., experimentally induced) states. Recent years have seen an increasing interest in understanding the complexity of these temporal variations, for example with respect to developmental changes in brain function or between-person differences in healthy and clinical populations. However, the psychometric reliability of brain signal variability and complexity measures-which is an important precondition for robust individual differences as well as longitudinal research-is not yet sufficiently studied. We examined reliability (split-half correlations) and test-retest correlations for task-free (resting-state) BOLD fMRI as well as split-half correlations for seven functional task data sets from the Human Connectome Project to evaluate their reliability. We observed good to excellent split-half reliability for temporal variability measures derived from rest and task fMRI activation time series (standard deviation, mean absolute successive difference, mean squared successive difference), and moderate test-retest correlations for the same variability measures under rest conditions. Brain signal complexity estimates (several entropy and dimensionality measures) showed moderate to good reliabilities under both, rest and task activation conditions. We calculated the same measures also for time-resolved (dynamic) functional connectivity time series and observed moderate to good reliabilities for variability measures, but poor reliabilities for complexity measures derived from functional connectivity time series. Global (i.e., mean across cortical regions) measures tended to show higher reliability than region-specific variability or complexity estimates. Larger subcortical regions showed similar reliability as cortical regions, but small regions showed lower reliability, especially for complexity measures. Lastly, we also show that reliability scores are only minorly dependent on differences in scan length and replicate our results across different parcellation and denoising strategies. These results suggest that the variability and complexity of BOLD activation time series are robust measures well-suited for individual differences research. Temporal variability of global functional connectivity over time provides an important novel approach to robustly quantifying the dynamics of brain function. PRACTITIONER POINTS: Variability and complexity measures of BOLD activation show good split-half reliability and moderate test-retest reliability. Measures of variability of global functional connectivity over time can robustly quantify neural dynamics. Length of fMRI data has only a minor effect on reliability.

Uncovering individual differences in fine-scale dynamics of functional connectivity.

Functional connectivity (FC) profiles contain subject-specific features that are conserved across time and have potential to capture brain-behavior relationships. Most prior work has focused on spatial features (nodes and systems) of these FC fingerprints, computed over entire imaging sessions. We propose a method for temporally filtering FC, which allows selecting specific moments in time while also maintaining the spatial pattern of node-based activity. To this end, we leverage a recently proposed decomposition of FC into edge time series (eTS). We systematically analyze functional magnetic resonance imaging frames to define features that enhance identifiability across multiple fingerprinting metrics, similarity metrics, and data sets. Results show that these metrics characteristically vary with eTS cofluctuation amplitude, similarity of frames within a run, transition velocity, and expression of functional systems. We further show that data-driven optimization of features that maximize fingerprinting metrics isolates multiple spatial patterns of system expression at specific moments in time. Selecting just 10% of the data can yield stronger fingerprints than are obtained from the full data set. Our findings support the idea that FC fingerprints are differentially expressed across time and suggest that multiple distinct fingerprints can be identified when spatial and temporal characteristics are considered simultaneously.

Co-evolving dynamics and topology in a coupled oscillator model of resting brain function.

Dynamic models of ongoing BOLD fMRI brain dynamics and models of communication strategies have been two important approaches to understanding how brain network structure constrains function. However, dynamic models have yet to widely incorporate one of the most important insights from communication models: the brain may not use all of its connections in the same way or at the same time. Here we present a variation of a phase delayed Kuramoto coupled oscillator model that dynamically limits communication between nodes on each time step. An active subgraph of the empirically derived anatomical brain network is chosen in accordance with the local dynamic state on every time step, thus coupling dynamics and network structure in a novel way. We analyze this model with respect to its fit to empirical time-averaged functional connectivity, finding that, with the addition of only one parameter, it significantly outperforms standard Kuramoto models with phase delays. We also perform analyses on the novel time series of active edges it produces, demonstrating a slowly evolving topology moving through intermittent episodes of integration and segregation. We hope to demonstrate that the exploration of novel modeling mechanisms and the investigation of dynamics of networks in addition to dynamics on networks may advance our understanding of the relationship between brain structure and function.

Few temporally distributed brain connectivity states predict human cognitive abilities.

Human functional brain connectivity can be temporally decomposed into states of high and low cofluctuation, defined as coactivation of brain regions over time. Rare states of particularly high cofluctuation have been shown to reflect fundamentals of intrinsic functional network architecture and to be highly subject-specific. However, it is unclear whether such network-defining states also contribute to individual variations in cognitive abilities - which strongly rely on the interactions among distributed brain regions. By introducing CMEP, a new eigenvector-based prediction framework, we show that as few as 16 temporally separated time frames (< 1.5% of 10 min resting-state fMRI) can significantly predict individual differences in intelligence (N = 263, p < .001). Against previous expectations, individual's network-defining time frames of particularly high cofluctuation do not predict intelligence. Multiple functional brain networks contribute to the prediction, and all results replicate in an independent sample (N = 831). Our results suggest that although fundamentals of person-specific functional connectomes can be derived from few time frames of highest connectivity, temporally distributed information is necessary to extract information about cognitive abilities. This information is not restricted to specific connectivity states, like network-defining high-cofluctuation states, but rather reflected across the entire length of the brain connectivity time series.

Multitask brain network reconfiguration is inversely associated with human intelligence.

Intelligence describes the general cognitive ability level of a person. It is one of the most fundamental concepts in psychological science and is crucial for the effective adaption of behavior to varying environmental demands. Changing external task demands have been shown to induce reconfiguration of functional brain networks. However, whether neural reconfiguration between different tasks is associated with intelligence has not yet been investigated. We used functional magnetic resonance imaging data from 812 subjects to show that higher scores of general intelligence are related to less brain network reconfiguration between resting state and seven different task states as well as to network reconfiguration between tasks. This association holds for all functional brain networks except the motor system and replicates in two independent samples (n = 138 and n = 184). Our findings suggest that the intrinsic network architecture of individuals with higher intelligence scores is closer to the network architecture as required by various cognitive demands. Multitask brain network reconfiguration may, therefore, represent a neural reflection of the behavioral positive manifold - the essence of the concept of general intelligence. Finally, our results support neural efficiency theories of cognitive ability and reveal insights into human intelligence as an emergent property from a distributed multitask brain network.

High-amplitude cofluctuations in cortical activity drive functional connectivity.

Resting-state functional connectivity is used throughout neuroscience to study brain organization and to generate biomarkers of development, disease, and cognition. The processes that give rise to correlated activity are, however, poorly understood. Here we decompose resting-state functional connectivity using a temporal unwrapping procedure to assess the contributions of moment-to-moment activity cofluctuations to the overall connectivity pattern. This approach temporally resolves functional connectivity at a timescale of single frames, which enables us to make direct comparisons of cofluctuations of network organization with fluctuations in the blood oxygen level-dependent (BOLD) time series. We show that surprisingly, only a small fraction of frames exhibiting the strongest cofluctuation amplitude are required to explain a significant fraction of variance in the overall pattern of connection weights as well as the network's modular structure. These frames coincide with frames of high BOLD activity amplitude, corresponding to activity patterns that are remarkably consistent across individuals and identify fluctuations in default mode and control network activity as the primary driver of resting-state functional connectivity. Finally, we demonstrate that cofluctuation amplitude synchronizes across subjects during movie watching and that high-amplitude frames carry detailed information about individual subjects (whereas low-amplitude frames carry little). Our approach reveals fine-scale temporal structure of resting-state functional connectivity and discloses that frame-wise contributions vary across time. These observations illuminate the relation of brain activity to functional connectivity and open a number of directions for future research.

Individualized event structure drives individual differences in whole-brain functional connectivity.

Resting-state functional connectivity is typically modeled as the correlation structure of whole-brain regional activity. It is studied widely, both to gain insight into the brain's intrinsic organization but also to develop markers sensitive to changes in an individual's cognitive, clinical, and developmental state. Despite this, the origins and drivers of functional connectivity, especially at the level of densely sampled individuals, remain elusive. Here, we leverage novel methodology to decompose functional connectivity into its precise framewise contributions. Using two dense sampling datasets, we investigate the origins of individualized functional connectivity, focusing specifically on the role of brain network "events" - short-lived and peaked patterns of high-amplitude cofluctuations. Here, we develop a statistical test to identify events in empirical recordings. We show that the patterns of cofluctuation expressed during events are repeated across multiple scans of the same individual and represent idiosyncratic variants of template patterns that are expressed at the group level. Lastly, we propose a simple model of functional connectivity based on event cofluctuations, demonstrating that group-averaged cofluctuations are suboptimal for explaining participant-specific connectivity. Our work complements recent studies implicating brief instants of high-amplitude cofluctuations as the primary drivers of static, whole-brain functional connectivity. Our work also extends those studies, demonstrating that cofluctuations during events are individualized, positing a dynamic basis for functional connectivity.

Multi-modal and multi-subject modular organization of human brain networks.

The human brain is a complex network of anatomically interconnected brain areas. Spontaneous neural activity is constrained by this architecture, giving rise to patterns of statistical dependencies between the activity of remote neural elements. The non-trivial relationship between structural and functional connectivity poses many unsolved challenges about cognition, disease, development, learning and aging. While numerous studies have focused on statistical relationships between edge weights in anatomical and functional networks, less is known about dependencies between their modules and communities. In this work, we investigate and characterize the relationship between anatomical and functional modular organization of the human brain, developing a novel multi-layer framework that expands the classical concept of multi-layer modularity. By simultaneously mapping anatomical and functional networks estimated from different subjects into communities, this approach allows us to carry out a multi-subject and multi-modal analysis of the brain's modular organization. Here, we investigate the relationship between anatomical and functional modules during resting state, finding unique and shared structures. The proposed framework constitutes a methodological advance in the context of multi-layer network analysis and paves the way to further investigate the relationship between structural and functional network organization in clinical cohorts, during cognitively demanding tasks, and in developmental or lifespan studies.

Cortico-subcortical interactions in overlapping communities of edge functional connectivity.

Both cortical and subcortical regions can be functionally organized into networks. Regions of the basal ganglia are extensively interconnected with the cortex via reciprocal connections that relay and modulate cortical function. Here we employ an edge-centric approach, which computes co-fluctuations among region pairs in a network to investigate the role and interaction of subcortical regions with cortical systems. By clustering edges into communities, we show that cortical systems and subcortical regions couple via multiple edge communities, with hippocampus and amygdala having a distinct pattern from striatum and thalamus. We show that the edge community structure of cortical networks is highly similar to one obtained from cortical nodes when the subcortex is present in the network. Additionally, we show that the edge community profile of both cortical and subcortical nodes can be estimates solely from cortico-subcortical interactions. Finally, we used a motif analysis focusing on edge community triads where a subcortical region coupled to two cortical regions and found that two community triads where one community couples the subcortex to the cortex were overrepresented. In summary, our results show organized coupling of the subcortex to the cortex that may play a role in cortical organization of primary sensorimotor/attention and heteromodal systems and puts forth the motif analysis of edge community triads as a promising method for investigation of communication patterns in networks.

Effects of diffusion signal modeling and segmentation approaches on subthalamic nucleus parcellation.

The subthalamic nucleus (STN) is commonly used as a surgical target for deep brain stimulation in movement disorders such as Parkinson's Disease. Tractography-derived connectivity-based parcellation (CBP) has been recently proposed as a suitable tool for non-invasive in vivo identification and pre-operative targeting of specific functional territories within the human STN. However, a well-established, accurate and reproducible protocol for STN parcellation is still lacking. The present work aims at testing the effects of different tractography-based approaches for the reconstruction of STN functional territories. We reconstructed functional territories of the STN on the high-quality dataset of 100 unrelated healthy subjects and on the test-retest dataset of the Human Connectome Project (HCP) repository. Connectivity-based parcellation was performed with a hypothesis-driven approach according to cortico-subthalamic connectivity, after dividing cortical areas into three groups: associative, limbic and sensorimotor. Four parcellation pipelines were compared, combining different signal modeling techniques (single-fiber vs multi-fiber) and different parcellation approaches (winner takes all parcellation vs fiber density thresholding). We tested these procedures on STN regions of interest obtained from three different, commonly employed, subcortical atlases. We evaluated the pipelines both in terms of between-subject similarity, assessed on the cohort of 100 unrelated healthy subjects, and of within-subject similarity, using a second cohort of 44 subjects with available test-retest data. We found that each parcellation provides converging results in terms of location of the identified parcels, but with significative variations in size and shape. All pipelines obtained very high within-subject similarity, with tensor-based approaches outperforming multi-fiber pipelines. On the other hand, higher between-subject similarity was found with multi-fiber signal modeling techniques combined with fiber density thresholding. We suggest that a fine-tuning of tractography-based parcellation may lead to higher reproducibility and aid the development of an optimized surgical targeting protocol.

Mapping individual differences across brain network structure to function and behavior with connectome embedding.

The connectome, a comprehensive map of the brain's anatomical connections, is often summarized as a matrix comprising all dyadic connections among pairs of brain regions. This representation cannot capture higher-order relations within the brain graph. Connectome embedding (CE) addresses this limitation by creating compact vectorized representations of brain nodes capturing their context in the global network topology. Here, nodes "context" is defined as random walks on the brain graph and as such, represents a generative model of diffusive communication around nodes. Applied to group-averaged structural connectivity, CE was previously shown to capture relations between inter-hemispheric homologous brain regions and uncover putative missing edges from the network reconstruction. Here we extend this framework to explore individual differences with a novel embedding alignment approach. We test this approach in two lifespan datasets (NKI: n = 542; Cam-CAN: n = 601) that include diffusion-weighted imaging, resting-state fMRI, demographics and behavioral measures. We demonstrate that modeling functional connectivity with CE substantially improves structural to functional connectivity mapping both at the group and subject level. Furthermore, age-related differences in this structure-function mapping, are preserved and enhanced. Importantly, CE captures individual differences by out-of-sample prediction of age and intelligence. The resulting predictive accuracy was higher compared to using structural connectivity and functional connectivity. We attribute these findings to the capacity of the CE to incorporate aspects of both anatomy (the structural graph) and function (diffusive communication). Our novel approach allows mapping individual differences in the connectome through structure to function and behavior.

Subject identification using edge-centric functional connectivity.

Group-level studies do not capture individual differences in network organization, an important prerequisite for understanding neural substrates shaping behavior and for developing interventions in clinical conditions. Recent studies have employed 'fingerprinting' analyses on functional connectivity to identify subjects' idiosyncratic features. Here, we develop a complementary approach based on an edge-centric model of functional connectivity, which focuses on the co-fluctuations of edges. We first show whole-brain edge functional connectivity (eFC) to be a robust substrate that improves identifiability over nodal FC (nFC) across different datasets and parcellations. Next, we characterize subjects' identifiability at different spatial scales, from single nodes to the level of functional systems and clusters using k-means clustering. Across spatial scales, we find that heteromodal brain regions exhibit consistently greater identifiability than unimodal, sensorimotor, and limbic regions. Lastly, we show that identifiability can be further improved by reconstructing eFC using specific subsets of its principal components. In summary, our results highlight the utility of the edge-centric network model for capturing meaningful subject-specific features and sets the stage for future investigations into individual differences using edge-centric models.

Modularity maximization as a flexible and generic framework for brain network exploratory analysis.

The modular structure of brain networks supports specialized information processing, complex dynamics, and cost-efficient spatial embedding. Inter-individual variation in modular structure has been linked to differences in performance, disease, and development. There exist many data-driven methods for detecting and comparing modular structure, the most popular of which is modularity maximization. Although modularity maximization is a general framework that can be modified and reparamaterized to address domain-specific research questions, its application to neuroscientific datasets has, thus far, been narrow. Here, we highlight several strategies in which the "out-of-the-box" version of modularity maximization can be extended to address questions specific to neuroscience. First, we present approaches for detecting "space-independent" modules and for applying modularity maximization to signed matrices. Next, we show that the modularity maximization frame is well-suited for detecting task- and condition-specific modules. Finally, we highlight the role of multi-layer models in detecting and tracking modules across time, tasks, subjects, and modalities. In summary, modularity maximization is a flexible and general framework that can be adapted to detect modular structure resulting from a wide range of hypotheses. This article highlights multiple frontiers for future research and applications.

Aging relates to a disproportionately weaker functional architecture of brain networks during rest and task states.

Functional connectivity - the co-activation of brain regions - forms the basis of the brain's functional architecture. Often measured during resting-state (i.e., in a task-free setting), patterns of functional connectivity within and between brain networks change with age. These patterns are of interest to aging researchers because age differences in resting-state connectivity relate to older adults' relative cognitive declines. Less is known about age differences in large-scale brain networks during directed tasks. Recent work in younger adults has shown that patterns of functional connectivity are highly correlated between rest and task states. Whether this finding extends to older adults remains largely unexplored. To this end, we assessed younger and older adults' functional connectivity across the whole brain using fMRI while participants underwent resting-state or completed directed tasks (e.g., a reasoning judgement task). Resting-state and task functional connectivity were less strongly correlated in older as compared to younger adults. This age-dependent difference could be attributed to significantly lower consistency in network organization between rest and task states among older adults. Older adults had less distinct or segregated networks during resting-state. This more diffuse pattern of organization was exacerbated during directed tasks. Finally, the default mode network, often implicated in neurocognitive aging, contributed strongly to this pattern. These findings establish that age differences in functional connectivity are state-dependent, providing greater insight into the mechanisms by which aging may lead to cognitive declines.