Generic aspects of complexity in brain imaging data and other biological systems

Hyperbolic multivariate feature learning in higher-order heterogeneous networks for drug-disease prediction

New drug discovery has always been a costly, time-consuming process with a high failure rate. Repurposing existing drugs offers a valuable alternative and reduces the risks associated with developing new drugs. Various experimental methods have been employed to facilitate drug repositioning; however, associations prediction between drugs and diseases through biological experiments is both expensive and time-consuming. Consequently, it is imperative to develop efficient and highly precise computational methods for predicting these associations. Based on this, we propose a drug-disease associations prediction method based on Hyperbolic Multivariate feature Learning in High-order Heterogeneous Networks for Drug-Disease Prediction, called H3ML. Our approach begins by mining high-order information from protein-disease and drug-protein networks to construct high-order heterogeneous networks. Subsequently, we employ multivariate feature learning to create hyperbolic representations, and then enhance the features of the heterogeneous network. Finally, we utilize a hyperbolic graph attention network in the hyperbolic space to aggregate neighbor information and perform the final prediction task. In addition, we evaluate the performance of H3ML by comparing it with some state-of-the-art methods across different datasets. The case study further validate the effectiveness of H3ML. Our implementation will be publicly available at: https://github.com/jianruichen/H-3ML.

Can ephapticity contribute to brain complexity?

The inquiry into the origin of brain complexity remains a pivotal question in neuroscience. While synaptic stimuli are acknowledged as significant, their efficacy often falls short in elucidating the extensive interconnections of the brain and nuanced levels of cognitive integration. Recent advances in neuroscience have brought the mechanisms underlying the generation of highly intricate dynamics, emergent patterns, and sophisticated oscillatory signals into question. Within this context, our study, in alignment with current research, postulates the hypothesis that ephaptic communication, in addition to synaptic mediation's, may emerge as a prime candidate for unraveling optimal brain complexity. Ephaptic communication, hitherto little studied, refers to direct interactions of the electric field between adjacent neurons, without the mediation of traditional synapses (electrical or chemical). We propose that these electric field couplings may provide an additional layer of connectivity that facilitates the formation of complex patterns and emergent dynamics in the brain. In this investigation, we conducted a comparative analysis between two types of networks utilizing the Quadratic Integrate-and-Fire Ephaptic model (QIF-E): (I) a small-world synaptic network (ephaptic-off) and (II) a mixed composite network comprising a small-world synaptic network with the addition of an ephaptic network (ephaptic-on). Utilizing the Multiscale Entropy methodology, we conducted an in-depth analysis of the responses generated by both network configurations, with complexity assessed by integrating across all temporal scales. Our findings demonstrate that ephaptic coupling enhances complexity under specific topological conditions, considering variables such as time, spatial scales, and synaptic intensity. These results offer fresh insights into the dynamics of communication within the nervous system and underscore the fundamental role of ephapticity in regulating complex brain functions.

Psychoactive substances for the treatment of neuropsychiatric disorders

In the contemporary landscape of psychiatric medicine, critical advancements have been noted in the utilization of psychoactive substances such as hallucinogens, 3,4-methylenedioxymethamphetamine (MDMA), and ketamine for the treatment of severe mental health disorders. This review provides a detailed evaluation of these substances, focusing on their mechanisms of action and the profound clinical outcomes observed in controlled environments. Hallucinogens like lysergic acid diethylamide and psilocybin primarily target the 5-HT2A receptor agonist-2 (5-HT2AR), inducing substantial perceptual and cognitive shifts that facilitate deep psychological introspection and significant therapeutic advances, particularly in patients suffering from depression and anxiety disorders. MDMA, influencing multiple neurotransmitter systems including 5-Hydroxytryptamine (5-HT), dopamine, and norepinephrine, has been demonstrated to effectively alleviate symptoms of post-traumatic stress disorder, enhancing patients' emotional engagement and resilience during psychotherapy. Meanwhile, ketamine, a glutamate receptor antagonist, rapidly alleviates depressive symptoms, offering a lifeline for individuals with treatment-resistant depression through its fast-acting antidepressant properties. The integration of these substances into psychiatric practice has shown promising results, fundamentally changing the therapeutic landscape for patients unresponsive to traditional treatment modalities. However, the potent effects of these agents also necessitate a cautious approach in clinical application, ensuring careful dosage control, monitoring, and risk management to prevent potential abuse and mitigate adverse effects.

Valorization of Pig Brains for Prime Quality Oil: A Comparative Evaluation of Organic-Solvent-Based and Solvent-Free Extractions

Pig processing industries have produced large quantities of by-products, which have either been discarded or used to make low-value products. This study aimed to provide recommendations for manufacturing edible oil from pig brains, thereby increasing the value of pork by-products. The experiment compared non-solvent extraction methods, specifically wet rendering and aqueous saline, to a standard solvent extraction method, the Bligh and Dyer method, for extracting oil from pig brains. The yield, color, fatty acid profile, a number of lipid classes, and lipid stability against lipolysis and oxidation of the pig brain oil were comprehensively compared, and the results revealed that these parameters varied depending on the extraction method. The wet rendering process provided the highest extracted oil yield (~13%), followed by the Bligh and Dyer method (~7%) and the aqueous saline method (~2.5%). The Bligh and Dyer method and wet rendering techniques produced a translucent yellow oil; however, an opaque light-brown-red oil was found in the aqueous saline method. The Bligh and Dyer method yielded the oil with the highest phospholipid, cholesterol, carotenoid, tocopherol, and free fatty acid contents (p < 0.05). Although the Bligh and Dyer method recovered the most unsaturated fatty acids, it also recovered more trans-fatty acids. Aqueous saline and wet rendering procedures yielded oil with low FFA levels (<1 g/100 g). The PV of the oil extracted using all methods was <1 meq/kg; however, the Bligh and Dyer method had a significant TBARS content (7.85 mg MDA equivalent/kg) compared to aqueous saline (1.75 mg MDA equivalent/kg) and wet rendering (1.14 mg MDA equivalent/kg) (p < 0.05). FTIR spectra of the pig brain oil revealed the presence of multiple components in varying quantities, as determined by chemical analysis experiments. Given the higher yield and lipid stability and the lower cholesterol and trans-fatty acid content, wet rendering can be regarded as a simple and environmentally friendly method for safely extracting quality edible oil from pig brains, which may play an important role in obtaining financial benefits, nutrition, the zero-waste approach, and increasing the utilization of by-products in the meat industry.

Screen time, brain network development and socio-emotional competence in childhood: moderation of associations by parent-child reading

BACKGROUND

Screen time in infancy is linked to changes in social-emotional development but the pathway underlying this association remains unknown. We aim to provide mechanistic insights into this association using brain network topology and to examine the potential role of parent-child reading in mitigating the effects of screen time.

METHODS

We examined the association of screen time on brain network topology using linear regression analysis and tested if the network topology mediated the association between screen time and later socio-emotional competence. Lastly, we tested if parent-child reading time was a moderator of the link between screen time and brain network topology.

RESULTS

Infant screen time was significantly associated with the emotion processing-cognitive control network integration (p = 0.005). This network integration also significantly mediated the association between screen time and both measures of socio-emotional competence (BRIEF-2 Emotion Regulation Index, p = 0.04; SEARS total score, p = 0.04). Parent-child reading time significantly moderated the association between screen time and emotion processing-cognitive control network integration (β = -0.640, p = 0.005).

CONCLUSION

Our study identified emotion processing-cognitive control network integration as a plausible biological pathway linking screen time in infancy and later socio-emotional competence. We also provided novel evidence for the role of parent-child reading in moderating the association between screen time and topological brain restructuring in early childhood.

Weighting the structural connectome: Exploring its impact on network properties and predicting cognitive performance in the human brain

