Next Article in Journal
The Affect Misattribution in the Interpretation of Ambiguous Stimuli in Terms of Warmth vs. Competence: Behavioral Phenomenon and Its Neural Correlates
Next Article in Special Issue
Effect of Tennis Expertise on Motion-in-Depth Perception at Different Speeds: An Event-Related Potential Study
Previous Article in Journal
Hyperprolactinemia Associated with Attentional Processing and Interference Control Impairments in Patients with Prolactinomas
Previous Article in Special Issue
Research on Top Archer’s EEG Microstates and Source Analysis in Different States
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Functional Connectivity as an Index of Brain Changes Following a Unicycle Intervention: A Graph-Theoretical Network Analysis

Institute of Psychology, University of Graz, 8010 Graz, Austria
*
Author to whom correspondence should be addressed.
Brain Sci. 2022, 12(8), 1092; https://doi.org/10.3390/brainsci12081092
Submission received: 8 July 2022 / Revised: 24 July 2022 / Accepted: 2 August 2022 / Published: 17 August 2022
(This article belongs to the Special Issue Physical Exercise-Driven Brain Plasticity)

Abstract

:
Grey matter volume reductions in the right superior temporal gyrus (rSTG) were observed in young adults who learned to ride a unicycle. As these decreases were correlated with the acquired ability in unicycling, the authors interpreted the change as a brain tissue reorganization to increase postural control’s automated and efficient coordination. The current study aims to further corroborate this interpretation by looking at changes in the functional brain network in the very same sample of participants. For this reason, we applied graph theory, a mathematics field used to study network structure functionality. Four global and two local graph-theoretical parameters were calculated to measure whole brain and rSTG specific changes in functional network activity following the three-week-unicycle training. Findings revealed that the Local Efficiency of the rSTG was significantly elevated after the intervention indicating an increase in fault tolerance of the rSTG, possibly reflecting decentralisation of rSTG specific functions to neighbouring nodes. Thus, the increased Local Efficiency may indicate increased task efficiency by decentralising the processing of functions crucial for balance.

1. Introduction

Recent research indicates that balance interventions result in various changes in structural and functional characteristics of the brain. This has mainly been shown in research dealing with the ageing brain. For example, the observation that balance-related abilities decrease with advancing age leading to increased risks of injury from falling, has been primarily caused by age-related brain changes [1]. With advancing age, balance control is pushed into a more cortically focused, attention-driven control [1,2,3,4].
This coincides with the fact that older people exhibit more brain activation even when just imagining standing or walking [5] and during an actual balance task [2]. Balance interventions can help older adults to improve their balance performance and decrease their risk of injury. Moreover, a decrease in brain activation is also associated with changes in balance tasks [6], indicating a direct relationship between brain activation and behavioural changes.
While showing more activation when performing balance tasks, the activation decreases in older adults (but not in younger adults) when increasing the cognitive load while balancing. This was shown by instructing participants to simultaneously perform a balance game and a serial subtraction task. The balancing game adapted to individual skills, and older adults performed an easier subtraction task to equalize the relative cognitive load. Still, the decrease in neuronal activation when performing both tasks at the same time was only shown in older adults [7]. The authors interpret this as age-related deterioration causing downgrading of neural activity.
Sehm and co-workers [8] found that a balance intervention caused grey matter changes in the hippocampus in healthy patients. In addition, performance improvements in participants who have Parkinson’s disease were correlated with grey matter changes in the bilateral anterior cingulate cortex, left inferior parietal cortex, right anterior precuneus, left ventral premotor cortex and left middle temporal gyrus. Moreover, the association between memory and balance interventions in healthy participants was later shown to persist in the behavioural dimension, in increased memory performance following a balance intervention [9]. The intervention had participants perform multiple different balance tasks twice a week for a total of 12 weeks. The training led to increases in balance performance as well as in memory and spatial cognition.
In [10] participants had to learn to ride a unicycle over a period of three weeks. Prior to and after the intervention structural and functional brain scans were obtained. The structural brain findings revealed a decrease in grey matter volume in the right superior temporal gyrus (rSTG) and a smaller cluster of reduction in the left parahippocampus. The rSTG is involved in different visuospatial processes such as spatial awareness, body positioning, the space-dependent motion of one’s body and more [11,12,13,14]. Similarly, the parahippocampus is associated with visuospatial processing along with cognitive processes [15]. Five weeks later, an additional follow-up scan showed a re-increase in the grey matter of the rSTG. The amount of grey matter reduction was also correlated with the unicycle-riding performance.
Additionally, this study also found an increase in fractional anisotropy in the right forceps major of the corpus callosum and the right corticospinal tract, along with an increase in the cortical thickness of the left superior part of the left precentral gyrus between the pre-test and the follow-up assessment [10].
The decrease in grey matter volume in the rSTG is somehow puzzling and difficult to interpret since, intuitively, an increase is commonly more likely associated with performance gains. However, as grey matter changes can be caused by several different underlying neuronal changes [16], the interpretation of the underlying cellular mechanisms resulting in the rSTG’s grey matter volume reduction is challenging. Weber et al. [10] preliminarily concluded that the observed GM decreases could probably hint at a reorganization of the brain tissue to support more automated and efficient coordination of postural control.
To further investigate this interpretation, we reanalysed the [10] dataset using graph-theoretical parameters of the functional resting-state scan. Graph-theoretic parameters quantify the qualities of networks on a local or global scale. Among other things, they can capture: (a) how efficiently a network is structured for communication (efficiency); (b) to what extent it is split into modules (modularity); as well as (c) the importance of a node to the function of the whole network (centrality). Using two local (rSTG) and four global (whole brain) parameters, we dissected regional changes coinciding with the grey matter reduction in the rSTG, while also investigating the effects of the unicycle intervention on the functional network spanning the whole brain.

