Endogenous multifractal brain dynamics are modulated by age, cholinergic blockade and cognitive performance

Pallidal multifractal complexity is a new potential physiomarker of dystonia

OBJECTIVE

Low-frequency 4-12 Hz pallidal oscillations are being considered as potential physiomarkers for dystonia. We suggest investigating the multifractal properties of pallidal activity as an additional marker.

METHODS

We employed local field potentials (LFP) recordings from 23 patients with dystonia who were undergoing deep brain stimulation (DBS) surgery to explore the connection between disease severity and the multifractal characteristics of pallidal activity. Furthermore, we performed an analysis of LFP recordings from four patients, following the externalization of DBS lead electrodes, to investigate the impact of DBS and neck muscle vibration on multifractal parameters.

RESULTS

Greater dystonia severity exhibited a correlation with a narrower multifractal spectrum width but higher multifractal spectral asymmetry. Both GPi DBS and muscle vibration in dystonia patients expanded the multifractal spectrum width while restoring multifractal spectral symmetry. Notably, the threshold peak intensities for an increase in multifractal spectrum width substantially overlapped with the optimal volume of tissue activated. A broader multifractal spectrum during DBS corresponded to more favorable clinical outcomes.

CONCLUSIONS

Multifractal properties of pallidal neuronal activity serve as indicators of neural dysfunction in dystonia.

SIGNIFICANCE

These findings suggest the potential of utilizing multifractal characteristics as predictive factors for the DBS outcome in dystonia.

Postural control in gymnasts: anisotropic fractal scaling reveals proprioceptive reintegration in vestibular perturbation

Dexterous postural control subtly complements movement variability with sensory correlations at many scales. The expressive poise of gymnasts exemplifies this lyrical punctuation of release with constraint, from coarse grain to fine scales. Dexterous postural control upon a 2D support surface might collapse the variation of center of pressure (CoP) to a relatively 1D orientation-a direction often oriented towards the focal point of a visual task. Sensory corrections in dexterous postural control might manifest in temporal correlations, specifically as fractional Brownian motions whose differences are more and less correlated with fractional Gaussian noises (fGns) with progressively larger and smaller Hurst exponent H. Traditional empirical work examines this arrangement of lower-dimensional compression of CoP along two orthogonal axes, anteroposterior (AP) and mediolateral (ML). Eyes-open and face-forward orientations cultivate greater variability along AP than ML axes, and the orthogonal distribution of spatial variability has so far gone hand in hand with an orthogonal distribution of H, for example, larger in AP and lower in ML. However, perturbing the orientation of task focus might destabilize the postural synergy away from its 1D distribution and homogenize the temporal correlations across the 2D support surface, resulting in narrower angles between the directions of the largest and smallest H. We used oriented fractal scaling component analysis (OFSCA) to investigate whether sensory corrections in postural control might thus become suborthogonal. OFSCA models raw 2D CoP trajectory by decomposing it in all directions along the 2D support surface and fits the directions with the largest and smallest H. We studied a sample of gymnasts in eyes-open and face-forward quiet posture, and results from OFSCA confirm that such posture exhibits the classic orthogonal distribution of temporal correlations. Head-turning resulted in a simultaneous decrease in this angle Δθ, which promptly reversed once gymnasts reoriented their heads forward. However, when vision was absent, there was only a discernible negative trend in Δθ, indicating a shift in the angle's direction but not a statistically significant one. Thus, the narrowing of Δθ may signify an adaptive strategy in postural control. The swift recovery of Δθ upon returning to a forward-facing posture suggests that the temporary reduction is specific to head-turning and does not impose a lasting burden on postural control. Turning the head reduced the angle between these two orientations, facilitating the release of postural degrees of freedom towards a more uniform spread of the CoP across both dimensions of the support surface. The innovative aspect of this work is that it shows how fractality might serve as a control parameter of adaptive mechanisms of dexterous postural control.

Multifractal spectral features enhance classification of anomalous diffusion

Anomalous diffusion processes, characterized by their nonstandard scaling of the mean-squared displacement, pose a unique challenge in classification and characterization. In a previous study [Mangalam et al., Phys. Rev. Res. 5, 023144 (2023)2643-156410.1103/PhysRevResearch.5.023144], we established a comprehensive framework for understanding anomalous diffusion using multifractal formalism. The present study delves into the potential of multifractal spectral features for effectively distinguishing anomalous diffusion trajectories from five widely used models: fractional Brownian motion, scaled Brownian motion, continuous-time random walk, annealed transient time motion, and Lévy walk. We generate extensive datasets comprising 10^{6} trajectories from these five anomalous diffusion models and extract multiple multifractal spectra from each trajectory to accomplish this. Our investigation entails a thorough analysis of neural network performance, encompassing features derived from varying numbers of spectra. We also explore the integration of multifractal spectra into traditional feature datasets, enabling us to assess their impact comprehensively. To ensure a statistically meaningful comparison, we categorize features into concept groups and train neural networks using features from each designated group. Notably, several feature groups demonstrate similar levels of accuracy, with the highest performance observed in groups utilizing moving-window characteristics and p varation features. Multifractal spectral features, particularly those derived from three spectra involving different timescales and cutoffs, closely follow, highlighting their robust discriminatory potential. Remarkably, a neural network exclusively trained on features from a single multifractal spectrum exhibits commendable performance, surpassing other feature groups. In summary, our findings underscore the diverse and potent efficacy of multifractal spectral features in enhancing the predictive capacity of machine learning to classify anomalous diffusion processes.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

A Self-Similarity Logic May Shape the Organization of the Nervous System

