Power spectrum scale invariance quantifies limbic dysregulation in trait anxious adults using fMRI: adapting methods optimized for characterizing autonomic dysregulation to neural dynamic time series

Resting-State EEG Alterations of Practice-Related Spectral Activity and Connectivity Patterns in Depression

BACKGROUND

Depression presents with altered energy regulation and neural plasticity. Previous electroencephalography (EEG) studies showed that practice in learning tasks increases power in beta range (13-30 Hz) in healthy subjects but not in those with impaired plasticity. Here, we ascertain whether depression presents with alterations of spectral activity and connectivity before and after a learning task.

METHODS

We used publicly available resting-state EEG recordings (64 electrodes) from 122 subjects. Based on Beck Depression Inventory (BDI) scores, they were assigned to either a high BDI (hBDI, BDI > 13, N = 46) or a control (CTL, BDI < 7, N = 75) group. We analyzed spectral activity, theta-beta, and theta-gamma phase-amplitude coupling (PAC) of EEG recorded at rest before and after a learning task.

RESULTS

At baseline, compared to CTL, hBDI exhibited greater power in beta over fronto-parietal regions and in gamma over the right parieto-occipital area. At post task, power increased in all frequency ranges only in CTL. Theta-beta and theta-gamma PAC were greater in hBDI at baseline but not after the task.

CONCLUSIONS

The lack of substantial post-task growth of beta power in depressed subjects likely represents power saturation due to greater baseline values. We speculate that inhibitory/excitatory imbalance, altered plasticity mechanisms, and energy dysregulation present in depression may contribute to this phenomenon.

Frequency-dependent alterations of global signal topography in patients with major depressive disorder

BACKGROUND

Major depressive disorder (MDD) is associated not only with disorders in multiple brain networks but also with frequency-specific brain activities. The abnormality of spatiotemporal networks in patients with MDD remains largely unclear.

METHODS

We investigated the alterations of the global spatiotemporal network in MDD patients using a large-sample multicenter resting-state functional magnetic resonance imaging dataset. The spatiotemporal characteristics were measured by the variability of global signal (GS) and its correlation with local signals (GSCORR) at multiple frequency bands. The association between these indicators and clinical scores was further assessed.

RESULTS

The GS fluctuations were reduced in patients with MDD across the full frequency range (0-0.1852 Hz). The GSCORR was also reduced in the MDD group, especially in the relatively higher frequency range (0.0728-0.1852 Hz). Interestingly, these indicators showed positive correlations with depressive scores in the MDD group and relative negative correlations in the control group.

CONCLUSION

The GS and its spatiotemporal effects on local signals were weakened in patients with MDD, which may impair inter-regional synchronization and related functions. Patients with severe depression may use the compensatory mechanism to make up for the functional impairments.

EEG spectral slope: A reliable indicator for continuous evaluation of consciousness levels during propofol anesthesia

The level of consciousness undergoes continuous alterations during anesthesia. Prior to the onset of propofol-induced complete unconsciousness, degraded levels of behavioral responsiveness can be observed. However, a reliable index to monitor altered consciousness levels during anesthesia has not been sufficiently investigated. In this study, we obtained 60-channel EEG data from 24 healthy participants during an ultra-slow propofol infusion protocol starting with an initial concentration of 1 μg/ml and a stepwise increase of 0.2 μg/ml in concentration. Consecutive auditory stimuli were delivered every 5 to 6 s, and the response time to the stimuli was used to assess the responsiveness levels. We calculated the spectral slope in a time-resolved manner by extracting 5-second EEG segments at each auditory stimulus and estimated their correlation with the corresponding response time. Our results demonstrated that during slow propofol infusion, the response time to external stimuli increased, while the EEG spectral slope, fitted at 15-45 Hz, became steeper, and a significant negative correlation was observed between them. Moreover, the spectral slope further steepened at deeper anesthetic levels and became flatter during anesthesia recovery. We verified these findings using an external dataset. Additionally, we found that the spectral slope of frontal electrodes over the prefrontal lobe had the best performance in predicting the response time. Overall, this study used a time-resolved analysis to suggest that the EEG spectral slope could reliably track continuously altered consciousness levels during propofol anesthesia. Furthermore, the frontal spectral slope may be a promising index for clinical monitoring of anesthesia depth.

Hierarchical network structure as the source of hierarchical dynamics (power-law frequency spectra) in living and non-living systems: How state-trait continua (body plans, personalities) emerge from first principles in biophysics

Living systems are hierarchical control systems that display a small world network structure. In such structures, many smaller clusters are nested within fewer larger ones, producing a fractal-like structure with a 'power-law' cluster size distribution (a mereology). Just like their structure, the dynamics of living systems shows fractal-like qualities: the timeseries of inner message passing and overt behavior contain high frequencies or 'states' (treble) that are nested within lower frequencies or 'traits' (bass), producing a power-law frequency spectrum that is known as a 'state-trait continuum' in the behavioral sciences. Here, we argue that the power-law dynamics of living systems results from their power-law network structure: organisms 'vertically encode' the deep spatiotemporal structure of their (anticipated) environments, to the effect that many small clusters near the base of the hierarchy produce high frequency signal changes and fewer larger clusters at its top produce ultra-low frequencies. Such ultra-low frequencies exert a tonic regulatory pressure that produces morphological as well as behavioral traits (i.e., body plans and personalities). Nested-modular structure causes higher frequencies to be embedded within lower frequencies, producing a power-law state-trait continuum. At the heart of such dynamics lies the need for efficient energy dissipation through networks of coupled oscillators, which also governs the dynamics of non-living systems (e.q., earthquakes, stock market fluctuations). Since hierarchical structure produces hierarchical dynamics, the development and collapse of hierarchical structure (e.g., during maturation and disease) should leave specific traces in system dynamics (shifts in lower frequencies, i.e. morphological and behavioral traits) that may serve as early warning signs to system failure. The applications of this idea range from (bio)physics and phylogenesis to ontogenesis and clinical medicine.

Persistent post-COVID headache is associated with suppression of scale-free functional brain dynamics in non-hospitalized individuals

INTRODUCTION

Post-acute coronavirus disease 2019 (COVID-19) syndrome (PACS) is a growing concern, with headache being a particularly debilitating symptom with high prevalence. The long-term effects of COVID-19 and post-COVID headache on brain function remain poorly understood, particularly among non-hospitalized individuals. This study focused on the power-law scaling behavior of functional brain dynamics, indexed by the Hurst exponent (H). This measure is suppressed during physiological and psychological distress and was thus hypothesized to be reduced in individuals with post-COVID syndrome, with greatest reductions among those with persistent headache.

METHODS

Resting-state blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging data were collected for 57 individuals who had COVID-19 (32 with no headache, 14 with ongoing headache, 11 recovered) and 17 controls who had cold and flu-like symptoms but  tested negative for COVID-19. Individuals were assessed an average of 4-5 months after COVID testing, in a cross-sectional, observational study design.

RESULTS

