Functional brain networks are shaped and constrained by the underlying structural network. However, functional networks are not merely a one-to-one reflection of the structural network. Several theories have been put forward to understand the relationship between structural and functional networks. However, it remains unclear how these theories can be unified. Two existing recent theories state that 1) functional networks can be explained by all possible walks in the structural network, which we will refer to as the series expansion approach, and 2) functional networks can be explained by a weighted combination of the eigenmodes of the structural network, the so-called eigenmode approach. To elucidate the unique or common explanatory power of these approaches to estimate functional networks from the structural network, we analysed the relationship between these two existing views. Using linear algebra, we first show that the eigenmode approach can be written in terms of the series expansion approach, i.e., walks on the structural network associated with different hop counts correspond to different weightings of the eigenvectors of this network. Second, we provide explicit expressions for the coefficients for both the eigenmode and series expansion approach. These theoretical results were verified by empirical data from Diffusion Tensor Imaging (DTI) and functional Magnetic Resonance Imaging (fMRI), demonstrating a strong correlation between the mappings based on both approaches. Third, we analytically and empirically demonstrate that the fit of the eigenmode approach to measured functional data is always at least as good as the fit of the series expansion approach, and that errors in the structural data lead to large errors of the estimated coefficients for the series expansion approach. Therefore, we argue that the eigenmode approach should be preferred over the series expansion approach. Results hold for eigenmodes of the weighted adjacency matrices as well as eigenmodes of the graph Laplacian. ​Taken together, these results provide an important step towards unification of existing theories regarding the structure-function relationships in brain networks.

@en

There is an increasing awareness of the advantages of multi-modal neuroimaging. Networks obtained from different modalities are usually treated in isolation, which is however contradictory to accumulating evidence that these networks show non-trivial interdependencies. Even networks obtained from a single modality, such as frequency-band specific functional networks measured from magnetoencephalography (MEG) are often treated independently. Here, we discuss how a multilayer network framework allows for integration of multiple networks into a single network description and how graph metrics can be applied to quantify multilayer network organisation for group comparison. We analyse how well-known biases for single layer networks, such as effects of group differences in link density and/or average connectivity, influence multilayer networks, and we compare four schemes that aim to correct for such biases: the minimum spanning tree (MST), effective graph resistance cost minimisation, efficiency cost optimisation (ECO) and a normalisation scheme based on singular value decomposition (SVD). These schemes can be applied to the layers independently or to the multilayer network as a whole. For correction applied to whole multilayer networks, only the SVD showed sufficient bias correction. For correction applied to individual layers, three schemes (ECO, MST, SVD) could correct for biases. By using generative models as well as empirical MEG and functional magnetic resonance imaging (fMRI) data, we further demonstrated that all schemes were sensitive to identify network topology when the original networks were perturbed. In conclusion, uncorrected multilayer network analysis leads to biases. These biases may differ between centres and studies and could consequently lead to unreproducible results in a similar manner as for single layer networks. We therefore recommend using correction schemes prior to multilayer network analysis for group comparisons.

@en

Recent work has revealed frequency-dependent global patterns of information flow by a network analysis of magnetoencephalography data of the human brain. However, it is unknown which properties on a small subgraph-scale of those functional brain networks are dominant at different frequencies bands. Motifs are the building blocks of networks on this level and have previously been identified as important features for healthy and abnormal brain function. In this study, we present a network construction that enables us to search and analyze motifs in different frequency bands. We give evidence that the bi-directional two-hop path is the most important motif for the information flow in functional brain networks. A clustering based on this motif exposes a spatially coherent yet frequency-dependent sub-division between the posterior, occipital and frontal brain regions.@en

The interplay between structural connections and emerging information flow in the human brain remains an open research problem. A recent study observed global patterns of directional information flow in empirical data using the measure of transfer entropy. For higher frequency bands, the overall direction of information flow was from posterior to anterior regions whereas an anterior-to-posterior pattern was observed in lower frequency bands. In this study, we applied a simple Susceptible-Infected-Susceptible (SIS) epidemic spreading model on the human connectome with the aim to reveal the topological properties of the structural network that give rise to these global patterns. We found that direct structural connections induced higher transfer entropy between two brain regions and that transfer entropy decreased with increasing distance between nodes (in terms of hops in the structural network). Applying the SIS model, we were able to confirm the empirically observed opposite information flow patterns and posterior hubs in the structural network seem to play a dominant role in the network dynamics. For small time scales, when these hubs acted as strong receivers of information, the global pattern of information flow was in the posterior-to-anterior direction and in the opposite direction when they were strong senders. Our analysis suggests that these global patterns of directional information flow are the result of an unequal spatial distribution of the structural degree between posterior and anterior regions and their directions seem to be linked to different time scales of the spreading process.