Brain function does not emerge from isolated activity, but rather from the interactions and exchanges between neural elements that form a network known as the connectome. The human connectome consists of structural and functional aspects. The structural connectome (SC) represents the anatomical connections, and the functional connectome represents the resulting dynamics that emerge from this arrangement of structures. As there are different ways of weighting these connections, it is important to consider how such different approaches impact study conclusions. Here, we propose that different weighted connectomes result in varied network properties, and while neither superior the other, selection might affect interpretation and conclusions in different study cases. We present three different weighting models, namely, number of streamlines (NOS), fractional anisotropy (FA), and axon diameter distribution (ADD), to demonstrate these differences. The later, is extracted using recently published AxSI method and is first compared to commonly used weighting methods. Moreover, we explore the functional relevance of each weighted SC, using the Human Connectome Project (HCP) database. By analyzing intelligence-related data, we develop a predictive model for cognitive performance based on graph properties and the National Institutes of Health (NIH) toolbox. Results demonstrate that the ADD SC, combined with a functional subnetwork model, outperforms other models in estimating cognitive performance.

Complexity Measures for EEG Microstate Sequences: Concepts and Algorithms

EEG microstate sequence analysis quantifies properties of ongoing brain electrical activity which is known to exhibit complex dynamics across many time scales. In this report we review recent developments in quantifying microstate sequence complexity, we classify these approaches with regard to different complexity concepts, and we evaluate excess entropy as a yet unexplored quantity in microstate research. We determined the quantities entropy rate, excess entropy, Lempel-Ziv complexity (LZC), and Hurst exponents on Potts model data, a discrete statistical mechanics model with a temperature-controlled phase transition. We then applied the same techniques to EEG microstate sequences from wakefulness and non-REM sleep stages and used first-order Markov surrogate data to determine which time scales contributed to the different complexity measures. We demonstrate that entropy rate and LZC measure the Kolmogorov complexity (randomness) of microstate sequences, whereas excess entropy and Hurst exponents describe statistical complexity which attains its maximum at intermediate levels of randomness. We confirmed the equivalence of entropy rate and LZC when the LZ-76 algorithm is used, a result previously reported for neural spike train analysis (Amigó et al., Neural Comput 16:717-736, https://doi.org/10.1162/089976604322860677 , 2004). Surrogate data analyses prove that entropy-based quantities and LZC focus on short-range temporal correlations, whereas Hurst exponents include short and long time scales. Sleep data analysis reveals that deeper sleep stages are accompanied by a decrease in Kolmogorov complexity and an increase in statistical complexity. Microstate jump sequences, where duplicate states have been removed, show higher randomness, lower statistical complexity, and no long-range correlations. Regarding the practical use of these methods, we suggest that LZC can be used as an efficient entropy rate estimator that avoids the estimation of joint entropies, whereas entropy rate estimation via joint entropies has the advantage of providing excess entropy as the second parameter of the same linear fit. We conclude that metrics of statistical complexity are a useful addition to microstate analysis and address a complexity concept that is not yet covered by existing microstate algorithms while being actively explored in other areas of brain research.

Understanding divergence: Placing developmental neuroscience in its dynamic context

Neurodevelopment is not merely a process of brain maturation, but an adaptation to constraints unique to each individual and to the environments we co-create. However, our theoretical and methodological toolkits often ignore this reality. There is growing awareness that a shift is needed that allows us to study divergence of brain and behaviour across conventional categorical boundaries. However, we argue that in future our study of divergence must also incorporate the developmental dynamics that capture the emergence of those neurodevelopmental differences. This crucial step will require adjustments in study design and methodology. If our ultimate aim is to incorporate the developmental dynamics that capture how, and ultimately when, divergence takes place then we will need an analytic toolkit equal to these ambitions. We argue that the over reliance on group averages has been a conceptual dead-end with regard to the neurodevelopmental differences. This is in part because any individual differences and developmental dynamics are inevitably lost within the group average. Instead, analytic approaches which are themselves new, or simply newly applied within this context, may allow us to shift our theoretical and methodological frameworks from groups to individuals. Likewise, methods capable of modelling complex dynamic systems may allow us to understand the emergent dynamics only possible at the level of an interacting neural system.

New Physico-Chemical Analysis of Magnesium-Doped Hydroxyapatite in Dextran Matrix Nanocomposites

The new magnesium-doped hydroxyapatite in dextran matrix (10MgHApD) nanocomposites were synthesized using coprecipitation technique. A spherical morphology was observed by scanning electron microscopy (SEM). The X-ray diffraction (XRD) characterization results show hydroxyapatite hexagonal phase formation. The element map scanning during the EDS analysis revealed homogenous distribution of constituent elements of calcium, phosphor, oxygen and magnesium. The presence of dextran in the sample was revealed by Fourier transform infrared (FTIR) spectroscopy. The antimicrobial activity of the 10MgHAPD nanocomposites was assessed by in vitro assays using Staphylococcus aureus ATCC 25923, Pseudomonas aeruginosa ATCC 27853, Streptococcus mutans ATCC 25175, Porphyromonas gingivalis ATCC 33277 and Candida albicans ATCC 10231 microbial strains. The results of the antimicrobial assays highlighted that the 10MgHApD nanocomposites presented excellent antimicrobial activity against all the tested microorganisms and for all the tested time intervals. Furthermore, the biocompatibility assays determined that the 10MgHApD nanocomposites did not exhibit any toxicity towards Human gingival fibroblast (HGF-1) cells.

Brain Network Analysis: A Review on Multivariate Analytical Methods

Despite the explosive growth of neuroimaging studies aimed at analyzing the brain as a complex system, critical methodological gaps remain to be addressed. Most tools currently used for analyzing network data of the brain are univariate in nature and are based on assumptions borne out of previous techniques not directly related to the big and complex data of the brain. Although graph-based methods have shown great promise, the development of principled multivariate models to address inherent limitations of graph-based methods, such as their dependence on network size and degree distributions, and to allow assessing the effects of multiple phenotypes on the brain and simulating brain networks has largely lagged behind. Although some studies have been made in developing multivariate frameworks to fill this gap, in the absence of a "gold-standard" method or guidelines, choosing the most appropriate method for each study can be another critical challenge for investigators in this multidisciplinary field. Here, we briefly introduce important multivariate methods for brain network analyses in two main categories: data-driven and model-based methods. We discuss whether/how such methods are suited for examining connectivity (edge-level), topology (system-level), or both. This review will aid in choosing an appropriate multivariate method with respect to variables such as network type, number of subjects and brain regions included, and the interest in connectivity, topology, or both. This review is aimed to be accessible to investigators from different backgrounds, with a focus on applications in brain network studies, though the methods may be applicable in other areas too. Impact statement As the U.S. National Institute of Health notes, the rich biomedical data can greatly improve our knowledge of human health if new analytical tools are developed, and their applications are broadly disseminated. A major challenge in analyzing the brain as a complex system is about developing parsimonious multivariate methods, and particularly choosing the most appropriate one among the existing methods with respect to the study variables in this multidisciplinary field. This study provides a review on the most important multivariate methods to aid in helping the most appropriate ones with respect to the desired variables for each study.

Age Related Functional Connectivity Signature Extraction Using Energy-Based Machine Learning Techniques

The study of brain connectivity plays an important role in understanding the functional organizations of the brain. It also helps to identify connectivity signatures that can be used for evaluating neural disorders and monitoring treatment efficacy. In this work, age-related changes in brain connectivity are studied to obtain aging signatures based on various modeling techniques. These include an energy-based machine learning technique to identify brain network interaction differences between two age groups with a large (30 years) age gap between them. Disconnectivity graphs and activation maps of the seven prominent resting-state networks (RSN) were obtained from functional MRI data of old and young adult subjects. Two-sample t-tests were performed on the local minimums with Bonferroni correction to control the family-wise error rate. These local minimums are connectivity states showing not only which brain regions but also how strong they are working together. They work as aging signatures that can be used to differentiate young and old groups. We found that the attention network's connectivity signature is a state with all the regions working together and young subjects have a stronger average connectivity among these regions. We have also found a common pattern between young and old subjects where the left and right brain regions of the frontal network are sometimes working separately instead of together. In summary, in this work, we combined machine learning and statistical approaches to extract connectivity signatures, which can be utilized to distinguish aging brains and monitor possible treatment efficacy.

Tissue-Engineered Models of the Human Brain: State-of-the-Art Analysis and Challenges

Neurological disorders affect billions of people across the world, making the discovery of effective treatments an important challenge. The evaluation of drug efficacy is further complicated because of the lack of in vitro models able to reproduce the complexity of the human brain structure and functions. Some limitations of 2D preclinical models of the human brain have been overcome by the use of 3D cultures such as cell spheroids, organoids and organs-on-chip. However, one of the most promising approaches for mimicking not only cell structure, but also brain architecture, is currently represented by tissue-engineered brain models. Both conventional (particularly electrospinning and salt leaching) and unconventional (particularly bioprinting) techniques have been exploited, making use of natural polymers or combinations between natural and synthetic polymers. Moreover, the use of induced pluripotent stem cells (iPSCs) has allowed the co-culture of different human brain cells (neurons, astrocytes, oligodendrocytes, microglia), helping towards approaching the central nervous system complexity. In this review article, we explain the importance of in vitro brain modeling, and present the main in vitro brain models developed to date, with a special focus on the most recent advancements in tissue-engineered brain models making use of iPSCs. Finally, we critically discuss achievements, main challenges and future perspectives.

Monofractal analysis of functional magnetic resonance imaging: An introductory review

The following review will aid readers in providing an overview of scale-free dynamics and monofractal analysis, as well as its applications and potential in functional magnetic resonance imaging (fMRI) neuroscience and clinical research. Like natural phenomena such as the growth of a tree or crashing ocean waves, the brain expresses scale-invariant, or fractal, patterns in neural signals that can be measured. While neural phenomena may represent both monofractal and multifractal processes and can be quantified with many different interrelated parameters, this review will focus on monofractal analysis using the Hurst exponent (H). Monofractal analysis of fMRI data is an advanced analysis technique that measures the complexity of brain signaling by quantifying its degree of scale-invariance. As such, the H value of the blood oxygenation level-dependent (BOLD) signal specifies how the degree of correlation in the signal may mediate brain functions. This review presents a brief overview of the theory of fMRI monofractal analysis followed by notable findings in the field. Through highlighting the advantages and challenges of the technique, the article provides insight into how to best conduct fMRI fractal analysis and properly interpret the findings with physiological relevance. Furthermore, we identify the future directions necessary for its progression towards impactful functional neuroscience discoveries and widespread clinical use. Ultimately, this presenting review aims to build a foundation of knowledge among readers to facilitate greater understanding, discussion, and use of this unique yet powerful imaging analysis technique.

Decoding functional brain networks through graph measures in infancy: The case of emotional faces

Graph measures represent an optimal way to investigate neural networks' organization, yet their application is still limited in developmental samples. To uncover the organization of 7-month-old infants' functional brain networks during an emotional perception task, we combined a decoding technique (i.e., Principal Component Regression) to graph metrics computation. Nodes' Within Module Degree Z Score (WMDZ) was computed as a measure of modular organization, and we decoded networks' functional organizations across EEG alpha and theta bands in response to static and dynamic facial expressions of emotions. We found that infants' brain topological activity differentiates between static and dynamic emotional faces due to the involvement of visual streams and sensorimotor areas, as often observed in adults. Moreover, network invariances point toward an already present rudimental network structure tuned to face processing already at 7-months of age. Overall, our results affirm the fruitfulness of the application of graph measures in developmental samples, due to their flexibility and the wealth of information they provide over infants' networks functional organization.

Neuronal Electrical Ongoing Activity as Cortical Areas Signature: An Insight from MNI Intracerebral Recording Atlas

The time course of the neuronal activity in the brain network, the neurodynamics, reflects the structure and functionality of the generating neuronal pools. Here, using the intracranial stereo-electroencephalographic (sEEG) recordings of the public Montreal Neurological Institute (MNI) atlas, we investigated the neurodynamics of primary motor (M1), somatosensory (S1) and auditory (A1) cortices measuring power spectral densities (PSD) and Higuchi fractal dimension (HFD) in the same subject (M1 vs. S1 in 16 subjects, M1 vs. A1 in 9, S1 vs. A1 in 6). We observed specific spectral features in M1, which prevailed above beta band, S1 in the alpha band, and A1 in the delta band. M1 HFD was higher than S1, both higher than A1. A clear distinction of neurodynamics properties of specific primary cortices supports the efforts in cortical parceling based on this expression of the local cytoarchitecture and connectivity. In this perspective, we selected within the MNI intracortical database a first set of primary motor, somatosensory and auditory cortices' representatives to query in recognizing ongoing patterns of neuronal communication. Potential clinical impact stands primarily in exploiting such exchange patterns to enhance the efficacy of neuromodulation intervention to cure symptoms secondary to neuronal activity unbalances.

Internet-like brain hierarchical network model: Alzheimer's disease study as an example

BACKGROUND AND OBJECTIVE

The modular structure and hierarchy are important topological characteristics in real complex networks such as brain networks on temporal scale. However, there are few studies investigating the hierarchical structure at the spatial scale of brain networks, the application of which still remains to be further studied.

METHODS

In this study, a novel model of brain hierarchical network based on the hierarchical characteristic of Internet topology is proposed for the first time, which is called Internet-like brain hierarchical network (IBHN). In this model, the whole brain network is partitioned into multiple levels: brain wide area network (Brain-WAN), brain metropolitan network (Brain-MAN), and brain local area network (Brain-LAN). A Brain-MAN is formed by the interconnection of multiple Brain-LANs, and the interconnection of multiple Brain-MANs forms a Brain-WAN. A multivariate analysis method is employed to measure overall functional connectivity between two brain networks at the same network level rather than detecting the change of each node pair's functional connection. Furthermore, we demonstrate the utility of IBHN model with application to a practical case-control study involving 64 patients with Alzheimer's disease and 75 healthy controls.

RESULTS

The proposed model identified enhanced functional connectivity (P-value<0.05) at Brain-WAN level and reduced functional connectivity (P-value=0.004) at Brain-LAN level of Alzheimer's disease patients, which can be used as a multi-dimension functional reference for AD's diagnosis.

CONCLUSIONS

This study not only provides insight into AD pathophysiology, but also further proves the effectiveness of the proposed IBHN model. In addition, the IBHN model makes it possible to explore the brain's functional organization from multiple dimensions and offers a multi-level perspective for the research of complex brain network.

Sensorimotor Activity and Network Connectivity to Dynamic and Static Emotional Faces in 7-Month-Old Infants

The present study investigated whether, as in adults, 7-month-old infants’ sensorimotor brain areas are recruited in response to the observation of emotional facial expressions. Activity of the sensorimotor cortex, as indexed by µ rhythm suppression, was recorded using electroencephalography (EEG) while infants observed neutral, angry, and happy facial expressions either in a static (N = 19) or dynamic (N = 19) condition. Graph theory analysis was used to investigate to which extent neural activity was functionally localized in specific cortical areas. Happy facial expressions elicited greater sensorimotor activation compared to angry faces in the dynamic experimental condition, while no difference was found between the three expressions in the static condition. Results also revealed that happy but not angry nor neutral expressions elicited a significant right-lateralized activation in the dynamic condition. Furthermore, dynamic emotional faces generated more efficient processing as they elicited higher global efficiency and lower networks’ diameter compared to static faces. Overall, current results suggest that, contrarily to neutral and angry faces, happy expressions elicit sensorimotor activity at 7 months and dynamic emotional faces are more efficiently processed by functional brain networks. Finally, current data provide evidence of the existence of a right-lateralized activity for the processing of happy facial expressions.

Separating scale-free and oscillatory components of neural activity in schizophrenia

INTRODUCTION

Alterations in narrow-band spectral power of electroencephalography (EEG) recordings are commonly reported in patients with schizophrenia (SZ). It is well established however that electrophysiological signals comprise a broadband scale-free (or fractal) component generated by mechanisms different from those producing oscillatory neural activity. Despite this known feature, it has not yet been investigated if spectral abnormalities found in SZ could be attributed to scale-free or oscillatory brain function.

METHODS

In this study, we analyzed resting-state EEG recordings of 14 SZ patients and 14 healthy controls. Scale-free and oscillatory components of the power spectral density (PSD) were separated, and band-limited power (BLP) of the original (mixed) PSD, as well as its fractal and oscillatory components, was estimated in five frequency bands. The scaling property of the fractal component was characterized by its spectral exponent in two distinct frequency ranges (1-13 and 13-30 Hz).

RESULTS

Analysis of the mixed PSD revealed a decrease of BLP in the delta band in SZ over the central regions; however, this difference could be attributed almost exclusively to a shift of power toward higher frequencies in the fractal component. Broadband neural activity expressed a true bimodal nature in all except frontal regions. Furthermore, both low- and high-range spectral exponents exhibited a characteristic topology over the cortex in both groups.

CONCLUSION

Our results imply strong functional significance of scale-free neural activity in SZ and suggest that abnormalities in PSD may emerge from alterations of the fractal and not only the oscillatory components of neural activity.

A Complex Systems Perspective on Neuroimaging Studies of Behavior and Its Disorders

The study of complex systems deals with emergent behavior that arises as a result of nonlinear spatiotemporal interactions between a large number of components both within the system, as well as between the system and its environment. There is a strong case to be made that neural systems as well as their emergent behavior and disorders can be studied within the framework of complexity science. In particular, the field of neuroimaging has begun to apply both theoretical and experimental procedures originating in complexity science—usually in parallel with traditional methodologies. Here, we illustrate the basic properties that characterize complex systems and evaluate how they relate to what we have learned about brain structure and function from neuroimaging experiments. We then argue in favor of adopting a complex systems-based methodology in the study of neuroimaging, alongside appropriate experimental paradigms, and with minimal influences from noncomplex system approaches. Our exposition includes a review of the fundamental mathematical concepts, combined with practical examples and a compilation of results from the literature.

Decreased connection density and modularity of functional brain networks during n-back working memory paradigm

INTRODUCTION

Investigating how the brain adapts to increased mental workload through large-scale functional reorganization appears as an important research question. Functional connectivity (FC) aims at capturing how disparate regions of the brain dynamically interact, while graph theory provides tools for the topological characterization of the reconstructed functional networks. Although numerous studies investigated how FC is altered in response to increased working memory (WM) demand, current results are still contradictory as few studies confirmed the robustness of these findings in a low-density setting.

METHODS

In this study, we utilized the n-back WM paradigm, in which subjects were presented stimuli (single digits) sequentially, and their task was to decide for each given stimulus if it matched the one presented n-times earlier. Electroencephalography recordings were performed under a control (0-back) and two task conditions of varying difficulty (2- and 3-back). We captured the characteristic connectivity patterns for each difficulty level by performing FC analysis and described the reconstructed functional networks with various graph theoretical measures.

RESULTS

We found a substantial decrease in FC when transitioning from the 0- to the 2- or 3-back conditions, however, no differences relating to task difficulty were identified. The observed changes in brain network topology could be attributed to the dissociation of two (frontal and occipitotemporal) functional modules that were only present during the control condition. Furthermore, behavioral and performance measures showed both positive and negative correlations to connectivity indices, although only in the higher frequency bands.

CONCLUSION

The marked decrease in FC may be due to temporarily abandoned connections that are redundant or irrelevant in solving the specific task. Our results indicate that FC analysis is a robust tool for investigating the response of the brain to increased cognitive workload.

Variance and Scale-Free Properties of Resting-State Blood Oxygenation Level-Dependent Signal After Fear Memory Acquisition and Extinction

Recently, the dynamic properties of brain activity rather than its stationary values have attracted more interest in clinical applications. It has been shown that brain signals exhibit scale-free dynamics or long-range temporal correlations (LRTC) that differ between rest and cognitive tasks in healthy controls and clinical groups. Little is known about how fear-inducing tasks may influence dispersion and the LRTC of subsequent resting-state brain activity. In this study, we aimed to explore the changes in the variance and scale-free properties of the brain’s blood oxygenation level-dependent (BOLD) signal during the resting-state sessions before and after fear learning and fear memory extinction. During a 1-h break between magnetic resonance imaging (MRI) scanning, 23 healthy, right-handed volunteers experienced a fear extinction procedure, followed by Pavlovian fear conditioning that included partial reinforcement using mild electrical stimulation. We extracted the average time course of the BOLD signal from 245 regions of interest (ROIs) taken from the resting-state functional atlas. The variance of the BOLD signal and the Hurst exponent (H), which reflects the scale-free dynamic, were compared in the resting states before and after fear learning and fear memory extinction. After fear extinction, six ROIs showed a difference in H at the uncorrected level of significance, including areas associated with fear processing. H decreased during fear extinction but then became higher than before fear learning, specifically in areas related to the fear extinction network (FEN). However, activity in the other ROIs restored the H to its initial level. The variance of the BOLD signal in six ROIs demonstrated a significant increase from initial rest to the post-task rest. A limited number of ROIs showed changes in both H and variance. Our results imply that the variability and scale-free properties of the BOLD signal might serve as additional indicators of changes in spontaneous brain activity related to recent experience.

Multifractal and Entropy-Based Analysis of Delta Band Neural Activity Reveals Altered Functional Connectivity Dynamics in Schizophrenia

Dynamic functional connectivity (DFC) was established in the past decade as a potent approach to reveal non-trivial, time-varying properties of neural interactions – such as their multifractality or information content –, that otherwise remain hidden from conventional static methods. Several neuropsychiatric disorders were shown to be associated with altered DFC, with schizophrenia (SZ) being one of the most intensely studied among such conditions. Here we analyzed resting-state electroencephalography recordings of 14 SZ patients and 14 age- and gender-matched healthy controls (HC). We reconstructed dynamic functional networks from delta band (0.5–4 Hz) neural activity and captured their spatiotemporal dynamics in various global network topological measures. The acquired network measure time series were made subject to dynamic analyses including multifractal analysis and entropy estimation. Besides group-level comparisons, we built a classifier to explore the potential of DFC features in classifying individual cases. We found stronger delta-band connectivity, as well as increased variance of DFC in SZ patients. Surrogate data testing verified the true multifractal nature of DFC in SZ, with patients expressing stronger long-range autocorrelation and degree of multifractality when compared to controls. Entropy analysis indicated reduced temporal complexity of DFC in SZ. When using these indices as features, an overall cross-validation accuracy surpassing 89% could be achieved in classifying individual cases. Our results imply that dynamic features of DFC such as its multifractal properties and entropy are potent markers of altered neural dynamics in SZ and carry significant potential not only in better understanding its pathophysiology but also in improving its diagnosis. The proposed framework is readily applicable for neuropsychiatric disorders other than schizophrenia.

Additional Insights into the Relationship Between Brain Network Architecture and Susceptibility and Resilience to the Psychiatric Sequelae of Childhood Maltreatment

Childhood maltreatment is associated with increased risk for psychiatric and substance use disorders. However, some maltreated individuals appear resilient to these consequences while manifesting the same array of brain changes as maltreated individuals with psychopathology. Hence, a critical issue has been to identify compensatory brain alterations in these resilient individuals. We recently reported that maltreatment is associated with a more vulnerable structural brain network architecture. Resilient individuals have this same vulnerability but appeared to be able to effectively compensate due to reduced nodal efficiency (ability of a node to influence the global network) in 9 specific brain regions that moderate the relationship between maltreatment and psychopathology. Following up we now report that network vulnerability increases progressively during late adolescence to plateau at about 21 years of age, which may help to explain age of onset of psychopathology. Further, we found that network vulnerability was most significantly affected by parental verbal abuse between 16-18 years of age and number of types of maltreatment during childhood. Asymptomatic individuals with no history of psychopathology had more prominent alterations in nodal efficiency than asymptomatic individuals with prior history, who specifically showed reduced nodal efficiency in right amygdala and right subcallosal gyrus. Experiencing inadequate financial sufficiency during childhood increased risk of susceptibility versus resilience by 2.98-fold (95% CI 1.49-5.97, p = 0.002) after adjusting for differences in exposure to maltreatment. Interestingly, adequate-to-higher financial sufficiency appeared to be protective and was associated with reduced nodal efficiency in the right postcentral gyrus and subcallosal gyrus 'resilience' nodes.

Connectivity differences between Gulf War Illness (GWI) phenotypes during a test of attention

One quarter of veterans returning from the 1990–1991 Persian Gulf War have developed Gulf War Illness (GWI) with chronic pain, fatigue, cognitive and gastrointestinal dysfunction. Exertion leads to characteristic, delayed onset exacerbations that are not relieved by sleep. We have modeled exertional exhaustion by comparing magnetic resonance images from before and after submaximal exercise. One third of the 27 GWI participants had brain stem atrophy and developed postural tachycardia after exercise (START: Stress Test Activated Reversible Tachycardia). The remainder activated basal ganglia and anterior insulae during a cognitive task (STOPP: Stress Test Originated Phantom Perception). Here, the role of attention in cognitive dysfunction was assessed by seed region correlations during a simple 0-back stimulus matching task (“see a letter, push a button”) performed before exercise. Analysis was analogous to resting state, but different from psychophysiological interactions (PPI). The patterns of correlations between nodes in task and default networks were significantly different for START (n = 9), STOPP (n = 18) and control (n = 8) subjects. Edges shared by the 3 groups may represent co-activation caused by the 0-back task. Controls had a task network of right dorsolateral and left ventrolateral prefrontal cortex, dorsal anterior cingulate cortex, posterior insulae and frontal eye fields (dorsal attention network). START had a large task module centered on the dorsal anterior cingulate cortex with direct links to basal ganglia, anterior insulae, and right dorsolateral prefrontal cortex nodes, and through dorsal attention network (intraparietal sulci and frontal eye fields) nodes to a default module. STOPP had 2 task submodules of basal ganglia–anterior insulae, and dorsolateral prefrontal executive control regions. Dorsal attention and posterior insulae nodes were embedded in the default module and were distant from the task networks. These three unique connectivity patterns during an attention task support the concept of Gulf War Disease with recognizable, objective patterns of cognitive dysfunction.

How Big Data and High-performance Computing Drive Brain Science

Brain science accelerates the study of intelligence and behavior, contributes fundamental insights into human cognition, and offers prospective treatments for brain disease. Faced with the challenges posed by imaging technologies and deep learning computational models, big data and high-performance computing (HPC) play essential roles in studying brain function, brain diseases, and large-scale brain models or connectomes. We review the driving forces behind big data and HPC methods applied to brain science, including deep learning, powerful data analysis capabilities, and computational performance solutions, each of which can be used to improve diagnostic accuracy and research output. This work reinforces predictions that big data and HPC will continue to improve brain science by making ultrahigh-performance analysis possible, by improving data standardization and sharing, and by providing new neuromorphic insights.

Increased segregation of functional networks in developing brains

A growing literature conceptualises typical brain development from a network perspective. However, largely due to technical and methodological challenges inherent in paediatric functional neuroimaging, there remains an important gap in our knowledge regarding the typical development of functional brain networks in "preschool" childhood (i.e., children younger than 6 years of age). In this study, we recorded brain oscillatory activity using age-appropriate magnetoencephalography in 24 children, including 14 preschool children aged from 4 to 6 years and 10 school children aged from 7 to 12 years. We compared the topology of the resting-state brain networks in these children, estimated using minimum spanning tree (MST) constructed from phase synchrony between beamformer-reconstructed time-series, with that of 24 adults. Our results show that during childhood the MST topology shifts from a star-like (centralised) toward a more line-like (de-centralised) configuration, indicating the functional brain networks become increasingly segregated. In addition, the increasing global network segregation is frequency-independent and accompanied by decreases in centrality (or connectedness) of cortical regions with age, especially in areas of the default mode network. We propose a heuristic MST model of "network space", which posits a clear developmental trajectory for the emergence of complex brain networks. Our results not only revealed topological reorganisation of functional networks across multiple temporal and spatial scales in childhood, but also fill a gap in the literature regarding neurophysiological mechanisms of functional brain maturation during the preschool years of childhood.

EEG-Based Prediction of Cognitive Load in Intelligence Tests

Measuring and assessing the cognitive load associated with different tasks is crucial for many applications, from the design of instructional materials to monitoring the mental well-being of aircraft pilots. The goal of this paper is to utilize EEG to infer the cognitive workload of subjects during intelligence tests. We chose the well established advanced progressive matrices test, an ideal framework because it presents problems at increasing levels of difficulty and has been rigorously validated in past experiments. We train classic machine learning models using basic EEG measures as well as measures of network connectivity and signal complexity. Our findings demonstrate that cognitive load can be well predicted using these features, even for a low number of channels. We show that by creating an individually tuned neural network for each subject, we can improve prediction compared to a generic model and that such models are robust to decreasing the number of available channels as well.

Multi-scale detection of hierarchical community architecture in structural and functional brain networks

Community detection algorithms have been widely used to study the organization of complex networks like the brain. These techniques provide a partition of brain regions (or nodes) into clusters (or communities), where nodes within a community are densely interconnected with one another. In their simplest application, community detection algorithms are agnostic to the presence of community hierarchies: clusters embedded within clusters of other clusters. To address this limitation, we exercise a multi-scale extension of a common community detection technique, and we apply the tool to synthetic graphs and to graphs derived from human neuroimaging data, including structural and functional imaging data. Our multi-scale community detection algorithm links a graph to copies of itself across neighboring topological scales, thereby becoming sensitive to conserved community organization across levels of the hierarchy. We demonstrate that this method is sensitive to topological inhomogeneities of the graph's hierarchy by providing a local measure of community stability and inter-scale reliability across topological scales. We compare the brain's structural and functional network architectures, and we demonstrate that structural graphs display a more prominent hierarchical community organization than functional graphs. Finally, we build an explicitly multimodal multiplex graph that combines both structural and functional connectivity in a single model, and we identify the topological scales where resting state functional connectivity and underlying structural connectivity show similar versus unique hierarchical community architecture. Together, our results demonstrate the advantages of the multi-scale community detection algorithm in studying hierarchical community structure in brain graphs, and they illustrate its utility in modeling multimodal neuroimaging data.

Susceptibility or Resilience to Maltreatment Can Be Explained by Specific Differences in Brain Network Architecture

BACKGROUND

Childhood maltreatment is a major risk factor for psychopathology. However, some maltreated individuals appear remarkably resilient to the psychiatric effects while manifesting the same array of brain abnormalities as maltreated individuals with psychopathology. Hence, a critical aim is to identify compensatory brain alterations that enable resilient individuals to maintain mental well-being despite alterations in stress-susceptible regions.

METHODS

Network models were constructed from diffusion tensor imaging and tractography in physically healthy unmedicated 18- to 25-year-old participants (N = 342, n = 192 maltreated) to develop network-based explanatory models.

RESULTS

First, we determined that susceptible and resilient individuals had the same alterations in global fiber stream network architecture using two different definitions of resilience: 1) no lifetime history of Axis I or II disorders, and 2) no clinically significant symptoms of anxiety, depression, anger-hostility, or somatization. Second, we confirmed an a priori hypothesis that right amygdala nodal efficiency was lower in asymptomatic resilient than in susceptible participants or control subjects. Third, we identified eight other nodes with reduced nodal efficiency in resilient individuals and showed that nodal efficiency moderated the relationship between maltreatment and psychopathology. Fourth, we found that models based on global network architecture and nodal efficiency could delineate group membership (control, susceptible, resilient) with 75%, 82%, and 80% cross-validated accuracy.