No significant differences in H values were found between non-headache COVID-19 and control groups., while those with ongoing headache had significantly reduced H values, and those who had recovered from headache had elevated H values, relative to non-headache groups. Effects were greatest in temporal, sensorimotor, and insular brain regions. Reduced H in these regions was also associated with decreased BOLD activity and local functional connectivity.

CONCLUSIONS

These findings provide new insights into the neurophysiological mechanisms that underlie persistent post-COVID headache, with reduced BOLD scaling as a potential biomarker that is specific to this debilitating condition.

Empirical mode decomposition of local field potential data from optogenetic experiments

Introduction

This study investigated the effects of cocaine administration and parvalbumin-type interneuron stimulation on local field potentials (LFPs) recorded in vivo from the medial prefrontal cortex (mPFC) of six mice using optogenetic tools.

Methods

The local network was subject to a brief 10 ms laser pulse, and the response was recorded for 2 s over 100 trials for each of the six subjects who showed stable coupling between the mPFC and the optrode. Due to the strong non-stationary and nonlinearity of the LFP, we used the adaptive, data-driven, Empirical Mode Decomposition (EMD) method to decompose the signal into orthogonal Intrinsic Mode Functions (IMFs).

Results

Through trial and error, we found that seven is the optimum number of orthogonal IMFs that overlaps with known frequency bands of brain activity. We found that the Index of Orthogonality (IO) of IMF amplitudes was close to zero. The Index of Energy Conservation (IEC) for each decomposition was close to unity, as expected for orthogonal decompositions. We found that the power density distribution vs. frequency follows a power law with an average scaling exponent of ~1.4 over the entire range of IMF frequencies 2-2,000 Hz.

Discussion

The scaling exponent is slightly smaller for cocaine than the control, suggesting that neural activity avalanches under cocaine have longer life spans and sizes.

Multiple-Resampling Cross-Spectral Analysis: An Unbiased Tool for Estimating Fractal Connectivity With an Application to Neurophysiological Signals

Investigating scale-free (i.e., fractal) functional connectivity in the brain has recently attracted increasing attention. Although numerous methods have been developed to assess the fractal nature of functional coupling, these typically ignore that neurophysiological signals are assemblies of broadband, arrhythmic activities as well as oscillatory activities at characteristic frequencies such as the alpha waves. While contribution of such rhythmic components may bias estimates of fractal connectivity, they are also likely to represent neural activity and coupling emerging from distinct mechanisms. Irregular-resampling auto-spectral analysis (IRASA) was recently introduced as a tool to separate fractal and oscillatory components in the power spectrum of neurophysiological signals by statistically summarizing the power spectra obtained when resampling the original signal by several non-integer factors. Here we introduce multiple-resampling cross-spectral analysis (MRCSA) as an extension of IRASA from the univariate to the bivariate case, namely, to separate the fractal component of the cross-spectrum between two simultaneously recorded neural signals by applying the same principle. MRCSA does not only provide a theoretically unbiased estimate of the fractal cross-spectrum (and thus its spectral exponent) but also allows for computing the proportion of scale-free coupling between brain regions. As a demonstration, we apply MRCSA to human electroencephalographic recordings obtained in a word generation paradigm. We show that the cross-spectral exponent as well as the proportion of fractal coupling increases almost uniformly over the cortex during the rest-task transition, likely reflecting neural desynchronization. Our results indicate that MRCSA can be a valuable tool for scale-free connectivity studies in characterizing various cognitive states, while it also can be generalized to other applications outside the field of neuroscience.

ENIGMA-anxiety working group: Rationale for and organization of large-scale neuroimaging studies of anxiety disorders

Anxiety disorders are highly prevalent and disabling but seem particularly tractable to investigation with translational neuroscience methodologies. Neuroimaging has informed our understanding of the neurobiology of anxiety disorders, but research has been limited by small sample sizes and low statistical power, as well as heterogenous imaging methodology. The ENIGMA-Anxiety Working Group has brought together researchers from around the world, in a harmonized and coordinated effort to address these challenges and generate more robust and reproducible findings. This paper elaborates on the concepts and methods informing the work of the working group to date, and describes the initial approach of the four subgroups studying generalized anxiety disorder, panic disorder, social anxiety disorder, and specific phobia. At present, the ENIGMA-Anxiety database contains information about more than 100 unique samples, from 16 countries and 59 institutes. Future directions include examining additional imaging modalities, integrating imaging and genetic data, and collaborating with other ENIGMA working groups. The ENIGMA consortium creates synergy at the intersection of global mental health and clinical neuroscience, and the ENIGMA-Anxiety Working Group extends the promise of this approach to neuroimaging research on anxiety disorders.

Comparing Autoregressive and Network Features for Classification of Depression and Anxiety

Autocorrelation in functional MRI (fMRI) time series has been studied for decades, mostly considered as noise in the time series which is removed via prewhitening with an autoregressive model. Recent results suggest that the coefficients of an autoregressive model t to fMRI data may provide an indicator of underlying brain activity, suggesting that prewhitening could be removing important diagnostic information. This paper explores the explanatory value of these autoregressive features extracted from fMRI by considering the use of these features in a classification task. As a point of comparison, functional network based features are extracted from the same data and used in the same classification task. We find that in most cases, network based features provide better classification accuracy. However, using principal component analysis to combine network based features and autoregressive features for classification based on a support vector machine provides improved classification accuracy compared to single features or network features, suggesting that when properly combined there may be additional information to be gained from autoregressive features.

Trait anxiety predicts amygdalar responses during direct processing of threat-related pictures

Previous studies on the associations between trait anxiety and amygdalar responses to threat stimuli have resulted in mixed findings, possibly due to sample characteristics, specific tasks, and analytical methods. The present functional magnetic resonance imaging (fMRI) study aimed to investigate linear or non-linear associations between trait anxiety and amygdalar responses in a sample of participants with low, medium, and high trait anxiety scores. During scanning, participants were presented with threat-related or neutral pictures and had either to solve an emotional task or an emotional-unrelated distraction task. Results showed that only during the explicit task trait anxiety was associated with right amygdalar responses to threat-related pictures as compared to neutral pictures. The best model was a cubic model with increased amygdala responses for very low and medium trait anxiety values but decreased amygdala activation for very high trait anxiety values. The findings imply a non-linear relation between trait anxiety and amygdala activation depending on task conditions.

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.

The Continuum Between Temperament and Mental Illness as Dynamical Phases and Transitions

