Information processing architecture of functionally defined clusters in the macaque cortex

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Early-in-life isoflurane exposure alters resting-state functional connectivity in juvenile non-human primates

BACKGROUND

Clinical studies suggest that anaesthesia exposure early in life affects neurobehavioural development. We designed a non-human primate (NHP) study to evaluate cognitive, behavioural, and brain functional and structural alterations after isoflurane exposure during infancy. These NHPs displayed decreased close social behaviour and increased astrogliosis in specific brain regions, most notably in the amygdala. Here we hypothesise that resting-state functional connectivity MRI can detect alterations in connectivity of brain areas that relate to these social behaviours and astrogliosis.

METHODS

Imaging was performed in 2-yr-old NHPs under light anaesthesia, after early-in-life (postnatal days 6-12) exposure to 5 h of isoflurane either one or three times, or to room air. Brain images were segmented into 82 regions of interest; the amygdala and the posterior cingulate cortex were chosen for a seed-based resting-state functional connectivity MRI analysis.

RESULTS

We found differences between groups in resting-state functional connectivity of the amygdala and the auditory cortices, medial premotor cortex, and posterior cingulate cortex. There were also alterations in resting-state functional connectivity between the posterior cingulate cortex and secondary auditory, polar prefrontal, and temporal cortices, and the anterior insula. Relationships were identified between resting-state functional connectivity alterations and the decrease in close social behaviour and increased astrogliosis.

CONCLUSIONS

Early-in-life anaesthesia exposure in NHPs is associated with resting-state functional connectivity alterations of the amygdala and the posterior cingulate cortex with other brain regions, evident at the juvenile age of 2 yr. These changes in resting-state functional connectivity correlate with the decrease in close social behaviour and increased astrogliosis. Using resting-state functional connectivity MRI to study the neuronal underpinnings of early-in-life anaesthesia-induced behavioural alterations could facilitate development of a biomarker for anaesthesia-induced developmental neurotoxicity.

Abnormal and Changing Information Interaction in Adults with Attention-Deficit/Hyperactivity Disorder Based on Network Motifs

Network motif analysis approaches provide insights into the complexity of the brain's functional network. In recent years, attention-deficit/hyperactivity disorder (ADHD) has been reported to result in abnormal information interactions in macro- and micro-scale functional networks. However, most existing studies remain limited due to potentially ignoring meso-scale topology information. To address this gap, we aimed to investigate functional motif patterns in ADHD to unravel the underlying information flow and analyze motif-based node roles to characterize the different information interaction methods for identifying the abnormal and changing lesion sites of ADHD. The results showed that the interaction functions of the right hippocampus and the right amygdala were significantly increased, which could lead patients to develop mood disorders. The information interaction of the bilateral thalamus changed, influencing and modifying behavioral results. Notably, the capability of receiving information in the left inferior temporal and the right lingual gyrus decreased, which may cause difficulties for patients in processing visual information in a timely manner, resulting in inattention. This study revealed abnormal and changing information interactions based on network motifs, providing important evidence for understanding information interactions at the meso-scale level in ADHD patients.

An open resource combining multi-contrast MRI and microscopy in the macaque brain

Understanding brain structure and function often requires combining data across different modalities and scales to link microscale cellular structures to macroscale features of whole brain organisation. Here we introduce the BigMac dataset, a resource combining in vivo MRI, extensive postmortem MRI and multi-contrast microscopy for multimodal characterisation of a single whole macaque brain. The data spans modalities (MRI and microscopy), tissue states (in vivo and postmortem), and four orders of spatial magnitude, from microscopy images with micrometre or sub-micrometre resolution, to MRI signals on the order of millimetres. Crucially, the MRI and microscopy images are carefully co-registered together to facilitate quantitative multimodal analyses. Here we detail the acquisition, curation, and first release of the data, that together make BigMac a unique, openly-disseminated resource available to researchers worldwide. Further, we demonstrate example analyses and opportunities afforded by the data, including improvement of connectivity estimates from ultra-high angular resolution diffusion MRI, neuroanatomical insight provided by polarised light imaging and myelin-stained histology, and the joint analysis of MRI and microscopy data for reconstruction of the microscopy-inspired connectome. All data and code are made openly available.

Explaining cocktail party effect and McGurk effect with a spiking neural network improved by Motif-topology

Network architectures and learning principles have been critical in developing complex cognitive capabilities in artificial neural networks (ANNs). Spiking neural networks (SNNs) are a subset of ANNs that incorporate additional biological features such as dynamic spiking neurons, biologically specified architectures, and efficient and useful paradigms. Here we focus more on network architectures in SNNs, such as the meta operator called 3-node network motifs, which is borrowed from the biological network. We proposed a Motif-topology improved SNN (M-SNN), which is further verified efficient in explaining key cognitive phenomenon such as the cocktail party effect (a typical noise-robust speech-recognition task) and McGurk effect (a typical multi-sensory integration task). For M-SNN, the Motif topology is obtained by integrating the spatial and temporal motifs. These spatial and temporal motifs are first generated from the pre-training of spatial (e.g., MNIST) and temporal (e.g., TIDigits) datasets, respectively, and then applied to the previously introduced two cognitive effect tasks. The experimental results showed a lower computational cost and higher accuracy and a better explanation of some key phenomena of these two effects, such as new concept generation and anti-background noise. This mesoscale network motifs topology has much room for the future.

Structure–Function Decoupling: A Novel Perspective for Understanding the Radiation-Induced Brain Injury in Patients With Nasopharyngeal Carcinoma

Radiation-induced functional and structural brain alterations are well documented in patients with nasopharyngeal carcinoma (NPC), followed by radiotherapy (RT); however, alterations in structure–function coupling remain largely unknown. Herein, we aimed to assess radiation-induced structure–function decoupling and its importance in predicting radiation encephalopathy (RE). We included 62 patients with NPC (22 patients in the pre-RT cohort, 18 patients in the post-RT-RE_+ve cohort, and 22 patients in the post-RT-RE_–ve cohort). A metric of regional homogeneity (ReHo)/voxel-based morphometry (VBM) was used to detect radiation-induced structure–function decoupling, which was then used as a feature to construct a predictive model for RE. Compared with the pre-RT group, patients in the post-RT group (which included post-RT-RE_+ve and post-RT-RE_–ve) showed higher ReHo/VBM coupling values in the substantia nigra (SN), the putamen, and the bilateral thalamus and lower values in the brain stem, the cerebellum, the bilateral medial temporal lobes (MTLs), the bilateral insula, the right precentral and postcentral gyri, the medial prefrontal cortex (MPFC), and the left inferior parietal lobule (IPL). In the post-RT group, negative correlations were observed between maximum dosage of RT (MDRT) to the ipsilateral temporal lobe and ReHo/VBM values in the ipsilateral middle temporal gyrus (MTG). Moreover, structure–function decoupling in the bilateral superior temporal gyrus (STG), the bilateral precentral and postcentral gyri, the paracentral lobules, the right precuneus and IPL, and the right MPFC exhibited excellent predictive performance (accuracy = 88.0%) in identifying patients likely to develop RE. These findings show that ReHo/VBM may be a novel effective imaging metric that reflects the neural mechanism underlying RE in patients with NPC.