2. Materials and Methods

2.1. Participants

The data used for the present graph-theoretical analyses were obtained in the study by [10]. A total of 23 right-handed participants (13 female) aged from 20 to 51 (M = 30.42, SD = 9.14) learned to ride a unicycle over three weeks. None of them had any prior experience in unicycle riding, and none had any prior diagnosed medical, neurological, or psychiatric illness. They also signed an informed consent agreement prior to their first MRI session.

2.2. Study Design and Procedure

A pre-post design was used to assess brain changes caused by learning unicycling. Before their first unicycle training session, all participants underwent an MRI session to assess structural data as well as functional resting-state data. The post-scan was obtained following a three-week unicycle training intervention. The scans were acquired as close to the last training session as possible. At the end of the study, participants were assessed regarding their unicycling proficiency by independent experts. For additional information about the dataset and study design, please see [10].

2.3. MRI Data Acquisition

MRI data were acquired using a 3T Magnetom Skyra scanner (Siemens Healthineers, Erlangen, Germany) using a 32-channel head coil. For structural images a high resolution (0.7 mm isotropic) T1-weighted MPRage-Sequence was used (TR = 2400 ms, TI = 1000 ms, TE = 2.32 ms, matrix = 320 × 320, FOV = 224 mm, 192 slices, thickness 0.7 mm, no gap, no PAT, FA = 8°) [17]. For functional data acquisition, 260 volumes (6 min 4 s) were recorded while participants were in a passive resting state. This was achieved using a T2-weighted multiband SMS-EPI-Sequence (Voxelsize = 2.5 mm isotropic, TR = 1.4 s, TE = 35 ms, matrix = 90 × 90, FoV = 225, FA = 65°).

2.4. Preprocessing

Results included in this manuscript come from preprocessing performed using fMRIPrep 20.1.1 (RRID:SCR_016216) [18,19], which is based on Nipype 1.5.0 (RRID:SCR_002502) [20,21]. Many internal operations of fMRIPrep use Nilearn 0.6.2 (RRID:SCR_001362) [22], mainly within the functional processing workflow.

2.4.1. Anatomical Data Preprocessing

The T1-weighted images were corrected for intensity non-uniformity (INU) with N4BiasFieldCorrection [23], distributed with ANTs 2.2.0 (RRID:SCR_004757) [24]. After skull-stripping, brain tissue segmentation of cerebrospinal fluid (CSF), white matter (WM) and grey matter (GM) was performed on the brain-extracted T1w using fast (FSL 5.0.9, RRID:SCR_002823) [25]. Volume-based spatial normalization to the MNI152NLin2009cAsym standard spaces was performed through nonlinear registration with antsRegistration (ANTs 2.2.0).

2.4.2. Functional Data Preprocessing

First, a reference volume and its skull-stripped version were generated using a custom methodology of fMRIPrep. Head-motion parameters with respect to the BOLD reference are estimated before any spatiotemporal filtering using mcflirt (FSL 5.0.9) [26]. BOLD runs were slice-time corrected using 3dTshift from AFNI 20160207 (RRID:SCR_005927) [27]. A deformation field to correct for susceptibility distortions was estimated based on fMRIPrep’s fieldmap-less approach. A corrected EPI (echo-planar imaging) reference was calculated for a more accurate co-registration with the anatomical reference based on the estimated susceptibility distortion. The BOLD reference was then co-registered to the T1w reference using bbregister [28]. Co-registration was configured with six degrees of freedom. Finally, the BOLD time series were resampled into the MNI152NLin2009cAsym standard space. Several confounding time series were calculated based on the preprocessed BOLD, including framewise displacement (FD) and three region-wise global signals. FD was computed using two formulations following Power (absolute sum of relative motions, [29]) and Jenkinson (relative root mean square displacement between affines, [26]). The three global signals were extracted within the CSF, the WM, and the whole-brain masks. Six motion parameters, framewise displacement, and the mean signal from CSF were regressed out of the data. Finally, smoothing with an FWHM of 2mm was performed on the data before network construction.

2.5. Graph Theoretical Analysis

Network construction and parameter calculation were performed using the Matlab Toolbox GRETNA v2.0 [30]. Here, we present a brief overview of our chosen nodes, edges and thresholds.

2.5.1. Node Definition

Defining the nodes was achieved by partitioning the brain into 116 separate regions using the Automated Automatical Labeling atlas (AAL-116) [31]. The rSTG, which we were trying to analyze (see section Graph Theoretical Parameters), is one clearly defined node within this atlas. This enabled the possibility of qualitatively comparing structural changes in the rSTG with changes in the node’s graph parameters covering the rSTG.

2.5.2. Edges

Edges were calculated using Pearson correlations. Correlating the mean time-courses of the fMRI-data of every single node with every other node established the edge weights between all 116 nodes. The weights were not binarized, and only positive edge weights were considered leading to a weighted undirected network.

