Persistent homology reveals robustness loss in inhaled substance abuse rs-fMRI networks.

Analyzing functional brain activity through functional magnetic resonance imaging (fMRI) is commonly done using tools from graph theory for the analysis of the correlation matrices. A drawback of these methods is that the networks must be restricted to values of the weights of the edges within certain thresholds and there is no consensus about the best choice of such thresholds. Topological data analysis (TDA) is a recently-developed tool in algebraic topology which allows us to analyze networks through combinatorial spaces obtained from them, with the advantage that all the possible thresholds can be considered at once. In this paper we applied TDA, in particular persistent homology, to study correlation matrices from rs-fMRI, and through statistical analysis, we detected significant differences between the topological structures of adolescents with inhaled substance abuse disorder (ISAD) and healthy controls. We interpreted the topological differences as indicative of a loss of robustness in the functional brain networks of the ISAD population.

Directed brain interactions over time: A resting-state EEG comparison between schizophrenia and healthy individuals.

Understanding the neurophysiological mechanisms of schizophrenia (SZ) is one of the challenges of neuroscience. Many anatomical and functional studies have pointed to problems in brain connectivity in SZ individuals. However, little is known about the relationships between specific brain regions and impairments in brain connectivity in SZ individuals. Herein we propose a new approach using time-varying graphs and the motif synchronization method to build dynamic brain functional networks (BFNs). Dynamic BFNs were constructed from resting-state electroencephalography (rs-EEG) of 14 schizophrenia (SZ) individuals and 14 healthy controls (HCs). BFNs were evaluated based on the percentage of synchronization importance between a pair of regions (considering external and internal interactions) over time. We found differences in the directed interaction between brain regions in SZ individuals compared to the control group. Our method revealed low bilaterally directed interactions between the temporal lobes in SZ individuals compared to HCs, indicating a potential link between altered brain connectivity and the characteristic symptoms of schizophrenia. From a clinical perspective, these results shed light on developing new therapeutic approaches targeting these specific neural interactions that are altered in individuals with SZ. This knowledge allows the application of better interventions focused on restoring or compensating for interrupted connectivity patterns.

Greater Neighborhood Disadvantage Is Associated with Alterations in Fetal Functional Brain Network Structure.

OBJECTIVE: To determine the association between neighborhood disadvantage (ND) and functional brain development of in utero fetuses. STUDY DESIGN: We conducted an observational study using Social Vulnerability Index (SVI) scores to assess the impact of ND on a prospectively recruited sample of healthy pregnant women from Washington, DC. Using 79 functional magnetic resonance imaging scans from 68 healthy pregnancies at a mean gestational age of 33.12 weeks, we characterized the overall functional brain network structure using a graph metric approach. We used linear mixed effects models to assess the relationship between SVI and gestational age on 5 graph metrics, adjusting for multiple scans. RESULTS: Exposure to greater ND was associated with less well integrated functional brain networks, as observed by longer characteristic path lengths and diminished global efficiency (GE), as well as diminished small world propensity (SWP). Across gestational ages, however, the association between SVI and network integration diminished to a negligible relationship in the third trimester. Conversely, SWP was significant across pregnancy, but the relationship changed such that there was a negative association with SWP earlier in the second trimester that inverted around the transition to the third trimester to a positive association. CONCLUSIONS: These data directly connect ND and altered functional brain maturation in fetuses. Our results suggest that, even before birth, proximity to environmental stressors in the wider neighborhood environment are associated with altered brain development.

Differences between Alzheimer's disease and mild cognitive impairment using brain networks from magnetic resonance texture analysis.

Several studies have aimed at identifying biomarkers in the initial phases of Alzheimer's disease (AD). Conversely, texture features, such as those from gray-level co-occurrence matrices (GLCMs), have highlighted important information from several types of medical images. More recently, texture-based brain networks have been shown to provide useful information in characterizing healthy individuals. However, no studies have yet explored the use of this type of network in the context of AD. This work aimed to employ texture brain networks to investigate the distinction between groups of patients with amnestic mild cognitive impairment (aMCI) and mild dementia due to AD, and a group of healthy subjects. Magnetic resonance (MR) images from the three groups acquired at two instances were used. Images were segmented and GLCM texture parameters were calculated for each region. Structural brain networks were generated using regions as nodes and the similarity among texture parameters as links, and graph theory was used to compute five network measures. An ANCOVA was performed for each network measure to assess statistical differences between groups. The thalamus showed significant differences between aMCI and AD patients for four network measures for the right hemisphere and one network measure for the left hemisphere. There were also significant differences between controls and AD patients for the left hippocampus, right superior parietal lobule, and right thalamus-one network measure each. These findings represent changes in the texture of these regions which can be associated with the cortical volume and thickness atrophies reported in the literature for AD. The texture networks showed potential to differentiate between aMCI and AD patients, as well as between controls and AD patients, offering a new tool to help understand these conditions and eventually aid early intervention and personalized treatment, thereby improving patient outcomes and advancing AD research.