From the morphological point of view, the nervous system exhibits a fractal, self-similar geometry at various levels of observations, from single cells up to cell networks. From the functional point of view, it is characterized by a hierarchical organization in which self-similar structures (networks) of different miniaturizations are nested within each other. In particular, neuronal networks, interconnected to form neuronal systems, are formed by neurons, which operate thanks to their molecular networks, mainly having proteins as components that via protein-protein interactions can be assembled in multimeric complexes working as micro-devices. On this basis, the term "self-similarity logic" was introduced to describe a nested organization where, at the various levels, almost the same rules (logic) to perform operations are used. Self-similarity and self-similarity logic both appear to be intimately linked to the biophysical evidence for the nervous system being a pattern-forming system that can flexibly switch from one coherent state to another. Thus, they can represent the key concepts to describe its complexity and its concerted, holistic behavior.

Morphology and Fractal-Based Classifications of Neurons and Microglia in Two and Three Dimensions

Microglia and neurons live physically intertwined, intimately related structurally and functionally in a dynamic relationship in which microglia change continuously over a much shorter timescale than do neurons. Although microglia may unwind and depart from the neurons they attend under certain circumstances, in general, together both contribute to the fractal topology of the brain that defines its computational capabilities. Both neuronal and microglial morphologies are well-described using fractal analysis complementary to more traditional measures. For neurons, the fractal dimension has proved valuable for classifying dendritic branching and other neuronal features relevant to pathology and development. For microglia, fractal geometry has substantially contributed to classifying functional categories, where, in general, the more pathological the biological status, the lower the fractal dimension for individual cells, with some exceptions, including hyper-ramification. This chapter provides a review of the intimate relationships between neurons and microglia, by introducing 2D and 3D fractal analysis methodology and its applications in neuron-microglia function in health and disease.

Fractals in Neuroanatomy and Basic Neurosciences: An Overview

The introduction of fractal geometry to the neurosciences has been a major paradigm shift over the last decades as it has helped overcome approximations and limitations that occur when Euclidean and reductionist approaches are used to analyze neurons or the entire brain. Fractal geometry allows for quantitative analysis and description of the geometric complexity of the brain, from its single units to the neuronal networks.As illustrated in the second section of this book, fractal analysis provides a quantitative tool for the study of the morphology of brain cells (i.e., neurons and microglia) and its components (e.g., dendritic trees, synapses), as well as the brain structure itself (cortex, functional modules, neuronal networks). The self-similar logic which generates and shapes the different hierarchical systems of the brain and even some structures related to its "container," that is, the cranial sutures on the skull, is widely discussed in the following chapters, with a link between the applications of fractal analysis to the neuroanatomy and basic neurosciences to the clinical applications discussed in the third section.

Improvements in task performance after practice are associated with scale-free dynamics of brain activity

Although practicing a task generally benefits later performance on that same task, there are individual differences in practice effects. One avenue to model such differences comes from research showing that brain networks extract functional advantages from operating in the vicinity of criticality, a state in which brain network activity is more scale-free. We hypothesized that higher scale-free signal from fMRI data, measured with the Hurst exponent (H), indicates closer proximity to critical states. We tested whether individuals with higher H during repeated task performance would show greater practice effects. In Study 1, participants performed a dual-n-back task (DNB) twice during MRI (n = 56). In Study 2, we used two runs of n-back task (NBK) data from the Human Connectome Project sample (n = 599). In Study 3, participants performed a word completion task (CAST) across six runs (n = 44). In all three studies, multivariate analysis was used to test whether higher H was related to greater practice-related performance improvement. Supporting our hypothesis, we found patterns of higher H that reliably correlated with greater performance improvement across participants in all three studies. However, the predictive brain regions were distinct, suggesting that the specific spatial H↑ patterns are not task-general.

A pattern of cognitive resource disruptions in childhood psychopathology

The Hurst exponent (H) isolated in fractal analyses of neuroimaging time series is implicated broadly in cognition. Within this literature, H is associated with multiple mental disorders, suggesting that H is transdimensionally associated with psychopathology. Here, we unify these results and demonstrate a pattern of decreased H with increased general psychopathology and attention-deficit/hyperactivity factor scores during a working memory task in 1,839 children. This pattern predicts current and future cognitive performance in children and some psychopathology in 703 adults. This pattern also defines psychological and functional axes associating psychopathology with an imbalance in resource allocation between fronto-parietal and sensorimotor regions, driven by reduced resource allocation to fronto-parietal regions. This suggests the hypothesis that impaired working memory function in psychopathology follows from a reduced cognitive resource pool and a reduction in resources allocated to the task at hand.

The time scales of irreversibility in spontaneous brain activity are altered in obsessive compulsive disorder

We study how obsessive-compulsive disorder (OCD) affects the complexity and time-reversal symmetry-breaking (irreversibility) of the brain resting-state activity as measured by magnetoencephalography (MEG). Comparing MEG recordings from OCD patients and age/sex matched control subjects, we find that irreversibility is more concentrated at faster time scales and more uniformly distributed across different channels of the same hemisphere in OCD patients than in control subjects. Furthermore, the interhemispheric asymmetry between homologous areas of OCD patients and controls is also markedly different. Some of these differences were reduced by 1-year of Kundalini Yoga meditation treatment. Taken together, these results suggest that OCD alters the dynamic attractor of the brain's resting state and hint at a possible novel neurophysiological characterization of this psychiatric disorder and how this therapy can possibly modulate brain function.

Brain entropy, fractal dimensions and predictability: A review of complexity measures for EEG in healthy and neuropsychiatric populations