The full range of biopsychosocial complexity is mind-boggling, spanning a vast range of spatiotemporal scales with complicated vertical, horizontal, and diagonal feedback interactions between contributing systems. It is unlikely that such complexity can be dealt with by a single model. One approach is to focus on a narrower range of phenomena which involve fewer systems but still cover the range of spatiotemporal scales. The suggestion is to focus on the relationship between temperament in healthy individuals and mental illness, which have been conjectured to lie along a continuum of neurobehavioral regulation involving neurochemical regulatory systems (e.g., monoamine and acetylcholine, opiate receptors, neuropeptides, oxytocin), and cortical regulatory systems (e.g., prefrontal, limbic). Temperament and mental illness are quintessentially dynamical phenomena, and need to be addressed in dynamical terms. A meteorological metaphor suggests similarities between temperament and chronic mental illness and climate, between individual behaviors and weather, and acute mental illness and frontal weather events. The transition from normative temperament to chronic mental illness is analogous to climate change. This leads to the conjecture that temperament and chronic mental illness describe distinct, high level, dynamical phases. This suggests approaching biopsychosocial complexity through the study of dynamical phases, their order and control parameters, and their phase transitions. Unlike transitions in physical systems, these biopsychosocial phase transitions involve information and semiotics. The application of complex adaptive dynamical systems theory has led to a host of markers including geometrical markers (periodicity, intermittency, recurrence, chaos) and analytical markers such as fluctuation spectroscopy, scaling, entropy, recurrence time. Clinically accessible biomarkers, in particular heart rate variability and activity markers have been suggested to distinguish these dynamical phases and to signal the presence of transitional states. A particular formal model of these dynamical phases will be presented based upon the process algebra, which has been used to model information flow in complex systems. In particular it describes the dual influences of energy and information on the dynamics of complex systems. The process algebra model is well-suited for dealing with the particular dynamical features of the continuum, which include transience, contextuality, and emergence. These dynamical phases will be described using the process algebra model and implications for clinical practice will be discussed.

Scaling Invariance of Sports Sex Gap

The controversy over the evolution of sex gap in sports stems from the reported that women’s performance will 1 day overtake men’s in the journal Nature. After debate, the recent studies suggest that the sports sex gap has been stable for a long time, due to insurmountable physiological differences. To find a mathematical model that accurately describes this stable gap, we analyze the best annual records of men and women in 25 events from 1992 to 2017, and find that power-law relationship could be acted as the best choice, with an R-squares as high as 0.999 (p ≤ 0.001). Then, based on the power law model, we use the records of men in 2018 to predict the performance of women in that year and compare them with real records. The results show that the deviation rate of the predicted value is only about 2.08%. As a conclusion, it could be said that there is a constant sex gap in sports, and the records of men and women evolve in parallel. This finding could serve as another quantitative rule in biology.

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.

Neurovisceral integration in the executive control network: A resting state analysis

Neurovisceral integration models emphasize the role of frontal lobes in cognitive, behavioral, and emotional regulation. Two candidate hubs for the regulation of cardio-autonomic control, anxiety, and executive attention are the dorsolateral prefrontal cortex (DLPFC) and middle frontal gyrus (MFG). Two-hundred and seventy-one adults (62.9 % female) aged 18-85 years were selected from the NKI-Rockland Sample. Resting state functional imaging data was preprocessed, and seeds extracted from bilateral DLPFC and MFG to test 4 regression models predicting connectivity with high frequency HRV (HF-HRV), trait anxiety (TA), and reaction time on an executive attention task. After controlling for age, sex, body mass index and head motion, the right DLPFC-MFG seed pair provided strongest support for neurovisceral integration indexed by HF-HRV, low TA and shorter reaction time on the attention network task. This hemispheric effect may underlie the inhibitory role of right PFC in the regulation of cardio-autonomic function, emotion, and executive attention.

Perspectives on the use of electrostimulation with the device “VEB”® in the management of disorders related to COVID-19

Background. One of the symptoms of COVID-19 is the so-called "cytokine storm". Its pathogenesis is that the initial release by lymphocytes and macrophages of proinflammatory cytokines in the classical immune response to SARS-CoV-2 is significantly enhanced and maintained due to excessive adrenergic stimulation of the immune cells. The proinflammatory adrenergic mechanism of the "cytokine storm" can be offset by the activation of the anti-inflammatory cholinergic mechanism by non-invasive stimulation of the vagus nerve. In 2015, a generator for electrotherapy and stimulation oh human nerve centers was created, called “VEB-1”®. Preliminary observation of volunteers revealed a modulating effect of a four-day course of electrical stimulation on the parameters of electroencephalogram, metabolism, as well as gas-discharge visualization (GDV). We hypothesized that changes in EEG parameters may be accompanied by a vagotonic shift of the sympatho-vagus balance, favorable for calming the “cytokine storm”. The main purpose of this study was to find out. In addition, concomitant changes in EEG, immunity, GDV, etc. due to the use of the devices "VEB-1"® and recently designed "VEB-2" had to be detected. Material and research methods. The object of observation were 18 volunteers: 11 women 33-62 y and 7 men 29-62 y (Mean±SD: 51±12 y) without clinical diagnose but with dysfunction of neuro-endocrine-immune complex and metabolism. In the morning registered HRV (“CardioLab+HRV”, “KhAI-Medica”, Kharkiv, UA), EEG (“NeuroCom Standard”, “KhAI-Medica”, Kharkiv, UA), kirlianogram by the method of GDV (“GDV Chamber”, “Biotechprogress”, SPb, RF), electroconductivity of skin in three pairs of points of acupuncture (“Medissa”), electrokinetic index of buccal epithelium ("Biotest", Kharkiv State University), as well as some parameters of immunity and metabolism. After the initial testing, an electrical stimulation session was performed with a “VEB-1”® or a “VEB-2” devices. The next morning after completing the four-day course, retesting was performed. Results. The effects of electrical stimulation can be divided into the following networks. Regarding EEG, this is a leveling of right-hand lateralization and normalizing decrease in the increased of the amplitude of the θ-rhythm and its spectral power density (SPD) at the loci F3, F7, F8, T3, T4, T6, P3, O1 and O2; further increase of SPD of δ-rhythm in loci F3, F4, T6, P3 and O1 as well as further decrease of SPD F4-α; reversion of the increased level of entropy in loci Fp1, F4, C3 and P3 to the lowered level. Regarding HRV, it is a vagotonic shift of sympatho-vagus balance due to a decrease in elevated levels of sympathetic tone markers and an increase in decreased levels of vagus tone markers, but without normalization. Neurotropic effects are accompanied by favorable changes in a number of immune parameters and a tendency to decrease the level of C-Reactive Protein. Regarding GDV, it is almost complete normalization of the initially increased GDI Area in the frontal projection and third Chakra Energy; normalizing decrease in the initially increased Energy of second and seventh Chakras; normalizing right-hand shift of more or less pronounced left-sided Asymmetry of first and third Chakra. These effects should be clearly interpreted as physiologically beneficial. The effects on these parameters are almost equally pronounced in people of both sexes when using both devices. Conclusion. Vagotonic and immunotropic effects of our device give us a reason to offer it for further research on the leveling of “cytokine storm” in patients with COVID-19.

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 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.

Dynamic changes in rhythmic and arrhythmic neural signatures in the subthalamic nucleus induced by anaesthesia and tracheal intubation

BACKGROUND

Subcortical structures, including the basal ganglia, have been proposed to be crucial for arousal, consciousness, and behavioural responsiveness. How the basal ganglia contribute to the loss and recovery of consciousness during anaesthesia has, however, not yet been well characterised.

METHODS