Resting state networks in patients with acute disorders of consciousness after severe traumatic brain injury.

OBJECTIVES: This study aims to describe resting state networks (RSN) in patients with disorders of consciousness (DOC)s after acute severe traumatic brain injury (TBI). METHODS: Adult patients with TBI with a GCS score <8 who remained in a coma, minimally conscious state (MCS), or unresponsive wakefulness syndrome (UWS), between 2017 and 2020 were included. Blood-oxygen-level dependent imaging was performed to compare their RSN with 10 healthy volunteers. RESULTS: Of a total of 293 patients evaluated, only 13 patients were included according to inclusion criteria: 7 in coma (54%), 2 in MCS (15%), and 4 (31%) had an UWS. RSN analysis showed that the default mode network (DMN) was present and symmetric in 6 patients (46%), absent in 1 (8%), and asymmetric in 6 (46%). The executive control network (ECN) was present in all patients but was asymmetric in 3 (23%). The right ECN was absent in 2 patients (15%) and the left ECN in 1 (7%). The medial visual network was present in 11 (85%) patients. Finally, the cerebellar network was symmetric in 8 patients (62%), asymmetric in 1 (8%), and absent in 4 (30%). CONCLUSIONS: A substantial impairment in activation of RSN is demonstrated in patients with DOC after severe TBI in comparison with healthy subjects. Three patterns of activation were found: normal/complete activation, 2) asymmetric activation or partially absent, and 3) absent activation.

Functional and Effective Connectivity Underlying Semantic Verbal Fluency.

Semantic verbal fluency (SVF) impairment is present in several neurological disorders. Although activation in SVF-related areas has been reported, how these regions are connected and their functional roles in the network remain divergent. We assessed SVF static and dynamic functional connectivity (FC) and effective connectivity in healthy participants using functional magnetic resonance imaging. We observed activation in the inferior frontal (IFG), middle temporal (pMTG) and angular gyri (AG), anterior cingulate (AC), insular cortex, and regions of the superior, middle, and medial frontal gyri (SFG, MFG, MidFG). Our static FC analysis showed a highly interconnected task and resting state network. Increased connectivity of AC with the pMTG and AG was observed for the task. The dynamic FC analysis provided circuits with connections similarly modulated across time and regions related to category identification, language comprehension, word selection and recovery, word generation, inhibition of speaking, speech planning, and articulatory planning of orofacial movements. Finally, the effective connectivity analysis provided a network that best explained our data, starting at the AG and going to the pMTG, from which there was a division between the ventral and dorsal streams. The SFG and MFG regions were connected and modulated by the MidFG, while the inferior regions formed the ventral stream. Therefore, we successfully assessed the SVF network, exploring regions associated with the entire processing, from category identification to word generation. The methodological approach can be helpful for further investigation of the SVF network in neurological disorders.

Neurocognition and brain functional connectivity in a non-clinical population-based sample with psychotic experiences.

We characterized the neurocognitive profile of communed-based individuals and unaffected siblings of patients with psychosis from Brazil reporting psychotic experiences (PEs). We also analyzed associations between PEs and the intra and inter-functional connectivity (FC) in the Default Mode Network (DMN), the Fronto-Parietal Network (FPN) and the Salience Network (SN) measured by functional magnetic resonance imaging. The combined sample of communed-based individuals and unaffected siblings of patients with psychosis comprised 417 (neurocognition) and 85 (FC) volunteers who were divided as having low (<75th percentile) and high (≥75th percentile) PEs (positive, negative, and depressive dimensions) assessed by the Community Assessment of Psychic Experiences. The neurocognitive profile and the estimated current brief intellectual quotient (IQ) were assessed using the digit symbol (processing speed), arithmetic (working memory), block design (visual learning) and information (verbal learning) subtests of Wechsler Adult Intelligence Scale-third edition. Logistic regression models were performed for neurocognitive analysis. For neuroimaging, we used the CONN toolbox to assess FC between the specified regions, and ROI-to-ROI analysis. In the combined sample, high PEs (all dimensions) were related to lower processing speed performance. High negative PEs were related to poor visual learning performance and lower IQ, while high depressive PEs were associated with poor working memory performance. Those with high negative PEs presented FPN hypoconnectivity between the right and left lateral prefrontal cortex. There were no associations between PEs and the DMN and SN FC. Brazilian individuals with high PEs showed neurocognitive impairments like those living in wealthier countries. Hypoconnectivity in the FPN in a community sample with high PEs is coherent with the hypothesis of functional dysconnectivity in schizophrenia.