There has been an increasing trend towards the use of complexity analysis in quantifying neural activity measured by electroencephalography (EEG) signals. On top of revealing complex neuronal processes of the brain that may not be possible with linear approaches, EEG complexity measures have also demonstrated their potential as biomarkers of psychopathology such as depression and schizophrenia. Unfortunately, the opacity of algorithms and descriptions originating from mathematical concepts have made it difficult to understand what complexity is and how to draw consistent conclusions when applied within psychology and neuropsychiatry research. In this review, we provide an overview and entry-level explanation of existing EEG complexity measures, which can be broadly categorized as measures of predictability and regularity. We then synthesize complexity findings across different areas of psychological science, namely, in consciousness research, mood and anxiety disorders, schizophrenia, neurodevelopmental and neurodegenerative disorders, as well as changes across the lifespan, while addressing some theoretical and methodological issues underlying the discrepancies in the data. Finally, we present important considerations when choosing and interpreting these metrics.

Turing's cascade instability supports the coordination of the mind, brain, and behavior

Turing inspired a computer metaphor of the mind and brain that has been handy and has spawned decades of empirical investigation, but he did much more and offered behavioral and cognitive sciences another metaphor-that of the cascade. The time has come to confront Turing's cascading instability, which suggests a geometrical framework driven by power laws and can be studied using multifractal formalism and multiscale probability density function analysis. Here, we review a rapidly growing body of scientific investigations revealing signatures of cascade instability and their consequences for a perceiving, acting, and thinking organism. We review work related to executive functioning (planning to act), postural control (bodily poise for turning plans into action), and effortful perception (action to gather information in a single modality and action to blend multimodal information). We also review findings on neuronal avalanches in the brain, specifically about neural participation in body-wide cascades. Turing's cascade instability blends the mind, brain, and behavior across space and time scales and provides an alternative to the dominant computer metaphor.

Fractal Structure of Brain Electrical Activity of Patients With Mental Disorders

This work was aimed at a comparative analysis of the degree of multifractality of electroencephalographic time series obtained from a group of healthy subjects and from patients with mental disorders. We analyzed long-term records of patients with paranoid schizophrenia and patients with depression. To evaluate the properties of multifractal scaling of various electroencephalographic time series, the method of maximum modulus of the wavelet transform and multifractal analysis of fluctuations without a trend were used. The stability of the width and position of the singularity spectrum for each of the test groups was revealed, and a relationship was established between the correlation and anticorrelation dynamics of successive values of the electroencephalographic time series and the type of mental disorders. It was shown that the main differences between the multifractal properties of brain activity in normal and pathological conditions lie in the different width of the multifractality spectrum and its location associated with the correlated or anticorrelated dynamics of the values of successive time series. It was found that the schizophrenia group is characterized by a greater degree of multifractality compared to the depression group. Thus, the degree of multifractality can be included in a set of tests for differential diagnosis and research of mental disorders.

Brain response to color stimuli: an EEG study with nonlinear approach

Color perception is a major guiding factor in the evolutionary process of human civilization, but most of the neurological background of the same are yet unknown. This work attempts to address this area with an EEG based neuro-cognitive study on response of brain to different color stimuli. With respect to a Grey baseline seven colors of the VIBGYOR were shown to 16 participants with normal color vision and corresponding EEG signals from different lobes (Frontal, Occipital & Parietal) were recorded. In an attempt to quantify the brain response while watching these colors, the corresponding EEG signals were analysed using two of the latest state of the art non-linear techniques (MFDFA and MFDXA) of dealing complex time series. MFDFA revealed that for all the participants the spectral width, and hence the complexity of the EEG signals, reaches a maximum while viewing color Blue, followed by colors Red and Green in all the brain lobes. MFDXA, on the other hand, suggests a lower degree of inter and intra lobe correlation while watching the VIBGYOR colors compared to baseline Grey, hinting towards a post processing of visual information. We hope that along with the novelty of methodologies, the unique outcomes of this study may leave a long term impact in the domain of color perception research.

Empirical Evidence for Extended Cognitive Systems

We present an empirically supported theoretical and methodological framework for quantifying the system‐level properties of person‐plus‐tool interactions in order to answer the question: “Are person‐plus‐tool‐systems extended cognitive systems?” Nineteen participants provided perceptual judgments regarding their ability to pass through apertures of various widths while using visual information, blindfolded wielding a rod, or blindfolded wielding an Enactive Torch—a vibrotactile sensory‐substitution device for detecting distance. Monofractal, multifractal, and recurrence quantification analyses were conducted to assess features of person‐plus‐tool movement dynamics. Trials where people utilized the rod or Enactive Torch demonstrated stable “self‐similarity,” or indices of healthy and adaptive single systems, regardless of aperture width, trial order, features of the participants’ judgments, and participant characteristics. Enactive Torch trials exhibited a somewhat greater range of dynamic fluctuations than the rod trials, as well as less movement recurrence, suggesting that the Enactive Torch allowed for more exploratory movements. Findings provide support for the notion that person‐plus‐tool systems can be classified as extended cognitive systems and a framework for quantifying system‐level properties of these systems. Implications concerning future research on extended cognition are discussed.

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.

Why Brain Criticality Is Clinically Relevant: A Scoping Review

The past 25 years have seen a strong increase in the number of publications related to criticality in different areas of neuroscience. The potential of criticality to explain various brain properties, including optimal information processing, has made it an increasingly exciting area of investigation for neuroscientists. Recent reviews on this topic, sometimes termed brain criticality, make brief mention of clinical applications of these findings to several neurological disorders such as epilepsy, neurodegenerative disease, and neonatal hypoxia. Other clinicallyrelevant domains - including anesthesia, sleep medicine, developmental-behavioral pediatrics, and psychiatry - are seldom discussed in review papers of brain criticality. Thorough assessments of these application areas and their relevance for clinicians have also yet to be published. In this scoping review, studies of brain criticality involving human data of all ages are evaluated for their current and future clinical relevance. To make the results of these studies understandable to a more clinical audience, a review of the key concepts behind criticality (e.g., phase transitions, long-range temporal correlation, self-organized criticality, power laws, branching processes) precedes the discussion of human clinical studies. Open questions and forthcoming areas of investigation are also considered.