CONCLUSIONS

Together these findings suggest that sparse fiber networks with increased small-worldness following maltreatment render individuals vulnerable to psychopathology if abnormalities occur in specific nodes, but that decreased ability of certain nodes to propagate information throughout the network mitigates the effects and leads to resilience.

Structural covariance and cortical reorganisation in schizophrenia: a MRI-based morphometric study

BACKGROUND

In patients with schizophrenia, distributed abnormalities are observed in grey matter volume. A recent hypothesis posits that these distributed changes are indicative of a plastic reorganisation process occurring in response to a functional defect in neuronal information transmission. We investigated the structural covariance across various brain regions in early-stage schizophrenia to determine if indeed the observed patterns of volumetric loss conform to a coordinated pattern of structural reorganisation.

METHODS

Structural magnetic resonance imaging scans were obtained from 40 healthy adults and 41 age, gender and parental socioeconomic status matched patients with schizophrenia. Volumes of grey matter tissue were estimated at the regional level across 90 atlas-based parcellations. Group-level structural covariance was studied using a graph theoretical framework.

RESULTS

Patients had distributed reduction in grey matter volume, with high degree of localised covariance (clustering) compared with controls. Patients with schizophrenia had reduced centrality of anterior cingulate and insula but increased centrality of the fusiform cortex, compared with controls. Simulating targeted removal of highly central nodes resulted in significant loss of the overall covariance patterns in patients compared with controls.