Intra and inter-individual variability in functional connectomes of patients with First Episode of Psychosis.

Patients with Schizophrenia may show different clinical presentations, not only regarding inter-individual comparisons but also in one specific subject over time. In fMRI studies, functional connectomes have been shown to carry valuable individual level information, which can be associated with cognitive and behavioral variables. Moreover, functional connectomes have been used to identify subjects within a group, as if they were fingerprints. For the particular case of Schizophrenia, it has been shown that there is reduced connectome stability as well as higher inter-individual variability. Here, we studied inter and intra-individual heterogeneity by exploring functional connectomes' variability and related it with clinical variables (PANSS Total scores and antipsychotic's doses). Our sample consisted of 30 patients with First Episode of Psychosis and 32 Healthy Controls, with a test-retest approach of two resting-state fMRI scanning sessions. In our patients' group, we found increased deviation from healthy functional connectomes and increased intragroup inter-subject variability, which was positively correlated to symptoms' levels in six subnetworks (visual, somatomotor, dorsal attention, ventral attention, frontoparietal and DMN). Moreover, changes in symptom severity were positively related to changes in deviation from healthy functional connectomes. Regarding intra-subject variability, we were unable to replicate previous findings of reduced connectome stability (i.e., increased intra-subject variability), but we found a trend suggesting that result. Our findings highlight the relevance of variability characterization in Schizophrenia, and they can be related to evidence of Schizophrenia patients having a noisy functional connectome.

High-order functional redundancy in ageing explained via alterations in the connectome in a whole-brain model.

The human brain generates a rich repertoire of spatio-temporal activity patterns, which support a wide variety of motor and cognitive functions. These patterns of activity change with age in a multi-factorial manner. One of these factors is the variations in the brain's connectomics that occurs along the lifespan. However, the precise relationship between high-order functional interactions and connnectomics, as well as their variations with age are largely unknown, in part due to the absence of mechanistic models that can efficiently map brain connnectomics to functional connectivity in aging. To investigate this issue, we have built a neurobiologically-realistic whole-brain computational model using both anatomical and functional MRI data from 161 participants ranging from 10 to 80 years old. We show that the differences in high-order functional interactions between age groups can be largely explained by variations in the connectome. Based on this finding, we propose a simple neurodegeneration model that is representative of normal physiological aging. As such, when applied to connectomes of young participant it reproduces the age-variations that occur in the high-order structure of the functional data. Overall, these results begin to disentangle the mechanisms by which structural changes in the connectome lead to functional differences in the ageing brain. Our model can also serve as a starting point for modeling more complex forms of pathological ageing or cognitive deficits.

Stationary EEG pattern relates to large-scale resting state networks - An EEG-fMRI study connecting brain networks across time-scales.

Relating brain dynamics acting on time scales that differ by at least an order of magnitude is a fundamental issue in brain research. The same is true for the observation of stable dynamical structures in otherwise highly non-stationary signals. The present study addresses both problems by the analysis of simultaneous resting state EEG-fMRI recordings of 53 patients with epilepsy. Confirming previous findings, we observe a generic and temporally stable average correlation pattern in EEG recordings. We design a predictor for the General Linear Model describing fluctuations around the stationary EEG correlation pattern and detect resting state networks in fMRI data. The acquired statistical maps are contrasted to several surrogate tests and compared with maps derived by spatial Independent Component Analysis of the fMRI data. By means of the proposed EEG-predictor we observe core nodes of known fMRI resting state networks with high specificity in the default mode, the executive control and the salience network. Our results suggest that both, the stationary EEG pattern as well as resting state fMRI networks are different expressions of the same brain activity. This activity is interpreted as the dynamics on (or close to) a stable attractor in phase space that is necessary to maintain the brain in an efficient operational mode. We discuss that this interpretation is congruent with the theoretical framework of complex systems as well as with the brain's energy balance.