2.5.3. Thresholds

When constructing a functional network, noise and especially the method selection for edge calculation can overestimate the number of edges [32]. To take this into account, we used proportional thresholding. To this end, a percentage defines how many edges are used when constructing the network. This is done based on edge weights, meaning all edges that fall in the defined percentage of highest weights are used in the network, and all others are set to zero [33]. However, this leads to the additional problem of choosing an appropriate percentage.
Because a percentage can be chosen arbitrarily, the validity of the functional network parameter potentially values suffers. To counteract that, we calculated every participant’s network and corresponding parameters multiple times using multiple thresholds before determining the area under the curve (AUC) of the individual parameters and comparing group averages [33]. Thus, we studied a wide range of thresholds while not inflating the multiple comparison problem. For this study, 18 thresholds ranged from the biggest 6% of edges to the biggest 40% of edges (2% steps).

2.5.4. Graph Theoretical Parameters

Graph theoretical parameters can be categorized into the network concepts they are measuring. We chose network concepts we found to be the most important and studied parameter of each category [34]. The global concepts we chose were network integration, network segregation, small-world characteristic and centrality. These concepts were measured by the respective parameters: Global Efficiency, Modularity, Small-Worldness and Hubness. Additionally, we looked at two local parameters representing node integration and node segregation, respectively, to analyse the right superior temporal gyrus (rSTG). Those were measured by Nodal Efficiency and Local Efficiency, respectively. We did not include parameters from all categories because the small-world property is not measurable on a local scale, and local centrality measures are highly correlated with local integration measures. The rSTG was used as the goal region for local analysis because there are no systematic studies on how changes in local functional graph-theoretic parameters correspond to changes in grey matter. As [10] have previously shown a decrease in grey matter volume in the rSTG of the same dataset, this analysis may help explain those finding by looking at the corresponding functional changes in parameters. Global Efficiency, Modularity, Small-Worldness, Nodal Efficiency and Local Efficiency were calculated as explained in [34].

Hubness

Hubness can be calculated in multiple ways [35]. Here, we defined the Degree Centrality value for every node in every threshold from each participant. Then we combined every single value from both pre and post-intervention data and calculated the mean and SD of the distribution. Lastly, we defined Hubs as nodes with a Degree Centrality value of at least one SD over the mean [36]. These categorizations were then used to compare whether the number of hubs changed between pre and post-intervention. Thus, this did not require the use of AUCs.

2.6. Statistical Analysis

Using a dependent t-test with the R package “rstatix” [37] all six parameters, Global Efficiency, Modularity, Small-Worldness, Hubness, Nodal Efficiency and Local Efficiency, were tested for group differences. Additionally, the AUCs for the Modularity parameter were calculated using the R package “DescTools” [38] because the deployed Matlab toolbox GRETNA did not output AUCs for the Modularity parameter.

2.7. Stability of Findings

Because there are no transparent best practices for atlas and threshold choices in functional network construction, we exploratively calculated the same parameters using different atlases and thresholds. Our focus was on exploring the stability of global parameters based on atlas size and network construction method. Thus, we used two spatially more refined graded atlases that were conventionalized based on functional data. However, this makes a reliable comparison of local parameter values between atlases impossible because there is no node fully representing the rSTG. Further, Nodes that are partially represented by voxels in the rSTG are not bound to have all their respective voxels in close vicinity of the rSTG. To explore whether a node from a functional atlas that is close to our original rSTG node has enough crossover to replicate the findings, we chose to use the node closest to our original rSTG node, based on MNI coordinates, as a representative.
For the alternative atlases, we used the Dosenbach-160 [39] and the Power-264 [40], dividing the brain into 160 and 264 segments, respectively. As for the thresholds, we used the original 18 thresholds and used the lower half (6% to 22% with 2% steps) and the higher half (24% to 40% with 2% steps) individually for all three atlases.

3. Results

3.1. Descriptive Statistics

First, we looked at the Small-Worldness measure of each constructed network to validate whether our choice in thresholds led to realistic network structures. Networks with a Small-Worldness value higher than one are considered networks with small-world properties [35]. In our analysis, only ten (1.2%) of the 828 constructed networks (23 persons × 2-time points × 18 thresholds) did not achieve a value of at least 1, showing that the chosen network construction parameters led to an overall reliable network structure.

3.2. Intervention Effects

All subsequently presented effects are calculated from parameters based on the AAL atlas network [31] and are Bonferroni corrected if not noted otherwise. Group differences are also graphed in Figure 1, visualizing mean parameter values and standard deviation across thresholds for every parameter except Hubness. The curves are a visual representation of the calculated AUCs used for group comparison. They show that some differences in pre- and post-parameter values are dependent on the used threshold while others seem to be stable over the whole threshold range.

3.2.1. Local Effects

Learning to ride a unicycle caused a significant change in Local Efficiency in the rSTG (t(22) = 3.51, p = 0.002, d = 0.886). However, Nodal Efficiency did not show a significant difference between pre- and post-intervention (t(22) = 0.56, p = 0.582, d = 0.361) (Figure 1).

3.2.2. Global Effects