Twelve patients with advanced Parkinson's disease, who were undergoing deep brain stimulation (DBS) electrode implantation in the subthalamic nucleus (STN), were included in this study. Local field potentials (LFPs) were recorded from the DBS electrodes and EEG was recorded from the scalp during induction of general anaesthesia (with propofol and sufentanil) and during tracheal intubation. Neural signatures of loss of consciousness and of the expected arousal during intubation were sought in the STN and EEG recordings.

RESULTS

Propofol-sufentanil anaesthesia resulted in power increases in delta, theta, and alpha frequencies, and broadband power decreases in higher frequencies in both STN and frontal cortical areas. This was accompanied by increased STN-frontal cortical coherence only in the alpha frequency band (119 [68]%; P=0.0049). We observed temporal activity changes in STN after tracheal intubation, including power increases in high-beta (22-40 Hz) frequency (98 [123]%; P=0.0064) and changes in the power-law exponent in the power spectra at lower frequencies (2-80 Hz), which were not observed in the frontal cortex. During anaesthesia, the dynamic changes in the high-gamma power in STN LFPs correlated with the power-law exponent in the power spectra at lower frequencies (2-80 Hz).

CONCLUSIONS

Apart from similar activity changes in both STN and cortex associated with anaesthesia-induced unresponsiveness, we observed specific neuronal activity changes in the STN in response to the anaesthesia and tracheal intubation. We also show that the power-law exponent in the power spectra in the STN was modulated by tracheal intubation in anaesthesia. Our results support the hypothesis that subcortical nuclei may play an important role in the loss and return of responsiveness.

Uncovering complex central autonomic networks at rest: a functional magnetic resonance imaging study on complex cardiovascular oscillations

This study aims to uncover brain areas that are functionally linked to complex cardiovascular oscillations in resting-state conditions. Multi-session functional magnetic resonance imaging (fMRI) and cardiovascular data were gathered from 34 healthy volunteers recruited within the human connectome project (the '100-unrelated subjects' release). Group-wise multi-level fMRI analyses in conjunction with complex instantaneous heartbeat correlates (entropy and Lyapunov exponent) revealed the existence of a specialized brain network, i.e. a complex central autonomic network (CCAN), reflecting what we refer to as complex autonomic control of the heart. Our results reveal CCAN areas comprised the paracingulate and cingulate gyri, temporal gyrus, frontal orbital cortex, planum temporale, temporal fusiform, superior and middle frontal gyri, lateral occipital cortex, angular gyrus, precuneous cortex, frontal pole, intracalcarine and supracalcarine cortices, parahippocampal gyrus and left hippocampus. The CCAN visible at rest does not include the insular cortex, thalamus, putamen, amygdala and right caudate, which are classical CAN regions peculiar to sympatho-vagal control. Our results also suggest that the CCAN is mainly involved in complex vagal control mechanisms, with possible links with emotional processing networks.

The Similarity Between Chinese Five-Pattern and Eysenck's Personality Traits: Evidence From Theory and Resting-State fMRI

Chinese five-pattern and Eysenck’s personality traits are two types of personality theories based on different cultural backgrounds. The former is an indigenous theory, and the latter is a cross-cultural theory. In order to verify the relationship between two different personality traits from theory and neuropsychology, the current study recruited 170 healthy adults to calculate their five-Pattern Personality Inventory (FPPI) and Eysenck Personality Questionnaire-Revised (EPQ) scales and to scan their brains using functional magnetic resonance imaging (fMRI). Then, we performed stepwise-regression analysis and mediation-effect analysis to explore the association between brain regional homogeneity (ReHo) and two types of personality traits. The results showed that the ReHo of the right superior temporal gyrus (STG) positively correlated with TaiYang traits for FPPI and that there was a significant linear relationship with extraversion and neuroticism for EPQ. Besides, the ReHo of the right medial prefrontal cortex (mPFC) positively correlated with TaiYin for FPPI, and it also showed a significant linear relationship with neuroticism for EPQ. Furthermore, we found that extroversion and neuroticism partially mediated the relationship between five-pattern personality traits and the regional brain function, based on the mediation-effect analysis. Our findings suggest that Chinese five-pattern personality traits have a close relationship with Eysenck’s personality traits and that both may be engaged in similar neurobiological mechanisms in common brain regions to some extent. Hence, these findings first reveal a relationship between Chinese traditional personality traits and Western Eysenck’s personality traits in terms of both theoretical and neurobiological contexts.

Analysis of brain connectivity during nitrous oxide sedation using graph theory

Nitrous oxide, the least potent inhalation anesthetic, is widely used for conscious sedation. Recently, it has been reported that the occurrence of anesthetic-induced loss of consciousness decreases the interconnection between brain regions, resulting in brain network changes. However, few studies have investigated these changes in conscious sedation using nitrous oxide. Therefore, the present study aimed to use graph theory to analyze changes in brain networks during nitrous oxide sedation. Participants were 20 healthy volunteers (10 men and 10 women, 20–40 years old) with no history of systemic disease. We acquired electroencephalogram (EEG) recordings of 32 channels during baseline, nitrous oxide inhalation sedation, and recovery. EEG epochs from the baseline and the sedation state (50% nitrous oxide) were extracted and analyzed with the network connection parameters of graph theory. Analysis of 1/f dynamics, revealed a steeper slope while in the sedation state than during the baseline. Network connectivity parameters showed significant differences between the baseline and sedation state, in delta, alpha1, alpha2, and beta2 frequency bands. The most pronounced differences in functional distance during nitrous oxide sedation were observed in the alpha1 and alpha2 frequency bands. Change in 1/f dynamics indicates that changes in brain network systems occur during nitrous oxide administration. Changes in network parameters imply that nitrous oxide interferes with the efficiency of information integration in the frequency bands important for cognitive processes and attention tasks. Alteration of brain network during nitrous oxide administration may be associated to the sedative mechanism of nitrous oxide.

Association between Scale-Free Brain Dynamics and Behavioral Performance: Functional MRI Study in Resting State and Face Processing Task

The scale-free dynamics of human brain activity, characterized by an elaborate temporal structure with scale-free properties, can be quantified using the power-law exponent (PLE) as an index. Power laws are well documented in nature in general, particularly in the brain. Some previous fMRI studies have demonstrated a lower PLE during cognitive-task-evoked activity than during resting state activity. However, PLE modulation during cognitive-task-evoked activity and its relationship with an associated behavior remain unclear. In this functional fMRI study in the resting state and face processing + control task, we investigated PLE during both the resting state and task-evoked activities, as well as its relationship with behavior measured using mean reaction time (mRT) during the task. We found that (1) face discrimination-induced BOLD signal changes in the medial prefrontal cortex (mPFC), posterior cingulate cortex (PCC), amygdala, and fusiform face area; (2) PLE significantly decreased during task-evoked activity specifically in mPFC compared with resting state activity; (3) most importantly, in mPFC, mRT significantly negatively correlated with both resting state PLE and the resting-task PLE difference. These results may lead to a better understanding of the associations between task performance parameters (e.g., mRT) and the scale-free dynamics of spontaneous and task-evoked brain activities.