A deep learning based approach identifies regions more relevant than resting-state networks to the prediction of general intelligence from resting-state fMRI.

Prediction of cognitive ability latent factors such as general intelligence from neuroimaging has elucidated questions pertaining to their neural origins. However, predicting general intelligence from functional connectivity limit hypotheses to that specific domain, being agnostic to time-distributed features and dynamics. We used an ensemble of recurrent neural networks to circumvent this limitation, bypassing feature extraction, to predict general intelligence from resting-state functional magnetic resonance imaging regional signals of a large sample (n = 873) of Human Connectome Project adult subjects. Ablating common resting-state networks (RSNs) and measuring degradation in performance, we show that model reliance can be mostly explained by network size. Using our approach based on the temporal variance of saliencies, that is, gradients of outputs with regards to inputs, we identify a candidate set of networks that more reliably affect performance in the prediction of general intelligence than similarly sized RSNs. Our approach allows us to further test the effect of local alterations on data and the expected changes in derived metrics such as functional connectivity and instantaneous innovations.

Long-term stability of the cortical volumetric profile and the functional human connectome throughout childhood and adolescence.

There is compelling evidence showing that between-subject variability in several functional and structural brain features is sufficient for unique identification in adults. However, individuation of brain functional connectomes depends on the stabilization of neurodevelopmental processes during childhood and adolescence. Here, we aimed to (1) evaluate the intra-subject functional connectome stability over time for the whole brain and for large scale functional networks and (2) determine the long-term identification accuracy or 'fingerprinting' for the cortical volumetric profile and the functional connectome. For these purposes, we analysed a longitudinal cohort of 239 children and adolescents scanned in two sessions with an interval of approximately 3 years (age range 6-15 years at baseline and 9-18 years at follow-up). Corroborating previous results using short between-scan intervals in children and adolescents, we observed a moderate identification accuracy (38%) for the whole functional profile. In contrast, identification accuracy using cortical volumetric profile was 95%. Among the large-scale networks, the default-mode (26.8%), the frontoparietal (23.4%) and the dorsal-attention (27.6%) networks were the most discriminative. Our results provide further evidence for a protracted development of specific individual structural and functional connectivity profiles.

Cerebral Structure and Function in Stroke-free Patients with Atrial Fibrillation.

OBJECTIVES: Atrial fibrillation (AF) is associated with high risk of dementia and brain atrophy in stroke-free patients, but the mechanisms underlying this association remain unclear. We aimed to examine the brain volume and connectivity of paramount cognitive brain networks in stroke-free patients with AF without dementia. MATERIALS AND METHODS: Twenty-six stroke-free patients with AF and 26 age and sex-matched subjects without AF were submitted to a 3-tesla brain structural and functional MRI. An extensive clinical evaluation excluded stroke, dementia, low cardiac output, carotid stenosis and metabolic diseases without optimal therapy. We used CHA2DS2-VASc score to classify the cardiovascular risk factor burden and a broad neuropsychological battery to assess the cognitive performance. Voxel based morphometry analysis of. structural MRI defined whole-brain gray and white matter volumes. Finally, we used eco-plannar MRI images to compare the differences of functional connectivity of 7 large-scale resting-state networks between AF patients and controls. RESULTS: Taking into account the history of hypertension and heart failure, AF was associated to volume decrease of the right basal frontal lobe and right inferior cerebellum. Decreased connectivity of the ventral Default Mode Network (vDMN) was observed in the AF group. No disruption of connectivity was observed in the executive, visuospatial and salience networks. CONCLUSION: Individuals with AF without stroke or dementia have subtle reduction of gray and white matter, restricted to frontal areas and cerebellum. These patients show decreased vDMN connectivity, without other large-scale brain network disruption.

Cholinergic neuromodulation of inhibitory interneurons facilitates functional integration in whole-brain models.