Scale-free amplitude modulation of low-frequency fluctuations in episodic migraine

Arrhythmic fluctuations in neural activity occur at many levels of the nervous system. Such activity does not have a characteristic temporal periodicity but can exhibit statistical similarities, most commonly power-law scaling behavior, which is indicative of scale-free dynamics. The recurrence of scaling laws across many different systems and its manifestation in behavior has prompted a search for unifying principles in human brain function. With this in mind, a focused search for abnormities in scale-free dynamics is of considerable clinical relevance to migraine and other clinical pain disorders. Here, we examined the scale-free properties of the resting-state functional magnetic resonance imaging (fMRI) signal in the broadband frequency range known to be related to spontaneous neural activity (0.01-0.1 Hz). In a large cohort of episodic migraine patients (N = 40), we observed that the strength of long-range temporal correlations in the fMRI signal (captured by the scaling exponent α) was significantly higher in the sensorimotor network compared with healthy controls. Increases in the scaling exponent were positively correlated with fMRI signal variance and negatively correlated with the patient's self-reported headache intensity. These changes in the fMRI signal suggest that the temporal structure of amplitude fluctuations carries valuable information about the dynamic state of the underlying neuronal networks and ensuing sensory impairments in migraine. The demonstrated scaling laws pose a novel quantitative approach for examining clinically relevant interindividual variability in migraine and other pain disorders.

Global Effects of Focal Brain Tumors on Functional Complexity and Network Robustness: A Prospective Cohort Study

BACKGROUND

Neurosurgical management of brain tumors has entered a paradigm of supramarginal resections that demands thorough understanding of peritumoral functional effects. Historically, the effects of tumors have been believed to be local, and long-range effects have not been considered.

OBJECTIVE

To test the hypothesis that tumors affect the brain globally, producing long-range gradients in cortical function.

METHODS

Resting-state functional magnetic resonance imaging (fMRI) data were acquired from 11 participants with glioblastoma and split into discovery and validation datasets in a single-center prospective cohort study. Fractal complexity was computed with a wavelet-based estimator of the Hurst exponent. Distance-related effects of the tumors were tested with a tumor mask-dilation technique and parcellation of the underlying Hurst maps.

RESULTS

Fractal complexity demonstrates a penumbra of suppression in the peritumoral region. At a global level, as distance from the tumor increases, this initial suppression is balanced by a subsequent overactivity before finally normalizing. These effects were best fit by a quadratic model and were consistent across different network construction pipelines. The Hurst exponent was correlated with graph theory measures of centrality including network robustness, but graph theory measures did not demonstrate distance-dependent effects.

CONCLUSION

This work provides evidence supporting the theory that focal brain tumors produce long-range gradients in function. Consequently, the effects of focal lesions need to be interpreted in terms of the global changes on functional complexity and network architecture rather than purely in terms of functional localization. Determining whether peritumoral changes represent potential plasticity may facilitate extended resection of tumors without functional cost.

Scale-free functional brain dynamics during recovery from sport-related concussion

Studies using blood‐oxygenation‐level‐dependent functional magnetic resonance imaging (BOLD fMRI) have characterized how the resting brain is affected by concussion. The literature to date, however, has largely focused on measuring changes in the spatial organization of functional brain networks. In the present study, changes in the temporal dynamics of BOLD signals are examined throughout concussion recovery using scaling (or fractal) analysis. Imaging data were collected for 228 university‐level athletes, 61 with concussion and 167 athletic controls. Concussed athletes were scanned at the acute phase of injury (1–7 days postinjury), the subacute phase (8–14 days postinjury), medical clearance to return to sport (RTS), 1 month post‐RTS and 1 year post‐RTS. The wavelet leader multifractal approach was used to assess scaling (c₁) and multifractal (c₂) behavior. Significant longitudinal changes were identified for c₁, which was lowest at acute injury, became significantly elevated at RTS, and returned near control levels by 1 year post‐RTS. No longitudinal changes were identified for c₂. Secondary analyses showed that clinical measures of acute symptom severity and time to RTP were related to longitudinal changes in c₁. Athletes with both higher symptoms and prolonged recovery had elevated c₁ values at RTS, while athletes with higher symptoms but rapid recovery had reduced c₁ at acute injury. This study provides the first evidence for long‐term recovery of BOLD scale‐free brain dynamics after a concussion.

Connections, Tracts, Fractals, and the Rest: A Working Guide to Network and Connectivity Studies in Neurosurgery

Brain mapping and connectomics can probe networks that span the entire brain, producing a diverse range of outputs for probing specific clinically relevant questions. The potential for understanding the effect of focal lesions on brain function, cognition, and plasticity abounds, any one of which would likely yield more effective and safer neurosurgical strategies. However, the possibilities of advanced magnetic resonance imaging and connectomics have been somewhat underused in neurosurgery, arising from actual or perceived difficulties in either application or analysis. The present review builds on previous work describing the theoretical attractions of connectomics to deliberate on the practical details of performing high-quality connectomics studies in neurosurgery. First, the data and methods involved in deriving connectomics models will be considered, specifically for the purpose of determining the nature of inferences that can be made subsequently. Next, a selection of key analysis methods will be explored using practical examples that illustrate their effective implementation and the insights that can be gleaned. The principles of study design will be introduced, including analysis tips and methods for making efficient use of available resources. Finally, a review of the best research practices for neuroimaging studies will be discussed, including principles of open access data sharing, study preregistration, and methods for improving replicability. Ultimately, we hope readers will be better placed to appraise the current connectomics studies in neurosurgery and empowered to develop their own high-quality studies, both of which are key steps in realizing the true potential of connectomics and advanced neuroimaging analyses in general.

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.

1/f electrophysiological spectra in resting and drug-induced states can be explained by the dynamics of multiple oscillatory relaxation processes