CONCLUSION

Regional volumetric deficits in schizophrenia are not a result of random, mutually independent processes. Our observations support the occurrence of a spatially interconnected reorganisation with the systematic de-escalation of conventional 'hub' regions. This raises the question of whether the morphological architecture in schizophrenia is primed for compensatory functions, albeit with a high risk of inefficiency.

The power of language: Functional brain network topology of deaf and hearing in relation to sign language experience

Prolonged auditory sensory deprivation leads to brain reorganization. This is indicated by functional enhancement in remaining sensory systems and known as cross-modal plasticity. In this study we investigated differences in functional brain network topology between deaf and hearing individuals. We also studied altered functional network responses between deaf and hearing individuals with a recording paradigm containing an eyes-closed and eyes-open condition. Electroencephalography activity was recorded in a group of sign language-trained deaf (N = 71) and hearing people (N = 122) living in rural Africa. Functional brain networks were constructed from the functional connectivity between fourteen electrodes distributed over the scalp. Functional connectivity was quantified with the phase lag index based on bandpass filtered epochs of brain signal. We studied the functional connectivity between the auditory, somatosensory and visual cortex and performed whole-brain minimum spanning tree analysis to capture network backbone characteristics. Functional connectivity between different regions involved in sensory information processing tended to be stronger in deaf people during the eyes-closed condition in both the alpha and beta frequency band. Furthermore, we found differences in functional backbone topology between deaf and hearing individuals. The backbone topology altered during transition from the eyes-closed to eyes-open condition irrespective of deafness, but was more pronounced in deaf individuals. The transition of backbone strength was different between individuals with congenital, pre-lingual or post-lingual deafness. Functional backbone characteristics correlated with the experience of sign language. Overall, our study revealed more insights in functional network reorganization caused by auditory deprivation and cross-modal plasticity. It further supports the idea of a brain plasticity potential in deaf and hearing people. The association between network organization and acquired sign language experience reflects the ability of ongoing brain adaptation in people with hearing disabilities.