After correcting for multiple comparisons, Global Efficiency did not significantly change between pre- and post-intervention (t(22) = 1.24, p = 0.226, d = 0.461). It is worth mentioning that without correction, the change in Global Efficiency was significant (tuncorrected(22) = 2.21, puncorrected = 0.038). Changes in Modularity (t(22) = 0.00, p = 1.000, d = −0.122), Small-Worldness (t(22) = 0.00, p = 1.000, d = −0.176) and Hubness (t(22) = 0.82, p = 0.422, d = 0.397) were not significant either (Figure 1).

3.3. Stability of Findings

The stability of the findings is qualitatively expressed by the difference in t-values a parameter has across different thresholds and atlases. As seen in Table 1, most of the global parameters (Global Efficiency, Small-Worldness and Modularity) showed relatively stable results across different thresholds and moderately stable results across atlases. For example, Global Efficiency t-values across AAL-116 thresholds lay between 1.93 and 2.28. As for t-values of different atlases using the total threshold, the values were between 1.64 and 2.21, which we classified as close given the big differences in atlas construction.
However, Hubness and the two local parameters (Local Efficiency and Nodal Efficiency) show large differences between atlases and, in the case of Hubness, the thresholds (AAL-116, Hubness t-value range: 0.53–2.02). See Discussion for further information and Table 1 for a list of all results.

4. Discussion

This study is among the first to look at changes in functional brain networks following a complex balance intervention. We have shown that local network changes coincide with the reported grey matter changes in the rSTG [10]. Specifically, the Local Efficiency was higher after the intervention, indicating that the rSTG’s neighbourhood nodes (nodes of the network directly connected to the rSTG) increased their efficiency in communicating with each other. This can also be interpreted as an increase in fault tolerance because it indicates how efficiently the neighbourhood of the rSTG can communicate in the case that the rSTG is disrupted [42,43]. Interestingly, none of the global parameters significantly changed following the intervention.
The presented findings show parallels to a study by [44], which investigated the effects of aerobic training on the brain network. Similar to this study, the training significantly modulated only local parameters. These findings could indicate that local parameters are better suited to identify exercise-driven brain changes, possible because these changes are primarily driven by functional specialization [45]; of course, it may also well be the case that a dataset involving larger samples of participants is needed to detect more global brain changes.

4.1. Implications

First and foremost, combining this study’s findings with the results reported in [10] gives us new insights into the changes that may occur in the rSTG when learning to ride a unicycle and possibly following balance interventions as a whole. Of course, this study does not answer which micro-scale changes occur when we measure changes in grey matter. There are multiple candidates; most likely, there may be combinations of all types of cellular changes [16]. However, our findings indicate that the reduced grey matter volume may be caused, or at least be associated with increases in the interconnectedness of topologically neighbouring nodes, which [43] was also interpreted as a local increase in fault tolerance of a node. This interpretation may provide insight into why the grey matter decreased.
The heightened interconnectedness in the resting state condition may indicate the decentralisation of functions prior mainly managed by the rSTG, such as spatial awareness. It is not uncommon for structurally different parts to be able to fulfil the same function [46,47]. Thus, a higher priority of the rSTG functions stimulated through unicycle training may have led to the decentralisation, indicated by the increase in interconnectedness, which in turn led to a decrease in the grey matter of the rSTG. While only speculative, testing this hypothesis is possible [47] and could significantly increase our understanding of neuroplasticity and brain network architecture.
An additional implication is that graph-theoretic network analysis is a viable tool to study the effects balance interventions have on the brain. This study also supports the promising notion that graph-theoretical parameters may be used as markers of optimal intervention choice [48].
Future studies should focus on using a mixture of neuroimaging tools and transcranial magnetic stimulation (TMS) to validate our interpretation of the functional changes [47]. The utilization of alternative balance training modalities will also show whether the findings are unicycle training specific or generalisable to all balance interventions.

4.2. Stability of Findings

The comparison of results between the different atlas and threshold choices is not standard practice, but given the lack of consensus in choosing these parameters, we made great efforts to demonstrate the reliability of our findings using this method. Papers that looked at how different analytic choices affect graph-theoretical parameters primarily focused on preprocessing [49] or structural networks [50]. That being said, previous studies took a similar route to test the reliability of their findings [51].
The findings showed that Global Efficiency, Small-Worldness and Modularity measures were relatively stable across different atlases and thresholds, while Hubness and the local parameter measures were more affected by the specific choices. As mentioned in the Methods section, this was expected for the local measurements because the nodes used to represent the rSTG do not actually cover the rSTG anatomically, as the alternative atlases were not constructed based on anatomical data. Therefore, this not only shows the importance of node choice in local network analyses but also validates our choice to use an anatomical atlas for the analysis.
The Hubness is the only parameter that shows variability in the results depending on atlas and threshold choices. This is most likely caused by the complex manner the hubs are calculated. For example, calculating a cutoff point for what defines a hub is especially affected by the number of regions the atlas has and the number of thresholds used. This increases the discrepancy between the calculated thresholds far more than for the other global parameters relying on AUCs. The question of how to define hubs reliably is not a new problem [36,52], but this study underlines the importance of cautious calculation and interpretation of results when working with Hubness.

4.3. Methodological Limitations