Neurophysiological recordings are dominated by arhythmical activity whose spectra can be characterised by power-law functions, and on this basis are often referred to as reflecting scale-free brain dynamics (1/f^β). Relatively little is known regarding the neural generators and temporal dynamics of this arhythmical behaviour compared to rhythmical behaviour. Here we used Irregularly Resampled AutoSpectral Analysis (IRASA) to quantify β, in both the high (5-100 Hz, β_hf) and low frequency bands (0.1-2.5 Hz, β_lf) in MEG/EEG/ECoG recordings and to separate arhythmical from rhythmical modes of activity, such as, alpha rhythms. In MEG/EEG/ECoG data, we demonstrate that oscillatory alpha power dynamically correlates over time with β_hf and similarly, participants with higher rhythmical alpha power have higher β_hf. In a series of pharmaco-MEG investigations using the GABA reuptake inhibitor tiagabine, the glutamatergic AMPA receptor antagonist perampanel, the NMDA receptor antagonist ketamine and the mixed partial serotonergic agonist LSD, a variety of effects on both β_hf and β_lf were observed. Additionally, strong modulations of β_hf were seen in monkey ECoG data during general anaesthesia using propofol and ketamine. We develop and test a unifying model which can explain, the 1/f nature of electrophysiological spectra, their dynamic interaction with oscillatory rhythms as well as the sensitivity of 1/f activity to drug interventions by considering electrophysiological spectra as being generated by a collection of stochastically perturbed damped oscillators having a distribution of relaxation rates.

The informational entropy endowed in cortical oscillations

A two-dimensional shadow may encompass more information than its corresponding three-dimensional object. Indeed, if we rotate the object, we achieve a pool of observed shadows from different angulations, gradients, shapes and variable length contours that make it possible for us to increase our available information. Starting from this simple observation, we show how informational entropies might turn out to be useful in the evaluation of scale-free dynamics in the brain. Indeed, brain activity exhibits a scale-free distribution that leads to the variations in the power law exponent typical of different functional neurophysiological states. Here we show that modifications in scaling slope are associated with variations in Rényi entropy, a generalization of Shannon informational entropy. From a three-dimensional object's perspective, by changing its orientation (standing for the cortical scale-free exponent), we detect different two-dimensional shadows from different perception angles (standing for Rényi entropy in different brain areas). We show how, starting from known values of Rényi entropy (easily detectable in brain fMRIs or EEG traces), it is feasible to calculate the scaling slope in a given moment and in a given brain area. Because changes in scale-free cortical dynamics modify brain activity, this issue points towards novel approaches to mind reading and description of the forces required for transcranial stimulation.

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.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.

The suppression of scale-free fMRI brain dynamics across three different sources of effort: aging, task novelty and task difficulty

There is growing evidence that fluctuations in brain activity may exhibit scale-free (“fractal”) dynamics. Scale-free signals follow a spectral-power curve of the form P(f ) ∝ f^−β, where spectral power decreases in a power-law fashion with increasing frequency. In this study, we demonstrated that fractal scaling of BOLD fMRI signal is consistently suppressed for different sources of cognitive effort. Decreases in the Hurst exponent (H), which quantifies scale-free signal, was related to three different sources of cognitive effort/task engagement: 1) task difficulty, 2) task novelty, and 3) aging effects. These results were consistently observed across multiple datasets and task paradigms. We also demonstrated that estimates of H are robust across a range of time-window sizes. H was also compared to alternative metrics of BOLD variability (SD_BOLD) and global connectivity (Gconn), with effort-related decreases in H producing similar decreases in SD_BOLD and Gconn. These results indicate a potential global brain phenomenon that unites research from different fields and indicates that fractal scaling may be a highly sensitive metric for indexing cognitive effort/task engagement.

Multi-Dimensional Dynamics of Human Electromagnetic Brain Activity

Magnetoencephalography (MEG) and electroencephalography (EEG) are invaluable neuroscientific tools for unveiling human neural dynamics in three dimensions (space, time, and frequency), which are associated with a wide variety of perceptions, cognition, and actions. MEG/EEG also provides different categories of neuronal indices including activity magnitude, connectivity, and network properties along the three dimensions. In the last 20 years, interest has increased in inter-regional connectivity and complex network properties assessed by various sophisticated scientific analyses. We herein review the definition, computation, short history, and pros and cons of connectivity and complex network (graph-theory) analyses applied to MEG/EEG signals. We briefly describe recent developments in source reconstruction algorithms essential for source-space connectivity and network analyses. Furthermore, we discuss a relatively novel approach used in MEG/EEG studies to examine the complex dynamics represented by human brain activity. The correct and effective use of these neuronal metrics provides a new insight into the multi-dimensional dynamics of the neural representations of various functions in the complex human brain.

Distinguishing cognitive state with multifractal complexity of hippocampal interspike interval sequences

Fractality, represented as self-similar repeating patterns, is ubiquitous in nature and the brain. Dynamic patterns of hippocampal spike trains are known to exhibit multifractal properties during working memory processing; however, it is unclear whether the multifractal properties inherent to hippocampal spike trains reflect active cognitive processing. To examine this possibility, hippocampal neuronal ensembles were recorded from rats before, during and after a spatial working memory task following administration of tetrahydrocannabinol (THC), a memory-impairing component of cannabis. Multifractal detrended fluctuation analysis was performed on hippocampal interspike interval sequences to determine characteristics of monofractal long-range temporal correlations (LRTCs), quantified by the Hurst exponent, and the degree/magnitude of multifractal complexity, quantified by the width of the singularity spectrum. Our results demonstrate that multifractal firing patterns of hippocampal spike trains are a marker of functional memory processing, as they are more complex during the working memory task and significantly reduced following administration of memory impairing THC doses. Conversely, LRTCs are largest during resting state recordings, therefore reflecting different information compared to multifractality. In order to deepen conceptual understanding of multifractal complexity and LRTCs, these measures were compared to classical methods using hippocampal frequency content and firing variability measures. These results showed that LRTCs, multifractality, and theta rhythm represent independent processes, while delta rhythm correlated with multifractality. Taken together, these results provide a novel perspective on memory function by demonstrating that the multifractal nature of spike trains reflects hippocampal microcircuit activity that can be used to detect and quantify cognitive, physiological, and pathological states.