A set of hub neurons and non-local connectivity features support global brain dynamics in C. elegans

The wiring architecture of neuronal networks is assumed to be a strong determinant of their dynamical computations. An ongoing effort in neuroscience is therefore to generate comprehensive synapse-resolution connectomes alongside brain-wide activity maps. However, the structure-function relationship, i.e., how the anatomical connectome and neuronal dynamics relate to each other on a global scale, remains unsolved. Systematically, comparing graph features in the C. elegans connectome with correlations in nervous system-wide neuronal dynamics, we found that few local connectivity motifs and mostly other non-local features such as triplet motifs and input similarities can predict functional relationships between neurons. Surprisingly, quantities such as connection strength and amount of common inputs do not improve these predictions, suggesting that the network's topology is sufficient. We demonstrate that hub neurons in the connectome are key to these relevant graph features. Consistently, inhibition of multiple hub neurons specifically disrupts brain-wide correlations. Thus, we propose that a set of hub neurons and non-local connectivity features provide an anatomical substrate for global brain dynamics.

A Riemannian approach to predicting brain function from the structural connectome

Ongoing brain function is largely determined by the underlying wiring of the brain, but the specific rules governing this relationship remain unknown. Emerging literature has suggested that functional interactions between brain regions emerge from the structural connections through mono- as well as polysynaptic mechanisms. Here, we propose a novel approach based on diffusion maps and Riemannian optimization to emulate this dynamic mechanism in the form of random walks on the structural connectome and predict functional interactions as a weighted combination of these random walks. Our proposed approach was evaluated in two different cohorts of healthy adults (Human Connectome Project, HCP; Microstructure-Informed Connectomics, MICs). Our approach outperformed existing approaches and showed that performance plateaus approximately around the third random walk. At macroscale, we found that the largest number of walks was required in nodes of the default mode and frontoparietal networks, underscoring an increasing relevance of polysynaptic communication mechanisms in transmodal cortical networks compared to primary and unimodal systems.

Structure-function coupling within the reward network in preschool children predicts executive functioning in later childhood

Early differences in reward behavior have been linked to executive functioning development. The nucleus accumbens (NAc) and orbitofrontal cortex (OFC) are activated by reward-related tasks and identified as key nodes of the brain circuit that underlie reward processing. We aimed to investigate the relation between NAc-OFC structural and functional connectivity in preschool children, as well as associations with future reward sensitivity and executive function. We showed that NAc-OFC structural and functional connectivity were not significantly associated in preschool children, but both independently predicted sensitivity to reward in males in a left-lateralized manner. Moreover, significant NAc-OFC structure-function coupling was only found in individuals who performed poorly on executive function tasks in later childhood, but not in the middle- and high-performing groups. As structure-function coupling is proposed to measure functional specialization, this finding suggests premature functional specialization within the reward network, which may impede dynamic communication with other regions, affects executive function development. Our study also highlights the utility of multimodal imaging data integration when studying the effects of reward network functional flexibility in the preschool age, a critical period in brain and executive function development.

•

Functional connectivity is not tethered to structural connectivity in preschool age.

•

Higher degree of SC-FC coupling reflects lower plasticity in early childhood.

•

Gender differences in reward sensitivity were present as early as in preschool age.

•

Early reward network SC-FC coupling affects later executive function.

Personalized Connectome-Based Modeling in Patients with Semi-Acute Phase TBI: Relationship to Acute Neuroimaging and 6 Month Follow-Up

Following traumatic brain injury (TBI), cognitive impairments manifest through interactions between microscopic and macroscopic changes. On the microscale, a neurometabolic cascade alters neurotransmission, while on the macroscale diffuse axonal injury impacts the integrity of long-range connections. Large-scale brain network modeling allows us to make predictions across these spatial scales by integrating neuroimaging data with biophysically based models to investigate how microscale changes invisible to conventional neuroimaging influence large-scale brain dynamics. To this end, we analyzed structural and functional neuroimaging data from a well characterized sample of 44 adult TBI patients recruited from a regional trauma center, scanned at 1–2 weeks postinjury, and with follow-up behavioral outcome assessed 6 months later. Thirty-six age-matched healthy adults served as comparison participants. Using The Virtual Brain, we fit simulations of whole-brain resting-state functional MRI to the empirical static and dynamic functional connectivity of each participant. Multivariate partial least squares (PLS) analysis showed that patients with acute traumatic intracranial lesions had lower cortical regional inhibitory connection strengths than comparison participants, while patients without acute lesions did not differ from the comparison group. Further multivariate PLS analyses found correlations between lower semiacute regional inhibitory connection strengths and more symptoms and lower cognitive performance at a 6 month follow-up. Critically, patients without acute lesions drove this relationship, suggesting clinical relevance of regional inhibitory connection strengths even when traumatic intracranial lesions were not present. Our results suggest that large-scale connectome-based models may be sensitive to pathophysiological changes in semi-acute phase TBI patients and predictive of their chronic outcomes.

Prefrontal connectomics: from anatomy to human imaging

The fundamental importance of prefrontal cortical connectivity to information processing and, therefore, disorders of cognition, emotion, and behavior has been recognized for decades. Anatomic tracing studies in animals have formed the basis for delineating the direct monosynaptic connectivity, from cells of origin, through axon trajectories, to synaptic terminals. Advances in neuroimaging combined with network science have taken the lead in developing complex wiring diagrams or connectomes of the human brain. A key question is how well these magnetic resonance imaging (MRI)-derived networks and hubs reflect the anatomic "hard wiring" first proposed to underlie the distribution of information for large-scale network interactions. In this review, we address this challenge by focusing on what is known about monosynaptic prefrontal cortical connections in non-human primates and how this compares to MRI-derived measurements of network organization in humans. First, we outline the anatomic cortical connections and pathways for each prefrontal cortex (PFC) region. We then review the available MRI-based techniques for indirectly measuring structural and functional connectivity, and introduce graph theoretical methods for analysis of hubs, modules, and topologically integrative features of the connectome. Finally, we bring these two approaches together, using specific examples, to demonstrate how monosynaptic connections, demonstrated by tract-tracing studies, can directly inform understanding of the composition of PFC nodes and hubs, and the edges or pathways that connect PFC to cortical and subcortical areas.

Differential structure-function network coupling in the inattentive and combined types of attention deficit hyperactivity disorder