Segregation and integration are two fundamental principles of brain structural and functional organization. Neuroimaging studies have shown that the brain transits between different functionally segregated and integrated states, and neuromodulatory systems have been proposed as key to facilitate these transitions. Although whole-brain computational models have reproduced this neuromodulatory effect, the role of local inhibitory circuits and their cholinergic modulation has not been studied. In this article, we consider a Jansen & Rit whole-brain model in a network interconnected using a human connectome, and study the influence of the cholinergic and noradrenergic neuromodulatory systems on the segregation/integration balance. In our model, we introduce a local inhibitory feedback as a plausible biophysical mechanism that enables the integration of whole-brain activity, and that interacts with the other neuromodulatory influences to facilitate the transition between different functional segregation/integration regimes in the brain.

Effects of ovariectomy on the inputs from the medial nucleus of the amygdala to the ventromedial nucleus of the hypothalamus in young adult rats.

During puberty, sexual hormones induce crucial changes in neural circuit organization, leading to significant sexual dimorphism in adult behaviours. The ventrolateral division of the ventromedial nucleus of the hypothalamus (VMHvl) is the major neural site controlling the receptive component of female sexual behaviour, which is dependent on ovarian hormones. The inputs to the VMHvl, originating from the medial nucleus of the amygdala (MeA), transmit essential information to trigger such behaviour. In this study, we investigated the projection pattern of the MeA to the VMHvl in ovariectomized rats at early puberty. Six-week-old Sprague-Dawley rats were ovariectomized (OVX) and, upon reaching 90 days of age, were subjected to iontophoretic injections of the neuronal anterograde tracer Phaseolus vulgaris leucoagglutinin into the MeA. Projections from the MeA to the VMHvl and to other structures included in the neural circuit responsible for female sexual behaviour were analysed in the Control and OVX groups. The results of the semi-quantitative analysis showed that peripubertal ovariectomy reduced the density of intra-amygdalar fibres. The stereological estimates, however, failed to find changes in the organization of the terminal fields of nerve fibres from the MeA to the VMHvl in the adult. The present data show that ovariectomized rats during the peripubertal phase did not undergo significant changes in MeA fibres reaching the VMHvl; however, they suggest a possible effect of ovariectomy on MeA connectivity under amygdalar subnuclei.

Development of the brain functional connectome follows puberty-dependent nonlinear trajectories.

Adolescence is a developmental period that dramatically impacts body and behavior, with pubertal hormones playing an important role not only in the morphological changes in the body but also in brain structure and function. Understanding brain development during adolescence has become a priority in neuroscience because it coincides with the onset of many psychiatric and behavioral disorders. However, little is known about how puberty influences the brain functional connectome. In this study, taking a longitudinal human sample of typically developing children and adolescents (of both sexes), we demonstrate that the development of the brain functional connectome better fits pubertal status than chronological age. In particular, centrality, segregation, efficiency, and integration of the brain functional connectome increase after the onset of the pubertal markers. We found that these effects are stronger in attention and task control networks. Lastly, after controlling for this effect, we showed that functional connectivity between these networks is related to better performance in cognitive flexibility. This study points out the importance of considering longitudinal nonlinear trends when exploring developmental trajectories, and emphasizes the impact of puberty on the functional organization of the brain in adolescence.

Anterior Insular Resting-State Functional Connectivity is Related to Cognitive Reserve in Multiple Sclerosis.

BACKGROUND AND PURPOSE: Cognitive dysfunction is common in multiple sclerosis (MS). The dorsal anterior insula (dAI) is a key hub of the salience network (SN) orchestrating access to critical cognitive brain regions. The aim of this study was to assess whole-brain dAI intrinsic functional connectivity (iFC) using resting-state functional MRI (rs-fMRI) in people with MS and healthy controls (HC) and test the relationship between cognitive reserve (CR) and dAI iFC in people with MS. METHODS: We studied 28 people with relapsing-remitting MS and 28 HC. CR index was quantified by combining premorbid IQ, leisure activities, and education level. For whole-brain iFC analyses, the bilateral dAI were used as seeds. Individual subject correlation maps were entered into general linear models for group comparison and to analyze the effect of CR index on dAI iFC, controlling for multiple comparisons. The correlation between CR index and iFC was assessed using a linear regression model. RESULTS: rs-fMRI analyses revealed a negative relationship between CR index and iFC within the left dAI and a left occipital cluster in people with MS including regions of the cuneus, superior occipital gyrus, and parieto-occipital sulcus. The regression analysis showed that people with MS and a higher CR index had a statistically significantly reduced iFC within the left dAI and the cluster. CONCLUSIONS: CR is relevant to functional connectivity within one of the main nodes of the SN, the dAI, and occipital regions in MS. These results have implications for how CR may modulate the susceptibility to cognitive dysfunction in MS.