Fractal and Multifractal Properties of Electrographic Recordings of Human Brain Activity: Toward Its Use as a Signal Feature for Machine Learning in Clinical Applications

The quantification of brain dynamics is essential to its understanding. However, the brain is a system operating on multiple time scales, and characterization of dynamics across time scales remains a challenge. One framework to study such dynamics is that of fractal geometry; and currently there exist several methods for the study of brain dynamics using fractal geometry. We aim to highlight some of the practical challenges of applying fractal geometry to brain dynamics—and as a putative feature for machine learning applications, and propose solutions to enable its wider use in neuroscience. Using intracranially recorded electroencephalogram (EEG) and simulated data, we compared monofractal and multifractal methods with regards to their sensitivity to signal variance. We found that both monofractal and multifractal properties correlate closely with signal variance, thus not being a useful feature of the signal. However, after applying an epoch-wise standardization procedure to the signal, we found that multifractal measures could offer non-redundant information compared to signal variance, power (in different frequency bands) and other established EEG signal measures. We also compared different multifractal estimation methods to each other in terms of reliability, and we found that the Chhabra-Jensen algorithm performed best. Finally, we investigated the impact of sampling frequency and epoch length on the estimation of multifractal properties. Using epileptic seizures as an example event in the EEG, we show that there may be an optimal time scale (i.e., combination of sampling frequency and epoch length) for detecting temporal changes in multifractal properties around seizures. The practical issues we highlighted and our suggested solutions should help in developing robust methods for the application of fractal geometry in EEG signals. Our analyses and observations also aid the theoretical understanding of the multifractal properties of the brain and might provide grounds for new discoveries in the study of brain signals. These could be crucial for the understanding of neurological function and for the developments of new treatments.