The heterogeneous presentation of inattentive and hyperactive-impulsive core symptoms in attention deficit hyperactivity disorder (ADHD) warrants further investigation into brain network connectivity as a basis for subtype divisions in this prevalent disorder. With diffusion and resting-state functional magnetic resonance imaging data from the Healthy Brain Network database, we analyzed both structural and functional network efficiency and structure-functional network (SC-FC) coupling at the default mode (DMN), executive control (ECN), and salience (SAN) intrinsic networks in 201 children diagnosed with the inattentive subtype (ADHD-I), the combined subtype (ADHD-C), and typically developing children (TDC) to characterize ADHD symptoms relative to TDC and to test differences between ADHD subtypes. Relative to TDC, children with ADHD had lower structural connectivity and network efficiency in the DMN, without significant group differences in functional networks. Children with ADHD-C had higher SC-FC coupling, a finding consistent with diminished cognitive flexibility, for all subnetworks compared to TDC. The ADHD-C group also demonstrated increased SC-FC coupling in the DMN compared to the ADHD-I group. The correlation between SC-FC coupling and hyperactivity scores was negative in the ADHD-I, but not in the ADHD-C group. The current study suggests that ADHD-C and ADHD-I may differ with respect to their underlying neuronal connectivity and that the added dimensionality of hyperactivity may not explain this distinction.

Mapping thalamocortical functional connectivity with large-scale brain networks in patients with first-episode psychosis

Abnormal thalamocortical networks involving specific thalamic nuclei have been implicated in schizophrenia pathophysiology. While comparable topography of anatomical and functional connectivity abnormalities has been reported in patients across illness stages, previous functional studies have been confined to anatomical pathways of thalamocortical networks. To address this issue, we incorporated large-scale brain network dynamics into examining thalamocortical functional connectivity. Forty patients with first-episode psychosis and forty healthy controls underwent T1-weighted and resting-state functional magnetic resonance imaging. Independent component analysis of voxelwise thalamic functional connectivity maps parcellated the cortex into thalamus-related networks, and thalamic subdivisions associated with these networks were delineated. Functional connectivity of (1) networks with the thalamus and (2) thalamic subdivision seeds were examined. In patients, functional connectivity of the salience network with the thalamus was decreased and localized to the ventrolateral (VL) and ventroposterior (VP) thalamus, while that of a network comprising the cerebellum, temporal and parietal regions was increased and localized to the mediodorsal (MD) thalamus. In patients, thalamic subdivision encompassing the VL and VP thalamus demonstrated hypoconnectivity and that encompassing the MD and pulvinar regions demonstrated hyperconnectivity. Our results extend the implications of disrupted thalamocortical networks involving specific thalamic nuclei to dysfunctional large-scale brain network dynamics in schizophrenia pathophysiology.

Heritability and interindividual variability of regional structure-function coupling

White matter structural connections are likely to support flow of functional activation or functional connectivity. While the relationship between structural and functional connectivity profiles, here called SC-FC coupling, has been studied on a whole-brain, global level, few studies have investigated this relationship at a regional scale. Here we quantify regional SC-FC coupling in healthy young adults using diffusion-weighted MRI and resting-state functional MRI data from the Human Connectome Project and study how SC-FC coupling may be heritable and varies between individuals. We show that regional SC-FC coupling strength varies widely across brain regions, but was strongest in highly structurally connected visual and subcortical areas. We also show interindividual regional differences based on age, sex and composite cognitive scores, and that SC-FC coupling was highly heritable within certain networks. These results suggest regional structure-function coupling is an idiosyncratic feature of brain organisation that may be influenced by genetic factors.

Virtual Connectomic Datasets in Alzheimer's Disease and Aging Using Whole-Brain Network Dynamics Modelling

Large neuroimaging datasets, including information about structural connectivity (SC) and functional connectivity (FC), play an increasingly important role in clinical research, where they guide the design of algorithms for automated stratification, diagnosis or prediction. A major obstacle is, however, the problem of missing features [e.g., lack of concurrent DTI SC and resting-state functional magnetic resonance imaging (rsfMRI) FC measurements for many of the subjects]. We propose here to address the missing connectivity features problem by introducing strategies based on computational whole-brain network modeling. Using two datasets, the Alzheimer’s Disease Neuroimaging Initiative (ADNI) dataset and a healthy aging dataset, for proof-of-concept, we demonstrate the feasibility of virtual data completion (i.e., inferring “virtual FC” from empirical SC or “virtual SC” from empirical FC), by using self-consistent simulations of linear and nonlinear brain network models. Furthermore, by performing machine learning classification (to separate age classes or control from patient subjects), we show that algorithms trained on virtual connectomes achieve discrimination performance comparable to when trained on actual empirical data; similarly, algorithms trained on virtual connectomes can be used to successfully classify novel empirical connectomes. Completion algorithms can be combined and reiterated to generate realistic surrogate connectivity matrices in arbitrarily large number, opening the way to the generation of virtual connectomic datasets with network connectivity information comparable to the one of the original data.

Integrating Multimodal and Longitudinal Neuroimaging Data with Multi-Source Network Representation Learning

Uncovering the complex network of the brain is of great interest to the field of neuroimaging. Mining from these rich datasets, scientists try to unveil the fundamental biological mechanisms in the human brain. However, neuroimaging data collected for constructing brain networks is generally costly, and thus extracting useful information from a limited sample size of brain networks is demanding. Currently, there are two common trends in neuroimaging data collection that could be exploited to gain more information: 1) multimodal data, and 2) longitudinal data. It has been shown that these two types of data provide complementary information. Nonetheless, it is challenging to learn brain network representations that can simultaneously capture network properties from multimodal as well as longitudinal datasets. Here we propose a general fusion framework for multi-source learning of brain networks - multimodal brain network fusion with longitudinal coupling (MMLC). In our framework, three layers of information are considered, including cross-sectional similarity, multimodal coupling, and longitudinal consistency. Specifically, we jointly factorize multimodal networks and construct a rotation-based constraint to couple network variance across time. We also adopt the consensus factorization as the group consistent pattern. Using two publicly available brain imaging datasets, we demonstrate that MMLC may better predict psychometric scores than some other state-of-the-art brain network representation learning algorithms. Additionally, the discovered significant brain regions are synergistic with previous literature. Our new approach may boost statistical power and sheds new light on neuroimaging network biomarkers for future psychometric prediction research by integrating longitudinal and multimodal neuroimaging data.

Brain network motifs are markers of loss and recovery of consciousness