@en

Over the last two decades the field of network science has been evolving fast. Many useful applications in a wide variety of disciplines have been found. The application of network science to the brain initiated the interdisciplinary field of complex brain networks. On a macroscopic level, brain regions are taken as nodes in a network. The analysis of pairwise connections between the brain regions as links has provided a new perspective on many problems. The application of network science to neuroscience data helped, for example, to identify the disruptions due to different neurological disorders when comparing healthy and abnormal brain networks. In this dissertation, we focus on the macroscopic level of brain regions and analyze their pairwise connections from a network science perspective. We address different general research questions from network science and exploit their application possibilities towards brain networks. Due to different measurement techniques, one can construct many different representations of brain networks. We thereby distinguish between the structural and functional brain network. Structural brain networks map the anatomical connections between the regions, which we could interpret as the ’streets’ of the brain. On top of these streets, we can measure the traffic with techniques like e.g. magnetoencephalography (MEG) or functional Magnetic Resonance Imaging (fMRI) resulting in so-called functional brain networks. However, the relation between the structural and the functional brain networks is still insufficiently understood. The first main research question of this dissertation focuses on the functional network layer and tries to identify the most important links and motifs of these networks. For this purpose, we propose the union of shortest path trees (USPT) as a new sampling method extracting all the shortest paths of a network (Chapter 2 and 3). After constructing the USPT, we compare the individual functional brain networks of multiple sclerosis patients and healthy controls (Chapter 2). Furthermore, we generalize this sampling method and present a new ranking of all the links based on the USPT (Chapter 3). Regarding the higher-order building blocks of the functional brain networks, we analyze the so-called information flow motifs based on MEG data from different frequency bands (Chapter 4). After researching the local properties of the functional brain networks, we analyze the influence of the underlying structural connections on the emerging information flow. Thus, the second main research question concerns the relationship between the functional and the underlying structural connectivity. Specifically, we analyze which topological properties of the structural networks drive the functional interactions. First, this question is approached in a mathematical and straightforward manner by assuming that an analytic function between the two networks exists (Chapter 5). We investigate this mapping function and its reverse by evaluating empirical individual and group-averaged multimodal data sets. A second approach towards the structure-function relationship employs a simple model of activity spread. The epidemic spreading model is applied on the human connectome to investigate the global patterns of directional information flow in brain networks (Chapter 6). The main focus here lies on the pairwise measure of transfer entropy to investigate the influence of one brain region on another. We present the results for the local and global outcomes of the dynamic spreading process aiming to identify the driving structural properties behind the observed global patterns. @en

Network science has widely studied the properties of brain networks. Recent work has observed a global back-to-front pattern of information flow for higher frequency bands in magnetoencephalography data. However, the effective connectivity at a local level remains yet to be analyzed. On a local level, the building blocks of all networks are motifs. In this study, we exploit the measure of dPTE to analyze motifs of the estimated effective connectivity networks. We find that some 3- and 4-motifs, the bidirectional two-hop path and its extended 4-node versions, are significantly overexpressed in the analyzed networks in comparison with random networks. With a recently developed motif-based clustering algorithm we separate the effective connectivity network in two main clusters which reveal its higher-order organization with a strong information flow between posterior hubs and anterior regions.

@en

Neuronal oscillations exist across a broad frequency spectrum, and are thought to provide a mechanism of interaction between spatially separated brain regions. Since ongoing mental activity necessitates the simultaneous formation of multiple networks, it seems likely that the brain employs interactions within multiple frequency bands, as well as cross-frequency coupling, to support such networks. Here, we propose a multi-layer network framework that elucidates this pan-spectral picture of network interactions. Our network consists of multiple layers (frequency-band specific networks) that influence each other via inter-layer (cross-frequency) coupling. Applying this model to MEG resting-state data and using envelope correlations as connectivity metric, we demonstrate strong dependency between within layer structure and inter-layer coupling, indicating that networks obtained in different frequency bands do not act as independent entities. More specifically, our results suggest that frequency band specific networks are characterised by a common structure seen across all layers, superimposed by layer specific connectivity, and inter-layer coupling is most strongly associated with this common mode. Finally, using a biophysical model, we demonstrate that there are two regimes of multi-layer network behaviour; one in which different layers are independent and a second in which they operate highly dependent. Results suggest that the healthy human brain operates at the transition point between these regimes, allowing for integration and segregation between layers. Overall, our observations show that a complete picture of global brain network connectivity requires integration of connectivity patterns across the full frequency spectrum.