This study is mainly limited by the lack of methodological consensus. The most important is that atlases based on anatomical data are not considered optimal for network construction. Even more, pre-conceptualized atlases themselves are in some ways inferior to brain parcellations that were specifically produced for the respective sample [53,54]. On the other hand, we hypothesized about structural findings in the rSTG, making node definition only possible by using an anatomical parcellation. Thus, choosing a structural atlas was not avoidable.
There is also much room for discussion that there may not be an optimal atlas, optimal parameter values and optimal preprocessing steps for network construction in the first place [55,56]. Generally, this does not invalidate the findings but it emphasizes the importance of further research seeking more reliable scientific standards for brain network analyses and, respectively, graph theoretical network analysis.

Author Contributions

Conceptualization, U.R. and K.K.; methodology, U.R.; validation, A.F. and K.K.; formal analysis, U.R.; investigation, B.W. and K.K.; resources, B.W. and K.K.; data curation, K.K.; writing—original draft preparation, U.R.; writing—review and editing, A.F., B.W. and K.K.; visualization, U.R.; supervision, K.K.; project administration, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

Open Access Funding by the University of Graz.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee at the University of Graz (GZ: 39/76/63 ex 2015/16).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data will be available upon request from the corresponding author.

Acknowledgments

The authors acknowledge the financial support by the University of Graz.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AUCarea under the curve
CSFcerebrospinal fluid
GMgrey matter
WMwhite matter
FDframewise displacement
rSTGright Superior Temporal Gyrus