Exploring multifractal-based features for mild Alzheimer's disease classification

PURPOSE

Multifractal applications to resting state functional MRI (rs-fMRI) time series for diagnosing Alzheimer's disease (AD) are still limited. We aim to address two issues: (I) if and what multifractal features are sufficiently discriminative to detect AD from the healthy; (II) if AD classification could be further improved by combining multifractal features with traditional features in this field.

METHODS

Rs-fMRI data of 25 AD patients and 38 normal controls were analyzed. A set of multifractal features were systematically investigated. Traditional features in monofractal, linear, and network-based categories were also extracted for comparison and combination. Both support vector machines and multiple kernel learning (MKL) were used to perform classification with individual and combined features.

RESULTS

We identified a multifractal feature, Δf, which has the strongest discriminative power among all the features in our study. Moreover, we found that the classification accuracy could be significantly improved from 69% (by Δf only) to up to 76%, when nonsparse MKL is used to combine Δf with the monofractal feature, Hurst. Finally, we showed that incorporating other multifractal features, α(0), Δα and Pc, could also improve traditional-feature-based AD classification.

CONCLUSION

Our work demonstrated the potential usefulness of multifractal analysis for AD research, especially when combining with the traditional rs-fMRI features. It contributes to distinguishing AD from NC subjects. Magn Reson Med 76:259-269, 2016. © 2015 Wiley Periodicals, Inc.

How can we study reasoning in the brain?

The brain did not develop a dedicated device for reasoning. This fact bears dramatic consequences. While for perceptuo-motor functions neural activity is shaped by the input's statistical properties, and processing is carried out at high speed in hardwired spatially segregated modules, in reasoning, neural activity is driven by internal dynamics and processing times, stages, and functional brain geometry are largely unconstrained a priori. Here, it is shown that the complex properties of spontaneous activity, which can be ignored in a short-lived event-related world, become prominent at the long time scales of certain forms of reasoning. It is argued that the neural correlates of reasoning should in fact be defined in terms of non-trivial generic properties of spontaneous brain activity, and that this implies resorting to concepts, analytical tools, and ways of designing experiments that are as yet non-standard in cognitive neuroscience. The implications in terms of models of brain activity, shape of the neural correlates, methods of data analysis, observability of the phenomenon, and experimental designs are discussed.

Nonlinear analysis of saccade speed fluctuations during combined action and perception tasks

BACKGROUND

Saccades are rapid eye movements used to gather information about a scene which requires both action and perception. These are usually studied separately, so that how perception influences action is not well understood. In a dual task, where the subject looks at a target and reports a decision, subtle changes in the saccades might be caused by action-perception interactions. Studying saccades might provide insight into how brain pathways for action and for perception interact.

NEW METHOD

We applied two complementary methods, multifractal detrended fluctuation analysis and Lempel-Ziv complexity index to eye peak speed recorded in two experiments, a pure action task and a combined action-perception task.

RESULTS

Multifractality strength is significantly different in the two experiments, showing smaller values for dual decision task saccades compared to simple-task saccades. The normalized Lempel-Ziv complexity index behaves similarly i.e. is significantly smaller in the decision saccade task than in the simple task.

COMPARISON WITH EXISTING METHODS

Compared to the usual statistical and linear approaches, these analyses emphasize the character of the dynamics involved in the fluctuations and offer a sensitive tool for quantitative evaluation of the multifractal features and of the complexity measure in the saccades peak speeds when different brain circuits are involved.

CONCLUSION

Our results prove that the peak speed fluctuations have multifractal characteristics with lower magnitude for the multifractality strength and for the complexity index when two neural pathways are simultaneously activated, demonstrating the nonlinear interaction in the brain pathways for action and perception.

Interplay between functional connectivity and scale-free dynamics in intrinsic fMRI networks

Studies employing functional connectivity-type analyses have established that spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals are organized within large-scale brain networks. Meanwhile, fMRI signals have been shown to exhibit 1/f-type power spectra - a hallmark of scale-free dynamics. We studied the interplay between functional connectivity and scale-free dynamics in fMRI signals, utilizing the fractal connectivity framework - a multivariate extension of the univariate fractional Gaussian noise model, which relies on a wavelet formulation for robust parameter estimation. We applied this framework to fMRI data acquired from healthy young adults at rest and while performing a visual detection task. First, we found that scale-invariance existed beyond univariate dynamics, being present also in bivariate cross-temporal dynamics. Second, we observed that frequencies within the scale-free range do not contribute evenly to inter-regional connectivity, with a systematically stronger contribution of the lowest frequencies, both at rest and during task. Third, in addition to a decrease of the Hurst exponent and inter-regional correlations, task performance modified cross-temporal dynamics, inducing a larger contribution of the highest frequencies within the scale-free range to global correlation. Lastly, we found that across individuals, a weaker task modulation of the frequency contribution to inter-regional connectivity was associated with better task performance manifesting as shorter and less variable reaction times. These findings bring together two related fields that have hitherto been studied separately - resting-state networks and scale-free dynamics, and show that scale-free dynamics of human brain activity manifest in cross-regional interactions as well.

Functional significance of complex fluctuations in brain activity: from resting state to cognitive neuroscience

