Richard F. Betzel

The development of new technologies for mapping structural and functional brain connectivity has led to the creation of comprehensive network maps of neuronal circuits and systems. The architecture of these brain networks can be examined and analyzed with a large variety of graph theory tools. Methods for detecting modules, or network communities, are of particular interest because they uncover major building blocks or subnetworks that are particularly densely connected, often corresponding to specialized functional components. A large number of methods for community detection have become available and are now widely applied in network neuroscience. This article first surveys a number of these methods, with an emphasis on their advantages and shortcomings; then it summarizes major findings on the existence of modules in both structural and functional brain networks and briefly considers their potential functional roles in brain evolution, wiring minimization, and the emergence of functional specialization and complex dynamics.

At rest, the brain's sensorimotor and higher cognitive systems engage in organized patterns of correlated activity forming resting-state networks. An important empirical question is how functional connectivity and structural connectivity within and between resting-state networks change with age. In this study we use network modeling techniques to identify significant changes in network organization across the human lifespan. The results of this study demonstrate that whole-brain functional and structural connectivity both exhibit reorganization with age. On average, functional connections within resting-state networks weaken in magnitude while connections between resting-state networks tend to increase. These changes can be localized to a small subset of functional connections that exhibit systematic changes across the lifespan. Collectively, changes in functional connectivity are also manifest at a system-wide level, as components of the control, default mode, saliency/ventral attention, dorsal attention, and visual networks become less functionally cohesive, as evidenced by decreased component modularity. Paralleling this functional reorganization is a decrease in the density and weight of anatomical white-matter connections. Hub regions are particularly affected by these changes, and the capacity of those regions to communicate with other regions exhibits a lifelong pattern of decline. Finally, the relationship between functional connectivity and structural connectivity also appears to change with age; functional connectivity along multi-step structural paths tends to be stronger in older subjects than in younger subjects. Overall, our analysis points to age-related changes in inter-regional communication unfolding within and between resting-state networks.

The complex relationship between structural and functional connectivity, as measured by noninvasive imaging of the human brain, poses many unresolved challenges and open questions. Here, we apply analytic measures of network communication to the structural connectivity of the human brain and explore the capacity of these measures to predict resting-state functional connectivity across three independently acquired datasets. We focus on the layout of shortest paths across the network and on two communication measures--search information and path transitivity--which account for how these paths are embedded in the rest of the network. Search information is an existing measure of information needed to access or trace shortest paths; we introduce path transitivity to measure the density of local detours along the shortest path. We find that both search information and path transitivity predict the strength of functional connectivity among both connected and unconnected node pairs. They do so at levels that match or significantly exceed path length measures, Euclidean distance, as well as computational models of neural dynamics. This capacity suggests that dynamic couplings due to interactions among neural elements in brain networks are substantially influenced by the broader network context adjacent to the shortest communication pathways.