Fractal Analysis of BOLD Time Series in a Network Associated With Waiting Impulsivity

Fractal phenomena can be found in numerous scientific areas including neuroscience. Fractals are structures, in which the whole has the same shape as its parts. A specific structure known as pink noise (also called fractal or 1/f noise) is one key fractal manifestation, exhibits both stability and adaptability, and can be addressed via the Hurst exponent (H). FMRI studies using H on regional fMRI time courses used fractality as an important characteristic to unravel neural networks from artificial noise. In this fMRI-study, we examined 103 healthy male students at rest and while performing the 5-choice serial reaction time task. We addressed fractality in a network associated with waiting impulsivity using the adaptive fractal analysis (AFA) approach to determine H. We revealed the fractal nature of the impulsivity network. Furthermore, fractality was influenced by individual impulsivity in terms of decreasing fractality with higher impulsivity in regions of top-down control (left middle frontal gyrus) as well as reward processing (nucleus accumbens and anterior cingulate cortex). We conclude that fractality as determined via H is a promising marker to quantify deviations in network functions at an early stage and, thus, to be able to inform preventive interventions before the manifestation of a disorder.

On the nature and use of models in network neuroscience

Network theory provides an intuitively appealing framework for studying relationships among interconnected brain mechanisms and their relevance to behaviour. As the space of its applications grows, so does the diversity of meanings of the term network model. This diversity can cause confusion, complicate efforts to assess model validity and efficacy, and hamper interdisciplinary collaboration. In this Review, we examine the field of network neuroscience, focusing on organizing principles that can help overcome these challenges. First, we describe the fundamental goals in constructing network models. Second, we review the most common forms of network models, which can be described parsimoniously along the following three primary dimensions: from data representations to first-principles theory; from biophysical realism to functional phenomenology; and from elementary descriptions to coarse-grained approximations. Third, we draw on biology, philosophy and other disciplines to establish validation principles for these models. We close with a discussion of opportunities to bridge model types and point to exciting frontiers for future pursuits.

Impact of Healthy Aging on Multifractal Hemodynamic Fluctuations in the Human Prefrontal Cortex

Fluctuations in resting-state cerebral hemodynamics show scale-free behavior over two distinct scaling ranges. Changes in such bimodal (multi) fractal pattern give insight to altered cerebrovascular or neural function. Our main goal was to assess the distribution of local scale-free properties characterizing cerebral hemodynamics and to disentangle the influence of aging on these multifractal parameters. To this end, we obtained extended resting-state records (N = 2¹⁴) of oxyhemoglobin (HbO), deoxyhemoglobin (HbR) and total hemoglobin (HbT) concentration time series with continuous-wave near-infrared spectroscopy technology from the brain cortex. 52 healthy volunteers were enrolled in this study: 24 young (30.6 ± 8.2 years), and 28 elderly (60.5 ± 12.0 years) subjects. Using screening tests on power-law, multifractal noise, and shuffled data sets we evaluated the presence of true multifractal hemodynamics reflecting long-range correlation (LRC). Subsequently, scaling-range adaptive bimodal signal summation conversion (SSC) was performed based on standard deviation (σ) of signal windows across a range of temporal scales (s). Building on moments of different order (q) of the measure, σ(s), multifractal SSC yielded generalized Hurst exponent function, H(q), and singularity spectrum, D(h) separately for a fast and slow component (the latter dominating the highest temporal scales). Parameters were calculated reflecting the estimated measure at s = N (focus), degree of LRC [Hurst exponent, H(2) and maximal Hölder exponent, h_max] and measuring strength of multifractality [full-width-half-maximum of D(h) and ΔH₁₅ = H(−15)−H(15)]. Correlation-based signal improvement (CBSI) enhanced our signal in terms of interpreting changes due to neural activity or local/systemic hemodynamic influences. We characterized the HbO-HbR relationship with the aid of fractal scale-wise correlation coefficient, r_σ(s) and SSC-based multifractal covariance analysis. In the majority of subjects, cerebral hemodynamic fluctuations proved bimodal multifractal. In case of slow component of raw HbT, h_max, and Ĥ(2) were lower in the young group explained by a significantly increased r_σ(s) among elderly at high temporal scales. Regarding the fast component of CBSI-pretreated HbT and that of HbO-HbR covariance, h_max, and focus were decreased in the elderly group. These observations suggest an attenuation of neurovascular coupling reflected by a decreased autocorrelation of the neuronal component concomitant with an accompanying increased autocorrelation of the non-neuronal component in the elderly group.

Concurrent Changes of Brain Functional Connectivity and Motor Variability When Adapting to Task Constraints