Non-classical behavior of the default mode network regions during an information processing task.

The default mode network (DMN) efficient deactivation and suppressed functional connectivity (FC) during goal-directed tasks, which require attentional resources, have been considered essential to healthy brain cognition. However, recent studies have shown that DMN regions do not always show the expected behavior. Then, we aimed to investigate the functional activation and connectivity of DMN nodes in young, healthy controls during a goal-directed task. We used an adaptation of the symbol digit modalities test (SDMT) to evaluate the information processing speed (IPS). Twenty-four subjects (10 women, age: 29 ± 7 years) underwent two functional Magnetic Resonance Imaging experiments: one during resting-state and one during a block-designed SDMT paradigm. We superimposed the templates of the DMN on the group activation map and observed the reorganization of the network. For the posterior cingulate cortex (PCC) node of the DMN, which is spatially extensive, comprising the precuneus (dorsal portion) and the posterior cingulate gyrus (PCG, ventral portion), the extent of each region was different between conditions, suggesting different functional roles for them. Therefore, for the functional connectivity (FC) analysis, we split the DMN-PCC region into two regions: left precuneus (BA 7) and PCG. The left precuneus (BA 7) was positively correlated with the left lingual gyrus (BA 17), a task-positive region, and negatively associated with the DMN nodes when comparing task performance with the resting-state condition. The other DMN regions presented the classical antagonistic role during the attentional task. In conclusion, we found that the activation and functional connectivity of the DMN is, in general, suppressed during the information processing. However, the left precuneus BA 7 presented a context-dependent modulatory behavior, working as a transient in-between hub connecting the DMN to task-positive areas. Such findings support studies that show increased activation and excitatory functional connectivity of DMN portions during goal-directed tasks. Moreover, our results may contribute to defining more precise functional correlates of IPS deficits in a wide range of clinical and neurological diseases.

FFClust: Fast fiber clustering for large tractography datasets for a detailed study of brain connectivity.

Automated methods that can identify white matter bundles from large tractography datasets have several applications in neuroscience research. In these applications, clustering algorithms have shown to play an important role in the analysis and visualization of white matter structure, generating useful data which can be the basis for further studies. This work proposes FFClust, an efficient fiber clustering method for large tractography datasets containing millions of fibers. Resulting clusters describe the whole set of main white matter fascicles present on an individual brain. The method aims to identify compact and homogeneous clusters, which enables several applications. In individuals, the clusters can be used to study the local connectivity in pathological brains, while at population level, the processing and analysis of reproducible bundles, and other post-processing algorithms can be carried out to study the brain connectivity and create new white matter bundle atlases. The proposed method was evaluated in terms of quality and execution time performance versus the state-of-the-art clustering techniques used in the area. Results show that FFClust is effective in the creation of compact clusters, with a low intra-cluster distance, while keeping a good quality Davies-Bouldin index, which is a metric that quantifies the quality of clustering approaches. Furthermore, it is about 8.6 times faster than the most efficient state-of-the-art method for one million fibers dataset. In addition, we show that FFClust is able to correctly identify atlas bundles connecting different brain regions, as an example of application and the utility of compact clusters.

Transitions between human functional brain networks reveal complex, cost-efficient and behaviorally-relevant temporal paths.

Resting-state functional MRI activity is organized as a complex network. However, this coordinated brain activity changes with time, raising questions about its evolving temporal arrangement. Does the brain visit different configurations through time in a random or ordered way? Advances in this area depend on developing novel paradigms that would allow us to shed light on these issues. We here propose to study the temporal changes in the functional connectome by looking at transition graphs of network activity. Nodes of these graphs correspond to brief whole-brain connectivity patterns (or meta-states), and directed links to the temporal transition between consecutive meta-states. We applied this method to two datasets of healthy subjects (160 subjects and a replication sample of 54), and found that transition networks had several non-trivial properties, such as a heavy-tailed degree distribution, high clustering, and a modular organization. This organization was implemented at a low biological cost with a high cost-efficiency of the dynamics. Furthermore, characteristics of the subjects' transition graphs, including global efficiency, local efficiency and their transition cost, were correlated with cognition and motor functioning. All these results were replicated in both datasets. We conclude that time-varying functional connectivity patterns of the brain in health progress in time in a highly organized and complex order, which is related to behavior.