The emergence of network neuroscience allows researchers to quantify the link between the organizational features of neuronal networks and the spectrum of cortical functions.Current models indicate that structure and function are significantly correlated, but the correspondence is not perfect because function reflects complex multisynaptic interactions in structural networks.Function cannot be directly estimated from structure, but must be inferred by models of higher-order interactions. Statistical, communication, and biophysical models have been used to translate brain structure to brain function.Structure–function coupling is regionally heterogeneous and follows molecular, cytoarchitectonic, and functional hierarchies. Structure–function relationships are a fundamental principle of many naturally occurring systems. However, network neuroscience research suggests that there is an imperfect link between structural connectivity and functional connectivity in the brain. Here, we synthesize the current state of knowledge linking structure and function in macroscale brain networks and discuss the different types of models used to assess this relationship. We argue that current models do not include the requisite biological detail to completely predict function. Structural network reconstructions enriched with local molecular and cellular metadata, in concert with more nuanced representations of functions and properties, hold great potential for a truly multiscale understanding of the structure–function relationship. Structure–function relationships are a fundamental principle of many naturally occurring systems. However, network neuroscience research suggests that there is an imperfect link between structural connectivity and functional connectivity in the brain. Here, we synthesize the current state of knowledge linking structure and function in macroscale brain networks and discuss the different types of models used to assess this relationship. We argue that current models do not include the requisite biological detail to completely predict function. Structural network reconstructions enriched with local molecular and cellular metadata, in concert with more nuanced representations of functions and properties, hold great potential for a truly multiscale understanding of the structure–function relationship. The relationship between structure and function is a central concept in natural sciences and engineering. Consider how the conformation of a protein determines its chemical properties and, ultimately, its biological function. The folding of the protein into a 3D structure promotes interactions among amino acids, allowing the protein to chemically interact with other molecules and endowing it with function. Conversely, disruption of the protein’s structure results in loss of function. Tellingly, the protein is said to be denatured, highlighting the idea that changing its structure has fundamentally altered its natural function. The function of the nervous system is analogously shaped by the structure and arrangement of neurons and neuronal populations. The complex network of synaptic projections forms a hierarchy (see Glossary) of nested and increasingly polyfunctional neural circuits that support perception, cognition, and action. Modern imaging technology permits high-throughput reconstruction of neural circuits across spatiotemporal scales and across species (Box 1). Through extensive international data sharing efforts, increasingly detailed reconstructions of the nervous system’s connection patterns are available in humans and in multiple model organisms, including invertebrate [1.Chiang A-S. et al.Three-dimensional reconstruction of brain-wide wiring networks in Drosophila at single-cell resolution.Curr. Biol. 2011; 21: 1-11Abstract Full Text Full Text PDF PubMed Scopus (324) Google Scholar], avian [2.Shanahan M. et al.Large-scale network organization in the avian forebrain: a connectivity matrix and theoretical analysis.Front. Comput. Neurosci. 2013; 7: 89Crossref PubMed Scopus (111) Google Scholar], rodent [3.Oh S.W. et al.A mesoscale connectome of the mouse brain.Nature. 2014; 508: 207Crossref PubMed Scopus (777) Google Scholar,4.Bota M. et al.Architecture of the cerebral cortical association connectome underlying cognition.Proc. Natl. Acad. Sci. U.S.A. 2015; 112: E2093-E2101Crossref PubMed Scopus (90) Google Scholar, and primate species [5.Markov N.T. et al.A weighted and directed interareal connectivity matrix for macaque cerebral cortex.Cereb. Cortex. 2012; 24: 17-36Crossref PubMed Scopus (272) Google Scholar,6.Majka P. et al.Towards a comprehensive atlas of cortical connections in a primate brain: Mapping tracer injection studies of the common marmoset into a reference digital template.J. Com. Neurol. 2016; 524: 2161-2181Crossref PubMed Scopus (33) Google Scholar. These comprehensive wiring diagrams of the nervous system, termed structural connectivity (SC) networks or connectomes, represent the physical connections between neural elements [7.Sporns O. et al.The human connectome: a structural description of the human brain.PLoS Comput. Biol. 2005; 1: e42Crossref PubMed Scopus (1521) Google Scholar]. The emergence of network neuroscience offers an opportunity to quantify and articulate the link between the organizational features of neuronal networks and the spectrum of cortical functions. SC networks possess distinctive and nonrandom attributes, including high local clustering and short path length, characteristic features of a small-world architecture [8.Watts D.J. Strogatz S.H. Collective dynamics of small-world networks.Nature. 1998; 393: 440Crossref PubMed Scopus (0) Google Scholar]. Populations with similar functional properties tend to cluster together, forming specialized modules crosslinked by hub nodes with diverse connectional fingerprints [9.Young M.P. The organization of neural systems in the primate cerebral cortex.J. Roy. Soc. Lond. B. 1993; 252: 13-18Crossref PubMed Google Scholar,10.Kötter R. et al.Connectional characteristics of areas in Walker’s map of primate prefrontal cortex.Neurocomputing. 2001; 38: 741-746Crossref Scopus (20) Google Scholar. The hubs of the networks are disproportionately interconnected with each other, forming a putative core [11.Hagmann P. et al.Mapping the structural core of human cerebral cortex.PLoS Biol. 2008; 6: e159Crossref PubMed Scopus (2384) Google Scholar] or ‘rich club’ [12.van den Heuvel M.P. et al.High-cost, high-capacity backbone for global brain communication.Proc. Natl. Acad. Sci. U.S.A. 2012; 109: 11372-11377Crossref PubMed Scopus (361) Google Scholar], an architectural feature that potentially allows signals to be sampled and integrated from specialized modules [13.Zamora-López G. et al.Cortical hubs form a module for multisensory integration on top of the hierarchy of cortical networks.Front. Neuroinform. 2010; 4: 1PubMed Google Scholar]. Finally, brain networks are spatially embedded, with finite metabolic and material resources [14.Bullmore E. Sporns O. The economy of brain network organization.Nat. Rev. Neurosci. 2012; 13: 336Crossref PubMed Scopus (1282) Google Scholar], resulting in increased prevalence of shorter, low-cost connections [15.Horvát S. et al.Spatial embedding and wiring cost constrain the functional layout of the cortical network of rodents and primates.PLoS Biol. 2016; 14: e1002512Crossref PubMed Google Scholar,16.Roberts J.A. et al.The contribution of geometry to the human connectome.NeuroImage. 2016; 124: 379-393Crossref PubMed Scopus (60) Google Scholar. These organizational attributes have been replicated across a range of species and tracing techniques, suggesting common organizational principles across phylogeny [17.Van den Heuvel M.P. et al.Comparative connectomics.Trends. Cogn. Sci. 2016; 20: 345-361Abstract Full Text Full Text PDF PubMed Scopus (118) Google Scholar]. The architecture of SC networks imparts a distinct signature on neuronal coactivation patterns. Inter-regional projections promote signaling and synchrony among distant neuronal populations, giving rise to coherent neural dynamics, measured as regional time series of electromagnetic or hemodynamic neural activity. Systematic coactivation among pairs of regions can be used to map functional connectivity (FC) networks. Over the past decade, these dynamics are increasingly recorded without task instruction or stimulation; the resulting ‘intrinsic’ or ‘ resting-state’ FC is thought to reflect spontaneous neural activity [18.Biswal B. et al.Functional connectivity in the motor cortex of resting human brain using echo-planar MRI.Magn. Reson. Med. 1995; 34: 537-541Crossref PubMed Scopus (5503) Google Scholar]. Intrinsic FC patterns are highly organized [19.Damoiseaux J. et al.Consistent resting-state networks across healthy subjects.Proc. Natl. Acad. Sci. U.S.A. 2006; 103: 13848-13853Crossref PubMed Scopus (2665) Google Scholar, 20.Bellec P. et al.Multi-level bootstrap analysis of stable clusters in resting-state fMRI.NeuroImage. 2010; 51: 1126-1139Crossref PubMed Scopus (170) Google Scholar, 21.Thomas Yeo B. et al.The organization of the human cerebral cortex estimated by intrinsic functional connectivity.J. Neurophysiol. 2011; 106: 1125-1165Crossref PubMed Scopus (2040) Google Scholar], reproducible [22.Gordon E.M. et al.Precision functional mapping of individual human brains.Neuron. 2017; 95: 791-807Abstract Full Text Full Text PDF PubMed Scopus (140) Google Scholar,23.Noble S. et al.A decade of test-retest reliability of functional connectivity: a systematic review and meta-analysis.NeuroImage. 2019; : 116157Crossref PubMed Scopus (5) Google Scholar, and comparable with task-driven coactivation patterns [24.Smith S.M. et al.Correspondence of the brain’s functional architecture during activation and rest.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 13040-13045Crossref PubMed Scopus (2684) Google Scholar,25.Cole M.W. et al.Intrinsic and task-evoked network architectures of the human brain.Neuron. 2014; 83: 238-251Abstract Full Text Full Text PDF PubMed Scopus (516) Google Scholar. The persistent and reproducible nature of brain activity during rest makes resting-state FC an ideal starting point to study structure–function relationships [26.Honey C.J. et al.Can structure predict function in the human brain?.NeuroImage. 2010; 52: 766-776Crossref PubMed Scopus (291) Google Scholar,27.Damoiseaux J.S. Greicius M.D. Greater than the sum of its parts: a review of studies combining structural connectivity and resting-state functional connectivity.Brain Struct. Funct. 2009; 213: 525-533Crossref PubMed Google Scholar. Here we synthesize the current state of knowledge linking structure and function in macroscale brain networks. We first show that direct one-to-one links between structure and function are limited and inherently obscured by the networked nature of the brain. We survey modern quantitative methods that move away from direct correlations between structure and function by conceptualizing function as emerging from higher-order interactions among multiple neuronal populations, with a focus on strengths, limitations, and commonalities. We posit that the next steps in understanding network-level structure–function relationships must take into account regional heterogeneity by enriching network reconstructions with microscale attributes, including transcriptomic, cytoarchitectonic, and neuromodulatory information. We close by highlighting emerging theories that macroscale structure–function relationships are not uniform across the brain, but vary in parallel with cytoarchitectonic and representational hierarchies. Early studies emphasized correlations between structural and functional connection weights. Structural weights are correlated with functional weights [28.Honey C. et al.Predicting human resting-state functional connectivity from structural connectivity.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 2035-2040Crossref PubMed Scopus (1543) Google Scholar], and nodes that are central to structural networks also tend to be central in functional networks [28.Honey C. et al.Predicting human resting-state functional connectivity from structural connectivity.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 2035-2040Crossref PubMed Scopus (1543) Google Scholar]. Furthermore, structurally connected pairs of neural elements display greater FC than structurally unconnected pairs [28.Honey C. et al.Predicting human resting-state functional connectivity from structural connectivity.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 2035-2040Crossref PubMed Scopus (1543) Google Scholar,29.Shen K. et al.Information processing architecture of functionally defined clusters in the macaque cortex.J. Neurosci. 2012; 32: 17465-17476Crossref PubMed Scopus (63) Google Scholar (Figure 1A). More globally, many intrinsic functional networks, particularly the visual and somatomotor networks, are circumscribed by patterns of dense anatomical connectivity [29.Shen K. et al.Information processing architecture of functionally defined clusters in the macaque cortex.J. Neurosci. 2012; 32: 17465-17476Crossref PubMed Scopus (63) Google Scholar, 30.Van Den Heuvel M.P. et al.Functionally linked resting-state networks reflect the underlying structural connectivity architecture of the human brain.Hum. Brain Mapp. 2009; 30: 3127-3141Crossref PubMed Scopus (625) Google Scholar, 31.Alves P.N. et al.An improved neuroanatomical model of the default-mode network reconciles previous neuroimaging and neuropathological findings.Commun. Biol. 2019; 2: 1-14Crossref PubMed Scopus (0) Google Scholar]. While SC and FC are significantly correlated, the correspondence is not perfect. Even the best-case estimates place the correlation at R2 ≈ 0.5 [28.Honey C. et al.Predicting human resting-state functional connectivity from structural connectivity.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 2035-2040Crossref PubMed Scopus (1543) Google Scholar], which means that considerable variance (at least half) in functional connection weights is unexplained by a simple 1:1 correspondence with structure. The discrepancy widens in the case of functional connections between regions that are not structurally connected (Figure 1B ). A particularly salient example is the case of homotopic functional connections between corresponding structures in the two hemispheres. These are typically the strongest subset of functional connections [32.Mišić B. et al.The functional connectivity landscape of the human brain.PLoS One. 2014; 9: e111007Crossref PubMed Scopus (14) Google Scholar], but not all homotopic functional connections are supported by a direct underlying callosal projection [33.Shen K. et al.Stable long-range interhemispheric coordination is supported by direct anatomical projections.Proc. Natl. Acad. Sci. U.S.A. 2015; 112: 6473-6478Crossref PubMed Scopus (52) Google Scholar], and strong homotopic functional connections may be observed even in individuals with no callosal connections [34.Uddin L.Q. et al.Residual functional connectivity in the split-brain revealed with resting-state fMRI.Neuroreport. 2008; 19: 703Crossref PubMed Scopus (96) Google Scholar, 35.O’Reilly J.X. et al.Causal effect of disconnection lesions on interhemispheric functional connectivity in rhesus monkeys.Proc. Natl. Acad. Sci. U.S.A. 2013; 110: 13982-13987Crossref PubMed Scopus (106) Google Scholar, 36.Layden E.A. et al.Interhemispheric functional connectivity in the zebra finch brain, absent the corpus callosum in normal ontogeny.NeuroImage. 2019; Crossref PubMed Scopus (1) Google Scholar]. These examples illustrate the simple point that sustained communication via indirect anatomical connections may manifest as strong FC. The discordance between structure and function is particularly pronounced at the mesoscopic scale. Intrinsic networks commonly observed in resting-state FC and meta-analytic coactivation cannot be recovered from structural networks [37.Mišić B. et al.Cooperative and competitive spreading dynamics on the human connectome.Neuron. 2015; 86: 1518-1529Abstract Full Text Full Text PDF PubMed Google Scholar,38.Betzel R.F. et al.Diversity of meso-scale architecture in human and non-human connectomes.Nat. Commun. 2018; 9: 346Crossref PubMed Scopus (21) Google Scholar (Figure 1C). While intrinsic networks can be reproducibly defined using independent component analysis [19.Damoiseaux J. et al.Consistent resting-state networks across healthy subjects.Proc. Natl. Acad. Sci. U.S.A. 2006; 103: 13848-13853Crossref PubMed Scopus (2665) Google Scholar], community detection [39.Power J.D. et al.Functional network organization of the human brain.Neuron. 2011; 72: 665-678Abstract Full Text Full Text PDF PubMed Scopus (1499) Google Scholar], or data-driven clustering [20.Bellec P. et al.Multi-level bootstrap analysis of stable clusters in resting-state fMRI.NeuroImage. 2010; 51: 1126-1139Crossref PubMed Scopus (170) Google Scholar,21.Thomas Yeo B. et al.The organization of the human cerebral cortex estimated by intrinsic functional connectivity.J. Neurophysiol. 2011; 106: 1125-1165Crossref PubMed Scopus (2040) Google Scholar, both in resting state recordings and in meta-analytic coactivation [24.Smith S.M. et al.Correspondence of the brain’s functional architecture during activation and rest.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 13040-13045Crossref PubMed Scopus (2684) Google Scholar], application of comparable methods to diffusion-weighted SC or anatomical covariance networks yields networks that are more spatially contiguous [37.Mišić B. et al.Cooperative and competitive spreading dynamics on the human connectome.Neuron. 2015; 86: 1518-1529Abstract Full Text Full Text PDF PubMed Google Scholar,40.Betzel R.F. et al.The modular organization of human anatomical brain networks: accounting for the cost of wiring.Net. Neurosci. 2017; 1: 42-68Crossref Scopus (33) Google Scholar. For example, clustering or community detection methods typically fail to identify a default mode-like structural network, perhaps because not all parts of the network are anatomically inter-connected [26.Honey C.J. et al.Can structure predict function in the human brain?.NeuroImage. 2010; 52: 766-776Crossref PubMed Scopus (291) Google Scholar]. Structural and functional networks also show global organizational differences. For example, structural networks show evidence of extensive assortative mixing, whereby nodes with similar properties (e.g., degrees) are more likely to be connected, whereas the same is not true of functional networks [50.Lim S. et al.Discordant attributes of structural and functional brain connectivity in a two-layer multiplex network.Sci. Rep. 2019; 9: 2885Crossref PubMed Scopus (2) Google Scholar]. At the mesoscopic scale, communities or modules recovered from structural networks are assortative, while communities recovered from functional networks are disassortative [38.Betzel R.F. et al.Diversity of meso-scale architecture in human and non-human connectomes.Nat. Commun. 2018; 9: 346Crossref PubMed Scopus (21) Google Scholar]. In other words, in functional networks there is a pronounced affinity among nodes with dissimilar attributes. As a result, tuning community detection algorithms to be sensitive to disassortative structures improves the match between structural and functional modules [38.Betzel R.F. et al.Diversity of meso-scale architecture in human and non-human connectomes.Nat. Commun. 2018; 9: 346Crossref PubMed Scopus (21) Google Scholar]. Altogether, a rich body of work demonstrates discordance between SC and FC that spans multiple scales, from the embedding of individual nodes and edges to their global arrangement. Why the discrepancy between SC and FC? Functional interactions may arise via indirect structural connections, resulting in coherent time courses among regions that are two or more synapses removed from each other. In other words, the propensity of two regional time courses to correlate is driven not only by direct signaling between them, but also by the common inputs they receive from sensory organs and from the entire network [27.Damoiseaux J.S. Greicius M.D. Greater than the sum of its parts: a review of studies combining structural connectivity and resting-state functional connectivity.Brain Struct. Funct. 2009; 213: 525-533Crossref PubMed Google Scholar,51.Bettinardi R.G. et al.How structure sculpts function: unveiling the contribution of anatomical connectivity to the brain’s spontaneous correlation structure.Chaos. 2017; 27: 047409Crossref PubMed Scopus (12) Google Scholar. A corollary is that functional interactions are much less distance-dependent than structural connections. Anatomical wiring is subject to material, spatial, and metabolic constraints [14.Bullmore E. Sporns O. The economy of brain network organization.Nat. Rev. Neurosci. 2012; 13: 336Crossref PubMed Scopus (1282) Google Scholar]; these pressures manifest in reduced connection probability and connection weight with increasing spatial separation [15.Horvát S. et al.Spatial embedding and wiring cost constrain the functional layout of the cortical network of rodents and primates.PLoS Biol. 2016; 14: e1002512Crossref PubMed Google Scholar,16.Roberts J.A. et al.The contribution of geometry to the human connectome.NeuroImage. 2016; 124: 379-393Crossref PubMed Scopus (60) Google Scholar. Although similar distance-dependence is observed for FC, its effect is weaker, ensuring systematic differences between structural and functional configurations. In the next section we consider models that translate structure to function by conceptualizing function as an emergent property of multiple structural links. As we have seen so far, there exists a nontrivial link between SC and FC, but the two are not perfectly aligned. A number of models have emerged that embody this link, including statistical models [41.Mišić B. et al.Network-level structure-function relationships in human neocortex.Cereb. Cortex. 2016; 26: 3285-3296Crossref PubMed Scopus (153) Google Scholar,42.Messé A. et al.Relating structure and function in the human brain: relative contributions of anatomy, stationary dynamics, and non-stationarities.PLoS. Comput. Biol. 2014; 10: e1003530Crossref PubMed Scopus (90) Google Scholar, communication models [37.Mišić B. et al.Cooperative and competitive spreading dynamics on the human connectome.Neuron. 2015; 86: 1518-1529Abstract Full Text Full Text PDF PubMed Google Scholar,43.Graham D. Rockmore D. The packet switching brain.J. Cogn. Neurosci. 2011; 23: 267-276Crossref PubMed Scopus (18) Google Scholar, 44.Goñi J. et al.Resting-brain functional connectivity predicted by analytic measures of network communication.Proc. Natl. Acad. Sci. U.S.A. 2014; 111: 833-838Crossref PubMed Scopus (208) Google Scholar, 45.Crofts J.J. Higham D.J. A weighted communicability measure applied to complex brain networks.J. Roy. Soc. Interf. 2009; 6: 411-414Crossref PubMed Scopus (61) Google Scholar and biophysical models [46.Honey C.J. et al.Network structure of cerebral cortex shapes functional connectivity on multiple time scales.Proc. Natl. Acad. Sci. U.S.A. 2007; 104: 10240-10245Crossref PubMed Scopus (941) Google Scholar, 47.Breakspear M. Dynamic models of large-scale brain activity.Nat. Neurosci. 2017; 20: 340Crossref PubMed Scopus (147) Google Scholar, 48.Sanz-Leon P. et al.Mathematical framework for large-scale brain network modeling in The Virtual Brain.NeuroImage. 2015; 111: 385-430Crossref PubMed Google Scholar, 49.Deco G. et al.Key role of coupling, delay, and noise in resting brain fluctuations.Proc. Natl. Acad. Sci. U.S.A. 2009; 106: 10302-10307Crossref PubMed Scopus (372) Google Scholar]. Though different in their implementation and assumptions, the common idea is to emphasize collective, higher-order interactions among neural elements that transcends the strong local clustering and geometric dependence of dyadic structural relationships. Here we briefly review each of these analytic strategies, with a focus on their biological interpretation and predictive utility, and most importantly, what they teach us about the nature of the structure–function relationship in the brain. Perhaps the simplest way to link structure and function is statistically. Varying forms of reduced rank regression have emerged as particularly useful, including canonical correlation [52.Deligianni F. et al.NODDI and tensor-based microstructural indices as predictors of functional connectivity.PLoS One. 2016; 11: e0153404Crossref PubMed Scopus (13) Google Scholar] and partial least squares [41.Mišić B. et al.Network-level structure-function relationships in human neocortex.Cereb. Cortex. 2016; 26: 3285-3296Crossref PubMed Scopus (153) Google Scholar]. In these data-driven models the objective is to simultaneously identify weighted combinations of structural and functional connections that are maximally correlated across individuals [53.McIntosh A.R. Mišić B. Multivariate statistical analyses for neuroimaging data.Annu. Rev. Psychol. 2013; 64: 499-525Crossref PubMed Scopus (73) Google Scholar] (Figure 2). An appealing feature of such models is that they embody multiple structure–function modes. In other words, a particular structural configuration or subnetwork may give rise to distinct patterns of functional interactions [41.Mišić B. et al.Network-level structure-function relationships in human neocortex.Cereb. Cortex. 2016; 26: 3285-3296Crossref PubMed Scopus (153) Google Scholar]. Taking this idea further, artificial neural networks can be used to learn functional networks from structural networks. For example, a recent study used a variant of the word2vec algorithm to build a low-dimensional embedding representation of the connectome and used it to train a deep neural network to predict edge-wise FC [54.Rosenthal G. et al.Mapping higher-order relations between brain structure and function with embedded vector representations of connectomes.Nat. Commun. 2018; 9: 2178Crossref PubMed Scopus (3) Google Scholar]. Altogether, statistical models offer a data-driven way to associate combinations of structural and functional connections without assuming a specific mode of interaction among neuronal populations. Communication models emerging from network science and telecommunication engineering conceptualize functional interactions as the superposition of elementary signaling events on the underlying anatomical network [43.Graham D. Rockmore D. The packet switching brain.J. Cogn. Neurosci. 2011; 23: 267-276Crossref PubMed Scopus (18) Google Scholar,55.Avena-Koenigsberger A. et al.Communication dynamics in complex brain networks.Nat. Rev. Neurosci. 2018; 19: 17Crossref Scopus (92) Google Scholar (Figure 2). By explicitly formulating a model of inter-regional signaling, these models open two important questions, namely: how biologically realistic is the model, and how well does the model fit the attributes of the functional network? Early studies focused on centralized forms of communication, such as shortest path routing, whereby discrete signals travel via the shortest contiguous set of edges from a source node to a prespecified target node. More recently, attention has shifted to decentralized mechanisms where signals diffuse through the network [56.Mišić B. et al.A network convergence zone in the hippocampus.PLoS. Comput. Biol. 2014; 10: e1003982Crossref PubMed Google Scholar,57.Atasoy S. et al.Human brain networks function in connectome-specific harmonic waves.Nat. Commun. 2016; 7: 10340Crossref PubMed Google Scholar, often broadcast on multiple fronts [37.Mišić B. et al.Cooperative and competitive spreading dynamics on the human connectome.Neuron. 2015; 86: 1518-1529Abstract Full Text Full Text PDF PubMed Google Scholar,51.Bettinardi R.G. et al.How structure sculpts function: unveiling the contribution of anatomical connectivity to the brain’s spontaneous correlation structure.Chaos. 2017; 27: 047409Crossref PubMed Scopus (12) Google Scholar,58.Abdelnour F. et al.Network diffusion accurately models the relationship between structural and functional brain connectivity networks.NeuroImage. 2014; 90: 335-347Crossref PubMed Scopus (71) Google Scholar,59.Worrell J.C. et al.Optimized connectome architecture for sensory-motor integration.Net. Neurosci. 2017; 1: 415-430Crossref Google Scholar. Others have considered mechanisms that are neither fully centralized nor decentralized, including communication via path ensembles [45.Crofts J.J. Higham D.J. A weighted communicability measure applied to complex brain networks.J. Roy. Soc. Interf. 2009; 6: 411-414Crossref PubMed Scopus (61) Google Scholar,60.Avena-Koenigsberger A. et al.Path ensembles and a tradeoff between communication efficiency and resilience in the human connectome.Brain Struct. Funct. 2017; 222: 603-618Crossref PubMed Scopus (17) Google Scholar or multiplexed strategies involving multiple mechanisms [44.Goñi J. et al.Resting-brain functional connectivity predicted by analytic measures of network communication.Proc. Natl. Acad. Sci. U.S.A. 2014; 111: 833-838Crossref PubMed Scopus (208) Google Scholar,61.Avena-Koenigsberger A. et al.A spectrum of routing strategies for brain networks.PLoS Comput. Biol. 2019; 15: e1006833Crossref PubMed Scopus (5) Google Scholar, 62.Betzel R.F. et al.Structural, geometric and genetic factors predict interregional brain connectivity patterns probed by electrocorticography.Nat. Biomed. Eng. 2019; : 1PubMed Google Scholar, 63.Vazquez-Rodriguez B. et al.Gradients of structure-function tethering across neocortex.Proc. Natl. Acad. Sci. U.S.A. 2019; 116: 21219-21227Crossref PubMed Scopus (6) Google Scholar. An emerging consensus is that, given the strong relationship between geometric and topological proximity in the brain, it is possible for decentralized mechanisms to utilize the shortest path structure of the network, either by diffusion [44.Goñi J. et al.Re

The network architecture of the human brain has become a feature of increasing interest to the neuroscientific community, largely because of its potential to illuminate human cognition, its variation over development and aging, and its alteration in disease or injury. Traditional tools and approaches to study this architecture have largely focused on single scales—of topology, time, and space. Expanding beyond this narrow view, we focus this review on pertinent questions and novel methodological advances for the multi-scale brain. We separate our exposition into content related to multi-scale topological structure, multi-scale temporal structure, and multi-scale spatial structure. In each case, we recount empirical evidence for such structures, survey network-based methodological approaches to reveal these structures, and outline current frontiers and open questions. Although predominantly peppered with examples from human neuroimaging, we hope that this account will offer an accessible guide to any neuroscientist aiming to measure, characterize, and understand the full richness of the brain's multiscale network structure—irrespective of species, imaging modality, or spatial resolution.

The brain is a complex, multiscale dynamical system composed of many interacting regions. Knowledge of the spatiotemporal organization of these interactions is critical for establishing a solid understanding of the brain’s functional architecture and the relationship between neural dynamics and cognition in health and disease. The possibility of studying these dynamics through careful analysis of neuroimaging data has catalyzed substantial interest in methods that estimate time-resolved fluctuations in functional connectivity (often referred to as “dynamic” or time-varying functional connectivity; TVFC). At the same time, debates have emerged regarding the application of TVFC analyses to resting fMRI data, and about the statistical validity, physiological origins, and cognitive and behavioral relevance of resting TVFC. These and other unresolved issues complicate interpretation of resting TVFC findings and limit the insights that can be gained from this promising new research area. This article brings together scientists with a variety of perspectives on resting TVFC to review the current literature in light of these issues. We introduce core concepts, define key terms, summarize controversies and open questions, and present a forward-looking perspective on how resting TVFC analyses can be rigorously and productively applied to investigate a wide range of questions in cognitive and systems neuroscience.

<h2>Summary</h2> Increasingly detailed data on the network topology of neural circuits create a need for theoretical principles that explain how these networks shape neural communication. Here we use a model of cascade spreading to reveal architectural features of human brain networks that facilitate spreading. Using an anatomical brain network derived from high-resolution diffusion spectrum imaging (DSI), we investigate scenarios where perturbations initiated at seed nodes result in global cascades that interact either cooperatively or competitively. We find that hub regions and a backbone of pathways facilitate early spreading, while the shortest path structure of the connectome enables cooperative effects, accelerating the spread of cascades. Finally, competing cascades become integrated by converging on polysensory associative areas. These findings show that the organizational principles of brain networks shape global communication and facilitate integrative function.

Neurobiological abnormalities associated with psychiatric disorders do not map well to existing diagnostic categories. High co-morbidity suggests dimensional circuit-level abnormalities that cross diagnoses. Here we seek to identify brain-based dimensions of psychopathology using sparse canonical correlation analysis in a sample of 663 youths. This analysis reveals correlated patterns of functional connectivity and psychiatric symptoms. We find that four dimensions of psychopathology - mood, psychosis, fear, and externalizing behavior - are associated (r = 0.68-0.71) with distinct patterns of connectivity. Loss of network segregation between the default mode network and executive networks emerges as a common feature across all dimensions. Connectivity linked to mood and psychosis becomes more prominent with development, and sex differences are present for connectivity related to mood and fear. Critically, findings largely replicate in an independent dataset (n = 336). These results delineate connectivity-guided dimensions of psychopathology that cross clinical diagnostic categories, which could serve as a foundation for developing network-based biomarkers in psychiatry.

The human brain is organized into large-scale functional modules that have been shown to evolve in childhood and adolescence. However, it remains unknown whether the underlying white matter architecture is similarly refined during development, potentially allowing for improvements in executive function. In a sample of 882 participants (ages 8-22) who underwent diffusion imaging as part of the Philadelphia Neurodevelopmental Cohort, we demonstrate that structural network modules become more segregated with age, with weaker connections between modules and stronger connections within modules. Evolving modular topology facilitates global network efficiency and is driven by age-related strengthening of hub edges present both within and between modules. Critically, both modular segregation and network efficiency are associated with enhanced executive performance and mediate the improvement of executive functioning with age. Together, results delineate a process of structural network maturation that supports executive function in youth.

The protracted development of structural and functional brain connectivity within distributed association networks coincides with improvements in higher-order cognitive processes such as executive function. However, it remains unclear how white-matter architecture develops during youth to directly support coordinated neural activity. Here, we characterize the development of structure–function coupling using diffusion-weighted imaging and n -back functional MRI data in a sample of 727 individuals (ages 8 to 23 y). We found that spatial variability in structure–function coupling aligned with cortical hierarchies of functional specialization and evolutionary expansion. Furthermore, hierarchy-dependent age effects on structure–function coupling localized to transmodal cortex in both cross-sectional data and a subset of participants with longitudinal data ( n = 294). Moreover, structure–function coupling in rostrolateral prefrontal cortex was associated with executive performance and partially mediated age-related improvements in executive function. Together, these findings delineate a critical dimension of adolescent brain development, whereby the coupling between structural and functional connectivity remodels to support functional specialization and cognition.

The dynamics of spontaneous fluctuations in neural activity are shaped by underlying patterns of anatomical connectivity. While numerous studies have demonstrated edge-wise correspondence between structural and functional connections, much less is known about how large-scale coherent functional network patterns emerge from the topology of structural networks. In the present study, we deploy a multivariate statistical technique, partial least squares, to investigate the association between spatially extended structural networks and functional networks. We find multiple statistically robust patterns, reflecting reliable combinations of structural and functional subnetworks that are optimally associated with one another. Importantly, these patterns generally do not show a one-to-one correspondence between structural and functional edges, but are instead distributed and heterogeneous, with many functional relationships arising from nonoverlapping sets of anatomical connections. We also find that structural connections between high-degree hubs are disproportionately represented, suggesting that these connections are particularly important in establishing coherent functional networks. Altogether, these results demonstrate that the network organization of the cerebral cortex supports the emergence of diverse functional network configurations that often diverge from the underlying anatomical substrate.

The human connectome represents a network map of the brain's wiring diagram and the pattern into which its connections are organized is thought to play an important role in cognitive function. The generative rules that shape the topology of the human connectome remain incompletely understood. Earlier work in model organisms has suggested that wiring rules based on geometric relationships (distance) can account for many but likely not all topological features. Here we systematically explore a family of generative models of the human connectome that yield synthetic networks designed according to different wiring rules combining geometric and a broad range of topological factors. We find that a combination of geometric constraints with a homophilic attachment mechanism can create synthetic networks that closely match many topological characteristics of individual human connectomes, including features that were not included in the optimization of the generative model itself. We use these models to investigate a lifespan dataset and show that, with age, the model parameters undergo progressive changes, suggesting a rebalancing of the generative factors underlying the connectome across the lifespan.

Encoding brain regions and their connections as a network of nodes and edges captures many of the possible paths along which information can be transmitted as humans process and perform complex behaviors. Because cognitive processes involve large, distributed networks of brain areas, principled examinations of multi-node routes within larger connection patterns can offer fundamental insights into the complexities of brain function. Here, we investigate both densely connected groups of nodes that could perform local computations as well as larger patterns of interactions that would allow for parallel processing. Finding such structures necessitates that we move from considering exclusively pairwise interactions to capturing higher order relations, concepts naturally expressed in the language of algebraic topology. These tools can be used to study mesoscale network structures that arise from the arrangement of densely connected substructures called cliques in otherwise sparsely connected brain networks. We detect cliques (all-to-all connected sets of brain regions) in the average structural connectomes of 8 healthy adults scanned in triplicate and discover the presence of more large cliques than expected in null networks constructed via wiring minimization, providing architecture through which brain network can perform rapid, local processing. We then locate topological cavities of different dimensions, around which information may flow in either diverging or converging patterns. These cavities exist consistently across subjects, differ from those observed in null model networks, and - importantly - link regions of early and late evolutionary origin in long loops, underscoring their unique role in controlling brain function. These results offer a first demonstration that techniques from algebraic topology offer a novel perspective on structural connectomics, highlighting loop-like paths as crucial features in the human brain's structural architecture.

We investigate the relationship of resting-state fMRI functional connectivity estimated over long periods of time with time-varying functional connectivity estimated over shorter time intervals. We show that using Pearson's correlation to estimate functional connectivity implies that the range of fluctuations of functional connections over short time-scales is subject to statistical constraints imposed by their connectivity strength over longer scales. We present a method for estimating time-varying functional connectivity that is designed to mitigate this issue and allows us to identify episodes where functional connections are unexpectedly strong or weak. We apply this method to data recorded from N = 80 participants, and show that the number of unexpectedly strong/weak connections fluctuates over time, and that these variations coincide with intermittent periods of high and low modularity in time-varying functional connectivity. We also find that during periods of relative quiescence regions associated with default mode network tend to join communities with attentional, control, and primary sensory systems. In contrast, during periods where many connections are unexpectedly strong/weak, default mode regions dissociate and form distinct modules. Finally, we go on to show that, while all functional connections can at times manifest stronger (more positively correlated) or weaker (more negatively correlated) than expected, a small number of connections, mostly within the visual and somatomotor networks, do so a disproportional number of times. Our statistical approach allows the detection of functional connections that fluctuate more or less than expected based on their long-time averages and may be of use in future studies characterizing the spatio-temporal patterns of time-varying functional connectivity.

Abstract To meet ongoing cognitive demands, the human brain must seamlessly transition from one brain state to another, in the process drawing on different cognitive systems. How does the brain’s network of anatomical connections help facilitate such transitions? Which features of this network contribute to making one transition easy and another transition difficult? Here, we address these questions using network control theory. We calculate the optimal input signals to drive the brain to and from states dominated by different cognitive systems. The input signals allow us to assess the contributions made by different brain regions. We show that such contributions, which we measure as energy, are correlated with regions’ weighted degrees. We also show that the network communicability, a measure of direct and indirect connectedness between brain regions, predicts the extent to which brain regions compensate when input to another region is suppressed. Finally, we identify optimal states in which the brain should start (and finish) in order to minimize transition energy. We show that the optimal target states display high activity in hub regions, implicating the brain’s rich club. Furthermore, when rich club organization is destroyed, the energy cost associated with state transitions increases significantly, demonstrating that it is the richness of brain regions that makes them ideal targets.

Connectomics has enhanced our understanding of neurocognitive development and decline by the integration of network sciences into studies across different stages of the human life span. However, these studies commonly occurred independently, missing the opportunity to test integrated models of the dynamical brain organization across the entire life span. In this review article, we survey empirical findings in life-span connectomics and propose a generative framework for computationally modeling the connectome over the human life span. This framework highlights initial findings that across the life span, the human connectome gradually shifts from an 'anatomically driven' organization to one that is more 'topological'. Finally, we consider recent advances that are promising to provide an integrative and systems perspective of human brain plasticity as well as underscore the pitfalls and challenges.

Network neuroscience has relied on a node-centric network model in which cells, populations and regions are linked to one another via anatomical or functional connections. This model cannot account for interactions of edges with one another. In this study, we developed an edge-centric network model that generates constructs ‘edge time series’ and ‘edge functional connectivity’ (eFC). Using network analysis, we show that, at rest, eFC is consistent across datasets and reproducible within the same individual over multiple scan sessions. We demonstrate that clustering eFC yields communities of edges that naturally divide the brain into overlapping clusters, with regions in sensorimotor and attentional networks exhibiting the greatest levels of overlap. We show that eFC is systematically modulated by variation in sensory input. In future work, the edge-centric approach could be useful for identifying novel biomarkers of disease, characterizing individual variation and mapping the architecture of highly resolved neural circuits. The authors present an edge-centric model of brain connectivity. Edge networks are stable across datasets, and their structure can be modulated by sensory input. When clustered, edge networks yield pervasively overlapping functional modules.

Brain areas' functional repertoires are shaped by their incoming and outgoing structural connections. In empirically measured networks, most connections are short, reflecting spatial and energetic constraints. Nonetheless, a small number of connections span long distances, consistent with the notion that the functionality of these connections must outweigh their cost. While the precise function of long-distance connections is unknown, the leading hypothesis is that they act to reduce the topological distance between brain areas and increase the efficiency of interareal communication. However, this hypothesis implies a nonspecificity of long-distance connections that we contend is unlikely. Instead, we propose that long-distance connections serve to diversify brain areas' inputs and outputs, thereby promoting complex dynamics. Through analysis of five weighted interareal network datasets, we show that long-distance connections play only minor roles in reducing average interareal topological distance. In contrast, areas' long-distance and short-range neighbors exhibit marked differences in their connectivity profiles, suggesting that long-distance connections enhance dissimilarity between areal inputs and outputs. Next, we show that-in isolation-areas' long-distance connectivity profiles exhibit nonrandom levels of similarity, suggesting that the communication pathways formed by long connections exhibit redundancies that may serve to promote robustness. Finally, we use a linearization of Wilson-Cowan dynamics to simulate the covariance structure of neural activity and show that in the absence of long-distance connections a common measure of functional diversity decreases. Collectively, our findings suggest that long-distance connections are necessary for supporting diverse and complex brain dynamics.

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.

The complexity of neural dynamics stems in part from the complexity of the underlying anatomy. Yet how white matter structure constrains how the brain transitions from one cognitive state to another remains unknown. Here we address this question by drawing on recent advances in network control theory to model the underlying mechanisms of brain state transitions as elicited by the collective control of region sets. We find that previously identified attention and executive control systems are poised to affect a broad array of state transitions that cannot easily be classified by traditional engineering-based notions of control. This theoretical versatility comes with a vulnerability to injury. In patients with mild traumatic brain injury, we observe a loss of specificity in putative control processes, suggesting greater susceptibility to neurophysiological noise. These results offer fundamental insights into the mechanisms driving brain state transitions in healthy cognition and their alteration following injury.

The development of advanced neuroimaging techniques and their deployment in large cohorts has enabled an assessment of functional and structural brain network architecture at an unprecedented level of detail. Across many temporal and spatial scales, network neuroscience has emerged as a central focus of intellectual efforts, seeking meaningful descriptions of brain networks and explanatory sets of network features that underlie circuit function in health and dysfunction in disease. However, the tools of network science commonly deployed provide insight into brain function at a fundamentally descriptive level, often failing to identify (patho-)physiological mechanisms that link system-level phenomena to the multiple hierarchies of brain function. Here we describe recently developed techniques stemming from advances in complex systems and network science that have the potential to overcome this limitation, thereby contributing mechanistic insights into neuroanatomy, functional dynamics, and pathology. Finally, we build on the Research Domain Criteria framework, highlighting the notion that mental illnesses can be conceptualized as dysfunctions of neural circuitry present across conventional diagnostic boundaries, to sketch how network-based methods can be combined with pharmacological, intermediate phenotype, genetic, and magnetic stimulation studies to probe mechanisms of psychopathology.

Brain networks are expected to be modular. However, existing techniques for estimating a network's modules make it difficult to assess the influence of organizational principles such as wiring cost reduction on the detected modules. Here we present a modification of an existing module detection algorithm that allowed us to focus on connections that are unexpected under a cost-reduction wiring rule and to identify modules from among these connections. We applied this technique to anatomical brain networks and showed that the modules we detected differ from those detected using the standard technique. We demonstrated that these novel modules are spatially distributed, exhibit unique functional fingerprints, and overlap considerably with rich clubs, giving rise to an alternative and complementary interpretation of the functional roles of specific brain regions. Finally, we demonstrated that, using the modified module detection approach, we can detect modules in a developmental dataset that track normative patterns of maturation. Collectively, these findings support the hypothesis that brain networks are composed of modules and provide additional insight into the function of those modules.

Abstract The human brain is organized into large-scale functional networks that can flexibly reconfigure their connectivity patterns, supporting both rapid adaptive control and long-term learning processes. However, it has remained unclear how short-term network dynamics support the rapid transformation of instructions into fluent behaviour. Comparing fMRI data of a learning sample ( N =70) with a control sample ( N =67), we find that increasingly efficient task processing during short-term practice is associated with a reorganization of large-scale network interactions. Practice-related efficiency gains are facilitated by enhanced coupling between the cingulo-opercular network and the dorsal attention network. Simultaneously, short-term task automatization is accompanied by decreasing activation of the fronto-parietal network, indicating a release of high-level cognitive control, and a segregation of the default mode network from task-related networks. These findings suggest that short-term task automatization is enabled by the brain’s ability to rapidly reconfigure its large-scale network organization involving complementary integration and segregation processes.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images. QSIPrep is a software platform for processing of most diffusion MRI datasets and ensures that adequate workflows are used.

The brain's functional diversity is reflected in the meso-scale architecture of its connectome, i.e. its division into clusters and communities of topologically-related brain regions. The dominant view, and one that is reinforced by current analysis techniques, is that communities are strictly assortative and segregated from one another, purportedly for the purpose of carrying out specialized information processing. Such a view, however, precludes the possibility of non-assortative communities that could engender a richer functional repertoire by allowing for a more complex set of inter-community interactions. Here, we use weighted stochastic blockmodels to uncover the meso-scale architecture of \emph{Drosophila}, mouse, rat, macaque, and human connectomes. We confirm that while many communities are assortative, others form core-periphery and disassortative structures, which in the human better recapitulate observed patterns of functional connectivity and in the mouse better recapitulate observed patterns of gene co-expression than other community detection techniques. We define a set of network measures for quantifying the diversity of community types in which brain regions participate. Finally, we show that diversity is peaked in control and subcortical systems in humans, and that individual differences in diversity within those systems predicts cognitive performance on Stroop and Navon tasks. In summary, our report paints a more diverse portrait of connectome meso-scale structure and demonstrates its relevance for cognitive performance.

Large-scale structural brain networks encode white matter connectivity patterns among distributed brain areas. These connection patterns are believed to support cognitive processes and, when compromised, can lead to neurocognitive deficits and maladaptive behavior. A powerful approach for studying the organizing principles of brain networks is to construct group-representative networks from multisubject cohorts. Doing so amplifies signal to noise ratios and provides a clearer picture of brain network organization. Here, we show that current approaches for generating sparse group-representative networks overestimate the proportion of short-range connections present in a network and, as a result, fail to match subject-level networks along a wide range of network statistics. We present an alternative approach that preserves the connection-length distribution of individual subjects. We have used this method in previous papers to generate group-representative networks, though to date its performance has not been appropriately benchmarked and compared against other methods. As a result of this simple modification, the networks generated using this approach successfully recapitulate subject-level properties, outperforming similar approaches by better preserving features that promote integrative brain function rather than segregative. The method developed here holds promise for future studies investigating basic organizational principles and features of large-scale structural brain networks.

Abstract Dynamical brain state transitions are critical for flexible working memory but the network mechanisms are incompletely understood. Here, we show that working memory performance entails brain-wide switching between activity states using a combination of functional magnetic resonance imaging in healthy controls and individuals with schizophrenia, pharmacological fMRI, genetic analyses and network control theory. The stability of states relates to dopamine D1 receptor gene expression while state transitions are influenced by D2 receptor expression and pharmacological modulation. Individuals with schizophrenia show altered network control properties, including a more diverse energy landscape and decreased stability of working memory representations. Our results demonstrate the relevance of dopamine signaling for the steering of whole-brain network dynamics during working memory and link these processes to schizophrenia pathophysiology.

Abstract A growing number of studies have used stylized network models of communication to predict brain function from structure. Most have focused on a small set of models applied globally. Here, we compare a large number of models at both global and regional levels. We find that globally most predictors perform poorly. At the regional level, performance improves but heterogeneously, both in terms of variance explained and the optimal model. Next, we expose synergies among predictors by using pairs to jointly predict FC. Finally, we assess age-related differences in global and regional coupling across the human lifespan. We find global decreases in the magnitude of structure-function coupling with age. We find that these decreases are driven by reduced coupling in sensorimotor regions, while higher-order cognitive systems preserve local coupling with age. Our results describe patterns of structure-function coupling across the cortex and how this may change with age.

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.

Graph theoretical analysis has played a key role in characterizing global features of the topology of complex networks, describing diverse systems such as protein interactions, food webs, social relations and brain connectivity. How system elements communicate with each other depends not only on the structure of the network, but also on the nature of the system's dynamics which are constrained by the amount of knowledge and resources available for communication processes. Complementing widely used measures that capture efficiency under the assumption that communication preferentially follows shortest paths across the network ("routing"), we define analytic measures directed at characterizing network communication when signals flow in a random walk process ("diffusion"). The two dimensions of routing and diffusion efficiency define a morphospace for complex networks, with different network topologies characterized by different combinations of efficiency measures and thus occupying different regions of this space. We explore the relation of network topologies and efficiency measures by examining canonical network models, by evolving networks using a multi-objective optimization strategy, and by investigating real-world network data sets. Within the efficiency morphospace, specific aspects of network topology that differentially favor efficient communication for routing and diffusion processes are identified. Charting regions of the morphospace that are occupied by canonical, evolved or real networks allows inferences about the limits of communication efficiency imposed by connectivity and dynamics, as well as the underlying selection pressures that have shaped network topology.

Abstract Advances in neuroimaging have made it possible to reconstruct functional networks from the activity patterns of brain regions distributed across the cerebral cortex. Recent work has shown that flexible reconfiguration of human brain networks over short timescales supports cognitive flexibility and learning. However, modulating network flexibility to enhance learning requires an understanding of an as-yet unknown relationship between flexibility and brain state. Here, we investigate the relationship between network flexibility and affect, leveraging an unprecedented longitudinal data set. We demonstrate that indices associated with positive mood and surprise are both associated with network flexibility – positive mood portends a more flexible brain while increased levels of surprise portend a less flexible brain. In both cases, these relationships are driven predominantly by a subset of brain regions comprising the somatomotor system. Our results simultaneously suggest a network-level mechanism underlying learning deficits in mood disorders as well as a potential target – altering an individual’s mood or task novelty – to improve learning.

The human connectome has been widely studied over the past decade. A principal finding is that it can be decomposed into communities of densely interconnected brain regions. This result, however, may be limited methodologically. Past studies have often used a flawed modularity measure in order to infer the connectome's community structure. Also, these studies relied on the intuition that community structure is best defined in terms of a network's static topology as opposed to a more dynamical definition. In this report we used the partition stability framework, which defines communities in terms of a Markov process (random walk), to infer the connectome's multi-scale community structure. Comparing the community structure to observed resting-state functional connectivity revealed communities across a broad range of dynamical scales that were closely related to functional connectivity. This result suggests a mapping between communities in structural networks, models of communication processes, and brain function. It further suggests that communication in the brain is not limited to a single characteristic scale, leading us to posit a heuristic for scale-selective communication in the cerebral cortex.

Structural white matter connections are thought to facilitate integration of neural information across functionally segregated systems. Recent studies have demonstrated that changes in the balance between segregation and integration in brain networks can be tracked by time-resolved functional connectivity derived from resting-state functional magnetic resonance imaging (rs-fMRI) data and that fluctuations between segregated and integrated network states are related to human behavior. However, how these network states relate to structural connectivity is largely unknown. To obtain a better understanding of structural substrates for these network states, we investigated how the relationship between structural connectivity, derived from diffusion tractography, and functional connectivity, as measured by rs-fMRI, changes with fluctuations between segregated and integrated states in the human brain. We found that the similarity of edge weights between structural and functional connectivity was greater in the integrated state, especially at edges connecting the default mode and the dorsal attention networks. We also demonstrated that the similarity of network partitions, evaluated between structural and functional connectivity, increased and the density of direct structural connections within modules in functional networks was elevated during the integrated state. These results suggest that, when functional connectivity exhibited an integrated network topology, structural connectivity and functional connectivity were more closely linked to each other and direct structural connections mediated a larger proportion of neural communication within functional modules. Our findings point out the possibility of significant contributions of structural connections to integrative neural processes underlying human behavior.

Recent advances in brain imaging techniques, measurement approaches, and storage capacities have provided an unprecedented supply of high temporal resolution neural data. These data present a remarkable opportunity to gain a mechanistic understanding not just of circuit structure, but also of circuit dynamics, and its role in cognition and disease. Such understanding necessitates a description of the raw observations, and a delineation of computational models and mathematical theories that accurately capture fundamental principles behind the observations. Here we review recent advances in a range of modeling approaches that embrace the temporally-evolving interconnected structure of the brain and summarize that structure in a dynamic graph. We describe recent efforts to model dynamic patterns of connectivity, dynamic patterns of activity, and patterns of activity atop connectivity. In the context of these models, we review important considerations in statistical testing, including parametric and non-parametric approaches. Finally, we offer thoughts on careful and accurate interpretation of dynamic graph architecture, and outline important future directions for method development.

Network neuroscience is the emerging discipline concerned with investigating the complex patterns of interconnections found in neural systems, and identifying principles with which to understand them. Within this discipline, one particularly powerful approach is network generative modelling, in which wiring rules are algorithmically implemented to produce synthetic network architectures with the same properties as observed in empirical network data. Successful models can highlight the principles by which a network is organized and potentially uncover the mechanisms by which it grows and develops. Here, we review the prospects and promise of generative models for network neuroscience. We begin with a primer on network generative models, with a discussion of compressibility and predictability, and utility in intuiting mechanisms, followed by a short history on their use in network science, broadly. We then discuss generative models in practice and application, paying particular attention to the critical need for cross-validation. Next, we review generative models of biological neural networks, both at the cellular and large-scale level, and across a variety of species including Caenorhabditis elegans, Drosophila, mouse, rat, cat, macaque and human. We offer a careful treatment of a few relevant distinctions, including differences between generative models and null models, sufficiency and redundancy, inferring and claiming mechanism, and functional and structural connectivity. We close with a discussion of future directions, outlining exciting frontiers both in empirical data collection efforts as well as in method and theory development that, together, further the utility of the generative network modelling approach for network neuroscience.

The network organization of the human brain varies across individuals, changes with development and aging, and differs in disease. Discovering the major dimensions along which this variability is displayed remains a central goal of both neuroscience and clinical medicine. Such efforts can be usefully framed within the context of the brain’s modular network organization, which can be assessed quantitatively using computational techniques and extended for the purposes of multi-scale analysis, dimensionality reduction, and biomarker generation. Although the concept of modularity and its utility in describing brain network organization is clear, principled methods for comparing multi-scale communities across individuals and time are surprisingly lacking. Here, we present a method that uses multi-layer networks to simultaneously discover the modular structure of many subjects at once. This method builds upon the well-known multi-layer modularity maximization technique, and provides a viable and principled tool for studying differences in network communities across individuals and within individuals across time. We test this method on two datasets and identify consistent patterns of inter-subject community variability, demonstrating that this variability – which would be undetectable using past approaches – is associated with measures of cognitive performance. In general, the multi-layer, multi-subject framework proposed here represents an advance over current approaches by straighforwardly mapping community assignments across subjects and holds promise for future investigations of inter-subject community variation in clinical populations or as a result of task constraints.

The hippocampal formation is a key structure for memory function in the brain. The functional anatomy of the brain suggests that the hippocampus may be a convergence zone, as it receives polysensory input from distributed association areas throughout the neocortex. However, recent quantitative graph-theoretic analyses of the static large-scale connectome have failed to demonstrate the centrality of the hippocampus; in the context of the whole brain, the hippocampus is not among the most connected or reachable nodes. Here we show that when communication dynamics are taken into account, the hippocampus is a key hub in the connectome. Using a novel computational model, we demonstrate that large-scale brain network topology is organized to funnel and concentrate information flow in the hippocampus, supporting the long-standing hypothesis that this region acts as a critical convergence zone. Our results indicate that the functional capacity of the hippocampus is shaped by its embedding in the large-scale connectome.

Intrinsically driven neural activity generated at rest exhibits complex spatiotemporal dynamics characterized by patterns of synchronization across distant brain regions. Mounting evidence suggests that these patterns exhibit fluctuations and nonstationarity at multiple time scales. Resting-state electroencephalographic (EEG) recordings were examined in 12 young adults for changes in synchronization patterns on a fast time scale in the range of tens to hundreds of milliseconds. Results revealed that EEG dynamics continuously underwent rapid transitions between intermittently stable states. Numerous approximate recurrences of states were observed within single recording epochs, across different epochs separated by longer times, and between participants. For broadband (4-30 Hz) data, a majority of states could be grouped into three families, suggesting the existence of a limited repertoire of core states that is continually revisited and shared across participants. Our results document the existence of fast synchronization dynamics iterating amongst a small set of core networks in the resting brain, complementing earlier findings of nonstationary dynamics in electromagnetic recordings and transient EEG microstates.

A diverse set of white matter connections supports seamless transitions between cognitive states. However, it remains unclear how these connections guide the temporal progression of large-scale brain activity patterns in different cognitive states. Here, we analyze the brain's trajectories across a set of single time point activity patterns from functional magnetic resonance imaging data acquired during the resting state and an n-back working memory task. We find that specific temporal sequences of brain activity are modulated by cognitive load, associated with age, and related to task performance. Using diffusion-weighted imaging acquired from the same subjects, we apply tools from network control theory to show that linear spread of activity along white matter connections constrains the probabilities of these sequences at rest, while stimulus-driven visual inputs explain the sequences observed during the n-back task. Overall, these results elucidate the structural underpinnings of cognitively and developmentally relevant spatiotemporal brain dynamics.

High-level cognitive constructs, such as creativity and intelligence, entail complex and multiple processes, including cognitive control processes. Recent neurocognitive research on these constructs highlight the importance of dynamic interaction across neural network systems and the role of cognitive control processes in guiding such a dynamic interaction. How can we quantitatively examine the extent and ways in which cognitive control contributes to creativity and intelligence? To address this question, we apply a computational network control theory (NCT) approach to structural brain imaging data acquired via diffusion tensor imaging in a large sample of participants, to examine how NCT relates to individual differences in distinct measures of creative ability and intelligence. Recent application of this theory at the neural level is built on a model of brain dynamics, which mathematically models patterns of inter-region activity propagated along the structure of an underlying network. The strength of this approach is its ability to characterize the potential role of each brain region in regulating whole-brain network function based on its anatomical fingerprint and a simplified model of node dynamics. We find that intelligence is related to the ability to “drive” the brain system into easy to reach neural states by the right inferior parietal lobe and lower integration abilities in the left retrosplenial cortex. We also find that creativity is related to the ability to “drive” the brain system into difficult to reach states by the right dorsolateral prefrontal cortex (inferior frontal junction) and higher integration abilities in sensorimotor areas. Furthermore, we found that different facets of creativity—fluency, flexibility, and originality—relate to generally similar but not identical network controllability processes. We relate our findings to general theories on intelligence and creativity.

The network architecture of the human brain contributes in shaping neural activity, influencing cognitive and behavioral processes. The availability of neuroimaging data across the lifespan allows us to monitor how this architecture reorganizes, influenced by processes like learning, adaptation, maturation, and senescence. Changing patterns in brain connectivity can be analyzed with the tools of network science, which can be used to reveal organizational principles such as modular network topology. The identification of network modules is fundamental, as they parse the brain into coherent sub-systems and allow for both functional integration and segregation among different brain areas. In this work we examined the brain’s modular organization by developing an ensemble-based multilayer network approach, allowing us to link changes of structural connectivity patterns to development and aging. We show that modular structure exhibits both linear and nonlinear age-related trends. In the early and late lifespan, communities are more modular, and we track the origins of this high modularity to two different substrates in brain connectivity, linked to the number and the weights of the intra-clusters edges. We also demonstrate that aging leads to a progressive and increasing reconfiguration of modules and a redistribution across hemispheres. Finally, we identify those brain regions that most contribute to network reconfiguration and those that remain more stable across the lifespan.

Executive function develops during adolescence, yet it remains unknown how structural brain networks mature to facilitate activation of the fronto-parietal system, which is critical for executive function. In a sample of 946 human youths (ages 8-23y) who completed diffusion imaging, we capitalized upon recent advances in linear dynamical network control theory to calculate the energetic cost necessary to activate the fronto-parietal system through the control of multiple brain regions given existing structural network topology. We found that the energy required to activate the fronto-parietal system declined with development, and the pattern of regional energetic cost predicts unseen individuals’ brain maturity. Finally, energetic requirements of the cingulate cortex were negatively correlated with executive performance, and partially mediated the development of executive performance with age. Our results reveal a mechanism by which structural networks develop during adolescence to reduce the theoretical energetic costs of transitions to activation states necessary for executive function.

Abstract Functional connectivity (FC) describes the statistical dependence between neuronal populations or brain regions in resting-state fMRI studies and is commonly estimated as the Pearson correlation of time courses. Clustering or community detection reveals densely coupled sets of regions constituting resting-state networks or functional systems. These systems manifest most clearly when FC is sampled over longer epochs but appear to fluctuate on shorter timescales. Here, we propose a new approach to reveal temporal fluctuations in neuronal time series. Unwrapping FC signal correlations yields pairwise co-fluctuation time series, one for each node pair or edge, and allows tracking of fine-scale dynamics across the network. Co-fluctuations partition the network, at each time step, into exactly two communities. Sampled over time, the overlay of these bipartitions, a binary decomposition of the original time series, very closely approximates functional connectivity. Bipartitions exhibit characteristic spatiotemporal patterns that are reproducible across participants and imaging runs, capture individual differences, and disclose fine-scale temporal expression of functional systems. Our findings document that functional systems appear transiently and intermittently, and that FC results from the overlay of many variable instances of system expression. Potential applications of this decomposition of functional connectivity into a set of binary patterns are discussed.

The relationship between structural and functional connectivity in the brain is a key question in systems neuroscience. Modern accounts assume a single global structure-function relationship that persists over time. Here we study structure-function coupling from a dynamic perspective, and show that it is regionally heterogeneous. We use a temporal unwrapping procedure to identify moment-to-moment co-fluctuations in neural activity, and reconstruct time-resolved structure-function coupling patterns. We find that patterns of dynamic structure-function coupling are region-specific. We observe stable coupling in unimodal and transmodal cortex, and dynamic coupling in intermediate regions, particularly in insular cortex (salience network) and frontal eye fields (dorsal attention network). Finally, we show that the variability of a region's structure-function coupling is related to the distribution of its connection lengths. Collectively, our findings provide a way to study structure-function relationships from a dynamic perspective.

Abstract Many studies have shown that the human endocrine system modulates brain function, reporting associations between fluctuations in hormone concentrations and brain connectivity. However, how hormonal fluctuations impact fast changes in brain network organization over short timescales remains unknown. Here, we leverage a recently proposed framework for modeling co-fluctuations between the activity of pairs of brain regions at a framewise timescale. In previous studies we showed that time points corresponding to high-amplitude co-fluctuations disproportionately contributed to the time-averaged functional connectivity pattern and that these co-fluctuation patterns could be clustered into a low-dimensional set of recurring “states.” Here, we assessed the relationship between these network states and quotidian variation in hormone concentrations. Specifically, we were interested in whether the frequency with which network states occurred was related to hormone concentration. We addressed this question using a dense-sampling dataset (N = 1 brain). In this dataset, a single individual was sampled over the course of two endocrine states: a natural menstrual cycle and while the subject underwent selective progesterone suppression via oral hormonal contraceptives. During each cycle, the subject underwent 30 daily resting-state fMRI scans and blood draws. Our analysis of the imaging data revealed two repeating network states. We found that the frequency with which state 1 occurred in scan sessions was significantly correlated with follicle-stimulating and luteinizing hormone concentrations. We also constructed representative networks for each scan session using only “event frames”—those time points when an event was determined to have occurred. We found that the weights of specific subsets of functional connections were robustly correlated with fluctuations in the concentration of not only luteinizing and follicle-stimulating hormones, but also progesterone and estradiol.

Generative models of the human connectome enable in silico generation of brain networks based on probabilistic wiring rules. These wiring rules are governed by a small number of parameters that are typically fitted to individual connectomes and quantify the extent to which geometry and topology shape the generative process. A significant shortcoming of generative modeling in large cohort studies is that parameter estimation is computationally burdensome, and the accuracy and reliability of current estimation methods remain untested. Here, we propose a fast, reliable, and accurate parameter estimation method for connectome generative models that is scalable to large sample sizes. Our method achieves improved estimation accuracy and reliability and reduces computational cost by orders of magnitude, compared to established methods. We demonstrate an inherent tradeoff between accuracy, reliability, and computational expense in parameter estimation and provide recommendations for leveraging this tradeoff. To enable power analyses in future studies, we empirically approximate the minimum sample size required to detect between-group differences in generative model parameters. While we focus on the classic two-parameter generative model based on connection length and the topological matching index, our method can be generalized to other growth-based generative models. Our work provides a statistical and practical guide to parameter estimation for connectome generative models.

Converging evidence from activation, connectivity, and stimulation studies suggests that auditory brain networks are lateralized. Here we show that these findings can be at least partly explained by the asymmetric network embedding of the primary auditory cortices. Using diffusion-weighted imaging in 3 independent datasets, we investigate the propensity for left and right auditory cortex to communicate with other brain areas by quantifying the centrality of the auditory network across a spectrum of communication mechanisms, from shortest path communication to diffusive spreading. Across all datasets, we find that the right auditory cortex is better integrated in the connectome, facilitating more efficient communication with other areas, with much of the asymmetry driven by differences in communication pathways to the opposite hemisphere. Critically, the primacy of the right auditory cortex emerges only when communication is conceptualized as a diffusive process, taking advantage of more than just the topologically shortest paths in the network. Altogether, these results highlight how the network configuration and embedding of a particular region may contribute to its functional lateralization.

Modularity is an important topological attribute for functional brain networks. Recent human fMRI studies have reported that modularity of functional networks varies not only across individuals being related to demographics and cognitive performance, but also within individuals co-occurring with fluctuations in network properties of functional connectivity, estimated over short time intervals. However, characteristics of these time-resolved functional networks during periods of high and low modularity have remained largely unexplored. In this study we investigate basic spatiotemporal properties of time-resolved networks in the high and low modularity periods during rest, with a particular focus on their spatial connectivity patterns, temporal homogeneity and test-retest reliability. We show that spatial connectivity patterns of time-resolved networks in the high and low modularity periods are represented by increased and decreased dissociation of the default mode network module from task-positive network modules, respectively. We also find that the instances of time-resolved functional connectivity sampled from within the high (respectively, low) modularity period are relatively homogeneous (respectively, heterogeneous) over time, indicating that during the low modularity period the default mode network interacts with other networks in a variable manner. We confirmed that the occurrence of the high and low modularity periods varies across individuals with moderate inter-session test-retest reliability and that it is correlated with previously-reported individual differences in the modularity of functional connectivity estimated over longer timescales. Our findings illustrate how time-resolved functional networks are spatiotemporally organized during periods of high and low modularity, allowing one to trace individual differences in long-timescale modularity to the variable occurrence of network configurations at shorter timescales.

Brain networks are flexible and reconfigure over time to support ongoing cognitive processes. However, tracking statistically meaningful reconfigurations across time has proven difficult. This has to do largely with issues related to sampling variability, making instantaneous estimation of network organization difficult, along with increased reliance on task-free (cognitively unconstrained) experimental paradigms, limiting the ability to interpret the origin of changes in network structure over time. Here, we address these challenges using time-varying network analysis in conjunction with a naturalistic viewing paradigm. Specifically, we developed a measure of inter-subject network similarity and used this measure as a coincidence filter to identify synchronous fluctuations in network organization across individuals. Applied to movie-watching data, we found that periods of high inter-subject similarity coincided with reductions in network modularity and increased connectivity between cognitive systems. In contrast, low inter-subject similarity was associated with increased system segregation and more rest-like architectures. We then used a data-driven approach to uncover clusters of functional connections that follow similar trajectories over time and are more strongly correlated during movie-watching than at rest. Finally, we show that synchronous fluctuations in network architecture over time can be linked to a subset of features in the movie. Our findings link dynamic fluctuations in network integration and segregation to patterns of inter-subject similarity, and suggest that moment-to-moment fluctuations in functional connectivity reflect shared cognitive processing across individuals.

The wiring of vertebrate and invertebrate brains provides the anatomical skeleton for cognition and behavior. Connections among brain regions are characterized by heterogeneous strength that is parsimoniously described by the wiring cost and homophily principles. Moreover, brains exhibit a characteristic global network topology, including modules and hubs. However, the mechanisms resulting in the observed interregional wiring principles and network topology of brains are unknown. Here, with the aid of computational modeling, we demonstrate that a mechanism based on heterochronous and spatially ordered neurodevelopmental gradients, without the involvement of activity-dependent plasticity or axonal guidance cues, can reconstruct a large part of the wiring principles (on average, 83%) and global network topology (on average, 80%) of diverse adult brain connectomes, including fly and human connectomes. In sum, space and time are key components of a parsimonious, plausible neurodevelopmental mechanism of brain wiring with a potential universal scope, encompassing vertebrate and invertebrate brains.

Significance Brain regions engage in complex patterns of activation over time. Relating these patterns to neural processing is a central challenge in cognitive neuroscience. Recent work has identified brief intermittent bursts of brain-wide signal cofluctuations, called events, and shown that events drive functional connectivity. The origins of events are unclear. Here, we address this gap in knowledge by implementing computational models of neural oscillators coupled by anatomical connections derived from maps of the human cerebral cortex. Analysis of the emerging large-scale brain dynamics reveals brief episodes with high system-wide signal amplitudes. Simulated events closely correspond to those seen recently in empirical recordings. Notably, simulated events are significantly aligned with underlying structural modules, thus suggesting an important role of modular network organization.

Abstract Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.

Progress in scientific disciplines is accompanied by standardization of terminology. Network neuroscience, at the level of macroscale organization of the brain, is beginning to confront the challenges associated with developing a taxonomy of its fundamental explanatory constructs. The Workgroup for HArmonized Taxonomy of NETworks (WHATNET) was formed in 2020 as an Organization for Human Brain Mapping (OHBM)-endorsed best practices committee to provide recommendations on points of consensus, identify open questions, and highlight areas of ongoing debate in the service of moving the field toward standardized reporting of network neuroscience results. The committee conducted a survey to catalog current practices in large-scale brain network nomenclature. A few well-known network names (e.g., default mode network) dominated responses to the survey, and a number of illuminating points of disagreement emerged. We summarize survey results and provide initial considerations and recommendations from the workgroup. This perspective piece includes a selective review of challenges to this enterprise, including (1) network scale, resolution, and hierarchies; (2) interindividual variability of networks; (3) dynamics and nonstationarity of networks; (4) consideration of network affiliations of subcortical structures; and (5) consideration of multimodal information. We close with minimal reporting guidelines for the cognitive and network neuroscience communities to adopt.

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.

Graph theory has provided a key mathematical framework to analyse the architecture of human brain networks. This architecture embodies an inherently complex relationship between connection topology, the spatial arrangement of network elements, and the resulting network cost and functional performance. An exploration of these interacting factors and driving forces may reveal salient network features that are critically important for shaping and constraining the brain's topological organization and its evolvability. Several studies have pointed to an economic balance between network cost and network efficiency with networks organized in an 'economical' small-world favouring high communication efficiency at a low wiring cost. In this study, we define and explore a network morphospace in order to characterize different aspects of communication efficiency in human brain networks. Using a multi-objective evolutionary approach that approximates a Pareto-optimal set within the morphospace, we investigate the capacity of anatomical brain networks to evolve towards topologies that exhibit optimal information processing features while preserving network cost. This approach allows us to investigate network topologies that emerge under specific selection pressures, thus providing some insight into the selectional forces that may have shaped the network architecture of existing human brains.

Network science has begun to reveal the fundamental principles by which large-scale brain networks are organized, including geometric constraints, a balance between segregative and integrative features, and functionally flexible brain areas. However, it remains unknown whether whole-brain networks imaged at the cellular level are organized according to similar principles. Here, we analyze whole-brain functional networks reconstructed from calcium imaging data recorded in larval zebrafish. Our analyses reveal that functional connections are distance-dependent and that networks exhibit hierarchical modular structure and hubs that span module boundaries. We go on to show that spontaneous network structure places constraints on stimulus-evoked reconfigurations of connections and that networks are highly consistent across individuals. Our analyses reveal basic organizing principles of whole-brain functional brain networks at the mesoscale. Our overarching methodological framework provides a blueprint for studying correlated activity at the cellular level using a low-dimensional network representation. Our work forms a conceptual bridge between macro- and mesoscale network neuroscience and opens myriad paths for future studies to investigate network structure of nervous systems at the cellular level.

The human brain is composed of functionally specialized systems that support cognition. Recently, we proposed an edge-centric model for detecting overlapping communities. It remains unclear how these communities and brain systems are related. Here, we address this question using data from the Midnight Scan Club and show that all brain systems are linked via at least two edge communities. We then examine the diversity of edge communities within each system, finding that heteromodal systems are more diverse than sensory systems. Next, we cluster the entire cortex to reveal it according to the regions' edge-community profiles. We find that regions in heteromodal systems are more likely to form their own clusters. Finally, we show that edge communities are personalized. Our work reveals the pervasive overlap of edge communities across the cortex and their relationship with brain systems. Our work provides pathways for future research using edge-centric brain networks.

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.

Alterations in the structural connectome of schizophrenia patients have been widely characterized, but the mechanisms remain largely unknown. Generative network models have recently been introduced as a tool to test the biological underpinnings of altered brain network formation. We evaluated different generative network models in healthy controls (n=152), schizophrenia patients (n=66), and their unaffected first-degree relatives (n=32), and we identified spatial and topological factors contributing to network formation. We further investigated how these factors relate to cognition and to polygenic risk for schizophrenia. Our data show that among the four tested classes of generative network models, structural brain networks were optimally accounted for by a two-factor model combining spatial constraints and topological neighborhood structure. The same wiring model explained brain network formation across study groups. However, relatives and schizophrenia patients exhibited significantly lower spatial constraints and lower topological facilitation compared to healthy controls. Further exploratory analyses point to potential associations of the model parameter reflecting spatial constraints with the polygenic risk for schizophrenia and cognitive performance. Our results identify spatial constraints and local topological structure as two interrelated mechanisms contributing to regular brain network formation as well as altered connectomes in schizophrenia and healthy individuals at familial risk for schizophrenia. On an exploratory level, our data further point to the potential relevance of spatial constraints for the genetic risk for schizophrenia and general cognitive functioning, thereby encouraging future studies in following up on these observations to gain further insights into the biological basis and behavioral relevance of model parameters.

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.

Circadian rhythms (lasting approximately 24 h) control and entrain various physiological processes, ranging from neural activity and hormone secretion to sleep cycles and eating habits. Several studies have shown that time of day (TOD) is associated with human cognition and brain functions. In this study, utilizing a chronotype-based paradigm, we applied a graph theory approach on resting-state functional MRI (rs-fMRI) data to compare whole-brain functional network topology between morning and evening sessions and between morning-type (MT) and evening-type (ET) participants. Sixty-two individuals (31 MT and 31 ET) underwent two fMRI sessions, approximately 1 hour (morning) and 10 h (evening) after their wake-up time, according to their declared habitual sleep-wake pattern on a regular working day. In the global analysis, the findings revealed the effect of TOD on functional connectivity (FC) patterns, including increased small-worldness, assortativity, and synchronization across the day. However, we identified no significant differences based on chronotype categories. The study of the modular structure of the brain at mesoscale showed that functional networks tended to be more integrated with one another in the evening session than in the morning session. Local/regional changes were affected by both factors (i.e., TOD and chronotype), mostly in areas associated with somatomotor, attention, frontoparietal, and default networks. Furthermore, connectivity and hub analyses revealed that the somatomotor, ventral attention, and visual networks covered the most highly connected areas in the morning and evening sessions: the latter two were more active in the morning sessions, and the first was identified as being more active in the evening. Finally, we performed a correlation analysis to determine whether global and nodal measures were associated with subjective assessments across participants. Collectively, these findings contribute to an increased understanding of diurnal fluctuations in resting brain activity and highlight the role of TOD in future studies on brain function and the design of fMRI experiments.

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.

Most neuroimaging studies of post-stroke recovery rely on analyses derived from standard node-centric functional connectivity to map the distributed effects in stroke patients. Here, given the importance of nonlocal and diffuse damage, we use an edge-centric approach to functional connectivity in order to provide an alternative description of the effects of this disorder. These techniques allow for the rendering of metrics such as normalized entropy, which describes the diversity of edge communities at each node. Moreover, the approach enables the identification of high amplitude co-fluctuations in fMRI time series. We found that normalized entropy is associated with stroke lesion severity and continually increases across the time of patients' recovery. Furthermore, high amplitude co-fluctuations not only relate to the lesion severity but are also associated with patients' level of recovery. The current study is the first edge-centric application for a clinical population in a longitudinal dataset and demonstrates how a different perspective for functional data analysis can further characterize topographic modulations of brain dynamics.

Network models of communication, e.g. shortest paths, diffusion, navigation, have become useful tools for studying structure-function relationships in the brain. These models generate estimates of communication efficiency between all pairs of brain regions, which can then be linked to the correlation structure of recorded activity, i.e. functional connectivity (FC). At present, however, communication models have a number of limitations, including difficulty adjudicating between models and the absence of a generic framework for modeling multiple interacting communication policies at the regional level. Here, we present a framework that allows us to incorporate multiple region-specific policies and fit them to empirical estimates of FC. Briefly, we show that many communication policies, including shortest paths and greedy navigation, can be modeled as biased random walks, enabling these policies to be incorporated into the same multi-policy communication model alongside unbiased processes, e.g. diffusion. We show that these multi-policy models outperform existing communication measures while yielding neurobiologically interpretable regional preferences. Further, we show that these models explain the majority of variance in time-varying patterns of FC. Collectively, our framework represents an advance in network-based communication models and establishes a strong link between these patterns and FC. Our findings open up many new avenues for future inquiries and present a flexible framework for modeling anatomically-constrained communication.

The interaction between brain regions changes over time, which can be characterized using time-varying functional connectivity (tvFC). The common approach to estimate tvFC uses sliding windows and offers limited temporal resolution. An alternative method is to use the recently proposed edge-centric approach, which enables the tracking of moment-to-moment changes in co-fluctuation patterns between pairs of brain regions. Here, we first examined the dynamic features of edge time series and compared them to those in the sliding window tvFC (sw-tvFC). Then, we used edge time series to compare subjects with autism spectrum disorder (ASD) and healthy controls (CN). Our results indicate that relative to sw-tvFC, edge time series captured rapid and bursty network-level fluctuations that synchronize across subjects during movie-watching. The results from the second part of the study suggested that the magnitude of peak amplitude in the collective co-fluctuations of brain regions (estimated as root sum square (RSS) of edge time series) is similar in CN and ASD. However, the trough-to-trough duration in RSS signal is greater in ASD, compared to CN. Furthermore, an edge-wise comparison of high-amplitude co-fluctuations showed that the within-network edges exhibited greater magnitude fluctuations in CN. Our findings suggest that high-amplitude co-fluctuations captured by edge time series provide details about the disruption of functional brain dynamics that could potentially be used in developing new biomarkers of mental disorders.

A rapidly emerging application of network neuroscience in neuroimaging studies has provided useful tools to understand individual differences in intrinsic brain function by mapping spontaneous brain activity, namely intrinsic functional network neuroscience (ifNN). However, the variability of methodologies applied across the ifNN studies-with respect to node definition, edge construction, and graph measurements-makes it difficult to directly compare findings and also challenging for end users to select the optimal strategies for mapping individual differences in brain networks. Here, we aim to provide a benchmark for best ifNN practices by systematically comparing the measurement reliability of individual differences under different ifNN analytical strategies using the test-retest design of the Human Connectome Project. The results uncovered four essential principles to guide ifNN studies: (1) use a whole brain parcellation to define network nodes, including subcortical and cerebellar regions; (2) construct functional networks using spontaneous brain activity in multiple slow bands; and (3) optimize topological economy of networks at individual level; and (4) characterize information flow with specific metrics of integration and segregation. We built an interactive online resource of reliability assessments for future ifNN (https://ibraindata.com/research/ifNN).

Abstract Recent studies have shown that functional connectivity can be decomposed into its exact frame-wise contributions, revealing short-lived, infrequent, and high-amplitude time points referred to as “events.” Events contribute disproportionately to the time-averaged connectivity pattern, improve identifiability and brain-behavior associations, and differences in their expression have been linked to endogenous hormonal fluctuations and autism. Here, we explore the characteristics of events while subjects watch movies. Using two independently-acquired imaging datasets in which participants passively watched movies, we find that events synchronize across individuals and based on the level of synchronization, can be categorized into three distinct classes: those that synchronize at the boundaries between movies, those that synchronize during movies, and those that do not synchronize at all. We find that boundary events, compared to the other categories, exhibit greater amplitude, distinct co-fluctuation patterns, and temporal propagation. We show that underlying boundary events1 is a specific mode of co-fluctuation involving the activation of control and salience systems alongside the deactivation of visual systems. Events that synchronize during the movie, on the other hand, display a pattern of co-fluctuation that is time-locked to the movie stimulus. Finally, we found that subjects’ time-varying brain networks are most similar to one another during these synchronous events.

Patterns of neural activity underlie human cognition. Transitions between these patterns are orchestrated by the brain's network architecture. What are the mechanisms linking network structure to cognitively relevant activation patterns? Here we implement principles of network control to investigate how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic engine. We also systematically incorporate neurotransmitter receptor density maps (18 receptors and transporters) and disease-related cortical abnormality maps (11 neurodegenerative, psychiatric and neurodevelopmental diseases; N = 17 000 patients, N = 22 000 controls). Integrating large-scale multimodal neuroimaging data from functional MRI, diffusion tractography, cortical morphometry, and positron emission tomography, we simulate how anatomically-guided transitions between cognitive states can be reshaped by pharmacological or pathological perturbation. Our results provide a comprehensive look-up table charting how brain network organisation and chemoarchitecture interact to manifest different cognitive topographies. This computational framework establishes a principled foundation for systematically identifying novel ways to promote selective transitions between desired cognitive topographies.

Abstract Pang et al. (2023) present novel analyses demonstrating that brain dynamics can be understood as resulting from the excitation of geometric modes, derived from the shape of the brain. Notably, they demonstrate that linear combinations of geometric modes can reconstruct patterns of fMRI data more accurately, and with fewer dimensions, than comparable connectivity-derived modes. Equipped with these results, and underpinned by neural field theory, the authors contend that the geometry of the cortical surface provides a more parsimonious explanation of brain activity than structural brain connectivity. This claim runs counter to prevailing theories of information flow in the brain, which emphasize the role of long-distance axonal projections and fasciculated white matter in relaying signals between cortical regions (Honey et al. 2009; Deco et al. 2011; Seguin et al., 2023). While we acknowledge that cortical geometry plays an important role in shaping human brain function, we feel that the presented work falls short of establishing that the brain’s geometry is “a more fundamental constraint on dynamics than complex interregional connectivity” (Pang et al. 2023). Here, we provide 1) a brief critique of the paper’s framing and 2) evidence showing that their methodology lacks specificity to the brain’s orientation and shape. Ultimately, we recognize that the geometric mode approach is a powerful representational framework for brain dynamics analysis, but we also believe that there are key caveats to consider alongside the claims made in the manuscript.

Many real-world systems are composed of many individual components that interact with one another in a complex pattern to produce diverse behaviors. Understanding how to intervene in these systems to guide behaviors is critically important to facilitate new discoveries and therapies in systems biology and neuroscience. A promising approach to optimizing interventions in complex systems is network control theory, an emerging conceptual framework and associated mathematics to understand how targeted input to nodes in a network system can predictably alter system dynamics. While network control theory is currently being applied to real-world data, the practical performance of these measures on simple networks with pre-specified structure is not well understood. In this study, we benchmark measures of network controllability on canonical graph models, providing an intuition for how control strategy, graph topology, and edge weight distribution mutually depend on one another. Our numerical studies motivate future analytical efforts to gain a mechanistic understanding of the relationship between graph topology and control, as well as efforts to design networks with specific control profiles.

The intricate connectivity patterns of neural circuits support a wide repertoire of communication processes and functional interactions. Here we systematically investigate how neural signaling is constrained by anatomical connectivity in the mesoscale Drosophila (fruit fly) brain network. We use a spreading model that describes how local perturbations, such as external stimuli, trigger global signaling cascades that spread through the network. Through a series of simple biological scenarios we demonstrate that anatomical embedding potentiates sensory-motor integration. We find that signal spreading is faster from nodes associated with sensory transduction (sensors) to nodes associated with motor output (effectors). Signal propagation was accelerated if sensor nodes were activated simultaneously, suggesting a topologically mediated synergy among sensors. In addition, the organization of the network increases the likelihood of convergence of multiple cascades towards effector nodes, thereby facilitating integration prior to motor output. Moreover, effector nodes tend to coactivate more frequently than other pairs of nodes, suggesting an anatomically enhanced coordination of motor output. Altogether, our results show that the organization of the mesoscale Drosophila connectome imparts privileged, behaviorally relevant communication patterns among sensors and effectors, shaping their capacity to collectively integrate information.

This scientific commentary refers to ‘Neural, electrophysiological and anatomical basis of brain-network variability and its characteristic changes in mental disorders’ by Zhang et al. (doi:10.1093/aww143). The vastness of the brain’s dynamic repertoire is one of the remarkable features of brain function, making it possible to adapt rapidly and efficiently to external task demands, implement novel behaviours, and switch from one task to another. Variability in the neural dynamics is, nonetheless, constrained and displays heterogeneous topography—specific regions appear more or less variable over time, both in terms of their activity time courses (Garrett et al. , 2011) and also in terms of their interactions with other brain regions (their functional connectivity). In this issue of Brain , Zhang and co-workers present a novel method for characterizing the temporal variability of a region’s functional connectivity profile [estimated from blood oxygen level-dependent (BOLD) functional MRI], relating this variability to the region’s electrophysiology (measured with electroencephalography; EEG) and to its macroscale structural connectivity (white-matter pathways), and further demonstrating its potential utility as a neural marker for mental disorders (Zhang et al. , 2016). The past decade has witnessed a burgeoning interest in the functional network architecture of the human brain. Most of these earlier studies have adopted a ‘static’ point of view, wherein functional connections between regions are characterized over long time scales, obscuring faster dynamics. More recently, however, it has become apparent that over the course of seconds to minutes, the human brain displays network-wide reconfigurations both at rest (Zalesky et al. , 2014) and during task performance (Braun et al. , 2015). These findings …

Coordinated brain activity reflects underlying cognitive processes and can be modeled as a network of inter-regional functional connections. The most costly connections in the network are long-distance correlations that, in the absence of underlying structural connections, are maintained by sustained energetic inputs. Here, we present a spatial modeling approach that amplifies contributions made by long-distance functional connections to whole-brain network architecture, while simultaneously suppressing contributions made by short-range connections. We use this method to characterize the long-distance architecture of functional networks and to identify aspects of community and hub structure that are driven by long-distance correlations and that, we argue, are of greater functional significance. We find that based only on patterns of long-distance connectivity, primary sensory cortices occupy increasingly central positions and appear more "hub-like". Additionally, we show that the community structure of long-distance connections spans multiple topological levels and differs from the community structure detected in networks that include both short-range and long-distance connections. In summary, these findings highlight the complex relationship between the brain's physical layout and its functional architecture. The results presented here inform future analyses of community structure and network hubs in health, across development, and in the case of neuropsychiatric disorders.

Recent studies have provided insight into inter-individual differences in creative thinking, focusing on characterizations of distributed large-scale brain networks both at the local level of regions and their pairwise interactions and at the global level of the brain as a whole. However, it remains unclear how creative thinking relates to mesoscale network features, e.g. community and hub organization. We applied a data-driven approach to examine community and hub structure in resting-state functional imaging data from a large sample of participants, and how they relate to individual differences in creative thinking. First, we computed for every participant the co-assignment probability of brain regions to the same community. We found that greater capacity for creative thinking was related to increased and decreased co-assignment of medial-temporal and subcortical regions to the same community, respectively, suggesting that creative capacity may be reflected in inter-individual differences in the meso-scale organization of brain networks. We then used participant-specific communities to identify network hubs—nodes whose connections form bridges across the boundaries of different communities—quantified based on their participation coefficients. We found that increased hubness of DMN and medial-temporal regions were positively and negatively related with creative ability, respectively. These findings suggest that creative capacity may be reflected in inter-individual differences in community interactions of DMN and medial-temporal structures. Collectively, these results demonstrate the fruitfulness of investigating mesoscale brain network features in relation to creative thinking.

Abstract 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.

Network neuroscience has emphasized the connectional properties of neural elements – cells, populations, and regions. This has come at the expense of the anatomical and functional connections that link these elements to one another. A new perspective – namely one that emphasizes 'edges' – may prove fruitful in addressing outstanding questions in network neuroscience. We highlight one recently proposed 'edge-centric' method and review its current applications, merits, and limitations. We also seek to establish conceptual and mathematical links between this method and previously proposed approaches in the network science and neuroimaging literature. We conclude by presenting several avenues for future work to extend and refine existing edge-centric analysis.

Brain networks exhibit community structure that reconfigures during cognitively demanding tasks. Extant work has emphasized a single class of communities: those that are assortative, or internally dense and externally sparse. Other classes that may play key functional roles in brain function have largely been ignored, leading to an impoverished view in the best case and a mischaracterization in the worst case. Here, we leverage weighted stochastic blockmodeling, a community detection method capable of detecting diverse classes of communities, to study the community structure of functional brain networks while subjects either rest or perform cognitively demanding tasks. We find evidence that the resting brain is largely assortative, although higher order association areas exhibit non-assortative organization, forming cores and peripheries. Surprisingly, this assortative structure breaks down during tasks and is supplanted by core, periphery, and disassortative communities. Using measures derived from the community structure, we show that it is possible to classify an individual’s task state with an accuracy that is well above average. Finally, we show that inter-individual differences in the composition of assortative and non-assortative communities is correlated with subject performance on in-scanner cognitive tasks. These findings offer a new perspective on the community organization of functional brain networks and its relation to cognition.

<h2>Summary</h2> Microbial organisms of the human gut microbiome do not exist in isolation but form complex and diverse interactions to maintain health and reduce risk of disease development. The organization of the gut microbiome is assumed to be a singular assortative network, where interactions between operational taxonomic units (OTUs) can readily be clustered into segregated and distinct communities. Here, we leverage recent methodological advances in network modeling to assess whether communities in the human microbiome exhibit a single network structure or whether co-existing mesoscale network architectures are present. We found evidence for core-periphery structures in the microbiome, supported by strong, assortative community interactions. This complex architecture, coupled with previously reported functional roles of OTUs, provides a nuanced understanding of how the microbiome simultaneously promotes high microbial diversity and maintains functional redundancy.

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 co-fluctuations 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 co-fluctuations of network organization with fluctuations in the BOLD time series. We show that, surprisingly, only a small fraction of frames exhibiting the strongest co-fluctuation 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 co-fluctuation 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 new directions for future research.

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.

Connectomics and network neuroscience offer quantitative scientific frameworks for modeling and analyzing networks of structurally and functionally interacting neurons, neuronal populations, and macroscopic brain areas. This shift in perspective and emphasis on distributed brain function has provided fundamental insight into the role played by the brain’s network architecture in cognition, disease, development, and aging. In this chapter, we review the core concepts of human connectomics at the macroscale. From the construction of networks using functional and diffusion MRI data, to their subsequent analysis using methods from network neuroscience, this review highlights key findings, commonly used methodologies, and discusses several emerging frontiers in connectomics.

Progress in scientific disciplines is accompanied by standardization of terminology. Network neuroscience, at the level of macro-scale organization of the brain, is beginning to confront the challenges associated with developing a taxonomy of its fundamental explanatory constructs. The Workgroup for HArmonized Taxonomy of NETworks (WHATNET) was formed in 2020 as an Organization for Human Brain Mapping (OHBM)-endorsed best practices committee to provide recommendations on points of consensus, identify open questions, and highlight areas of ongoing debate in the service of moving the field towards standardized reporting of network neuroscience results. The committee conducted a survey to catalog current practices in large-scale brain network nomenclature. A few well-known network names (e.g., default mode network) dominated responses to the survey, and a number of illuminating points of disagreement emerged. We summarize survey results and provide initial considerations and recommendations from the workgroup. This perspective piece includes a selective review of challenges to this enterprise, including 1) network scale, resolution, and hierarchies; 2) inter-individual variability of networks; 3) dynamics and non-stationarity of networks; 4) consideration of network affiliations of subcortical structures; and 5) consideration of multi-modal information. We close with minimal reporting guidelines for the cognitive and network neuroscience communities to adopt.

Over the past three decades, research from the field of social neuroscience has identified a constellation of brain regions that relate to social cognition. Although these studies have provided important insights into the specific neural regions underlying social behavior, they may overlook the broader neural context in which those regions and the interactions between them are embedded. Network neuroscience is an emerging discipline that focuses on modeling and analyzing brain networks-collections of interacting neural elements. Because human cognition requires integrating information across multiple brain regions and systems, we argue that a novel social cognitive network neuroscience approach-which leverages methods from the field of network neuroscience and graph theory-can advance our understanding of how brain systems give rise to social behavior. This review provides an overview of the field of network neuroscience, discusses studies that have leveraged this approach to advance social neuroscience research, highlights the potential contributions of social cognitive network neuroscience to understanding social behavior and provides suggested tools and resources for conducting network neuroscience research.

Neural systems can be modeled as complex networks in which neural elements are represented as nodes linked to one another through structural or functional connections. The resulting network can be analyzed using mathematical tools from network science and graph theory to quantify the system's topological organization and to better understand its function. Here, we used two-photon calcium imaging to record spontaneous activity from the same set of cells in mouse auditory cortex over the course of several weeks. We reconstruct functional networks in which cells are linked to one another by edges weighted according to the correlation of their fluorescence traces. We show that the networks exhibit modular structure across multiple topological scales and that these multi-scale modules unfold as part of a hierarchy. We also show that, on average, network architecture becomes increasingly dissimilar over time, with similarity decaying monotonically with the distance (in time) between sessions. Finally, we show that a small fraction of cells maintain strongly-correlated activity over multiple days, forming a stable temporal core surrounded by a fluctuating and variable periphery. Our work indicates a framework for studying spontaneous activity measured by two-photon calcium imaging using computational methods and graphical models from network science. The methods are flexible and easily extended to additional datasets, opening the possibility of studying cellular level network organization of neural systems and how that organization is modulated by stimuli or altered in models of disease.

The brain is a complex dynamical system composed of many interacting sub-regions. Knowledge of how these interactions reconfigure over time is critical to a full understanding of the brain’s functional architecture, the neural basis of flexible cognition and behavior, and how neural systems are disrupted in psychiatric and neurological illness. The idea that we might be able to study neural and cognitive dynamics through analysis of neuroimaging data has catalyzed substantial interest in methods which seek to estimate moment-to-moment fluctuations in functional connectivity (often referred to as “dynamic” or time-varying connectivity; TVC). At the same time, debates have emerged regarding the application of TVC analyses to resting fMRI data, and about the statistical validity, physiological origins, and cognitive relevance of resting TVC. These and other unresolved issues complicate the interpretation of resting TVC findings and limit the insights which can be gained from this otherwise promising research area. This article reviews the current resting TVC literature in light of these issues. We introduce core concepts, define key terms, summarize current controversies and open questions, and present a forward-looking perspective on how resting TVC analyses can be rigorously applied to investigate a wide range of questions in cognitive and systems neuroscience.

ABSTRACT Diffusion-weighted magnetic resonance imaging (dMRI) has become the primary method for non-invasively studying the organization of white matter in the human brain. While many dMRI acquisition sequences have been developed, they all sample q-space in order to characterize water diffusion. Numerous software platforms have been developed for processing dMRI data, but most work on only a subset of sampling schemes or implement only parts of the processing workflow. Reproducible research and comparisons across dMRI methods are hindered by incompatible software, diverse file formats, and inconsistent naming conventions. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing upon a diverse set of software suites to capitalize upon their complementary strengths, QSIPrep automatically applies best practices for dMRI preprocessing, including denoising, distortion correction, head motion correction, coregistration, and spatial normalization. Throughout, QSIPrep provides both visual and quantitative measures of data quality as well as “glass-box” methods reporting. Taken together, these features facilitate easy implementation of best practices for processing of diffusion images while simultaneously ensuring reproducibility.

Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher-order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher-order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.

Previous studies have adopted an edge-centric framework to study fine-scale dynamics in human fMRI. To date, however, no studies have applied this same framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of “bursty” dynamics and events – brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

Abstract Edge time series decompose functional connectivity into its framewise contributions. Previous studies have focused on characterizing the properties of high-amplitude frames (time points when the global co-fluctuation amplitude takes on its largest value), including their cluster structure. Less is known about middle- and low-amplitude co-fluctuations (peaks in co-fluctuation time series but of lower amplitude). Here, we directly address those questions, using data from two dense-sampling studies: the MyConnectome project and Midnight Scan Club. We develop a hierarchical clustering algorithm to group peak co-fluctuations of all magnitudes into nested and multiscale clusters based on their pairwise concordance. At a coarse scale, we find evidence of three large clusters that, collectively, engage virtually all canonical brain systems. At finer scales, however, each cluster is dissolved, giving way to increasingly refined patterns of co-fluctuations involving specific sets of brain systems. We also find an increase in global co-fluctuation magnitude with hierarchical scale. Finally, we comment on the amount of data needed to estimate co-fluctuation pattern clusters and implications for brain-behavior studies. Collectively, the findings reported here fill several gaps in current knowledge concerning the heterogeneity and richness of co-fluctuation patterns as estimated with edge time series while providing some practical guidance for future studies.

Dynamical brain state transitions are critical for flexible working memory but the network mechanisms are incompletely understood. Here, we show that working memory entails brain-wide switching between activity states. The stability of states relates to dopamine D1 receptor gene expression while state transitions are influenced by D2 receptor expression and pharmacological modulation. Schizophrenia patients show altered network control properties, including a more diverse energy landscape and decreased stability of working memory representations.

Subjective mood is a psychophysiological property that depends on complex interactions among the central and peripheral nervous systems. How network interactions in the brain drive temporal fluctuations in mood is unknown. Here we investigate how functional network configuration relates to mood profiles in a single individual over the course of 1 year. Using data from the ‘MyConnectome Project’, we construct a comprehensive mapping between resting-state functional connectivity (FC) patterns and subjective mood scales using an associative multivariate technique (partial least squares). We report three principal findings. First, FC patterns reliably tracked daily fluctuations in mood. Second, positive mood was marked by an integrated architecture, with prominent interactions between canonical resting-state networks. Finally, one of the top-ranked nodes in mood-related network reconfiguration was the subgenual anterior cingulate cortex, an area commonly associated with mood regulation and dysregulation. Altogether, these results showcase the utility of highly sampled individual-focused data sets for affective neuroscience.

Many studies have shown that the human endocrine system modulates brain function, reporting associations between fluctuations in hormone concentrations and both brain activity and connectivity. However, how hormonal fluctuations impact fast changes in brain network structure over short timescales remains unknown. Here, we leverage “edge time series” analysis to investigate the relationship between high-amplitude network states and quotidian variation in sex steroid and gonadotropic hormones in a single individual sampled over the course of two endocrine states, across a natural menstrual cycle and under a hormonal regimen. We find that the frequency of high-amplitude network states are associated with follicle-stimulating and luteinizing hormone, but not the sex hormones estradiol and progesterone. Nevertheless, we show that scan-to-scan variation in the co-fluctuation patterns expressed during network states are robustly linked with the concentration of all four hormones, positing a network-level target of hormonal control. We conclude by speculating on the role of hormones in shaping ongoing brain dynamics.

ABSTRACT Edge time series decompose FC into its framewise contributions. Previous studies have focused on characterizing the properties of high-amplitude frames, including their cluster structure. Less is known about middle- and low-amplitude co-fluctuations. Here, we address those questions directly, using data from two dense-sampling studies: the MyConnectome project and Midnight Scan Club. We develop a hierarchical clustering algorithm to group peak co-fluctuations of all magnitudes into nested and multi-scale clusters based on their pairwise concordance. At a coarse scale, we find evidence of three large clusters that, collectively, engage virtually all canonical brain systems. At finer scales, however, each cluster is dissolved, giving way to increasingly refined patterns of co-fluctuations involving specific sets of brain systems. We also find an increase in global co-fluctuation magnitude with hierarchical scale. Finally, we comment on the amount of data needed to estimate co-fluctuation pattern clusters and implications for brain-behavior studies. Collectively, the findings reported here fill several gaps in current knowledge concerning the heterogeneity and richness of co-fluctuation patterns as estimated with edge time series while providing some practical guidance for future studies.

Network models of anatomical connections allow for the extraction of quantitative features describing brain organization, and their comparison across brains from different species. Such comparisons can inform our understanding of between-species differences in brain architecture and can be compared to existing taxonomies and phylogenies. Here we performed a quantitative comparative analysis using the MaMI database (Tel Aviv University), a collection of brain networks reconstructed from ex vivo diffusion MRI spanning 125 species and 12 taxonomic orders or superorders. We used a broad range of metrics to measure between-mammal distances and compare these estimates to the separation of species as derived from taxonomy and phylogeny. We found that within-taxonomy order network distances are significantly closer than between-taxonomy network distances, and this relation holds for several measures of network distance. Furthermore, to estimate the evolutionary divergence between species, we obtained phylogenetic distances across 10,000 plausible phylogenetic trees. The anatomical network distances were rank-correlated with phylogenetic distances 10,000 times, creating a distribution of coefficients that demonstrate significantly positive correlations between network and phylogenetic distances. Collectively, these analyses demonstrate species-level organization across scales and informational sources: we relate brain networks distances, derived from MRI, with evolutionary distances, derived from genotyping data.

Abstract It is commonplace in neuroscience to assume that if two tasks activate the same brain areas in the same way, then they are recruiting the same underlying networks. Yet computational theory has shown that the same pattern of activity can emerge from many different underlying network representations. Here we evaluated whether similarity in activation necessarily implies similarity in network architecture by comparing region-wise activation patterns and functional correlation profiles from a large sample of healthy subjects (N = 242). Participants performed two executive control tasks known to recruit nearly identical brain areas, the color-word Stroop task and the Multi-Source Interference Task (MSIT). Using a measure of instantaneous functional correlations, based on edge time series, we estimated the task-related networks that differed between incongruent and congruent conditions. We found that the two tasks were much more different in their network profiles than in their evoked activity patterns at different analytical levels, as well as for a wide range of methodological pipelines. Our results reject the notion that having the same activation patterns means two tasks engage the same underlying representations, suggesting that task representations should be independently evaluated at both node and edge (connectivity) levels.

Abstract Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of “bursty” dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

Background: Properties of functional connectivity (FC), such as network integration and segregation, are shown to be associated with various human behaviors. For example, Godwin et al. and Sun et al. found increased integration with attention allocation, whereas Cohen and D'Esposito and Shine et al. observed increased segregation with simple motor tasks. The current study investigated how viewing video clips with different valence and arousal influenced integration-segregation properties in task-based FC networks. Methods: We analyzed an open dataset collected by Kim et al. We performed a generalized psychophysiological interaction (gPPI) analysis paired with network analysis and community detection to investigate changes in brain network dynamics when people watched four types of videos that differed by affective valence (unpleasant or pleasant) and arousal (arousing or calm). Results: Results showed that unpleasant arousing videos produced greater FC deviation from the baseline (task-induced FC deviation [tiFCd]) and perturbed the brain into a more segregated state than other kinds of video. Increased segregation was only observed in association systems, not sensorimotor systems. Discussion: Unpleasant arousing content perturbed the brain to a functionally distinct state from the other three types of affective videos. We suggest that the change in brain state was related to people disengaging from the unpleasant arousing content or, alternatively, staying alert while exposed to unpleasant arousing stimuli. The study also added to our understanding of how combining task-based gPPI analysis with community detection methods and network segregation measures can advance our knowledge of the links between behavior and brain state changes. Impact statement Network integration and segregation is an important property of the human brain. We address the question of how affective stimuli influence brain dynamics from a functional connectivity (FC) network integration-segregation perspective. By conducting a whole-brain generalized psychophysiological interaction (gPPI) analysis paired with community detection methods, we found that highly aversive video content induced significant FC changes and perturbed the brain to a more segregated state.

Scientific discovery in connectomics relies on the use of network null models. To systematically evaluate the prominence of brain network features, empirical measures are compared against null statistics computed in randomized networks. Modern imaging and tracing technologies provide an increasingly rich repertoire of biologically meaningful edge weights. Despite the prevalence of weighted graph analysis in connectomics, randomization models that only preserve binary node degree remain most widely used. Here, to adapt network null models to weighted network inference, we propose a simulated annealing procedure for generating strength sequence-preserving randomized networks. This model outperforms other commonly used rewiring algorithms in preserving weighted degree (strength). We show that these results generalize to directed networks as well as a wide range of real-world networks, making them generically applicable in neuroscience and in other scientific disciplines. Furthermore, we introduce morphospace representation as a tool for the assessment of null network ensemble variability and feature preservation. Finally, we show how the choice of a network null model can yield fundamentally different inferences about established organizational features of the brain such as the rich-club phenomenon and lay out best practices for the use of rewiring algorithms in brain network inference. Collectively, this work provides a simple but powerful inferential method to meet the challenges of analyzing richly detailed next-generation connectomics datasets.

The human brain is never at "rest"; its activity is constantly fluctuating over time, transitioning from one brain state--a whole-brain pattern of activity--to another. Network control theory offers a framework for understanding the effort -- energy -- associated with these transitions. One branch of control theory that is especially useful in this context is "optimal control", in which input signals are used to selectively drive the brain into a target state. Typically, these inputs are introduced independently to the nodes of the network (each input signal is associated with exactly one node). Though convenient, this input strategy ignores the continuity of cerebral cortex -- geometrically, each region is connected to its spatial neighbors, allowing control signals, both exogenous and endogenous, to spread from their foci to nearby regions. Here, we adapt the network control model so that input signals have a spatial extent that decays exponentially from the input site. We show that this more realistic strategy takes advantage of spatial dependencies in structural connectivity and activity to reduce the energy (effort) associated with brain state transitions. We further leverage these dependencies to explore near-optimal control strategies such that, on a per-transition basis, the number of input signals required for a given control task is reduced, in some cases by two orders of magnitude. This approximation yields network-wide maps of input site density, which we compare to an existing database of functional, metabolic, genetic, and neurochemical maps, finding a close correspondence. Ultimately, not only do we propose a more efficient framework that is also more adherent to well-established brain organizational principles, but we also posit neurobiologically grounded bases for optimal control.