@en

The relationship between structural and functional brain networks is still highly debated. Most previous studies have used a single functional imaging modality to analyze this relationship. In this work, we use multimodal data, from functional MRI, magnetoencephalography, and diffusion tensor imaging, and assume that there exists a mapping between the connectivity matrices of the resting-state functional and structural networks. We investigate this mapping employing group averaged as well as individual data. We indeed find a significantly high goodness of fit level for this structure-function mapping. Our analysis suggests that a functional connection is shaped by all walks up to the diameter in the structural network in both modality cases. When analyzing the inverse mapping, from function to structure, longer walks in the functional network also seem to possess minor influence on the structural connection strengths. Even though similar overall properties for the structure-function mapping are found for different functional modalities, our results indicate that the structure-function relationship is modality dependent.

@en

Objective
An important problem in systems neuroscience is the relation between complex structural and functional brain networks. Here we use simulations of a simple dynamic process based upon the susceptible–infected–susceptible (SIS) model of infection dynamics on an empirical structural brain network to investigate the extent to which the functional interactions between any two brain areas depend upon (i) the presence of a direct structural connection; and (ii) the degree product of the two areas in the structural network.
Methods
For the structural brain network, we used a 78 × 78 matrix representing known anatomical connections between brain regions at the level of the AAL atlas (Gong et al., 2009). On this structural network we simulated brain dynamics using a model derived from the study of epidemic processes on networks. Analogous to the SIS model, each vertex/brain region could be in one of two states (inactive/active) with two parameters β and δ determining the transition probabilities. First, the phase transition between the fully inactive and partially active state was investigated as a function of β and δ. Second, the statistical interdependencies between time series of node states were determined (close to and far away from the critical state) with two measures: (i) functional connectivity based upon the correlation coefficient of integrated activation time series; and (ii) effective connectivity based upon conditional co-activation at different time intervals.
Results
We find a phase transition between an inactive and a partially active state for a critical ratio τ = β/δ of the transition rates in agreement with the theory of SIS models. Slightly above the critical threshold, node activity increases with degree, also in line with epidemic theory. The functional, but not the effective connectivity matrix closely resembled the underlying structural matrix. Both functional connectivity and, to a lesser extent, effective connectivity were higher for connected as compared to disconnected (i.e.: not directly connected) nodes. Effective connectivity scaled with the degree product. For functional connectivity, a weaker scaling relation was only observed for disconnected node pairs. For random networks with the same degree distribution as the original structural network, similar patterns were seen, but the scaling exponent was significantly decreased especially for effective connectivity.
Conclusions
Even with a very simple dynamical model it can be shown that functional relations between nodes of a realistic anatomical network display clear patterns if the system is studied near the critical transition. The detailed nature of these patterns depends on the properties of the functional or effective connectivity measure that is used. While the strength of functional interactions between any two nodes clearly depends upon the presence or absence of a direct connection, this study has shown that the degree product of the nodes also plays a large role in explaining interaction strength, especially for disconnected nodes and in combination with an effective connectivity measure. The influence of degree product on node interaction strength probably reflects the presence of large numbers of indirect connections.@en

Understanding the electrophysiological basis of resting state networks (RSNs) in the human brain is a critical step towards elucidating how inter-areal connectivity supports healthy brain function. In recent years, the relationship between RSNs (typically measured using haemodynamic signals) and electrophysiology has been explored using functional Magnetic Resonance Imaging (fMRI) and magnetoencephalography (MEG). Significant progress has been made, with similar spatial structure observable in both modalities. However, there is a pressing need to understand this relationship beyond simple visual similarity of RSN patterns. Here, we introduce a mathematical model to predict fMRI-based RSNs using MEG. Our unique model, based upon a multivariate Taylor series, incorporates both phase and amplitude based MEG connectivity metrics, as well as linear and non-linear interactions within and between neural oscillations measured in multiple frequency bands. We show that including non-linear interactions, multiple frequency bands and cross-frequency terms significantly improves fMRI network prediction. This shows that fMRI connectivity is not only the result of direct electrophysiological connections, but is also driven by the overlap of connectivity profiles between separate regions. Our results indicate that a complete understanding of the electrophysiological basis of RSNs goes beyond simple frequency-specific analysis, and further exploration of non-linear and cross-frequency interactions will shed new light on distributed network connectivity, and its perturbation in pathology.

@en