Motivation but not valence modulates neuroticism-dependent cingulate cortex and insula activity

Neuroticism has been found to specifically modulate amygdala activations during differential processing of valence and motivation while other brain networks yet are unexplored for associated effects. The main purpose of this study was to investigate whether neural mechanisms processing valence or motivation are prone to neuroticism in the salience network (SN), a network that is anchored in the anterior cingulate cortex (ACC) and the anterior insula. This study used functional magnetic resonance imaging (fMRI) and an approach/avoid emotional pictures task to investigate brain activations modulated by pictures' valence or motivational status between high and low neurotic individuals. We found that neuroticism-dependent SN and the parahippocampal-fusiform area activations were modulated by motivation but not valence. Valence in contrast interacted with neuroticism in the lateral orbitofrontal cortex. We suggested that neuroticism modulated valence and motivation processing, however, under the influence of the two distinct networks. Neuroticism modulated the motivation through the SN while it modulated the valence through the orbitofrontal networks.

Causal Neuro-immune Relationships at Patients with Chronic Pyelonephritis and Cholecystitis. Correlations between Parameters EEG, HRV and White Blood Cell Count

We aim to analyze in bounds KJ Tracey's immunological homunculus conception the relationships between parameters of electroencephalogram (EEG) and heart rate variability (HRV), on the one hand, and the parameters of bhite blood cell count, on the other hand.

METHODS

In basal conditions in 23 men, patients with chronic pyelonephritis and cholecystitis in remission, recorded EEG ("NeuroCom Standard", KhAI Medica, Ukraine) and HRV ("Cardiolab+VSR", KhAI Medica, Ukraine). In portion of blood counted up white blood cell count.

RESULTS

Revealed that canonical correlation between constellation EEG and HRV parameters form with blood level of leukocytes 0.92 (p<10-5), with relative content in white blood cell count stubnuclear neutrophiles 0.93 (p<10-5), segmentonucleary neutrophiles 0.89 (p<10-3), eosinophiles 0.87 (p=0.003), lymphocytes 0.77 (p<10-3) and with monocytes 0.75 (p=0.003).

CONCLUSION

Parameters of white blood cell count significantly modulated by electrical activity some structures of central and autonomic nervous systems.

From Anxious to Reckless: A Control Systems Approach Unifies Prefrontal-Limbic Regulation Across the Spectrum of Threat Detection

Here we provide an integrative review of basic control circuits, and introduce techniques by which their regulation can be quantitatively measured using human neuroimaging. We illustrate the utility of the control systems approach using four human neuroimaging threat detection studies (N = 226), to which we applied circuit-wide analyses in order to identify the key mechanism underlying individual variation. In so doing, we build upon the canonical prefrontal-limbic control system to integrate circuit-wide influence from the inferior frontal gyrus (IFG). These were incorporated into a computational control systems model constrained by neuroanatomy and designed to replicate our experimental data. In this model, the IFG acts as an informational set point, gating signals between the primary prefrontal-limbic negative feedback loop and its cortical information-gathering loop. Along the cortical route, if the sensory cortex provides sufficient information to make a threat assessment, the signal passes to the ventromedial prefrontal cortex (vmPFC), whose threat-detection threshold subsequently modulates amygdala outputs. However, if signal outputs from the sensory cortex do not provide sufficient information during the first pass, the signal loops back to the sensory cortex, with each cycle providing increasingly fine-grained processing of sensory data. Simulations replicate IFG (chaotic) dynamics experimentally observed at both ends at the threat-detection spectrum. As such, they identify distinct types of IFG disconnection from the circuit, with associated clinical outcomes. If IFG thresholds are too high, the IFG and sensory cortex cycle for too long; in the meantime the coarse-grained (excitatory) pathway will dominate, biasing ambiguous stimuli as false positives. On the other hand, if cortical IFG thresholds are too low, the inhibitory pathway will suppress the amygdala without cycling back to the sensory cortex for much-needed fine-grained sensory cortical data, biasing ambiguous stimuli as false negatives. Thus, the control systems model provides a consistent mechanism for IFG regulation, capable of producing results consistent with our data for the full spectrum of threat-detection: from fearful to optimal to reckless. More generally, it illustrates how quantitative characterization of circuit dynamics can be used to unify a fundamental dimension across psychiatric affective symptoms, with implications for populations that range from anxiety disorders to addiction.

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.

Clinically Anxious Individuals Show Disrupted Feedback between Inferior Frontal Gyrus and Prefrontal-Limbic Control Circuit

Clinical anxiety is associated with generalization of conditioned fear, in which innocuous stimuli elicit alarm. Using Pavlovian fear conditioning (electric shock), we quantify generalization as the degree to which subjects' neurobiological responses track perceptual similarity gradients to a conditioned stimulus. Previous studies show that the ventromedial prefrontal cortex (vmPFC) inversely and ventral tegmental area directly track the gradient of perceptual similarity to the conditioned stimulus in healthy individuals, whereas clinically anxious individuals fail to discriminate. Here, we extend this work by identifying specific functional roles within the prefrontal-limbic circuit. We analyzed fMRI time-series acquired from 57 human subjects during a fear generalization task using entropic measures of circuit-wide regulation and feedback (power spectrum scale invariance/autocorrelation), in combination with structural (diffusion MRI-probabilistic tractography) and functional (stochastic dynamic causal modeling) measures of prefrontal-limbic connectivity within the circuit. Group comparison and correlations with anxiety severity across 57 subjects revealed dysregulatory dynamic signatures within the inferior frontal gyrus (IFG), which our prior work has linked to impaired feedback within the circuit. Bayesian model selection then identified a fully connected prefrontal-limbic model comprising the IFG, vmPFC, and amygdala. Dysregulatory IFG dynamics were associated with weaker reciprocal excitatory connectivity between the IFG and the vmPFC. The vmPFC exhibited inhibitory influence on the amygdala. Our current results, combined with our previous work across a threat-perception spectrum of 137 subjects and a meta-analysis of 366 fMRI studies, dissociate distinct roles for three prefrontal-limbic regions, wherein the IFG provides evaluation of stimulus meaning, which then informs the vmPFC in inhibiting the amygdala.

SIGNIFICANCE STATEMENT

Affective neuroscience has generally treated prefrontal regions (orbitofrontal cortex, dorsolateral prefrontal cortex, inferior frontal gyrus, ventromedial prefrontal cortex) equivalently as inhibitory components of the prefrontal-limbic system. Yet research across the anxiety spectrum suggests that the inferior frontal gyrus may have a more complex role in emotion regulation, as this region shows abnormal function in disorders of both hyperarousal and hypoarousal. Using entropic measures of circuit-wide regulation and feedback, in combination with measures of structural and functional connectivity, we dissociate distinct roles for three prefrontal-limbic regions, wherein the inferior frontal gyrus provides evaluation of stimulus meaning, which then informs the ventromedial prefrontal cortex in inhibiting the amygdala. This reconfiguration coheres with studies of conceptual disambiguation also implicating the inferior frontal gyrus.