Behavioral studies have shown that human cognition is characterized by properties such as temporal scale invariance, heavy-tailed non-Gaussian distributions, and long-range correlations at long time scales, suggesting models of how (non observable) components of cognition interact. On the other hand, results from functional neuroimaging studies show that complex scaling and intermittency may be generic spatio-temporal properties of the brain at rest. Somehow surprisingly, though, hardly ever have the neural correlates of cognition been studied at time scales comparable to those at which cognition shows scaling properties. Here, we analyze the meanings of scaling properties and the significance of their task-related modulations for cognitive neuroscience. It is proposed that cognitive processes can be framed in terms of complex generic properties of brain activity at rest and, ultimately, of functional equations, limiting distributions, symmetries, and possibly universality classes characterizing them.

Fractal analysis of spontaneous fluctuations of the BOLD signal in the human brain networks

PURPOSE

To investigate what extent brain regions are continuously interacting during resting-state, independent component analyses (ICA) was applied to analyze resting-state functional MRI (RS-fMRI) data. According to the analyzed results, it was surprisingly found that low frequency fluctuations (LFFs), which belong to the 1/f signal (a signal with power spectrum whose power spectral density is inversely proportional to the frequency), have been classified into groups using ICA; furthermore, the spatial distributions of these groups within the brain were found to resemble the spatial distributions of different networks, which manifests that the signal characteristics of RS LFFs are distinct across networks. In our work, we applied the 1/f model in the fractal analyses to further investigate this distinction.

MATERIALS AND METHODS

Twenty healthy participants got involved in this study. They were scanned to acquire the RS-fMRI data. The acquired data were first processed with ICA to obtain the networks of the resting brain. Afterward, the blood-oxygenation level dependent (BOLD) signals of these networks were processed with the fractal analyses for obtaining the fractal parameter α.

RESULTS

α was found to significantly vary across networks, which reveals that the fractal characteristic of LFFs differs across networks. According to prior literatures, this difference could be brought by the discrepancy of hemodynamic response amplitude (HRA) between networks. Hence, in our work, we also performed the computational simulation to discover the relationship between α and HRA. Based on the simulation results, HRA is highly linear-correlated with the fractal characteristics of LFFs which is revealed by α.

CONCLUSION

Our results support that the origin of RS-fMRI signals contains arterial fluctuations. Hence, in addition to the commonly used method such as synchrony analysis and power spectral analysis, another approach, the fractal analysis, is suggested for acquiring the information of hemodynamic responses by means of RS-fMRI data.

Optimizing complexity measures for FMRI data: algorithm, artifact, and sensitivity

INTRODUCTION

Complexity in the brain has been well-documented at both neuronal and hemodynamic scales, with increasing evidence supporting its use in sensitively differentiating between mental states and disorders. However, application of complexity measures to fMRI time-series, which are short, sparse, and have low signal/noise, requires careful modality-specific optimization.

METHODS

HERE WE USE BOTH SIMULATED AND REAL DATA TO ADDRESS TWO FUNDAMENTAL ISSUES: choice of algorithm and degree/type of signal processing. Methods were evaluated with regard to resilience to acquisition artifacts common to fMRI as well as detection sensitivity. Detection sensitivity was quantified in terms of grey-white matter contrast and overlap with activation. We additionally investigated the variation of complexity with activation and emotional content, optimal task length, and the degree to which results scaled with scanner using the same paradigm with two 3T magnets made by different manufacturers. Methods for evaluating complexity were: power spectrum, structure function, wavelet decomposition, second derivative, rescaled range, Higuchi's estimate of fractal dimension, aggregated variance, and detrended fluctuation analysis. To permit direct comparison across methods, all results were normalized to Hurst exponents.

RESULTS

Power-spectrum, Higuchi's fractal dimension, and generalized Hurst exponent based estimates were most successful by all criteria; the poorest-performing measures were wavelet, detrended fluctuation analysis, aggregated variance, and rescaled range.

CONCLUSIONS

Functional MRI data have artifacts that interact with complexity calculations in nontrivially distinct ways compared to other physiological data (such as EKG, EEG) for which these measures are typically used. Our results clearly demonstrate that decisions regarding choice of algorithm, signal processing, time-series length, and scanner have a significant impact on the reliability and sensitivity of complexity estimates.

Moment-to-moment brain signal variability: a next frontier in human brain mapping?

Neuroscientists have long observed that brain activity is naturally variable from moment-to-moment, but neuroimaging research has largely ignored the potential importance of this phenomenon. An emerging research focus on within-person brain signal variability is providing novel insights, and offering highly predictive, complementary, and even orthogonal views of brain function in relation to human lifespan development, cognitive performance, and various clinical conditions. As a result, brain signal variability is evolving as a bona fide signal of interest, and should no longer be dismissed as meaningless noise when mapping the human brain.

Neuronal long-range temporal correlations and avalanche dynamics are correlated with behavioral scaling laws

Scale-free fluctuations are ubiquitous in behavioral performance and neuronal activity. In time scales from seconds to hundreds of seconds, psychophysical dynamics and the amplitude fluctuations of neuronal oscillations are governed by power-law-form long-range temporal correlations (LRTCs). In millisecond time scales, neuronal activity comprises cascade-like neuronal avalanches that exhibit power-law size and lifetime distributions. However, it remains unknown whether these neuronal scaling laws are correlated with those characterizing behavioral performance or whether neuronal LRTCs and avalanches are related. Here, we show that the neuronal scaling laws are strongly correlated both with each other and with behavioral scaling laws. We used source reconstructed magneto- and electroencephalographic recordings to characterize the dynamics of ongoing cortical activity. We found robust power-law scaling in neuronal LRTCs and avalanches in resting-state data and during the performance of audiovisual threshold stimulus detection tasks. The LRTC scaling exponents of the behavioral performance fluctuations were correlated with those of concurrent neuronal avalanches and LRTCs in anatomically identified brain systems. The behavioral exponents also were correlated with neuronal scaling laws derived from a resting-state condition and with a similar anatomical topography. Finally, despite the difference in time scales, the scaling exponents of neuronal LRTCs and avalanches were strongly correlated during both rest and task performance. Thus, long and short time-scale neuronal dynamics are related and functionally significant at the behavioral level. These data suggest that the temporal structures of human cognitive fluctuations and behavioral variability stem from the scaling laws of individual and intrinsic brain dynamics.