In behavioral neuroscience, the adaptability of humans facing different constraints has been addressed on one side at the brain level, where a variety of functional networks dynamically support the same performance, and on the other side at the behavioral level, where fractal properties in sensorimotor variables have been considered as a hallmark of adaptability. To bridge the gap between the two levels of observation, we have jointly investigated the changes of network connectivity in the sensorimotor cortex assessed by modularity analysis and the properties of motor variability assessed by multifractal analysis during a prolonged tapping task. Four groups of participants had to produce the same tapping performance while being deprived from 0, 1, 2, or 3 sensory feedbacks simultaneously (auditory and/or visual and/or tactile). Whereas tapping performance was not statistically different across groups, the number of brain networks involved and the degree of multifractality of the inter-tap interval series were significantly correlated, increasing as a function of feedback deprivation. Our findings provide first evidence that concomitant changes in brain modularity and multifractal properties characterize adaptations underlying unchanged performance. We discuss implications of our findings with respect to the degeneracy properties of complex systems, and the entanglement of adaptability and effective adaptation.

Predicting future learning from baseline network architecture

Human behavior and cognition result from a complex pattern of interactions between brain regions. The flexible reconfiguration of these patterns enables behavioral adaptation, such as the acquisition of a new motor skill. Yet, the degree to which these reconfigurations depend on the brain's baseline sensorimotor integration is far from understood. Here, we asked whether spontaneous fluctuations in sensorimotor networks at baseline were predictive of individual differences in future learning. We analyzed functional MRI data from 19 participants prior to six weeks of training on a new motor skill. We found that visual-motor connectivity was inversely related to learning rate: sensorimotor autonomy at baseline corresponded to faster learning in the future. Using three additional scans, we found that visual-motor connectivity at baseline is a relatively stable individual trait. These results suggest that individual differences in motor skill learning can be predicted from sensorimotor autonomy at baseline prior to task execution.

Temporal fractal analysis of the rs-BOLD signal identifies brain abnormalities in autism spectrum disorder

Background

Brain connectivity in autism spectrum disorders (ASD) has proven difficult to characterize due to the heterogeneous nature of the spectrum. Connectivity in the brain occurs in a complex, multilevel and multi-temporal manner, driving the fluctuations observed in local oxygen demand. These fluctuations can be characterized as fractals, as they auto-correlate at different time scales. In this study, we propose a model-free complexity analysis based on the fractal dimension of the rs-BOLD signal, acquired with magnetic resonance imaging. The fractal dimension can be interpreted as measure of signal complexity and connectivity. Previous studies have suggested that reduction in signal complexity can be associated with disease. Therefore, we hypothesized that a detectable difference in rs-BOLD signal complexity could be observed between ASD patients and Controls.

Methods and findings

Anatomical and functional data from fifty-five subjects with ASD (12.7 ± 2.4 y/o) and 55 age-matched (14.1 ± 3.1 y/o) healthy controls were accessed through the NITRC database and the ABIDE project. Subjects were scanned using a 3T GE Signa MRI and a 32-channel RF-coil. Axial FSPGR-3D images were used to prescribe rs-BOLD (TE/TR = 30/2000ms) where 300 time points were acquired. Motion correction was performed on the functional data and anatomical and functional images were aligned and spatially warped to the N27 standard brain atlas. Fractal analysis, performed on a grey matter mask, was done by estimating the Hurst exponent in the frequency domain using a power spectral density approach and refining the estimation in the time domain with de-trended fluctuation analysis and signal summation conversion methods. Voxel-wise fractal dimension (FD) was calculated for every subject in the control group and in the ASD group to create ROI-based Z-scores for the ASD patients. Voxel-wise validation of FD normality across controls was confirmed, and non-Gaussian voxels were eliminated from subsequent analysis. To maintain a 95% confidence level, only regions where Z-score values were at least 2 standard deviations away from the mean (i.e. where |Z| > 2.0) were included in the analysis. We found that the main regions, where signal complexity significantly decreased among ASD patients, were the amygdala (p = 0.001), the vermis (p = 0.02), the basal ganglia (p = 0.01) and the hippocampus (p = 0.02). No regions reported significant increase in signal complexity in this study. Our findings were correlated with ADIR and ADOS assessment tools, reporting the highest correlation with the ADOS metrics.

Conclusions

Brain connectivity is best modeled as a complex system. Therefore, a measure of complexity as the fractal dimension of fluctuations in brain oxygen demand and utilization could provide important information about connectivity issues in ASD. Moreover, this technique can be used in the characterization of a single subject, with respect to controls, without the need for group analysis. Our novel approach provides an ideal avenue for personalized diagnostics, thus providing unique patient specific assessment that could help in individualizing treatments.

Decomposing Multifractal Crossovers

Physiological processes-such as, the brain's resting-state electrical activity or hemodynamic fluctuations-exhibit scale-free temporal structuring. However, impacts common in biological systems such as, noise, multiple signal generators, or filtering by transport function, result in multimodal scaling that cannot be reliably assessed by standard analytical tools that assume unimodal scaling. Here, we present two methods to identify breakpoints or crossovers in multimodal multifractal scaling functions. These methods incorporate the robust iterative fitting approach of the focus-based multifractal formalism (FMF). The first approach (moment-wise scaling range adaptivity) allows for a breakpoint-based adaptive treatment that analyzes segregated scale-invariant ranges. The second method (scaling function decomposition method, SFD) is a crossover-based design aimed at decomposing signal constituents from multimodal scaling functions resulting from signal addition or co-sampling, such as, contamination by uncorrelated fractals. We demonstrated that these methods could handle multimodal, mono- or multifractal, and exact or empirical signals alike. Their precision was numerically characterized on ideal signals, and a robust performance was demonstrated on exemplary empirical signals capturing resting-state brain dynamics by near infrared spectroscopy (NIRS), electroencephalography (EEG), and blood oxygen level-dependent functional magnetic resonance imaging (fMRI-BOLD). The NIRS and fMRI-BOLD low-frequency fluctuations were dominated by a multifractal component over an underlying biologically relevant random noise, thus forming a bimodal signal. The crossover between the EEG signal components was found at the boundary between the δ and θ bands, suggesting an independent generator for the multifractal δ rhythm. The robust implementation of the SFD method should be regarded as essential in the seamless processing of large volumes of bimodal fMRI-BOLD imaging data for the topology of multifractal metrics free of the masking effect of the underlying random noise.

Emerging Frontiers of Neuroengineering: A Network Science of Brain Connectivity

Neuroengineering is faced with unique challenges in repairing or replacing complex neural systems that are composed of many interacting parts. These interactions form intricate patterns over large spatiotemporal scales and produce emergent behaviors that are difficult to predict from individual elements. Network science provides a particularly appropriate framework in which to study and intervene in such systems by treating neural elements (cells, volumes) as nodes in a graph and neural interactions (synapses, white matter tracts) as edges in that graph. Here, we review the emerging discipline of network neuroscience, which uses and develops tools from graph theory to better understand and manipulate neural systems from micro- to macroscales. We present examples of how human brain imaging data are being modeled with network analysis and underscore potential pitfalls. We then highlight current computational and theoretical frontiers and emphasize their utility in informing diagnosis and monitoring, brain-machine interfaces, and brain stimulation. A flexible and rapidly evolving enterprise, network neuroscience provides a set of powerful approaches and fundamental insights that are critical for the neuroengineer's tool kit.

Fluid and flexible minds: Intelligence reflects synchrony in the brain's intrinsic network architecture

Human intelligence has been conceptualized as a complex system of dissociable cognitive processes, yet studies investigating the neural basis of intelligence have typically emphasized the contributions of discrete brain regions or, more recently, of specific networks of functionally connected regions. Here we take a broader, systems perspective in order to investigate whether intelligence is an emergent property of synchrony within the brain’s intrinsic network architecture. Using a large sample of resting-state fMRI and cognitive data (n = 830), we report that the synchrony of functional interactions within and across distributed brain networks reliably predicts fluid and flexible intellectual functioning. By adopting a whole-brain, systems-level approach, we were able to reliably predict individual differences in human intelligence by characterizing features of the brain’s intrinsic network architecture. These findings hold promise for the eventual development of neural markers to predict changes in intellectual function that are associated with neurodevelopment, normal aging, and brain disease.