References

  1. Seidler, R.D.; Bernard, J.A.; Burutolu, T.B.; Fling, B.W.; Gordon, M.T.; Gwin, J.T.; Kwak, Y.; Lipps, D.B. Motor control and aging: Links to age-related brain structural, functional, and biochemical effects. Neurosci. Biobehav. Rev. 2010, 34, 721–733. [Google Scholar] [CrossRef] [PubMed]
  2. Baudry, S. Aging Changes the Contribution of Spinal and Corticospinal Pathways to Control Balance. Exerc. Sport Sci. Rev. 2016, 44, 104–109. [Google Scholar] [CrossRef] [PubMed]
  3. Boisgontier, M.P.; Beets, I.A.; Duysens, J.; Nieuwboer, A.; Krampe, R.T.; Swinnen, S.P. Age-related differences in attentional cost associated with postural dual tasks: Increased recruitment of generic cognitive resources in older adults. Neurosci. Biobehav. Rev. 2013, 37, 1824–1837. [Google Scholar] [CrossRef] [PubMed]
  4. Papegaaij, S.; Taube, W.; Baudry, S.; Otten, E.; Hortobágyi, T. Aging causes a reorganization of cortical and spinal control of posture. Front. Aging Neurosci. 2014, 6, 28. [Google Scholar] [CrossRef] [PubMed]
  5. Allali, G.; van der Meulen, M.; Beauchet, O.; Rieger, S.W.; Vuilleumier, P.; Assal, F. The Neural Basis of Age-Related Changes in Motor Imagery of Gait: An fMRI Study. J. Gerontol. Ser. A 2013, 69, 1389–1398. [Google Scholar] [CrossRef]
  6. Ruffieux, J.; Mouthon, A.; Keller, M.; Mouthon, M.; Annoni, J.M.; Taube, W. Balance Training Reduces Brain Activity during Motor Simulation of a Challenging Balance Task in Older Adults: An fMRI Study. Front. Behav. Neurosci. 2018, 12, 10. [Google Scholar] [CrossRef]
  7. De Rond, V.; Orcioli-Silva, D.; Dijkstra, B.W.; de Xivry, J.J.O.; Pantall, A.; Nieuwboer, A. Compromised Brain Activity With Age During a Game-Like Dynamic Balance Task: Single- vs. Dual-Task Performance. Front. Aging Neurosci. 2021, 13, 657308. [Google Scholar] [CrossRef]
  8. Sehm, B.; Taubert, M.; Conde, V.; Weise, D.; Classen, J.; Dukart, J.; Draganski, B.; Villringer, A.; Ragert, P. Structural brain plasticity in Parkinson’s disease induced by balance training. Neurobiol. Aging 2014, 35, 232–239. [Google Scholar] [CrossRef]
  9. Rogge, A.K.; Röder, B.; Zech, A.; Nagel, V.; Hollander, K.; Braumann, K.M.; Hötting, K. Balance training improves memory and spatial cognition in healthy adults. Sci. Rep. 2017, 7, 5661. [Google Scholar] [CrossRef]
  10. Weber, B.; Koschutnig, K.; Schwerdtfeger, A.; Rominger, C.; Papousek, I.; Weiss, E.M.; Tilp, M.; Fink, A. Learning Unicycling Evokes Manifold Changes in Gray and White Matter Networks Related to Motor and Cognitive Functions. Sci. Rep. 2019, 9, 4324. [Google Scholar] [CrossRef]
  11. Baier, B.; Suchan, J.; Karnath, H.O.; Dieterich, M. Neural correlates of disturbed perception of verticality. Neurology 2012, 78, 728–735. [Google Scholar] [CrossRef] [PubMed]
  12. Ellison, A. An exploration of the role of the superior temporal gyrus in visual search and spatial perception using TMS. Brain 2004, 127, 2307–2315. [Google Scholar] [CrossRef] [PubMed]
  13. Karnath, H.O.; Ferber, S.; Himmelbach, M. Spatial awareness is a function of the temporal not the posterior parietal lobe. Nature 2001, 411, 950–953. [Google Scholar] [CrossRef] [PubMed]
  14. Karnath, H.O. New insights into the functions of the superior temporal cortex. Nat. Rev. Neurosci. 2001, 2, 568–576. [Google Scholar] [CrossRef] [PubMed]
  15. Aminoff, E.M.; Kveraga, K.; Bar, M. The role of the parahippocampal cortex in cognition. Trends Cogn. Sci. 2013, 17, 379–390. [Google Scholar] [CrossRef]
  16. Zatorre, R.J.; Fields, R.D.; Johansen-Berg, H. Plasticity in gray and white: Neuroimaging changes in brain structure during learning. Nat. Neurosci. 2012, 15, 528–536. [Google Scholar] [CrossRef]
  17. Glasser, M.F.; Sotiropoulos, S.N.; Wilson, J.A.; Coalson, T.S.; Fischl, B.; Andersson, J.L.; Xu, J.; Jbabdi, S.; Webster, M.; Polimeni, J.R.; et al. The minimal preprocessing pipelines for the Human Connectome Project. NeuroImage 2013, 80, 105–124. [Google Scholar] [CrossRef]
  18. Esteban, O.; Markiewicz, C.J.; Goncalves, M.; DuPre, E.; Kent, J.D.; Salo, T.; Ciric, R.; Pinsard, B.; Blair, R.W.; Poldrack, R.A.; et al. fMRIPrep: A Robust Preprocessing Pipeline for Functional MRI. 2018. Available online: https://doi.org/10.5281/ZENODO.852659 (accessed on 2 June 2020).
  19. Esteban, O.; Markiewicz, C.J.; Blair, R.W.; Moodie, C.A.; Isik, A.I.; Erramuzpe, A.; Kent, J.D.; Goncalves, M.; DuPre, E.; Snyder, M.; et al. fMRIPrep: A robust preprocessing pipeline for functional MRI. Nat. Methods 2018, 16, 111–116. [Google Scholar] [CrossRef]
  20. Gorgolewski, K.; Burns, C.D.; Madison, C.; Clark, D.; Halchenko, Y.O.; Waskom, M.L.; Ghosh, S.S. Nipype: A Flexible, Lightweight and Extensible Neuroimaging Data Processing Framework in Python. Front. Neuroinformatics 2011, 5, 13. [Google Scholar] [CrossRef]
  21. Gorgolewski, K.J.; Esteban, O.; Markiewicz, C.J.; Burns, C.; Goncalves, M.; Jarecka, D.; Ziegler, E.; Berleant, S.; Ellis, D.G.; Pinsard, B.; et al. nipy/nipype: 1.5.1. 2018. Available online: https://doi.org/10.5281/ZENODO.596855 (accessed on 2 June 2020).
  22. Abraham, A.; Pedregosa, F.; Eickenberg, M.; Gervais, P.; Mueller, A.; Kossaifi, J.; Gramfort, A.; Thirion, B.; Varoquaux, G. Machine learning for neuroimaging with scikit-learn. Front. Neuroinform. 2014, 8, 14. [Google Scholar] [CrossRef]
  23. Tustison, N.J.; Avants, B.B.; Cook, P.A.; Zheng, Y.; Egan, A.; Yushkevich, P.A.; Gee, J.C. N4ITK: Improved N3 Bias Correction. IEEE Trans. Med. Imaging 2010, 29, 1310–1320. [Google Scholar] [CrossRef] [PubMed]
  24. Avants, B.; Epstein, C.; Grossman, M.; Gee, J. Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain. Med. Image Anal. 2008, 12, 26–41. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, Y.; Brady, M.; Smith, S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imaging 2001, 20, 45–57. [Google Scholar] [CrossRef] [PubMed]
  26. Jenkinson, M.; Bannister, P.; Brady, M.; Smith, S. Improved Optimization for the Robust and Accurate Linear Registration and Motion Correction of Brain Images. NeuroImage 2002, 17, 825–841. [Google Scholar] [CrossRef] [PubMed]
  27. Cox, R.; Hyde, J. Software tools for analysis and visualization of fMRI data. NMR Biomed. 1997, 10, 171–178. [Google Scholar] [CrossRef]
  28. Greve, D.N.; Fischl, B. Accurate and robust brain image alignment using boundary-based registration. NeuroImage 2009, 48, 63–72. [Google Scholar] [CrossRef] [PubMed]
  29. Power, J.D.; Mitra, A.; Laumann, T.O.; Snyder, A.Z.; Schlaggar, B.L.; Petersen, S.E. Methods to detect, characterize, and remove motion artifact in resting state fMRI. NeuroImage 2014, 84, 320–341. [Google Scholar] [CrossRef]
  30. Wang, J.; Wange, X.; Xia, M.; Liao, X.; Evans, A.; He, Y. GRETNA: A graph theoretical network analysis toolbox for imaging connectomics. Front. Hum. Neurosci. 2015, 9, 386. [Google Scholar] [CrossRef]
  31. Tzourio-Mazoyer, N.; Landeau, B.; Papathanassiou, D.; Crivello, F.; Etard, O.; Delcroix, N.; Mazoyer, B.; Joliot, M. Automated Anatomical Labeling of Activations in SPM Using a Macroscopic Anatomical Parcellation of the MNI MRI Single-Subject Brain. NeuroImage 2002, 15, 273–289. [Google Scholar] [CrossRef]
  32. Fornito, A.; Zalesky, A.; Breakspear, M. Graph analysis of the human connectome: Promise, progress, and pitfalls. NeuroImage 2013, 80, 426–444. [Google Scholar] [CrossRef]
  33. Zalesky, A.; Fornito, A.; Bullmore, E.T. Network-based statistic: Identifying differences in brain networks. NeuroImage 2010, 53, 1197–1207. [Google Scholar] [CrossRef] [PubMed]
  34. Rubinov, M.; Sporns, O. Complex network measures of brain connectivity: Uses and interpretations. NeuroImage 2010, 52, 1059–1069. [Google Scholar] [CrossRef] [PubMed]
  35. Fornito, A.; Zalesky, A.; Bullmore, E.T. Fundamentals of Brain Network Analysis; Elsevier: Amsterdam, The Netherlands, 2016; p. xvii, 476. [Google Scholar]
  36. van den Heuvel, M.P.; Stam, C.; Boersma, M.; Pol, H.H. Small-world and scale-free organization of voxel-based resting-state functional connectivity in the human brain. NeuroImage 2008, 43, 528–539. [Google Scholar] [CrossRef] [PubMed]
  37. Kassambara, A. rstatix: Pipe-Friendly Framework for Basic Statistical Tests, R Package Version 0.7.0. 2021. Available online: https://rpkgs.datanovia.com/rstatix/(accessed on 7 July 2022).
  38. Signorell, A.; Aho, K.; Alfons, A.; Anderegg, N.; Aragon, T.; Arachchige, C.; Arppe, A.; Baddeley, A.; Barton, K.; Bolker, B.; et al. DescTools: Tools for Descriptive Statistics, R Package Version 0.99.36. 2022. Available online: https://cran.r-project.org/web/packages/DescTools/index.html(accessed on 7 July 2022).
  39. Dosenbach, N.U.F.; Nardos, B.; Cohen, A.L.; Fair, D.A.; Power, J.D.; Church, J.A.; Nelson, S.M.; Wig, G.S.; Vogel, A.C.; Lessov-Schlaggar, C.N.; et al. Prediction of Individual Brain Maturity Using fMRI. Science 2010, 329, 1358–1361. [Google Scholar] [CrossRef] [PubMed]
  40. Power, J.D.; Barnes, K.A.; Snyder, A.Z.; Schlaggar, B.L.; Petersen, S.E. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. NeuroImage 2012, 59, 2142–2154. [Google Scholar] [CrossRef]
  41. Rey, D.; Neuhäuser, M. Wilcoxon-Signed-Rank Test. In International Encyclopedia of Statistical Science; Springer: Berlin/Heidelberg, Germany, 2011; pp. 1658–1659. [Google Scholar]
  42. Achard, S.; Bullmore, E. Efficiency and Cost of Economical Brain Functional Networks. PLoS Comput. Biol. 2007, 3, 1–10. [Google Scholar] [CrossRef]
  43. Latora, V.; Marchiori, M. Economic small-world behavior in weighted networks. Eur. Phys. J. B—Condens. Matter Complex Systems 2003, 32, 249–263. [Google Scholar] [CrossRef]
  44. Burdette, J.H.; Laurienti, P.J.; Espeland, M.A.; Morgan, A.; Telesford, Q.; Vechlekar, C.D.; Hayasaka, S.; Jennings, J.M.; Katula, J.A.; Kraft, R.A.; et al. Using network science to evaluate exercise-associated brain changes in older adults. Front. Aging Neurosci. 2010, 2, 23. [Google Scholar] [CrossRef]
  45. Passingham, R.E.; Rowe, J.B. Functional specialization. In A Short Guide to Brain Imaging; Oxford University Press: Oxford, UK, 2015; pp. 71–90. [Google Scholar] [CrossRef]
  46. Kitano, H. Biological robustness. Nat. Rev. Genet. 2004, 5, 826–837. [Google Scholar] [CrossRef]
  47. Noppeney, U.; Friston, K.J.; Price, C.J. Degenerate neuronal systems sustaining cognitive functions. J. Anat. 2004, 205, 433–442. [Google Scholar] [CrossRef]
  48. Gallen, C.L.; D’Esposito, M. Brain Modularity: A Biomarker of Intervention-related Plasticity. Trends Cogn. Sci. 2019, 23, 293–304. [Google Scholar] [CrossRef] [PubMed]
  49. Aurich, N.K.; Alves Filho, J.O.; Marques da Silva, A.M.; Franco, A.R. Evaluating the reliability of different preprocessing steps to estimate graph theoretical measures in resting state fMRI data. Front. Neurosci. 2015, 9, 48. [Google Scholar] [CrossRef] [PubMed]
  50. Mårtensson, G.; Pereira, J.B.; Mecocci, P.; Vellas, B.; Tsolaki, M.; Kłoszewska, I.; Soininen, H.; Lovestone, S.; Simmons, A.; Volpe, G.; et al. Stability of graph theoretical measures in structural brain networks in Alzheimer’s disease. Sci. Rep. 2018, 8, 11592. [Google Scholar] [CrossRef] [PubMed]
  51. Liu, Y.; Wang, H.; Duan, Y.; Huang, J.; Ren, Z.; Ye, J.; Dong, H.; Shi, F.; Li, K.; Wang, J. Functional Brain Network Alterations in Clinically Isolated Syndrome and Multiple Sclerosis: A Graph-based Connectome Study. Radiology 2017, 282, 534–541. [Google Scholar] [CrossRef] [PubMed]
  52. Sporns, O.; Honey, C.J.; Kötter, R. Identification and Classification of Hubs in Brain Networks. PLoS ONE 2007, 2, e1049. [Google Scholar] [CrossRef] [PubMed]
  53. Korhonen, O.; Saarimäki, H.; Glerean, E.; Sams, M.; Saramäki, J. Consistency of Regions of Interest as nodes of fMRI functional brain networks. Netw. Neurosci. 2017, 1, 254–274. [Google Scholar] [CrossRef]
  54. Stanley, M.L.; Moussa, M.N.; Paolini, B.M.; Lyday, R.G.; Burdette, J.H.; Laurienti, P.J. Defining nodes in complex brain networks. Front. Comput. Neurosci. 2013, 7. [Google Scholar] [CrossRef]
  55. Salehi, M.; Greene, A.S.; Karbasi, A.; Shen, X.; Scheinost, D.; Constable, R.T. There is no single functional atlas even for a single individual: Functional parcel definitions change with task. NeuroImage 2020, 208, 116366. [Google Scholar] [CrossRef]
  56. Van den Heuvel, M.P.; de Lange, S.C.; Zalesky, A.; Seguin, C.; Yeo, B.T.; Schmidt, R. Proportional thresholding in resting-state fMRI functional connectivity networks and consequences for patient-control connectome studies: Issues and recommendations. NeuroImage 2017, 152, 437–449. [Google Scholar] [CrossRef]