Power spectrum scale invariance as a neural marker of cocaine misuse and altered cognitive control

BACKGROUND

Magnetic resonance imaging (MRI) has highlighted the effects of chronic cocaine exposure on cerebral structures and functions, and implicated the prefrontal cortices in deficits of cognitive control. Recent investigations suggest power spectrum scale invariance (PSSI) of cerebral blood oxygenation level dependent (BOLD) signals as a neural marker of cerebral activity. We examined here how PSSI is altered in association with cocaine misuse and impaired cognitive control.

METHODS

Eighty-eight healthy (HC) and seventy-five age and gender matched cocaine dependent (CD) adults participated in functional MRI of a stop signal task (SST). BOLD images were preprocessed using standard procedures in SPM, including detrending, band-pass filtering (0.01-0.25 Hz), and correction for head motions. Voxel-wise PSSI measures were estimated by a linear fit of the power spectrum with a log-log scale. In group analyses, we examined differences in PSSI between HC and CD, and its association with clinical and behavioral variables using a multiple regression. A critical component of cognitive control is post-signal behavioral adjustment, which is compromised in cocaine dependence. Therefore, we examined the PSSI changes in association with post-signal slowing (PSS) in the SST.

RESULTS

Compared to HC, CD showed decreased PSS and PSSI in multiple frontoparietal regions. PSSI was positively correlated with PSS in HC in multiple regions, including the left inferior frontal gyrus (IFG) and right supramarginal gyrus (SMG), which showed reduced PSSI in CD.

CONCLUSIONS

These findings suggest disrupted connectivity dynamics in the fronto-parietal areas in association with post-signal behavioral adjustment in cocaine addicts. These new findings support PSSI as a neural marker of impaired cognitive control in cocaine addiction.

Graph theoretical analysis of brain connectivity in phantom sound perception

Tinnitus is a phantom sound commonly thought of to be produced by the brain related to auditory deafferentation. The current study applies concepts from graph theory to investigate the differences in lagged phase functional connectivity using the average resting state EEG of 311 tinnitus patients and 256 healthy controls. The primary finding of the study was a significant increase in connectivity in beta and gamma oscillations and a significant reduction in connectivity in the lower frequencies for the tinnitus group. There also seems to be parallel processing of long-distance information between delta, theta, alpha1 and gamma frequency bands that is significantly stronger in the tinnitus group. While the network reorganizes into a more regular topology in the low frequency carrier oscillations, development of a more random topology is witnessed in the high frequency oscillations. In summary, tinnitus can be regarded as a maladaptive ‘disconnection’ syndrome, which tries to both stabilize into a regular topology and broadcast the presence of a deafferentation-based bottom-up prediction error as a result of a top-down prediction.

The temporal structure of resting-state brain activity in the medial prefrontal cortex predicts self-consciousness

Recent studies have demonstrated an overlap between the neural substrate of resting-state activity and self-related processing in the cortical midline structures (CMS). However, the neural and psychological mechanisms mediating this so-called "rest-self overlap" remain unclear. To investigate the neural mechanisms, we estimated the temporal structure of spontaneous/resting-state activity, e.g. its long-range temporal correlations or self-affinity across time as indexed by the power-law exponent (PLE). The PLE was obtained in resting-state activity in the medial prefrontal cortex (MPFC) and the posterior cingulate cortex (PCC) in 47 healthy subjects by functional magnetic resonance imaging (fMRI). We performed correlation analyses of the PLE and Revised Self-Consciousness Scale (SCSR) scores, which enabled us to access different dimensions of self-consciousness and specified rest-self overlap in a psychological regard. The PLE in the MPFC's resting-state activity correlated with private self-consciousness scores from the SCSR. Conversely, we found no correlation between the PLE and the other subscales of the SCSR (public, social) or between other resting-state measures, including functional connectivity, and the SCSR subscales. This is the first evidence for the association between the scale-free dynamics of resting-state activity in the CMS and the private dimension of self-consciousness. This finding implies the relationship of especially the private dimension of self with the temporal structure of resting-state activity.

Using network dynamic fMRI for detection of epileptogenic foci

Background

Epilepsy is one of the most prevalent neurological disorders. It remains medically intractable for about one-third of patients with focal epilepsy, for whom precise localization of the epileptogenic zone responsible for seizure initiation may be critical for successful surgery. Existing fMRI literature points to widespread network disturbances in functional connectivity. Per previous scalp and intracranial EEG studies and consistent with excessive local synchronization during interictal discharges, we hypothesized that, relative to same regions in healthy controls, epileptogenic foci would exhibit less chaotic dynamics, identifiable via entropic analyses of resting state fMRI time series.

Methods

In order to first validate this hypothesis on a cohort of patients with known ground truth, here we test individuals with well-defined epileptogenic foci (left mesial temporal lobe epilepsy). We analyzed voxel-wise resting-state fMRI time-series using the autocorrelation function (ACF), an entropic measure of regulation and feedback, and performed follow-up seed-to-voxel functional connectivity analysis. Disruptions in connectivity of the region exhibiting abnormal dynamics were examined in relation to duration of epilepsy and patients’ cognitive performance using a delayed verbal memory recall task.

Results

ACF analysis revealed constrained (less chaotic) functional dynamics in left temporal lobe epilepsy patients, primarily localized to ipsilateral temporal pole, proximal to presumed focal points. Autocorrelation decay rates differentiated, with 100 % accuracy, between patients and healthy controls on a subject-by-subject basis within a leave-one-subject out classification framework. Regions identified via ACF analysis formed a less efficient network in patients, as compared to controls. Constrained dynamics were linked with locally increased and long-range decreased connectivity that, in turn, correlated significantly with impaired memory (local left temporal connectivity) and epilepsy duration (left temporal – posterior cingulate cortex connectivity).

Conclusions

Our current results suggest that data driven functional MRI methods that target network dynamics hold promise in providing clinically valuable tools for identification of epileptic regions.

Brief learning induces a memory bias for arousing-negative words: an fMRI study in high and low trait anxious persons

Persons suffering from anxiety disorders display facilitated processing of arousing and negative stimuli, such as negative words. This memory bias is reflected in better recall and increased amygdala activity in response to such stimuli. However, individual learning histories were not considered in most studies, a concern that we meet here. Thirty-four female persons (half with high-, half with low trait anxiety) participated in a criterion-based associative word-learning paradigm, in which neutral pseudowords were paired with aversive or neutral pictures, which should lead to a valence change for the negatively paired pseudowords. After learning, pseudowords were tested with fMRI to investigate differential brain activation of the amygdala evoked by the newly acquired valence. Explicit and implicit memory was assessed directly after training and in three follow-ups at 4-day intervals. The behavioral results demonstrate that associative word-learning leads to an explicit (but no implicit) memory bias for negatively linked pseudowords, relative to neutral ones, which confirms earlier studies. Bilateral amygdala activation underlines the behavioral effect: Higher trait anxiety is correlated with stronger amygdala activation for negatively linked pseudowords than for neutrally linked ones. Most interestingly, this effect is also present for negatively paired pseudowords that participants could not remember well. Moreover, neutrally paired pseudowords evoked higher amygdala reactivity than completely novel ones in highly anxious persons, which can be taken as evidence for generalization. These findings demonstrate that few word-learning trials generate a memory bias for emotional stimuli, indexed both behaviorally and neurophysiologically. Importantly, the typical memory bias for emotional stimuli and the generalization to neutral ones is larger in high anxious persons.