In our study, we aimed to understand how individual differences in intellectual functioning are reflected in the intrinsic network architecture of the human brain. We applied statistical methods, known as spectral decompositions, in order to identify individual differences in the synchronous patterns of spontaneous brain activity that reliably predict core aspects of human intelligence. The synchrony of brain activity at rest across multiple discrete neural networks demonstrated positive relationships with fluid intelligence. In contrast, global synchrony within the brain’s network architecture reliably, and inversely, predicted mental flexibility, a core facet of intellectual functioning. The multinetwork systems approach described here represents a methodological and conceptual extension of earlier efforts that related differences in intellectual ability to variations in specific brain regions, networks, or their interactions. Our findings suggest that the neural basis of complex, integrative cognitive functions can be most completely understood from the perspective of network neuroscience.

A network engineering perspective on probing and perturbing cognition with neurofeedback

Network science and engineering provide a flexible and generalizable tool set to describe and manipulate complex systems characterized by heterogeneous interaction patterns among component parts. While classically applied to social systems, these tools have recently proven to be particularly useful in the study of the brain. In this review, we describe the nascent use of these tools to understand human cognition, and we discuss their utility in informing the meaningful and predictable perturbation of cognition in combination with the emerging capabilities of neurofeedback. To blend these disparate strands of research, we build on emerging conceptualizations of how the brain functions (as a complex network) and how we can develop and target interventions or modulations (as a form of network control). We close with an outline of current frontiers that bridge neurofeedback, connectomics, and network control theory to better understand human cognition.

Applying fMRI complexity analyses to the single subject: a case study for proposed neurodiagnostics

Nonlinear dynamic tools have been statistically validated at the group level to identify subtle differences in system wide regulation of brain meso-circuits, often increasing clinical sensitivity over conventional analyses alone. We explored the feasibility of extracting information at the single-subject level, illustrating two pairs of healthy individuals with psychological differences in stress reactivity. We applied statistical and nonlinear dynamic tools to capture key characteristics of the prefrontal-limbic loop. We compared single subject results with statistical results for the larger group. We concluded that complexity analyses may identify important differences at the single-subject level, supporting their potential towards neurodiagnostic applications.

Childhood maltreatment is associated with alteration in global network fiber-tract architecture independent of history of depression and anxiety

Childhood maltreatment is a major risk factor for psychopathology. It is also associated with alterations in the network architecture of the brain, which we hypothesized may play a significant role in the development of psychopathology. In this study, we analyzed the global network architecture of physically healthy unmedicated 18-25 year old subjects (n=262) using diffusion tensor imaging (DTI) MRI and tractography. Anatomical networks were constructed from fiber streams interconnecting 90 cortical or subcortical regions for subjects with no-to-low (n=122) versus moderate-to-high (n=140) exposure to maltreatment. Graph theory analysis revealed lower degree, strength, global efficiency, and maximum Laplacian spectra, higher pathlength, small-worldness and Laplacian skewness, and less deviation from artificial networks in subjects with moderate-to-high exposure to maltreatment. On balance, local clustering was similar in both groups, but the different clusters were more strongly interconnected in the no-to-low exposure group. History of major depression, anxiety and attention deficit hyperactivity disorder did not have a significant impact on global network measures over and above the effect of maltreatment. Maltreatment is an important factor that needs to be taken into account in studies examining the relationship between network differences and psychopathology.

Fractal Analysis of Brain Blood Oxygenation Level Dependent (BOLD) Signals from Children with Mild Traumatic Brain Injury (mTBI)

BACKGROUND

Conventional imaging techniques are unable to detect abnormalities in the brain following mild traumatic brain injury (mTBI). Yet patients with mTBI typically show delayed response on neuropsychological evaluation. Because fractal geometry represents complexity, we explored its utility in measuring temporal fluctuations of brain resting state blood oxygen level dependent (rs-BOLD) signal. We hypothesized that there could be a detectable difference in rs-BOLD signal complexity between healthy subjects and mTBI patients based on previous studies that associated reduction in signal complexity with disease.

METHODS

Fifteen subjects (13.4 ± 2.3 y/o) and 56 age-matched (13.5 ± 2.34 y/o) healthy controls were scanned using a GE Discovery MR750 3T MRI and 32-channel RF-coil. Axial FSPGR-3D images were used to prescribe rs-BOLD (TE/TR = 35/2000ms), acquired over 6 minutes. Motion correction was performed and anatomical and functional images were aligned and spatially warped to the N27 standard atlas. Fractal analysis, performed on grey matter, was done by estimating the Hurst exponent using de-trended fluctuation analysis and signal summation conversion methods.

RESULTS AND CONCLUSIONS

Voxel-wise fractal dimension (FD) was calculated for every subject in the control group to generate mean and standard deviation maps for regional Z-score analysis. Voxel-wise validation of FD normality across controls was confirmed, and non-Gaussian voxels (3.05% over the brain) were eliminated from subsequent analysis. For each mTBI patient, regions where Z-score values were at least 2 standard deviations away from the mean (i.e. where |Z| > 2.0) were identified. In individual patients the frequently affected regions were amygdala (p = 0.02), vermis(p = 0.03), caudate head (p = 0.04), hippocampus(p = 0.03), and hypothalamus(p = 0.04), all previously reported as dysfunctional after mTBI, but based on group analysis. It is well known that the brain is best modeled as a complex system. Therefore a measure of complexity using rs-BOLD signal FD could provide an additional method to grade and monitor mTBI. Furthermore, this approach can be personalized thus providing unique patient specific assessment.

Resting-state theta band connectivity and graph analysis in generalized social anxiety disorder

BACKGROUND

Functional magnetic resonance imaging (fMRI) resting-state studies show generalized social anxiety disorder (gSAD) is associated with disturbances in networks involved in emotion regulation, emotion processing, and perceptual functions, suggesting a network framework is integral to elucidating the pathophysiology of gSAD. However, fMRI does not measure the fast dynamic interconnections of functional networks. Therefore, we examined whole-brain functional connectomics with electroencephalogram (EEG) during resting-state.

METHODS

Resting-state EEG data was recorded for 32 patients with gSAD and 32 demographically-matched healthy controls (HC). Sensor-level connectivity analysis was applied on EEG data by using Weighted Phase Lag Index (WPLI) and graph analysis based on WPLI was used to determine clustering coefficient and characteristic path length to estimate local integration and global segregation of networks.

RESULTS

WPLI results showed increased oscillatory midline coherence in the theta frequency band indicating higher connectivity in the gSAD relative to HC group during rest. Additionally, WPLI values positively correlated with state anxiety levels within the gSAD group but not the HC group. Our graph theory based connectomics analysis demonstrated increased clustering coefficient and decreased characteristic path length in theta-based whole brain functional organization in subjects with gSAD compared to HC.

CONCLUSIONS

Theta-dependent interconnectivity was associated with state anxiety in gSAD and an increase in information processing efficiency in gSAD (compared to controls). Results may represent enhanced baseline self-focused attention, which is consistent with cognitive models of gSAD and fMRI studies implicating emotion dysregulation and disturbances in task negative networks (e.g., default mode network) in gSAD.

Associations between autistic traits and fractional anisotropy values in white matter tracts in a nonclinical sample of young adults

Whereas a number of studies have examined relationships among brain activity, social cognitive skills, and autistic traits, fewer studies have evaluated whether structural connections among brain regions relate to these traits and skills. Uncinate fasciculus (UF) and inferior longitudinal fasciculus (ILF) are white matter tracts that may underpin the behavioral expression of these skills because they connect regions within or provide sensory information to brain areas implicated in social cognition, and structural differences in these tracts have been associated with autistic traits. We examined relationships among self-reported autistic traits, mentalizing, and water diffusivity in UF and ILF in a nonclinical sample of 24 young adults (mean age = 21.92 years, SD = 4.72 years; 15 women). We measured autistic traits using the Autism-Spectrum Quotient, and we measured mentalizing using the Dynamic Interactive Shapes Clips task. We used Tract-Based Spatial Statistics and randomize to examine relationships among fractional anisotropy (FA) values in bilateral ILF and UF, age, cognitive abilities, autistic traits, and mentalizing. Autistic traits were positively related to FA values in left ILF. No other relationships between FA values and other variables were significant. Results suggest that left ILF may be involved in the expression of autistic traits in individuals without clinical diagnoses.