Figure 1. Pre- and post-test distribution of parameters over the full range of used thresholds. The grey areas indicate the standard deviations. In the case of Hubness, the groups are represented in a boxplot. Enod = Nodal Efficiency; Eloc = Local Efficiency; Eglob = Global Efficiency; Q = Modularity; Sigma = Small-Worldness; Hub = number of Hubs; * = significant difference.
Figure 1. Pre- and post-test distribution of parameters over the full range of used thresholds. The grey areas indicate the standard deviations. In the case of Hubness, the groups are represented in a boxplot. Enod = Nodal Efficiency; Eloc = Local Efficiency; Eglob = Global Efficiency; Q = Modularity; Sigma = Small-Worldness; Hub = number of Hubs; * = significant difference.
Brainsci 12 01092 g001
Table 1. Overview of all uncorrected Results.
Table 1. Overview of all uncorrected Results.
TotalLowHigh
tptptp
AAL-116
          Eloc4.25<0.0014.38<0.0013.96<0.001
          Enod1.730.0971.730.0991.720.099
          Eglob2.210.0381.930.0672.280.033
          Mod−0.590.564−1.090.2860.060.956
          SW−0.840.409−1.060.3020.070.943
          Hub1.900.0702.020.0560.530.599
Dosenbach-160
          Eloc0.690.5000.290.7771.100.283
          Enod0.440.6680.500.6220.380.712
          Eglob1.830.0801.790.0881.510.146
          Mod0.030.9790.200.847−0.180.856
          SW0.010.9910.010.9920.030.978
          Hub0.390.7020.280.7830.440.666