Scale-free brain dynamics under physical and psychological distress: pre-treatment effects in women diagnosed with breast cancer

Stressful life events are related to negative outcomes, including physical and psychological manifestations of distress, and behavioral deficits. Patients diagnosed with breast cancer report impaired attention and working memory prior to adjuvant therapy, which may be induced by distress. In this article, we examine whether brain dynamics show systematic changes due to the distress associated with cancer diagnosis. We hypothesized that impaired working memory is associated with suppression of "long-memory" neuronal dynamics; we tested this by measuring scale-free ("fractal") brain dynamics, quantified by the Hurst exponent (H). Fractal scaling refers to signals that do not occur at a specific time-scale, possessing a spectral power curve P(f)∝ f(-β); they are "long-memory" processes, with significant autocorrelations. In a BOLD functional magnetic resonance imaging study, we scanned three groups during a working memory task: women scheduled to receive chemotherapy or radiotherapy and aged-matched controls. Surprisingly, patients' BOLD signal exhibited greater H with increasing intensity of anticipated treatment. However, an analysis of H and functional connectivity against self-reported measures of psychological distress (Worry, Anxiety, Depression) and physical distress (Fatigue, Sleep problems) revealed significant interactions. The modulation of (Worry, Anxiety) versus (Fatigue, Sleep Problems, Depression) showed the strongest effect, where higher worry and lower fatigue was related to reduced H in regions involved in visuospatial search, attention, and memory processing. This is also linked to decreased functional connectivity in these brain regions. Our results indicate that the distress associated with cancer diagnosis alters BOLD scaling, and H is a sensitive measure of the interaction between psychological versus physical distress.

Scale-free brain activity: past, present, and future

Brain activity observed at many spatiotemporal scales exhibits a 1/f-like power spectrum, including neuronal membrane potentials, neural field potentials, noninvasive electroencephalography (EEG), magnetoencephalography (MEG), and functional magnetic resonance imaging (fMRI) signals. A 1/f-like power spectrum is indicative of arrhythmic brain activity that does not contain a predominant temporal scale (hence, 'scale-free'). This characteristic of scale-free brain activity distinguishes it from brain oscillations. Although scale-free brain activity and brain oscillations coexist, our understanding of the former remains limited. Recent research has shed light on the spatiotemporal organization, functional significance, and potential generative mechanisms of scale-free brain activity, as well as its developmental and clinical relevance. A deeper understanding of this prevalent brain signal should provide new insights into, and analytical tools for, cognitive neuroscience.

Network connectivity modulates power spectrum scale invariance

Measures of complexity are sensitive in detecting disease, which has made them attractive candidates for diagnostic biomarkers; one complexity measure that has shown promise in fMRI is power spectrum scale invariance (PSSI). Even if scale-free features of neuroimaging turn out to be diagnostically useful, however, their underlying neurobiological basis is poorly understood. Using modeling and simulations of a schematic prefrontal-limbic meso-circuit, with excitatory and inhibitory networks of nodes, we present here a framework for how network density within a control system can affect the complexity of signal outputs. Our model demonstrates that scale-free behavior, similar to that observed in fMRI PSSI data, can be obtained for sufficiently large networks in a context as simple as a linear stochastic system of differential equations, although the scale-free range improves when introducing more realistic, nonlinear behavior in the system. PSSI values (reflective of complexity) vary as a function of both input type (excitatory, inhibitory) and input density (mean number of long-range connections, or strength), independent of their node-specific geometric distribution. Signals show pink noise (1/f) behavior when excitatory and inhibitory influences are balanced. As excitatory inputs are increased and decreased, signals shift towards white and brown noise, respectively. As inhibitory inputs are increased and decreased, signals shift towards brown and white noise, respectively. The results hold qualitatively at the hemodynamic scale, which we modeled by introducing a neurovascular component. Comparing hemodynamic simulation results to fMRI PSSI results from 96 individuals across a wide spectrum of anxiety-levels, we show how our model can generate concrete and testable hypotheses for understanding how connectivity affects regulation of meso-circuits in the brain.

Human gender differences in the perception of conspecific alarm chemosensory cues

It has previously been established that, in threatening situations, animals use alarm pheromones to communicate danger. There is emerging evidence of analogous chemosensory "stress" cues in humans. For this study, we collected alarm and exercise sweat from "donors," extracted it, pooled it and presented it to 16 unrelated "detector" subjects undergoing fMRI. The fMRI protocol consisted of four stimulus runs, with each combination of stimulus condition and donor gender represented four times. Because olfactory stimuli do not follow the canonical hemodynamic response, we used a model-free approach. We performed minimal preprocessing and worked directly with block-average time series and step-function estimates. We found that, while male stress sweat produced a comparably strong emotional response in both detector genders, female stress sweat produced a markedly stronger arousal in female than in male detectors. Our statistical tests pinpointed this gender-specificity to the right amygdala (strongest in the superficial nuclei). When comparing the olfactory bulb responses to the corresponding stimuli, we found no significant differences between male and female detectors. These imaging results complement existing behavioral evidence, by identifying whether gender differences in response to alarm chemosignals are initiated at the perceptual versus emotional level. Since we found no significant differences in the olfactory bulb (primary processing site for chemosensory signals in mammals), we infer that the specificity in responding to female fear is likely based on processing meaning, rather than strength, of chemosensory cues from each gender.

Small-world network properties in prefrontal cortex correlate with predictors of psychopathology risk in young children: a NIRS study

Near infrared spectroscopy (NIRS) is an emerging imaging technique that is relatively inexpensive, portable, and particularly well suited for collecting data in ecological settings. Therefore, it holds promise as a potential neurodiagnostic for young children. We set out to explore whether NIRS could be utilized in assessing the risk of developmental psychopathology in young children. A growing body of work indicates that temperament at young age is associated with vulnerability to psychopathology later on in life. In particular, it has been shown that low effortful control (EC), which includes the focusing and shifting of attention, inhibitory control, perceptual sensitivity, and a low threshold for pleasure, is linked to conditions such as anxiety, depression and attention deficit hyperactivity disorder (ADHD). Physiologically, EC has been linked to a control network spanning among other sites the prefrontal cortex. Several psychopathologies, such as depression and ADHD, have been shown to result in compromised small-world network properties. Therefore we set out to explore the relationship between EC and the small-world properties of PFC using NIRS. NIRS data were collected from 44 toddlers, ages 3-5, while watching naturalistic stimuli (movie clips). Derived complex network measures were then correlated to EC as derived from the Children's Behavior Questionnaire (CBQ). We found that reduced levels of EC were associated with compromised small-world properties of the prefrontal network. Our results suggest that the longitudinal NIRS studies of complex network properties in young children hold promise in furthering our understanding of developmental psychopathology.

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.

Brain areas controlling heart rate variability in tinnitus and tinnitus-related distress

BACKGROUND

Tinnitus is defined as an intrinsic sound perception that cannot be attributed to an external sound source. Distress in tinnitus patients is related to increased beta activity in the dorsal part of the anterior cingulate and the amount of distress correlates with network activity consisting of the amygdala-anterior cingulate cortex-insula-parahippocampus. Previous research also revealed that distress is associated to a higher sympathetic (OS) tone in tinnitus patients and tinnitus suppression to increased parasympathetic (PS) tone.

METHODOLOGY

The aim of the present study is to investigate the relationship between tinnitus distress and the autonomic nervous system and find out which cortical areas are involved in the autonomic nervous system influences in tinnitus distress by the use of source localized resting state electroencephalogram (EEG) recordings and electrocardiogram (ECG). Twenty-one tinnitus patients were included in this study.

CONCLUSIONS

The results indicate that the dorsal and subgenual anterior cingulate, as well as the left and right insula are important in the central control of heart rate variability in tinnitus patients. Whereas the sympathovagal balance is controlled by the subgenual and pregenual anterior cingulate cortex, the right insula controls sympathetic activity and the left insula the parasympathetic activity. The perceived distress in tinnitus patients seems to be sympathetically mediated.

Neural correlates of personality: an integrative review

This review examines the neural correlates of Gray's model (Gray and McNaughton, 2000; McNaughton and Corr, 2004), supplemented by a fourth dimension: constraint (Carver, 2005). The purpose of this review is to summarize findings from fMRI studies that tap on neural correlates of personality aspects in healthy subjects, in order to provide insight into the neural activity underlying human temperament. BAS-related personality traits were consistently reported to correlate positively to activity of the ventral and dorsal striatum and ventral PFC in response to positive stimuli. FFFS and BIS-related personality traits are positively correlated to activity in the amygdala in response to negative stimuli. There is limited evidence that constraint is associated with PFC and ACC activity. In conclusion, functional MRI research sheds some light on the specific neural networks underlying personality. It is clear that more sophisticated task paradigms are required, as well as personality questionnaires that effectively differentiate between BAS, FFFS, BIS, and constraint. Further research is proposed to potentially reveal new insight in the neural subsystems governing basic human behavior.

The redistribution of power: neurocardiac signaling, alcohol and gender

Human adaptability involves interconnected biological and psychological control processes that determine how successful we are in meeting internal and environmental challenges. Heart rate variability (HRV), the variability in consecutive R-wave to R-wave intervals (RRI) of the electrocardiogram, captures synergy between the brain and cardiovascular control systems that modulate adaptive responding. Here we introduce a qualitatively new dimension of adaptive change in HRV quantified as a redistribution of spectral power by applying the Wasserstein distance with exponent 1 metric (W(1)) to RRI spectral data. We further derived a new index, D, to specify the direction of spectral redistribution and clarify physiological interpretation. We examined gender differences in real time RRI spectral power response to alcohol, placebo and visual cue challenges. Adaptive changes were observed as changes in power of the various spectral frequency bands (i.e., standard frequency domain HRV indices) and, during both placebo and alcohol intoxication challenges, as changes in the structure (shape) of the RRI spectrum, with a redistribution towards lower frequency oscillations. The overall conclusions from the present study are that the RRI spectrum is capable of a fluid and highly flexible response, even when oscillations (and thus activity at the sinoatrial node) are pharmacologically suppressed, and that low frequency oscillations serve a crucial but less studied role in physical and mental health.

Scale-free properties of the functional magnetic resonance imaging signal during rest and task

It has been shown recently that a significant portion of brain electrical field potentials consists of scale-free dynamics. These scale-free brain dynamics contain complex spatiotemporal structures and are modulated by task performance. Here we show that the fMRI signal recorded from the human brain is also scale free; its power-law exponent differentiates between brain networks and correlates with fMRI signal variance and brain glucose metabolism. Importantly, in parallel to brain electrical field potentials, the variance and power-law exponent of the fMRI signal decrease during task activation, suggesting that the signal contains more long-range memory during rest and conversely is more efficient at online information processing during task. Remarkably, similar changes also occurred in task-deactivated brain regions, revealing the presence of an optimal dynamic range in the fMRI signal. The scale-free properties of the fMRI signal and brain electrical field potentials bespeak their respective stationarity and nonstationarity. This suggests that neurovascular coupling mechanism is likely to contain a transformation from nonstationarity to stationarity. In summary, our results demonstrate the functional relevance of scale-free properties of the fMRI signal and impose constraints on future models of neurovascular coupling.

A mouse model of high trait anxiety shows reduced heart rate variability that can be reversed by anxiolytic drug treatment

Increasing evidence suggests that specific physiological measures may serve as biomarkers for successful treatment to alleviate symptoms of pathological anxiety. Studies of autonomic function investigating parameters such as heart rate (HR), HR variability and blood pressure (BP) indicated that HR variability is consistently reduced in anxious patients, whereas HR and BP data show inconsistent results. Therefore, HR and HR variability were measured under various emotionally challenging conditions in a mouse model of high innate anxiety (high anxiety behaviour; HAB) vs. control normal anxiety-like behaviour (NAB) mice. Baseline HR, HR variability and activity did not differ between mouse lines. However, after cued Pavlovian fear conditioning, both elevated tachycardia and increased fear responses were observed in HAB mice compared to NAB mice upon re-exposure to the conditioning stimulus serving as the emotional stressor. When retention of conditioned fear was tested in the home cage, HAB mice again displayed higher fear responses than NAB mice, while the HR responses were similar. Conversely, in both experimental settings HAB mice consistently exhibited reduced HR variability. Repeated administration of the anxiolytic NK1 receptor antagonist L-822429 lowered the conditioned fear response and shifted HR dynamics in HAB mice to a more regular pattern, similar to that in NAB mice. Additional receiver-operating characteristic (ROC) analysis demonstrated the high specificity and sensitivity of HR variability to distinguish between normal and high anxiety trait. These findings indicate that assessment of autonomic response in addition to freezing might be a useful indicator of the efficacy of novel anxiolytic treatments.

The NIRS Analysis Package: noise reduction and statistical inference

Near infrared spectroscopy (NIRS) is a non-invasive optical imaging technique that can be used to measure cortical hemodynamic responses to specific stimuli or tasks. While analyses of NIRS data are normally adapted from established fMRI techniques, there are nevertheless substantial differences between the two modalities. Here, we investigate the impact of NIRS-specific noise; e.g., systemic (physiological), motion-related artifacts, and serial autocorrelations, upon the validity of statistical inference within the framework of the general linear model. We present a comprehensive framework for noise reduction and statistical inference, which is custom-tailored to the noise characteristics of NIRS. These methods have been implemented in a public domain Matlab toolbox, the NIRS Analysis Package (NAP). Finally, we validate NAP using both simulated and actual data, showing marked improvement in the detection power and reliability of NIRS.