Identifying multiplicative interactions between temporal scales of human movement variability

Conventional scaling analyses of human movement variability such as detrended fluctuation analyses assume that the movement variable can be decomposed into scale-dependent variation. However, these conventional scaling analyses are insensitive to multiplicative interactions within the movement variable. Multiplicative interactions refer to couplings between the scale-dependent variations across multiple scales that generate intermittent changes in the human movement variable. The mathematical concept for intermittent variability generated by multiplicative interactions is called multifractal variability. Multifractal variability is numerically defined by a spectrum of scaling exponents (i.e., a multifractal spectrum) that can be an important feature of coordinated movements and, consequently, relevant when aiming to identify movement disorders. In the current study, a new method is introduced based on detrended fluctuation analysis that can identify the multifractal spectrum from the temporal variation of local scaling exponents. The influence of multiplicative interactions on the local scaling exponents is tested by a Monte Carlo surrogate test. The methods are validated on multiplicative cascading processes with known multiplicative interactions. The application of the new methods is subsequently illustrated on an example of centre of pressure variations during quiet and relaxed standing. The results show that multiplicative interactions are present during periods with large movements of the center of gravity, where the movements of the centre of gravity and centre of pressure couple into coordinative structures. Further application and interpretation of the developed method for the study of human movement variability are discussed.

Introduction to multifractal detrended fluctuation analysis in matlab

Fractal structures are found in biomedical time series from a wide range of physiological phenomena. The multifractal spectrum identifies the deviations in fractal structure within time periods with large and small fluctuations. The present tutorial is an introduction to multifractal detrended fluctuation analysis (MFDFA) that estimates the multifractal spectrum of biomedical time series. The tutorial presents MFDFA step-by-step in an interactive Matlab session. All Matlab tools needed are available in Introduction to MFDFA folder at the website www.ntnu.edu/inm/geri/software. MFDFA are introduced in Matlab code boxes where the reader can employ pieces of, or the entire MFDFA to example time series. After introducing MFDFA, the tutorial discusses the best practice of MFDFA in biomedical signal processing. The main aim of the tutorial is to give the reader a simple self-sustained guide to the implementation of MFDFA and interpretation of the resulting multifractal spectra.

Correlating brain blood oxygenation level dependent (BOLD) fractal dimension mapping with magnetic resonance spectroscopy (MRS) in Alzheimer's disease

OBJECTIVES

To correlate temporal fractal structure of resting state blood oxygen level dependent (rsBOLD) functional magnetic resonance imaging (fMRI) with in vivo proton magnetic resonance spectroscopy ((1)H-MRS), in Alzheimer's disease (AD) and healthy age-matched normal controls (NC).

MATERIALS AND METHODS

High temporal resolution (4 Hz) rsBOLD signal and single voxel (left putamen) magnetic resonance spectroscopy data was acquired in 33 AD patients and 13 NC. The rsBOLD data was analyzed using two types of fractal dimension (FD) analysis based on relative dispersion and frequency power spectrum. Comparisons in FD were performed between AD and NC, and FD measures were correlated with (1)H-MRS findings.

RESULTS

Temporal fractal analysis of rsBOLD, was able to differentiate AD from NC subjects (P = 0.03). Low FD correlated with markers of AD severity including decreased concentrations of N-acetyl aspartate (R = 0.44, P = 0.015) and increased myoinositol (mI) (R = -0.45, P = 0.012).

CONCLUSION

Based on these results we suggest fractal analysis of rsBOLD could provide an early marker of AD.

Infra-slow fluctuations in electrophysiological recordings, blood-oxygenation-level-dependent signals, and psychophysical time series

Converging electrophysiological and neuroimaging data show that mammalian brain dynamics are governed by spontaneous modulations of neuronal activity levels in cortical and subcortical structures. The time scales of these fluctuations form a continuum from seconds to tens and hundreds of seconds corresponding to slow (0.1-1Hz), infra-slow (0.01-0.1Hz), and "ultradian" (<0.01Hz) frequency bands, respectively. We focus here on the spontaneous neuronal dynamics in the infra-slow frequency band, infra-slow fluctuations (ISFs), and explore their electrophysiological substrates and behavioral correlates. Although electrophysiological ISFs and the associated infra-slow modulations of fast (here, >1Hz) neuronal activities have been recognized on numerous occasions since late 50's, a resurgence in interest towards this frequency band has been driven by a discovery that ISFs in blood-oxygenation-level dependent (BOLD) signals are correlated among specific constellations of brain regions, which constitute intrinsic connectivity networks and define the dynamic architecture of spontaneous brain activity at large. Importantly, electrophysiological and BOLD signal ISFs are directly correlated both with ISFs in amplitudes of fast neuronal activities and with ISFs in behavioral performance. Moreover, both electrophysiological and neuroimaging data suggest that the apparently scale-free ISFs may arise from more local quasi-periodic infra-slow oscillations with a contribution of time-scale-specific cellular-level mechanisms. We conclude that ISFs in electrophysiological recordings, BOLD signals, neuronal activity levels, and behavioral time series are likely to reflect the same underlying phenomenon; a superstructure of interacting and transiently oscillatory ISFs that regulate both the integration within and decoupling between concurrently active neuronal communities.