Power-264
          Eloc0.040.9710.440.665*136**0.964*
          Enod0.380.7070.320.7490.450.657
          Eglob1.640.1151.700.1041.460.159
          Mod−0.100.920<0.010.997−0.240.813
          SW−0.410.690−0.400.690−0.380.710
          Hub0.770.4490.740.465*108.5**0.378*
Bold = p < 0.05; *Italic* = Wilcoxon-Signed-Rank Test; Total = 6–40% thresholds in 2% steps; Low = 6–22% thresholds in 2% steps; High = 24–40% thresholds in 2% steps; Eloc = Local Efficiency; Enod = Nodal Efficiency; Eglob = Global Efficiency; Mod = Modularity; SW = Small-Worldness; Hub = Hubness. Note. The Wilcoxon- Signed-Rank test [41] has been used in two of the presented comparisons. This was because, in both cases, the assumptions were not met. Therefore, their values are not directly comparable to the other tests and should be neglected when comparing outcomes.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Riedmann, U.; Fink, A.; Weber, B.; Koschutnig, K. Functional Connectivity as an Index of Brain Changes Following a Unicycle Intervention: A Graph-Theoretical Network Analysis. Brain Sci. 2022, 12, 1092. https://doi.org/10.3390/brainsci12081092

AMA Style

Riedmann U, Fink A, Weber B, Koschutnig K. Functional Connectivity as an Index of Brain Changes Following a Unicycle Intervention: A Graph-Theoretical Network Analysis. Brain Sciences. 2022; 12(8):1092. https://doi.org/10.3390/brainsci12081092

Chicago/Turabian Style

Riedmann, Uwe, Andreas Fink, Bernhard Weber, and Karl Koschutnig. 2022. "Functional Connectivity as an Index of Brain Changes Following a Unicycle Intervention: A Graph-Theoretical Network Analysis" Brain Sciences 12, no. 8: 1092. https://doi.org/10.3390/brainsci12081092

APA Style

Riedmann, U., Fink, A., Weber, B., & Koschutnig, K. (2022). Functional Connectivity as an Index of Brain Changes Following a Unicycle Intervention: A Graph-Theoretical Network Analysis. Brain Sciences, 12(8), 1092. https://doi.org/10.3390/brainsci12081092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop