This paper presents an integrative approach to investigating the functional network for spatial navigation in the human brain. This approach incorporates a large-scale neuroimaging meta-analytic database, resting-state functional magnetic resonance imaging, and network modeling and graph-theoretical techniques.
Spatial navigation is a complex function involving the integration and manipulation of multisensory information. Using different navigation tasks, many promising results have been achieved on the specific functions of various brain regions (e.g., hippocampus, entorhinal cortex, and parahippocampal place area). Recently, it has been suggested that a non-aggregate network process involving multiple interacting brain regions may better characterize the neural basis of this complex function. This paper presents an integrative approach for constructing and analyzing the functionally-specific network for spatial navigation in the human brain. Briefly, this integrative approach consists of three major steps: 1) to identify brain regions important for spatial navigation (nodes definition); 2) to estimate functional connectivity between each pair of these regions and construct the connectivity matrix (network construction); 3) to investigate the topological properties (e.g., modularity and small worldness) of the resulting network (network analysis). The presented approach, from a network perspective, could help us better understand how our brain supports flexible navigation in complex and dynamic environments, and the revealed topological properties of the network can also provide important biomarkers for guiding early identification and diagnosis of Alzheimer’s disease in clinical practice.
Functional specificity is a fundamental organization principle of the human brain, which plays a crucial role in shaping cognitive functions1. Abnormalities in the organization of functional specificity can reflect hallmark cognitive impairments and the associated pathological foundations of major brain disorders such as autism and Alzheimer's disease2,3. While conventional theories and research have tended to focus on single brain regions, such as the fusiform face area (FFA) for face recognition4 and parahippocampus place area (PPA)5 for scene processing, an increasing body of evidence suggests that complex cognitive functions, including spatial navigation and language, require coordinate activity across multiple brain regions6. Investigating the mechanisms underlying the interactions in support of complex cognitive functions is a critical scientific question that will help to shed light on the functional architecture and operation of the brain. Here, taking spatial navigation as an example, we present an integrative method for modeling the functional network for spatial navigation in the human brain.
Spatial navigation is a complex cognitive function, which involves the integration and manipulation of multiple cognitive components, such as visual-spatial coding, memory, and decision making7. With functional magnetic resonance imaging (fMRI), numerous studies have made significant advances in understanding the underlying cognitive processing and neural mechanisms. For instance, specific functions have been linked to different brain regions using various navigation tasks: scene processing is specifically associated with PPA, and transformation of navigation strategies is associated with the retrosplenial cortex (RSC)8,9. These studies provided important insights into the neural basis of spatial navigation. However, navigation is an internally dynamic and multimodal function, and the functions of single regions are not sufficient to explain large individual differences in spatial navigation10 that are commonly observed.
With the emergence of fMRI-based connectomics, researchers began to explore how some key brain regions interact with each other to support spatial navigation. For example, functional connectivity between the entorhinal and posterior cingulate cortices has been found to underpin navigation discrepancies in at-risk Alzheimer's disease11. In another study, we for the first time proposed a network approach by integrating connectome methods and almost all functionally relevant regions (nodes) for spatial navigation, and the results showed that topological properties of this network showed specific associations with navigation behaviors12. This study provides new insights into theories of how multiple brain regions interact with each other to support flexible navigation behaviors10,13.
The present work demonstrates an updated version of the integrative approach for modeling the functional network. Briefly, two updates were included: 1) While the nodes defined in the original study were identified based on an earlier and smaller database (55 studies with 2,765 activations, accessed in 2014), the present definition was based on the latest database (77 studies with 3,908 activations, accessed in 2022); 2) to increase functional homogeneity of each node, besides the original anatomical AAL (Anatomical Automatic Labeling) atlas14, we applied a new brain parcellation, which has a much finer resolution and higher functional homogeneity (see below). We expected that both updates would improve the modeling of the functional network. This updated protocol provides a detailed procedure for investigating the neural basis of spatial navigation from a network perspective and helps understand individual variations in navigation behaviors in health and disease. A similar procedure could also be used for network modeling for other cognitive constructs (e.g., language and memory).
NOTE: All the software used here is shown in the Table of Materials. The data used in this study for demonstration purposes were from the Human Connectome Project (HCP: http://www. humanconnectome.org)15. All experimental procedures were approved by the Institutional Review Board (IRB) at Washington University. Imaging data in the HCP dataset were acquired using a modified 3T Siemens Skyra scanner with a 32-channel head coil. Other image acquisition parameters are detailed in an earlier paper16. Minimal preprocessed data were downloaded for the demonstration, which had finished following preprocessing steps: gradient distortion correction, motion correction, field map preprocessing, spatial distortion correction, spatial normalization to the Montreal Neurological Institute (MNI) space, intensity normalization, and bias field removal. Resting-state fMRI data from researchers' projects can also be used.
1. Data preprocessing
Figure 1: Rs-fMRI preprocess and functional network connectivity estimation. The settings of preprocess (removing first 10 images, spatially smoothing with FWHM of 4 mm, linear temporally detrending, regressing out white matter signals, cerebrospinal fluid (CSF) signals, and head motion with 24 parameters, filtering the band of 0.01-0.1 HZ) and the static correlation with fisher' Z transformed. Abbreviations: Rs-fMRI = resting-state functional magnetic resonance imaging; FWHM = full width at half-maximum; CSF = cerebrospinal fluid. Please click here to view a larger version of this figure.
2. Network construction and analyses
NOTE: The general workflow for the construction and analyses of the navigation network are summarized into three main steps (Figure 2).
Figure 2: General workflow for the construction and analyses of the navigation network. (A) Choose navigation as the term to be searched in Neurosynth database. (B) A list of activation coordinates can be generated. (C) Run a meta-analysis using functions from the Neurosynth to get several brain maps. (D,E) By incorporating the meta-analytic map and a whole-brain parcellation atlas (AICHA), nodes (ROI) can be generated. (F) The construction of a navigation network using the resulting navigation nodes and their functional connectivity (Connectivity Estimation and Network Analysis). Abbreviations: ROI = region of interest; AICHA = atlas of intrinsic connectivity of homotopic areas. Please click here to view a larger version of this figure.
Figure 3: Network metrics analysis. This analysis defines the weighted positive networks with 10 thresholds. Calculate two global network metrics of small word and efficiency, four nodal network metrics of clustering coefficient, shortest path length, efficiency, and degree centrality. Please click here to view a larger version of this figure.
Figure 4: The calculation of average navigation networks. The averaged (functional) operation helps to calculate the average networks of all participants. Please click here to view a larger version of this figure.
The navigation networks
The present study identified 27 brain regions, which are associated with spatial navigation, by incorporating the latest meta-analysis neuroimaging database and the AICHA atlas. These regions consisted of the medial temporal and the parietal regions that have been commonly reported in navigation neuroimaging studies. The spatial distribution of these regions is shown in Figure 5A and Figure 5C. As a comparison, we also visualized an earlier definition of the spatial navigation regions in Figure 5B and Figure 5D. Twenty regions from the AAL atlas were included as a comparison. These two sets of regions showed a large overlap.
Figure 5: The modular average navigation networks. (A) Modularity of NaviNet_AICHA in REST1. (B) Modularity of NaviNet_AAL in REST1. (C,D) represent the modularity of NaviNet_AICHA and NaviNet_AAL in REST1 regressing out the whole brain signals respectively. Different node colors indicate different modules identified in each network. Two modules are shown in NaviNet_AICHA and NaviNet_AAL, which includes a medial temporal module and a parietal module. Abbreviations: AICHA = atlas of intrinsic connectivity of homotopic areas; AAL = anatomical automatic labeling. Please click here to view a larger version of this figure.
Interestingly, these two networks showed similar community distribution (Figure 5). Specifically, modularity and repeatability analyses showed two modules in both NaviNet_AICHA and NaviNet_AAL (one ventral module including the medial temporal regions and one dorsal module including the parietal regions) (Table 1). The ventral and dorsal modules were similar between NaviNet_AICHA and NaviNet_AAL, although the number of nodes was larger in the former, given the finer brain parcellation of the AICHAI atlas. These results were independent of the strategies used for dealing with global signals in the preprocessing procedure (Figure 5). In addition, similar community distributions were observed in the REST2 dataset (Supplementary Figure S1).
Similarity of the topological properties of two navigation networks
Next, we examined the similarity of each network measure between the two networks. The purpose of the similarity analyses was two-fold: (1) to evaluate the generalizability of the results when using different definition strategies and (2) to determine an optimal network threshold for the network analyses.
In general, five of the six metrics, except for the clustering coefficient, showed significant correlations between the two networks with the majority of the network sparsity thresholds used in the network analyses (Figure 6). The similarity values increased quickly with the sparsity threshold for all metrics, except for the mean node degree, which showed an excellent similarity value with either threshold. The small-world metric showed the highest similarity at a threshold between 0.30 and 0.40, where other metrics also showed the highest similarity. These results suggest that the network-level analyses could reflect stable individual differences independent of node definition choices, and that the sparsity threshold of 0.30-0.40 would result in better generalizability in navigation network analyses. See Supplementary Figure S2 for more similarities with REST2.
Figure 6: Similarity of topological properties of two networks. The results of REST 1 (A) without and (B) with regressing out whole brain signals are shown. The Pearson correlation coefficient on the y-axis indicates the similarity of topological properties of the two networks. The sparsity threshold ranges from 0.05 to 0.5. The asterisks in the figures indicate the level of significance, where *p < 0.05, ** p < 0.01, *** p < 0.001. Please click here to view a larger version of this figure.
The test-retest reliability
We also evaluated the test-retest reliability of the topological measures of the navigation networks. Various sparsity thresholds between 0.05 and 0.50 were used when calculating these network measures to remove potential spurious connectivity in the networks (see Supplementary Table S1 for details). Here, we mainly reported results with a threshold of 0.40, given the similarity results above. The majority of the network metrics showed fair to good reliability (ICC > 0.2) both in the network NaviNet_AICHA and NaviNet_AAL, while the NaviNet_AAL showed relatively higher reliability than the NaviNet_AICHA. In addition, we found that including global signal regression in fMRI data preprocessing could result in a higher reliability (Figure 7). The clustering coefficient, shortest path length, and small-world in the NaviNet_AAL network showed the highest test-retest reliability, while the clustering coefficient and small-world in NaviNet_AICHA also showed higher test-retest reliability than other measures. These results suggest that the clustering coefficient and small-world are the most reliable among these metrics.
Figure 7: Test-retest reliability of topological properties of navigation networks. (A) Test-retest reliability for data without regressing out the whole brain signals. (B) Test-retest reliability for data regressing out the whole brain signals. Abbreviations: ICC = intra-class correlation coefficient; Cc = clustering coefficient; Lp = shortest path length; Sw = small-worldness; Nd = mean nodal degree; Eg = global efficiency; Eloc = local efficiency; AICHA = atlas of intrinsic connectivity of homotopic areas; AAL = anatomical automatic labeling; cNGS = data without regressing out the whole brain signals; cWGS = regressing out the whole brain signals. Please click here to view a larger version of this figure.
module numbers | |||||
2 | 3 | 4 | 5 | ||
NaviNet_AICHA | cNGS | 1 | 1 | 0.96 | 0.67 |
cWGS | 1 | 0.96 | 0.78 | 0.89 | |
NaviNet_AAL | cNGS | 1 | 0.95 | 0.95 | 0.65 |
cWGS | 1 | 0.95 | 0.95 | 0.95 |
Table 1: The repeatability of module partition between REST 1 and REST 2. The first row indicates the number of modules. cNGS represents the rs-fMRI without regressing out the global signals and cWGS represents the rs-fMRI regressed out the global signals. A larger number indicates a higher repeatability, and two modules are chosen for NaviNet_AICHA and NaviNet_AAL in the present text. Abbreviations: AICHA = atlas of intrinsic connectivity of homotopic areas; AAL = anatomical automatic labeling.
Supplementary Figure S1: The modular average navigation networks in REST 2. (A,B) The modularity of NaviNet_AICHA and NaviNet_AAL in the REST 2 dataset without regressing out the whole brain signals. (C,D) The modularity of NaviNet_AICHA and NaviNet_AAL in the REST 2 dataset with the whole brain signals regression. Different node colors indicate different modules identified in each network. Both two modules are shown in NaviNet_AICHA and NaviNet_AAL. Please click here to download this File.
Supplementary Figure S2: Similarity of topological properties of two networks in REST 2. The results of REST 2 without/with regressing out whole brain signals are shown (A and B, respectively). The Pearson correlation coefficient on the y-axis indicates the similarity of topological properties of the two networks. The sparsity threshold ranges from 0.05 to 0.5. The asterisks in the figures indicate the level of significance, where * p < 0.05, ** p < 0.01, *** p < 0.001. Please click here to download this File.
Supplementary Table S1: Test-retest reliability of topological properties of navigation networks with different sparsity thresholds. The values indicate the intra-class correlation coefficient of NaviNet_AICHA and NaviNet_AAL with different sparsity thresholds. Abbreviations: Cc = clustering coefficient; Lp = shortest path length; Sw = small woldness; Nd = mean nodal degree; Eg = global efficiency; Eloc = local efficiency; cNGS = data without regressing out the whole brain signals; cWGS = regressing out the whole brain signals. Please click here to download this File.
Network neuroscience is expected to help in understanding how the brain network supports human cognitive functions32. This protocol demonstrates an integrative approach to studying the functional network for spatial navigation in the human brain, which can also inspire network modeling for other cognitive constructs (e.g., language).
This approach consisted of three main steps: node definition, network construction, and network analysis. While network construction and network analysis are the same as those in general network studies of the whole brain, node definition is the most critical step of this protocol. This step makes use of a large-scale meta-analysis of functional activation related to spatial navigation to localize the most important brain regions for navigation behaviors. Thus, we can model the functionally meaningful network, which helps understand the neural basis of the complex processing from a network perspective. Note that prefrontal regions were missing in the node definition results while an increasing number of navigation studies have suggested critical roles of these regions33. This could be due to the lack of activations of these regions in navigation-related studies within the database, which resulted in limited data for the meta-analysis. When more data are available for localizing these navigation-related prefrontal regions, it would be an interesting question to investigate their roles in the navigation network in future studies. Researchers could also apply this protocol to study other cognitive functions when localizing individual brain regions is possible. Field knowledge is needed for identifying regions of interest to ensure close associations with the specific function.
In this protocol, we focused on the spatial navigation network and showed a high cover of various brain regions reported in spatial navigation studies. Given the absence of a universally agreed-upon definition of brain regions that support spatial navigation, the demonstration used two sets of regions, one generated by incorporating the largest meta-analysis and the AICHA atlas, and the other with the AAL atlas. The network topological properties based on the two definitions generally showed high similarity, which supports the effectiveness of the functional-specific network modeling to a certain extent.
We noted that the similarity strength increased with the sparsity thresholds used in the network analyses, and the results suggested that a sparsity threshold of 0.30-0.40 would be a proper choice as all the network metrics showed the highest similarity with these thresholds. With such thresholds, the network metrics also showed fair to good test-retest reliability particularly for the shortest path length and small-worldness in the case where global signal regression was included in the data preprocessing. These results largely support the use of these metrics in studies of individual differences and related brain disorders.
Due to a lack of proper behavioral data, we could not present behavioral correlates of the network metrics with spatial navigation in this protocol. Based on a few previous studies on brain-behavior associations of functional connectivity metrics of navigation-related regions12,34, we expected that the network modeling with this protocol would show a specific association with spatial navigation. Large-scale samples with both neuroimaging and behavioral data are needed to further investigate these associations. In addition, while the test-retest reliability results were not very high, the strength was comparable with those reported by previous fMRI studies35.
Future studies can apply this protocol to better understand the neural basis of spatial navigation from a network perspective and explore its variations in humans. For instance, researchers can make use of this protocol to investigate the development and aging trajectory of navigation networks, and in clinical practice, the network properties provide important biomarkers for guiding early identification and diagnosis of brain disorders such as Alzheimer's disease. Besides, future studies could also apply a similar protocol to construct the network models for other cognitive constructs.
The authors have nothing to disclose.
Xiang-Zhen Kong was supported by the National Natural Science Foundation of China (32171031), STI 2030 – Major Project (2021ZD0200409), Fundamental Research Funds for the Central Universities (2021XZZX006), and Information Technology Center of Zhejiang University.
Brain connectivity toolbox (BCT) | Mikail Rubinov & Olaf Sporns | 2019 | The Brain Connectivity Toolbox (brain-connectivity-toolbox.net) is a MATLAB toolbox for complex-network (graph) analysis of structural and functional brain-connectivity data sets. |
GRETNA | Jinhui Wang et al. | 2 | GRETNA is a graph theoretical network analysis toolbox which allows researchers to perform comprehensive analysis on the topology of brain connectome by integrating the most of network measures studied in current neuroscience field. |
MATLAB | MathWorks | 2021a | MATLAB is a programming and numeric computing platform used by millions of engineers and scientists to analyze data, develop algorithms, and create models. |
Python | Guido van Rossum et al. | 3.8.6 | Python is a programming language that lets you work more quickly and integrate your systems more effectively. |
Statistical Parametric Mapping (SPM) | Karl Friston et.al | 12 | Statistical Parametric Mapping refers to the construction and assessment of spatially extended statistical processes used to test hypotheses about functional imaging data. |