Motifs are patterns of inter-connections between nodes of a network, and have been investigated as building blocks of directed networks. This study explored the re-organization of 3-node motifs during loss and recovery of consciousness. Nine healthy subjects underwent a 3-h anesthetic protocol while 128-channel electroencephalography (EEG) was recorded. In the alpha (8-13 Hz) band, 5-min epochs of EEG were extracted for: Baseline; Induction; Unconscious; 30-, 10- and 5-min pre-recovery of responsiveness; 30- and 180-min post-recovery of responsiveness. We constructed a functional brain network using the weighted and directed phase lag index, on which we calculated the frequency and topology of 3-node motifs. Three motifs (motifs 1, 2 and 5) were significantly present across participants and epochs, when compared to random networks (p < 0.05). The topology of motifs 1 and 5 changed significantly between responsive and unresponsive epochs (p-values < 0.01; Kendall's W = 0.664 (motif 1) and 0.529 (motif 5)). Motif 1 was constituted of long-range chain-like connections, while motif 5 was constituted of short-range, loop-like connections. Our results suggest that anesthetic-induced unconsciousness is associated with a topological re-organization of network motifs. As motif topological re-organization may precede (motif 5) or accompany (motif 1) the return of responsiveness, motifs could contribute to the understanding of the neural correlates of consciousness.

Network topology of the marmoset connectome

The brain is a complex network of interconnected and interacting neuronal populations. Global efforts to understand the emergence of behavior and the effect of perturbations depend on accurate reconstruction of white matter pathways, both in humans and in model organisms. An emerging animal model for next-generation applied neuroscience is the common marmoset (Callithrix jacchus). A recent open respository of retrograde and anterograde tract tracing presents an opportunity to systematically study the network architecture of the marmoset brain (Marmoset Brain Architecture Project; http://www.marmosetbrain.org). Here we comprehensively chart the topological organization of the mesoscale marmoset cortico-cortical connectome. The network possesses multiple nonrandom attributes that promote a balance between segregation and integration, including near-minimal path length, multiscale community structure, a connective core, a unique motif composition, and multiple cavities. Altogether, these structural attributes suggest a link between network architecture and function. Our findings are consistent with previous reports across a range of species, scales, and reconstruction technologies, suggesting a small set of organizational principles universal across phylogeny. Collectively, these results provide a foundation for future anatomical, functional, and behavioral studies in this model organism.

MECP2 Duplication Causes Aberrant GABA Pathways, Circuits and Behaviors in Transgenic Monkeys: Neural Mappings to Patients with Autism

MECP2 gain-of-function and loss-of-function in genetically engineered monkeys recapitulates typical phenotypes in patients with autism, yet where MECP2 mutation affects the monkey brain and whether/how it relates to autism pathology remain unknown. Here we report a combination of gene-circuit-behavior analyses including MECP2 coexpression network, locomotive and cognitive behaviors, and EEG and fMRI findings in 5 MECP2 overexpressed monkeys (Macaca fascicularis; 3 females) and 20 wild-type monkeys (Macaca fascicularis; 11 females). Whole-genome expression analysis revealed MECP2 coexpressed genes significantly enriched in GABA-related signaling pathways, whereby reduced β-synchronization within fronto-parieto-occipital networks was associated with abnormal locomotive behaviors. Meanwhile, MECP2-induced hyperconnectivity in prefrontal and cingulate networks accounted for regressive deficits in reversal learning tasks. Furthermore, we stratified a cohort of 49 patients with autism and 72 healthy controls of 1112 subjects using functional connectivity patterns, and identified dysconnectivity profiles similar to those in monkeys. By establishing a circuit-based construct link between genetically defined models and stratified patients, these results pave new avenues to deconstruct clinical heterogeneity and advance accurate diagnosis in psychiatric disorders.SIGNIFICANCE STATEMENT Autism spectrum disorder (ASD) is a complex disorder with co-occurring symptoms caused by multiple genetic variations and brain circuit abnormalities. To dissect the gene-circuit-behavior causal chain underlying ASD, animal models are established by manipulating causative genes such as MECP2 However, it is unknown whether such models have captured any circuit-level pathology in ASD patients, as demonstrated by human brain imaging studies. Here, we use transgenic macaques to examine the causal effect of MECP2 overexpression on gene coexpression, brain circuits, and behaviors. For the first time, we demonstrate that the circuit abnormalities linked to MECP2 and autism-like traits in the monkeys can be mapped to a homogeneous ASD subgroup, thereby offering a new strategy to deconstruct clinical heterogeneity in ASD.

Linking Structure and Function in Macroscale Brain Networks

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.

Distinct structure-function relationships across cortical regions and connectivity scales in the rat brain

An improved understanding of the structure-function relationship in the brain is necessary to know to what degree structural connectivity underpins abnormal functional connectivity seen in disorders. We integrated high-field resting-state fMRI-based functional connectivity with high-resolution macro-scale diffusion-based and meso-scale neuronal tracer-based structural connectivity, to obtain an accurate depiction of the structure-function relationship in the rat brain. Our main goal was to identify to what extent structural and functional connectivity strengths are correlated, macro- and meso-scopically, across the cortex. Correlation analyses revealed a positive correspondence between functional and macro-scale diffusion-based structural connectivity, but no significant correlation between functional connectivity and meso-scale neuronal tracer-based structural connectivity. Zooming in on individual connections, we found strong functional connectivity in two well-known resting-state networks: the sensorimotor and default mode network. Strong functional connectivity within these networks coincided with strong short-range intrahemispheric structural connectivity, but with weak heterotopic interhemispheric and long-range intrahemispheric structural connectivity. Our study indicates the importance of combining measures of connectivity at distinct hierarchical levels to accurately determine connectivity across networks in the healthy and diseased brain. Although characteristics of the applied techniques may affect where structural and functional networks (dis)agree, distinct structure-function relationships across the brain could also have a biological basis.

Development of structure-function coupling in human brain networks during youth

The human brain is organized into a hierarchy of functional systems that evolve in childhood and adolescence to support the dynamic control of attention and behavior. However, it remains unknown how developing white-matter architecture supports coordinated fluctuations in neural activity underlying cognition. We document marked remodeling of structure–function coupling in youth, which aligns with cortical hierarchies of functional specialization and evolutionary expansion. Further, we demonstrate that structure–function coupling in rostrolateral prefrontal cortex supports age-related improvements in executive ability. These findings have broad relevance for accounts of experience-dependent plasticity in healthy development and abnormal development associated with neuropsychiatric illness.

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.

A macaque connectome for large-scale network simulations in TheVirtualBrain

Models of large-scale brain networks that are informed by the underlying anatomical connectivity contribute to our understanding of the mapping between the structure of the brain and its dynamical function. Connectome-based modelling is a promising approach to a more comprehensive understanding of brain function across spatial and temporal scales, but it must be constrained by multi-scale empirical data from animal models. Here we describe the construction of a macaque (Macaca mulatta and Macaca fascicularis) connectome for whole-cortex simulations in TheVirtualBrain, an open-source simulation platform. We take advantage of available axonal tract-tracing datasets and enhance the existing connectome data using diffusion-based tractography in macaques. We illustrate the utility of the connectome as an extension of TheVirtualBrain by simulating resting-state BOLD-fMRI data and fitting it to empirical resting-state data.

Functional specialization of macaque premotor F5 subfields with respect to hand and mouth movements: A comparison of task and resting-state fMRI

Based on architectonic, tract-tracing or functional criteria, the rostral portion of ventral premotor cortex in the macaque monkey, also termed area F5, has been divided into several subfields. Cytoarchitectonical investigations suggest the existence of three subfields, F5c (convexity), F5p (posterior) and F5a (anterior). Electrophysiological investigations have suggested a gradual dorso-ventral transition from hand- to mouth-dominated motor fields, with F5p and ventral F5c strictly related to hand movements and mouth movements, respectively. The involvement of F5a in this respect, however, has received much less attention. Recently, data-driven resting-state fMRI approaches have also been used to examine the presence of distinct functional fields in macaque ventral premotor cortex. Although these studies have suggested several functional clusters in/near macaque F5, so far the parcellation schemes derived from these clustering methods do not completely retrieve the same level of F5 specialization as suggested by aforementioned invasive techniques. Here, using seed-based resting-state fMRI analyses, we examined the functional connectivity of different F5 seeds with key regions of the hand and face/mouth parieto-frontal-insular motor networks. In addition, we trained monkeys to perform either hand grasping or ingestive mouth movements in the scanner in order to compare resting-state with task-derived functional hand and mouth motor networks. In line with previous single-cell investigations, task-fMRI suggests involvement of F5p, dorsal F5c and F5a in the execution of hand grasping movements, while non-communicative mouth movements yielded particularly pronounced responses in ventral F5c. Corroborating with anatomical tracing data of macaque F5 subfields, seed-based resting-state fMRI suggests a transition from predominant functional correlations with the hand-motor network in F5p to mostly mouth-motor network functional correlations in ventral F5c. Dorsal F5c yielded robust functional correlations with both hand- and mouth-motor networks. In addition, the deepest part of the fundus of the inferior arcuate, corresponding to area 44, displayed a strikingly different functional connectivity profile compared to neighboring F5a, suggesting a different functional specialization for these two neighboring regions.

Exploring the limits of network topology estimation using diffusion-based tractography and tracer studies in the macaque cortex

Reconstructing the anatomical pathways of the brain to study the human connectome has become an important endeavour for understanding brain function and dynamics. Reconstruction of the cortico-cortical connectivity matrix in vivo often relies on noninvasive diffusion-weighted imaging (DWI) techniques but the extent to which they can accurately represent the topological characteristics of structural connectomes remains unknown. We addressed this question by constructing connectomes using DWI data collected from macaque monkeys in vivo and with data from published invasive tracer studies. We found the strength of fiber tracts was well estimated from DWI and topological properties like degree and modularity were captured by tractography-based connectomes. Rich-club/core-periphery type architecture could also be detected but the classification of hubs using betweenness centrality, participation coefficient and core-periphery identification techniques was inaccurate. Our findings indicate that certain aspects of cortical topology can be faithfully represented in noninvasively-obtained connectomes while other network analytic measures warrant cautionary interpretations.

A systematic review on the quantitative relationship between structural and functional network connectivity strength in mammalian brains

The mammalian brain is composed of densely connected and interacting regions, which form structural and functional networks. An improved understanding of the structure-function relation is crucial to understand the structural underpinnings of brain function and brain plasticity after injury. It is currently unclear how functional connectivity strength relates to structural connectivity strength. We obtained an overview of recent papers that report on correspondences between quantitative functional and structural connectivity measures in the mammalian brain. We included network studies in which functional connectivity was measured with resting-state fMRI, and structural connectivity with either diffusion-weighted MRI or neuronal tract tracers. Twenty-seven of the 28 included studies showed a positive structure-function relationship. Large inter-study variations were found comparing functional connectivity strength with either quantitative diffusion-based (correlation coefficient (r) ranges: 0.18-0.82) or neuronal tracer-based structural connectivity measures (r = 0.24-0.74). Two functional datasets demonstrated lower structure-function correlations with neuronal tracer-based (r = 0.22 and r = 0.30) than with diffusion-based measures (r = 0.49 and r = 0.65). The robust positive quantitative structure-function relationship supports the hypothesis that structural connectivity provides the hardware from which functional connectivity emerges. However, methodological differences between the included studies complicate the comparison across studies, which emphasize the need for validation and standardization in brain structure-function studies.

Macro-connectomics and microstructure predict dynamic plasticity patterns in the non-human primate brain

The brain displays a remarkable ability to adapt following injury by altering its connections through neural plasticity. Many of the biological mechanisms that underlie plasticity are known, but there is little knowledge as to when, or where in the brain plasticity will occur following injury. This knowledge could guide plasticity-promoting interventions and create a more accurate roadmap of the recovery process following injury. We causally investigated the time-course of plasticity after hippocampal lesions using multi-modal MRI in monkeys. We show that post-injury plasticity is highly dynamic, but also largely predictable on the basis of the functional connectivity of the lesioned region, gradients of cell densities across the cortex and the pre-lesion network structure of the brain. The ability to predict which brain areas will plastically adapt their functional connectivity following injury may allow us to decipher why some brain lesions lead to permanent loss of cognitive function, while others do not.

The brain has the ability to adapt after injury, a process known as plasticity. When one area sustains damage, for example following a car accident or stroke, other areas change their activity and structure to compensate. Understanding how this happens is critical to helping people recover from brain injuries. Certain factors may affect how well the brain can repair itself. These include how much the damaged area interacts with other areas, and which cell types different areas of the brain contain.

Froudist-Walsh et al. set out to determine how these factors influence recovery from brain injury in monkeys, whose brains are similar to our own. The monkeys had damage to a structure called the hippocampus. This part of the brain has a key role in memory, which is often impaired in patients with brain injuries. The hippocampus cannot repair itself because the brain has only a limited capacity to grow new neurons. Instead, the brain attempts to compensate for disruption to the hippocampus via changes in other, undamaged areas.

Using brain imaging, Froudist-Walsh et al. show that the types of changes that occur depend on how much time has passed since the injury. In the first three months, many areas of the brain change how much they coordinate their activity with other areas. Highly connected areas reduce their communication with other areas the most. In the long-term, the responses of brain areas depend more on which cell types they contain. Areas with more support cells known as “glia” – which supply nutrients and energy to neurons – are better able to adapt their connectivity up to a year after the injury.

These findings may ultimately benefit people who have suffered brain injuries after accidents or stroke. They suggest that stimulating intact brain areas may be helpful in the months immediately after an injury. By contrast, long-term therapy may need to focus more on structural repair. Future studies must build on these results to discover the best ways to induce successful recovery from brain injury.

A multi-scale layer-resolved spiking network model of resting-state dynamics in macaque visual cortical areas

Cortical activity has distinct features across scales, from the spiking statistics of individual cells to global resting-state networks. We here describe the first full-density multi-area spiking network model of cortex, using macaque visual cortex as a test system. The model represents each area by a microcircuit with area-specific architecture and features layer- and population-resolved connectivity between areas. Simulations reveal a structured asynchronous irregular ground state. In a metastable regime, the network reproduces spiking statistics from electrophysiological recordings and cortico-cortical interaction patterns in fMRI functional connectivity under resting-state conditions. Stable inter-area propagation is supported by cortico-cortical synapses that are moderately strong onto excitatory neurons and stronger onto inhibitory neurons. Causal interactions depend on both cortical structure and the dynamical state of populations. Activity propagates mainly in the feedback direction, similar to experimental results associated with visual imagery and sleep. The model unifies local and large-scale accounts of cortex, and clarifies how the detailed connectivity of cortex shapes its dynamics on multiple scales. Based on our simulations, we hypothesize that in the spontaneous condition the brain operates in a metastable regime where cortico-cortical projections target excitatory and inhibitory populations in a balanced manner that produces substantial inter-area interactions while maintaining global stability.

The mammalian cortex fulfills its complex tasks by operating on multiple temporal and spatial scales from single cells to entire areas comprising millions of cells. These multi-scale dynamics are supported by specific network structures at all levels of organization. Since models of cortex hitherto tend to concentrate on a single scale, little is known about how cortical structure shapes the multi-scale dynamics of the network. We here present dynamical simulations of a multi-area network model at neuronal and synaptic resolution with population-specific connectivity based on extensive experimental data which accounts for a wide range of dynamical phenomena. Our model elucidates relationships between local and global scales in cortex and provides a platform for future studies of cortical function.

Stage-dependent remodeling of projections to motor cortex in ALS mouse model revealed by a new variant retrograde-AAV9

Amyotrophic lateral sclerosis (ALS) is characterized by the progressive degeneration of motoneurons in the primary motor cortex (pMO) and in spinal cord. However, the pathogenic process involves multiple subnetworks in the brain and functional MRI studies demonstrate an increase in functional connectivity in areas connected to pMO despite the ongoing neurodegeneration. The extent and the structural basis of the motor subnetwork remodeling in experimentally tractable models remain unclear. We have developed a new retrograde AAV9 to quantitatively map the projections to pMO in the SOD1(G93A) ALS mouse model. We show an increase in the number of neurons projecting from somatosensory cortex to pMO at presymptomatic stages, followed by an increase in projections from thalamus, auditory cortex and contralateral MO (inputs from 20 other structures remains unchanged) as disease advances. The stage- and structure-dependent remodeling of projection to pMO in ALS may provide insights into the hyperconnectivity observed in ALS patients.

Frontoparietal Functional Connectivity in the Common Marmoset

In contrast to the well established macaque monkey, little is known about functional connectivity patterns of common marmoset monkey (Callithrix jacchus) that is poised to become the leading transgenic primate model. Here, we used resting-state ultra-high-field fMRI data collected from anesthetized marmosets and macaques along with awake human subjects, to examine and compare the brain's functional organization, with emphasis on the saccade system. Exploratory independent component analysis revealed eight resting-state networks in marmosets that greatly overlapped with corresponding macaque and human networks including a distributed frontoparietal network. Seed-region analyses of the superior colliculus (SC) showed homolog areas in macaques and marmosets. The marmoset SC displayed the strongest frontal functional connectivity with area 8aD at the border to area 6DR. Functional connectivity of this frontal region revealed a similar functional connectivity pattern as the frontal eye fields in macaques and humans. Furthermore, areas 8aD, 8aV, PG,TPO, TE2, and TE3 were identified as major hubs based on region-wise evaluation of betweeness centrality, suggesting that these cortical regions make up the functional core of the marmoset brain. The results support an evolutionarily preserved frontoparietal system and provide a starting point for invasive neurophysiological studies in the marmoset saccade and visual systems.

Delineating the Macroscale Areal Organization of the Macaque Cortex In Vivo

Complementing long-standing traditions centered on histology, fMRI approaches are rapidly maturing in delineating brain areal organization at the macroscale. The non-human primate (NHP) provides the opportunity to overcome critical barriers in translational research. Here, we establish the data requirements for achieving reproducible and internally valid parcellations in individuals. We demonstrate that functional boundaries serve as a functional fingerprint of the individual animals and can be achieved under anesthesia or awake conditions (rest, naturalistic viewing), though differences between awake and anesthetized states precluded the detection of individual differences across states. Comparison of awake and anesthetized states suggested a more nuanced picture of changes in connectivity for higher-order association areas, as well as visual and motor cortex. These results establish feasibility and data requirements for the generation of reproducible individual-specific parcellations in NHPs, provide insights into the impact of scan state, and motivate efforts toward harmonizing protocols.

Structure-function relationships during segregated and integrated network states of human brain functional connectivity

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.

Neural Reorganization Due to Neonatal Amygdala Lesions in the Rhesus Monkey: Changes in Morphology and Network Structure

It is generally believed that neural damage that occurs early in development is associated with greater adaptive capacity relative to similar damage in an older individual. However, few studies have surveyed whole brain changes following early focal damage. In this report, we employed multimodal magnetic resonance imaging analyses of adult rhesus macaque monkeys who had previously undergone bilateral, neurotoxic lesions of the amygdala at about 2 weeks of age. A deformation-based morphometric approach demonstrated reduction of the volumes of the anterior temporal lobe, anterior commissure, basal ganglia, and pulvinar in animals with early amygdala lesions compared to controls. In contrast, animals with early amygdala lesions had an enlarged cingulate cortex, medial superior frontal gyrus, and medial parietal cortex. Diffusion-weighted imaging tractography and network analysis were also used to compare connectivity patterns and higher-level measures of communication across the brain. Using the communicability metric, which integrates direct and indirect paths between regions, lesioned animals showed extensive degradation of network integrity in the temporal and orbitofrontal cortices. This work demonstrates both degenerative as well as progressive large-scale neural changes following long-term recovery from neonatal focal brain damage.

Contribution of diffusion, perfusion and functional MRI to the disconnection hypothesis in subcortical vascular cognitive impairment

Vascular cognitive impairment (VCI) describes all forms of cognitive impairment caused by any type of cerebrovascular disease. Early identification of VCI is quite difficult due to the lack of both sensitive and specific biomarkers. Extensive damage to the white matter tracts, which connect the cortical and subcortical regions, has been shown in subcortical VCI (SVCI), the most common subtype of VCI that is caused by small vessel disease. Two specific MRI sequences, including diffusion tensor imaging (DTI) and functional MRI (fMRI), have emerged as useful tools for identifying subtle white matter changes and the intrinsic connectivity between distinct cortical regions. This review describes the advantages of these two modalities in SVCI research and the current DTI and fMRI findings on SVCI. Using DTI technique, a variety of studies found that white matter microstructural damages in the anterior and superior areas are more specific to SVCI. Similarly, functional brain abnormalities detected by fMRI have also been mainly shown in anterior brain areas in SVCI. The characteristic distribution of brain abnormalities in SVCI interrupts the prefrontal-subcortical loop that results in cognitive impairments in particular domains, which further confirms the ‘disconnection syndrome’ hypothesis. In addition, another MRI technique, arterial spin labelling (ASL), has been used to describe the disconnection patterns in a variety of conditions by measuring cerebral blood flow. The role of the ASL technique in SVCI research is also assessed. Finally, the review proposes the application of multimodality fusion in the investigation of SVCI pathogenesis.

Development of large-scale functional networks from birth to adulthood: A guide to the neuroimaging literature

The development of human cognition results from the emergence of coordinated activity between distant brain areas. Network science, combined with non-invasive functional imaging, has generated unprecedented insights regarding the adult brain's functional organization, and promises to help elucidate the development of functional architectures supporting complex behavior. Here we review what is known about functional network development from birth until adulthood, particularly as understood through the use of resting-state functional connectivity MRI (rs-fcMRI). We attempt to synthesize rs-fcMRI findings with other functional imaging techniques, with macro-scale structural connectivity, and with knowledge regarding the development of micro-scale structure. We highlight a number of outstanding conceptual and technical barriers that need to be addressed, as well as previous developmental findings that may need to be revisited. Finally, we discuss key areas ripe for future research in order to (1) better characterize normative developmental trajectories, (2) link these trajectories to biologic mechanistic events, as well as component behaviors and (3) better understand the clinical implications and pathophysiological basis of aberrant network development.

The Rhesus Monkey Connectome Predicts Disrupted Functional Networks Resulting from Pharmacogenetic Inactivation of the Amygdala

Contemporary research suggests that the mammalian brain is a complex system, implying that damage to even a single functional area could have widespread consequences across the system. To test this hypothesis, we pharmacogenetically inactivated the rhesus monkey amygdala, a subcortical region with distributed and well-defined cortical connectivity. We then examined the impact of that perturbation on global network organization using resting-state functional connectivity MRI. Amygdala inactivation disrupted amygdalocortical communication and distributed corticocortical coupling across multiple functional brain systems. Altered coupling was explained using a graph-based analysis of experimentally established structural connectivity to simulate disconnection of the amygdala. Communication capacity via monosynaptic and polysynaptic pathways, in aggregate, largely accounted for the correlational structure of endogenous brain activity and many of the non-local changes that resulted from amygdala inactivation. These results highlight the structural basis of distributed neural activity and suggest a strategy for linking focal neuropathology to remote neurophysiological changes.

Fluctuations between high- and low-modularity topology in time-resolved functional connectivity

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.

Multilayer motif analysis of brain networks

In the last decade, network science has shed new light both on the structural (anatomical) and on the functional (correlations in the activity) connectivity among the different areas of the human brain. The analysis of brain networks has made possible to detect the central areas of a neural system and to identify its building blocks by looking at overabundant small subgraphs, known as motifs. However, network analysis of the brain has so far mainly focused on anatomical and functional networks as separate entities. The recently developed mathematical framework of multi-layer networks allows us to perform an analysis of the human brain where the structural and functional layers are considered together. In this work, we describe how to classify the subgraphs of a multiplex network, and we extend the motif analysis to networks with an arbitrary number of layers. We then extract multi-layer motifs in brain networks of healthy subjects by considering networks with two layers, anatomical and functional, respectively, obtained from diffusion and functional magnetic resonance imaging. Results indicate that subgraphs in which the presence of a physical connection between brain areas (links at the structural layer) coexists with a non-trivial positive correlation in their activities are statistically overabundant. Finally, we investigate the existence of a reinforcement mechanism between the two layers by looking at how the probability to find a link in one layer depends on the intensity of the connection in the other one. Showing that functional connectivity is non-trivially constrained by the underlying anatomical network, our work contributes to a better understanding of the interplay between the structure and function in the human brain.

Identifying topological motif patterns of human brain functional networks

Recent imaging connectome studies demonstrated that the human functional brain network follows an efficient small-world topology with cohesive functional modules and highly connected hubs. However, the functional motif patterns that represent the underlying information flow remain largely unknown. Here, we investigated motif patterns within directed human functional brain networks, which were derived from resting-state functional magnetic resonance imaging data with controlled confounding hemodynamic latencies. We found several significantly recurring motifs within the network, including the two-node reciprocal motif and five classes of three-node motifs. These recurring motifs were distributed in distinct patterns to support intra- and inter-module functional connectivity, which also promoted integration and segregation in network organization. Moreover, the significant participation of several functional hubs in the recurring motifs exhibited their critical role in global integration. Collectively, our findings highlight the basic architecture governing brain network organization and provide insight into the information flow mechanism underlying intrinsic brain activities. Hum Brain Mapp 38:2734-2750, 2017. © 2017 Wiley Periodicals, Inc.

Multiregional integration in the brain during resting-state fMRI activity

Data-driven models of functional magnetic resonance imaging (fMRI) activity can elucidate dependencies that involve the combination of multiple brain regions. Activity in some regions during resting-state fMRI can be predicted with high accuracy from the activities of other regions. However, it remains unclear in which regions activity depends on unique integration of multiple predictor regions. To address this question, sparse (parsimonious) models could serve to better determine key interregional dependencies by reducing false positives. We used resting-state fMRI data from 46 subjects, and for each region of interest (ROI) per subject we performed whole-brain recursive feature elimination (RFE) to select the minimal set of ROIs that best predicted activity in the modeled ROI. We quantified the dependence of activity on multiple predictor ROIs, by measuring the gain in prediction accuracy of models that incorporated multiple predictor ROIs compared to models that used a single predictor ROI. We identified regions that showed considerable evidence of multiregional integration and determined the key regions that contributed to their observed activity. Our models reveal fronto-parietal integration networks, little integration in primary sensory regions, as well as redundancy between some regions. Our study demonstrates the utility of whole-brain RFE to generate data-driven models with minimal sets of ROIs that predict activity with high accuracy. By determining the extent to which activity in each ROI depended on integration of signals from multiple ROIs, we find cortical integration networks during resting-state activity.

Models of fMRI activity can elucidate underlying dependencies that involve the combination of multiple brain regions. However, it remains unclear in which regions activity depends on unique integration of multiple predictor regions. To address this question, sparse (parsimonious) models could serve to better determine key interregional dependencies by reducing false positives. We used resting-state fMRI data, and for each brain region we performed whole-brain recursive feature elimination to select the minimal set of regions that best predicted activity in the region. We identified integrator regions by quantifying the gain in prediction accuracy of models that incorporated multiple predictor regions compared to single predictor region. Our study provides data-driven models that use minimal sets of regions to predict activity with high accuracy. By determining the extent to which activity in each region depended on integration of signals from multiple sources, we find cortical integration networks during resting-state activity.

Functional connectomics from a "big data" perspective

In the last decade, explosive growth regarding functional connectome studies has been observed. Accumulating knowledge has significantly contributed to our understanding of the brain's functional network architectures in health and disease. With the development of innovative neuroimaging techniques, the establishment of large brain datasets and the increasing accumulation of published findings, functional connectomic research has begun to move into the era of "big data", which generates unprecedented opportunities for discovery in brain science and simultaneously encounters various challenging issues, such as data acquisition, management and analyses. Big data on the functional connectome exhibits several critical features: high spatial and/or temporal precision, large sample sizes, long-term recording of brain activity, multidimensional biological variables (e.g., imaging, genetic, demographic, cognitive and clinic) and/or vast quantities of existing findings. We review studies regarding functional connectomics from a big data perspective, with a focus on recent methodological advances in state-of-the-art image acquisition (e.g., multiband imaging), analysis approaches and statistical strategies (e.g., graph theoretical analysis, dynamic network analysis, independent component analysis, multivariate pattern analysis and machine learning), as well as reliability and reproducibility validations. We highlight the novel findings in the application of functional connectomic big data to the exploration of the biological mechanisms of cognitive functions, normal development and aging and of neurological and psychiatric disorders. We advocate the urgent need to expand efforts directed at the methodological challenges and discuss the direction of applications in this field.

Mapping complementary features of cross-species structural connectivity to construct realistic "Virtual Brains"

Modern systems neuroscience increasingly leans on large-scale multi-lab neuroinformatics initiatives to provide necessary capacity for biologically realistic modeling of primate whole-brain activity. Here, we present a framework to assemble primate brain's biologically plausible anatomical backbone for such modeling initiatives. In this framework, structural connectivity is determined by adding complementary information from invasive macaque axonal tract tracing and non-invasive human diffusion tensor imaging. Both modalities are combined by means of available interspecies registration tools and a newly developed Bayesian probabilistic modeling approach to extract common connectivity evidence. We demonstrate how this novel framework is embedded in the whole-brain simulation platform called The Virtual Brain (TVB). Hum Brain Mapp 38:2080-2093, 2017. © 2017 Wiley Periodicals, Inc.

An Anatomical Interface between Memory and Oculomotor Systems

Visual behavior is guided by memories from prior experience and knowledge of the visual scene. The hippocampal system (HC), in particular, has been implicated in the guidance of saccades: Amnesic patients, following damage to the HC, exhibit selective deficits in their gaze patterns. However, the neural circuitry by which mnemonic representations influence the oculomotor system remains unknown. We used a data-driven, network-based approach on directed anatomical connectivity from the macaque brain to reveal an extensive set of polysnaptic pathways spanning the extrastriate, posterior parietal and prefrontal cortices that potentially mediate the exchange of information between the memory and visuo-oculomotor systems. We additionally show how the potential for directed information flow from the hippocampus to oculomotor control areas is exceptionally high. In particular, the dorsolateral pFC and FEF—regions known to be responsible for the cognitive control of saccades—are topologically well positioned to receive information from the hippocampus. Together with neuropsychological evidence of altered gaze patterns following damage to the hippocampus, our findings suggest that a reconsideration of hippocampal involvement in oculomotor guidance is needed.

Brain connectivity and auditory hallucinations: In search of novel noninvasive brain stimulation therapeutic approaches for schizophrenia

Auditory verbal hallucinations (AVH) are among the most characteristic symptoms of schizophrenia and have been linked to likely disturbances of structural and functional connectivity within frontal, temporal, parietal and subcortical networks involved in language and auditory functions. Resting-state functional magnetic resonance imaging (fMRI) has shown that alterations in the functional connectivity activity of the default-mode network (DMN) may also subtend hallucinations. Noninvasive neurostimulation techniques such as repetitive transcranial magnetic stimulation (rTMS) have the ability to modulate activity of targeted cortical sites and their associated networks, showing a high potential for modulating altered connectivity subtending schizophrenia. Notwithstanding, the clinical benefit of these approaches remains weak and variable. Further studies in the field should foster a better understanding concerning the status of networks subtending AVH and the neural impact of rTMS in relation with symptom improvement. Additionally, the identification and characterization of clinical biomarkers able to predict response to treatment would be a critical asset allowing better care for patients with schizophrenia.

Stability of Network Communities as a Function of Task Complexity

The analysis of the community architecture in functional brain networks has revealed important relations between specific behavioral patterns and characteristic features of the associated functional organization. Numerous studies have assessed changes in functional communities during different states of awareness, learning, information processing, and various behavioral patterns. The robustness of detected communities within a network has been an often-discussed topic in complex systems research. However, our knowledge regarding the intersubject stability of functional communities in the human brain while performing different tasks is still lacking. In this study, we examined the variability of functional communities in weighted undirected graphs based on fMRI recordings of healthy participants across three conditions: the resting state, syllable production as a simple vocal motor task, and meaningful speech production representing a complex behavioral pattern with cognitive involvement. On the basis of the constructed empirical networks, we simulated a large cohort of artificial graphs and performed a leave-one-out stability analysis to assess the sensitivity of communities in the group-averaged networks with respect to perturbations in the averaging cohort. We found that the stability of partitions derived from group-averaged networks depended on task complexity. The determined community architecture in mean networks reflected within-behavior network stability and between-behavior flexibility of the human functional connectome. The sensitivity of functional communities increased from rest to syllable production to speaking, which suggests that the approximation quality of the community structure in the average network to reflect individual per-participant partitions depends on task complexity.

Structural architecture supports functional organization in the human aging brain at a regionwise and network level

Functional interactions in the brain are constrained by the underlying anatomical architecture, and structural and functional networks share network features such as modularity. Accordingly, age-related changes of structural connectivity (SC) may be paralleled by changes in functional connectivity (FC). We provide a detailed qualitative and quantitative characterization of the SC-FC coupling in human aging as inferred from resting-state blood oxygen-level dependent functional magnetic resonance imaging and diffusion-weighted imaging in a sample of 47 adults with an age range of 18-82. We revealed that SC and FC decrease with age across most parts of the brain and there is a distinct age-dependency of regionwise SC-FC coupling and network-level SC-FC relations. A specific pattern of SC-FC coupling predicts age more reliably than does regionwise SC or FC alone (r = 0.73, 95% CI = [0.7093, 0.8522]). Hence, our data propose that regionwise SC-FC coupling can be used to characterize brain changes in aging. Hum Brain Mapp 37:2645-2661, 2016. © 2016 Wiley Periodicals, Inc.

Network-Level Structure-Function Relationships in Human Neocortex

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.