Gradual emergence of spontaneous correlated brain activity during fading of general anesthesia in rats: Evidences from fMRI and local field potentials

Chronic full-band recordings with graphene microtransistors as neural interfaces for discrimination of brain states

Brain states such as sleep, anesthesia, wakefulness, or coma are characterized by specific patterns of cortical activity dynamics, from local circuits to full-brain emergent properties. We previously demonstrated that full-spectrum signals, including the infraslow component (DC, direct current-coupled), can be recorded acutely in multiple sites using flexible arrays of graphene solution-gated field-effect transistors (gSGFETs). Here, we performed chronic implantation of 16-channel gSGFET arrays over the rat cerebral cortex and recorded full-band neuronal activity with two objectives: (1) to test the long-term stability of implanted devices; and (2) to investigate full-band activity during the transition across different levels of anesthesia. First, we demonstrate it is possible to record full-band signals with stability, fidelity, and spatiotemporal resolution for up to 5.5 months using chronic epicortical gSGFET implants. Second, brain states generated by progressive variation of levels of anesthesia could be identified as traditionally using the high-pass filtered (AC, alternating current-coupled) spectrogram: from synchronous slow oscillations in deep anesthesia through to asynchronous activity in the awake state. However, the DC signal introduced a highly significant improvement for brain-state discrimination: the DC band provided an almost linear information prediction of the depth of anesthesia, with about 85% precision, using a trained algorithm. This prediction rose to about 95% precision when the full-band (AC + DC) spectrogram was taken into account. We conclude that recording infraslow activity using gSGFET interfaces is superior for the identification of brain states, and further supports the preclinical and clinical use of graphene neural interfaces for long-term recordings of cortical activity.

M-current modulation of cortical slow oscillations: Network dynamics and computational modeling

The slow oscillation is a synchronized network activity expressed by the cortical network in slow wave sleep and under anesthesia. Waking up requires a transition from this synchronized brain state to a desynchronized one. Cholinergic innervation is critical for the transition from slow-wave-sleep to wakefulness, and muscarinic action is largely exerted through the muscarinic-sensitive potassium current (M-current) block. We investigated the dynamical impact of blocking the M-current on slow oscillations, both in cortical slices and in a cortical network computational model. Blocking M-current resulted in an elongation of Up states (by four times) and in a significant firing rate increase, reflecting an increased network excitability, albeit no epileptiform discharges occurred. These effects were replicated in a biophysical cortical model, where a parametric reduction of the M-current resulted in a progressive elongation of Up states and firing rate. All neurons, and not only those modeled with M-current, increased their firing rates due to network recurrency. Further increases in excitability induced even longer Up states, approaching the microarousals described in the transition towards wakefulness. Our results bridge an ionic current with network modulation, providing a mechanistic insight into network dynamics of awakening.

Assessing brain state and anesthesia level with two-photon calcium signals

Brain states, such as wake, sleep, or different depths of anesthesia are usually assessed using electrophysiological techniques, such as the local field potential (LFP) or the electroencephalogram (EEG), which are ideal signals for detecting activity patterns such as asynchronous or oscillatory activities. However, it is technically challenging to have these types of measures during calcium imaging recordings such as two-photon or wide-field techniques. Here, using simultaneous two-photon and LFP measurements, we demonstrate that despite the slower dynamics of the calcium signal, there is a high correlation between the LFP and two-photon signals taken from the neuropil outside neuronal somata. Moreover, we find the calcium signal to be systematically delayed from the LFP signal, and we use a model to show that the delay between the two signals is due to the physical distance between the recording sites. These results suggest that calcium signals alone can be used to detect activity patterns such as slow oscillations and ultimately assess the brain state and level of anesthesia.

Vascular and neuronal effects of general anesthesia on the brain: An fMRI study

BACKGROUND AND PURPOSE

A number of functional magnetic resonance imaging (fMRI) studies rely on application of anesthetic agents during scanning that can modulate and complicate interpretation of the measured hemodynamic blood oxygenation level-dependent (BOLD) response. The purpose of the present study was to investigate the effect of general anesthesia on two main components of BOLD signal including neuronal activity and vascular response.

METHODS

Breath-holding (BH) fMRI was conducted in wakefulness and under anesthesia states in 9 patients with drug-resistant epilepsy who needed to get scanned under anesthesia during laser interstitial thermal therapy. BOLD and BOLD cerebrovascular reactivity (BOLD-CVR) maps were compared using t-test between two states to assess the effect of anesthesia on neuronal activity and vascular factors (p < .05).

RESULTS

Overall, our findings revealed an increase in BOLD-CVR and decrease in BOLD response under anesthesia in several brain regions. The results proposed that the modulatory mechanism of anesthetics on neuronal and vascular components of BOLD signal may work in different ways.

CONCLUSION

This experiment for the first human study showed that anesthesia may play an important role in dissociation between neuronal and vascular responses contributed to hemodynamic BOLD signal using BH fMRI imaging that may assist the implication of general anesthesia and interpretation of outcomes in clinical setting.

The effect of general anesthesia on the test–retest reliability of resting-state fMRI metrics and optimization of scan length

Resting-state functional magnetic resonance imaging (rs-fMRI) has been known as a powerful tool in neuroscience. However, exploring the test–retest reliability of the metrics derived from the rs-fMRI BOLD signal is essential, particularly in the studies of patients with neurological disorders. Here, two factors, namely, the effect of anesthesia and scan length, have been estimated on the reliability of rs-fMRI measurements. A total of nine patients with drug-resistant epilepsy (DRE) requiring interstitial thermal therapy (LITT) were scanned in two states. The first scan was performed in an awake state before surgery on the same patient. The second scan was performed 2 weeks later under general anesthesia necessary for LITT surgery. At each state, two rs-fMRI sessions were obtained that each one lasted 15 min, and the effect of scan length was evaluated. Voxel-wise rs-fMRI metrics, including the amplitude of low-frequency fluctuation (ALFF), the fractional amplitude of low-frequency fluctuation (fALFF), functional connectivity (FC), and regional homogeneity (ReHo), were measured. Intraclass correlation coefficient (ICC) was calculated to estimate the reliability of the measurements in two states of awake and under anesthesia. Overall, it appeared that the reliability of rs-fMRI metrics improved under anesthesia. From the 15-min data, we found mean ICC values in awake state including 0.81, 0.51, 0.65, and 0.84 for ALFF, fALFF, FC, and ReHo, respectively, as well as 0.80, 0.59, 0.83, and 0.88 for ALFF, fALFF, FC, and ReHo, respectively, under anesthesia. Additionally, our findings revealed that reliability increases as the function of scan length. We showed that the optimized scan length to achieve less variability of rs-fMRI measurements was 3.1–7.5 min shorter in an anesthetized, compared to a wakeful state.

Selective blockade of rat brain T-type calcium channels provides insights on neurophysiological basis of arousal dependent resting state functional magnetic resonance imaging signals

A number of studies point to slow (0.1–2 Hz) brain rhythms as the basis for the resting-state functional magnetic resonance imaging (rsfMRI) signal. Slow waves exist in the absence of stimulation, propagate across the cortex, and are strongly modulated by vigilance similar to large portions of the rsfMRI signal. However, it is not clear if slow rhythms serve as the basis of all neural activity reflected in rsfMRI signals, or just the vigilance-dependent components. The rsfMRI data exhibit quasi-periodic patterns (QPPs) that appear to increase in strength with decreasing vigilance and propagate across the brain similar to slow rhythms. These QPPs can complicate the estimation of functional connectivity (FC) via rsfMRI, either by existing as unmodeled signal or by inducing additional wide-spread correlation between voxel-time courses of functionally connected brain regions. In this study, we examined the relationship between cortical slow rhythms and the rsfMRI signal, using a well-established pharmacological model of slow wave suppression. Suppression of cortical slow rhythms led to significant reduction in the amplitude of QPPs but increased rsfMRI measures of intrinsic FC in rats. The results suggest that cortical slow rhythms serve as the basis of only the vigilance-dependent components (e.g., QPPs) of rsfMRI signals. Further attenuation of these non-specific signals enhances delineation of brain functional networks.

Slow waves form expanding, memory-rich mesostates steered by local excitability in fading anesthesia

In the arousal process, the brain restores its integrative activity from the synchronized state of slow wave activity (SWA). The mechanisms underpinning this state transition remain, however, to be elucidated. Here we simultaneously probed neuronal assemblies throughout the whole cortex with micro-electrocorticographic recordings in mice. We investigated the progressive shaping of propagating SWA at different levels of isoflurane. We found a form of memory of the wavefront shapes at deep anesthesia, tightly alternating posterior-anterior-posterior patterns. At low isoflurane, metastable patterns propagated in more directions, reflecting an increased complexity. The wandering across these mesostates progressively increased its randomness, as predicted by simulations of a network of spiking neurons, and confirmed in our experimental data. The complexity increase is explained by the elevated excitability of local assemblies with no modifications of the network connectivity. These results shed new light on the functional reorganization of the cortical network as anesthesia fades out.

•

Complexity of isoflurane-induced slow waves reliably determines anesthesia level

•

In deep anesthesia, the propagation strictly alternates between front-back-front patterns

•

In light anesthesia, there is a continuum of directions and faster propagation

•

Local excitability underpins the cortical reorganization in fading anesthesia

Energy-Based Hierarchical Clustering of Cortical Slow Waves in Multi-Electrode Recordings

The recent development of novel multi-electrode recording technologies has revealed the existence of traveling patterns of cortical activity in many species and under different states of awareness. Among these, slow activation waves occurring under sleep and anesthesia have been widely investigated as they provide unique insights into network features such as excitability, connectivity, structure, and dynamics of the cerebral cortex. Such characterization is usually based on clustering methods which are constrained by a priori assumptions as to the number of clusters to be used or rely on wave-by-wave pattern reconstruction. Here, we introduce a new computational tool based on modal analysis of fluid flows which is robustly applied to multivariate electrophysiological data from cortical networks, namely the Energy-based Hierarchical Waves Clustering method (EHWC). EHWC is composed of three main steps: (1) detecting the occurrence of global waves; (2) reducing the data dimensionality via singular value decomposition; (3) clustering hierarchically the singled-out waves. The analysis does not require the single-channel contribution to the waves, which is a typical bottleneck in this kind of analysis due to the unavoidable intrinsic variability of locally recorded activity. For testing and validation, here we used in vivo extracellular recordings from mice cortex under three different levels of anesthesia. As a result, we found slow waves with an increasing number of propagation modes as the anesthesia level decreases, giving an estimate of the increasing complexity of network dynamics. This and other wave's features replicate and extend the findings from previous literature, paving the way to extend the same approach to non-invasive electrophysiological recordings like EEG and fMRI used clinically for the characterization of brain dynamics and clinical stratification in brain lesions.

Impact of GABAA and GABAB Inhibition on Cortical Dynamics and Perturbational Complexity during Synchronous and Desynchronized States

Quantitative estimations of spatiotemporal complexity of cortical activity patterns are used in the clinic as a measure of consciousness levels, but the cortical mechanisms involved are not fully understood. We used a version of the perturbational complexity index (PCI) adapted to multisite recordings from the ferret (either sex) cerebral cortex in vitro (sPCI) to investigate the role of GABAergic inhibition in cortical complexity. We studied two dynamical states: slow-wave activity (synchronous state) and desynchronized activity, that express low and high causal complexity respectively. Progressive blockade of GABAergic inhibition during both regimes revealed its impact on the emergent cortical activity and on sPCI. Gradual GABAA receptor blockade resulted in higher synchronization, being able to drive the network from a desynchronized to a synchronous state, with a progressive decrease of complexity (sPCI). Blocking GABAB receptors also resulted in a reduced sPCI, in particular when in a synchronous, slow wave state. Our findings demonstrate that physiological levels of inhibition contribute to the generation of dynamical richness and spatiotemporal complexity. However, if inhibition is diminished or enhanced, cortical complexity decreases. Using a computational model, we explored a larger parameter space in this relationship and demonstrate a link between excitatory/inhibitory balance and the complexity expressed by the cortical network.SIGNIFICANCE STATEMENT The spatiotemporal complexity of the activity expressed by the cerebral cortex is a highly revealing feature of the underlying network's state. Complexity varies with physiological brain states: it is higher during awake than during sleep states. But it also informs about pathologic states: in disorders of consciousness, complexity is lower in an unresponsive wakefulness syndrome than in a minimally conscious state. What are the network parameters that modulate complexity? Here we investigate how inhibition, mediated by either GABAA or GABAA receptors, influences cortical complexity. And we do this departing from two extreme functional states: a highly synchronous, slow-wave state, and a desynchronized one that mimics wakefulness. We find that there is an optimal level of inhibition in which complexity is highest.

Attractor competition enriches cortical dynamics during awakening from anesthesia

Slow oscillations (≲ 1 Hz), a hallmark of slow-wave sleep and deep anesthesia across species, arise from spatiotemporal patterns of activity whose complexity increases as wakefulness is approached and cognitive functions emerge. The arousal process constitutes an open window to the unknown mechanisms underlying the emergence of such dynamical richness in awake cortical networks. Here, we investigate the changes in network dynamics as anesthesia fades out in the rat visual cortex. Starting from deep anesthesia, slow oscillations gradually increase their frequency, eventually expressing maximum regularity. This stage is followed by the abrupt onset of an infra-slow (~0.2 Hz) alternation between sleep-like oscillations and activated states. A population rate model reproduces this transition driven by an increased excitability that brings it to periodically cross a critical point. Based on our model, dynamical richness emerges as a competition between two metastable attractor states, a conclusion strongly supported by the data.

Up and Down States During Slow Oscillations in Slow-Wave Sleep and Different Levels of Anesthesia

Slow oscillations are a pattern of synchronized network activity generated by the cerebral cortex. They consist of Up and Down states, which are periods of activity interspersed with periods of silence, respectively. However, even when this is a unique dynamic regime of transitions between Up and Down states, this pattern is not constant: there is a range of oscillatory frequencies (0.1-4 Hz), and the duration of Up vs. Down states during the cycles is variable. This opens many questions. Is there a constant relationship between the duration of Up and Down states? How much do they vary across conditions and oscillatory frequencies? Are there different sub regimes within the slow oscillations? To answer these questions, we aimed to explore a concrete aspect of slow oscillations, Up and Down state durations, across three conditions: deep anesthesia, light anesthesia, and slow-wave sleep (SWS), in the same chronically implanted rats. We found that light anesthesia and SWS have rather similar properties, occupying a small area of the Up and Down state duration space. Deeper levels of anesthesia occupy a larger region of this space, revealing that a large variety of Up and Down state durations can emerge within the slow oscillatory regime. In a network model, we investigated the network parameters that can explain the different points within our bifurcation diagram in which slow oscillations are expressed.

Muting, not fragmentation, of functional brain networks under general anesthesia

Changes in resting-state functional connectivity (rs-FC) under general anesthesia have been widely studied with the goal of identifying neural signatures of consciousness. This work has commonly revealed an apparent fragmentation of whole-brain network structure during unconsciousness, which has been interpreted as reflecting a break-down in connectivity and a disruption of the brain's ability to integrate information. Here we show, by studying rs-FC under varying depths of isoflurane-induced anesthesia in nonhuman primates, that this apparent fragmentation, rather than reflecting an actual change in network structure, can be simply explained as the result of a global reduction in FC. Specifically, by comparing the actual FC data to surrogate data sets that we derived to test competing hypotheses of how FC changes as a function of dose, we found that increases in whole-brain modularity and the number of network communities - considered hallmarks of fragmentation - are artifacts of constructing FC networks by thresholding based on correlation magnitude. Taken together, our findings suggest that deepening levels of unconsciousness are instead associated with the increasingly muted expression of functional networks, an observation that constrains current interpretations as to how anesthesia-induced FC changes map onto existing neurobiological theories of consciousness.

Understanding the Effects of Anesthesia on Cortical Electrophysiological Recordings: A Scoping Review

General anesthesia in animal experiments is an ethical must and is required for all the procedures that are likely to cause more than slight or momentary pain. As anesthetics are known to deeply affect experimental findings, including electrophysiological recordings of brain activity, understanding their mechanism of action is of paramount importance. It is widely recognized that the depth and type of anesthesia introduce significant bias in electrophysiological measurements by affecting the shape of both spontaneous and evoked signals, e.g., modifying their latency and relative amplitude. Therefore, for a given experimental protocol, it is relevant to identify the appropriate anesthetic, to minimize the impact on neuronal circuits and related signals under investigation. This review focuses on the effect of different anesthetics on cortical electrical recordings, examining their molecular mechanisms of action, their influence on neuronal microcircuits and, consequently, their impact on cortical measurements.

Systematic Review: Anesthetic Protocols and Management as Confounders in Rodent Blood Oxygen Level Dependent Functional Magnetic Resonance Imaging (BOLD fMRI)-Part B: Effects of Anesthetic Agents, Doses and Timing

Simple Summary

To understand brain function in rats and mice functional magnetic resonance imaging of the brain is used. With this type of “brain scan” regional changes in blood flow and oxygen consumption are measured as an indirect surrogate for activity of brain regions. Animals are often anesthetized for the experiments to prevent stress and blurred images due to movement. However, anesthesia may alter the measurements, as blood flow within the brain is differently affected by different anesthetics, and anesthetics also directly affect brain function. Consequently, results obtained under one anesthetic protocol may not be comparable with those obtained under another, and/or not representative for awake animals and humans. We have systematically searched the existing literature for studies analyzing the effects of different anesthesia methods or studies that compared anesthetized and awake animals. Most studies reported that anesthetic agents, doses and timing had an effect on functional magnetic resonance imaging results. To obtain results which promote our understanding of brain function, it is therefore essential that a standard for anesthetic protocols for functional magnetic resonance is defined and their impact is well characterized.

Abstract

In rodent models the use of functional magnetic resonance imaging (fMRI) under anesthesia is common. The anesthetic protocol might influence fMRI readouts either directly or via changes in physiological parameters. As long as those factors cannot be objectively quantified, the scientific validity of fMRI in rodents is impaired. In the present systematic review, literature analyzing in rats and mice the influence of anesthesia regimes and concurrent physiological functions on blood oxygen level dependent (BOLD) fMRI results was investigated. Studies from four databases that were searched were selected following pre-defined criteria. Two separate articles publish the results; the herewith presented article includes the analyses of 83 studies. Most studies found differences in BOLD fMRI readouts with different anesthesia drugs and dose rates, time points of imaging or when awake status was compared to anesthetized animals. To obtain scientifically valid, reproducible results from rodent fMRI studies, stable levels of anesthesia with agents suitable for the model under investigation as well as known and objectively quantifiable effects on readouts are, thus, mandatory. Further studies should establish dose ranges for standardized anesthetic protocols and determine time windows for imaging during which influence of anesthesia on readout is objectively quantifiable.

Resting-State functional networks of different topographic representations in the somatosensory cortex of macaque monkeys and humans

Information processing in the brain is mediated through a complex functional network architecture whose comprising nodes integrate and segregate themselves on different timescales. To gain an understanding of the network function it is imperative to identify and understand the network structure with respect to the underlying anatomical connectivity and the topographic organization. Here we show that the previously described resting-state network for the somatosensory area 3b comprises of distinct networks that are characteristic for different topographic representations. Seed-based resting-state functional connectivity analysis in macaque monkeys and humans using BOLD-fMRI signals from the face, the hand and rest of the medial somatosensory representations of area 3b revealed different correlation patterns. Both monkeys and humans have many similarities in the connectivity networks, although the networks are more complex in humans with many more nodes. In both the species face area network has the highest ipsilateral and contralateral connectivity, which included areas 3b and 4, and ventral premotor area. The area 3b hand network included ipsilateral hand representation in area 4. The emergent functional network structures largely reflect the known anatomical connectivity. Our results show that different body part representations in area 3b have independent functional networks perhaps reflecting differences in the behavioral use of different body parts. The results also show that large cortical areas if considered together, do not give a complete and accurate picture of the network architecture.

Brain networks of rats under anesthesia using resting-state fMRI: comparison with dead rats, random noise and generative models of networks

OBJECTIVE

Connectivity networks are crucial to understand the brain resting-state activity using functional magnetic resonance imaging (rs-fMRI). Alterations of these brain networks may highlight important findings concerning the resilience of the brain to different disorders. The focus of this paper is to evaluate the robustness of brain network estimations, discriminate them under anesthesia and compare them to generative models.

APPROACH

The extraction of brain functional connectivity (FC) networks is difficult and biased due to the properties of the data: low signal to noise ratio, high dimension low sample size. We propose to use wavelet correlations to assess FC between brain areas under anesthesia using four anesthetics (isoflurane, etomidate, medetomidine, urethane). The networks are then deduced from the functional connectivity matrices by applying statistical thresholds computed using the number of samples at a given scale of wavelet decomposition. Graph measures are extracted and extensive comparisons with generative models of structured networks are conducted.

MAIN RESULTS

The sample size and filtering are critical to obtain significant correlations values and thereby detect connections between regions. This is necessary to construct networks different from random ones as shown using rs-fMRI brain networks of dead rats. Brain networks under anesthesia on rats have topological features that are mixing small-world, scale-free and random networks. Betweenness centrality indicates that hubs are present in brain networks obtained from anesthetized rats but locations of these hubs are altered by anesthesia.

SIGNIFICANCE

Understanding the effects of anesthesia on brain areas is of particular importance in the context of animal research since animal models are commonly used to explore functions, evaluate lesions or illnesses, and test new drugs. More generally, results indicate that the use of correlations in the context of fMRI signals is robust but must be treated with caution. Solutions are proposed in order to control spurious correlations by setting them to zero.

Modulation of cortical slow oscillatory rhythm by GABAB receptors: an in vitro experimental and computational study

KEY POINTS

We confirm that GABAB receptors (GABAB -Rs) are involved in the termination of Up-states; their blockade consistently elongates Up-states. GABAB -Rs also modulate Down-states and the oscillatory cycle, thus having an impact on slow oscillation rhythm and its regularity. The most frequent effect of GABAB -R blockade is elongation of Down-states and subsequent decrease of oscillatory frequency, with an increased regularity. In a quarter of cases, GABAB -R blockade shortened Down-states and increased oscillatory frequency, changes that are independent of firing rates in Up-states. Our computer model provides mechanisms for the experimentally observed dynamics following blockade of GABAB -Rs, for Up/Down durations, oscillatory frequency and regularity. The time course of excitation, inhibition and adaptation can explain the observed dynamics of the network. This study brings novel insights into the role of GABAB -R-mediated slow inhibition on the slow oscillatory activity, which is considered the default activity pattern of the cortical network.

ABSTRACT

Slow wave oscillations (SWOs) dominate cortical activity during deep sleep, anaesthesia and in some brain lesions. SWOs are composed of periods of activity (Up states) interspersed with periods of silence (Down states). The rhythmicity expressed during SWOs integrates neuronal and connectivity properties of the network and is often altered under pathological conditions. Adaptation mechanisms as well as synaptic inhibition mediated by GABAB receptors (GABAB -Rs) have been proposed as mechanisms governing the termination of Up states. The interplay between these two mechanisms is not well understood, and the role of GABAB -Rs controlling the whole cycle of the SWO has not been described. Here we contribute to its understanding by combining in vitro experiments on spontaneously active cortical slices and computational techniques. GABAB -R blockade modified the whole SWO cycle, not only elongating Up states, but also affecting the subsequent Down state duration. Furthermore, while adaptation tends to yield a rather regular behaviour, we demonstrate that GABAB -R activation desynchronizes the SWOs. Interestingly, variability changes could be accomplished in two different ways: by either shortening or lengthening the duration of Down states. Even when the most common observation following GABAB -Rs blocking is the lengthening of Down states, both changes are expressed experimentally and also in numerical simulations. Our simulations suggest that the sluggishness of GABAB -Rs to follow the excitatory fluctuations of the cortical network can explain these different network dynamics modulated by GABAB -Rs.

Correlated activity favors synergistic processing in local cortical networks in vitro at synaptically relevant timescales

Neural information processing is widely understood to depend on correlations in neuronal activity. However, whether correlation is favorable or not is contentious. Here, we sought to determine how correlated activity and information processing are related in cortical circuits. Using recordings of hundreds of spiking neurons in organotypic cultures of mouse neocortex, we asked whether mutual information between neurons that feed into a common third neuron increased synergistic information processing by the receiving neuron. We found that mutual information and synergistic processing were positively related at synaptic timescales (0.05-14 ms), where mutual information values were low. This effect was mediated by the increase in information transmission-of which synergistic processing is a component-that resulted as mutual information grew. However, at extrasynaptic windows (up to 3,000 ms), where mutual information values were high, the relationship between mutual information and synergistic processing became negative. In this regime, greater mutual information resulted in a disproportionate increase in redundancy relative to information transmission. These results indicate that the emergence of synergistic processing from correlated activity differs according to timescale and correlation regime. In a low-correlation regime, synergistic processing increases with greater correlation, and in a high-correlation regime, synergistic processing decreases with greater correlation.

Dynamic causal modelling of phase-amplitude interactions

Models of coupled phase oscillators are used to describe a wide variety of phenomena in neuroimaging. These models typically rest on the premise that oscillator dynamics do not evolve beyond their respective limit cycles, and hence that interactions can be described purely in terms of phase differences. Whilst mathematically convenient, the restrictive nature of phase-only models can limit their explanatory power. We therefore propose a generalisation of dynamic causal modelling that incorporates both phase and amplitude. This allows for the separate quantifications of phase and amplitude contributions to the connectivity between neural regions. We show, using model-generated data and simulations of coupled pendula, that phase-amplitude models can describe strongly coupled systems more effectively than their phase-only counterparts. We relate our findings to four metrics commonly used in neuroimaging: the Kuramoto order parameter, cross-correlation, phase-lag index, and spectral entropy. We find that, with the exception of spectral entropy, the phase-amplitude model is able to capture all metrics more effectively than the phase-only model. We then demonstrate, using local field potential recordings in rodents and functional magnetic resonance imaging in macaque monkeys, that amplitudes in oscillator models play an important role in describing neural dynamics in anaesthetised brain states.

Bistability, Causality, and Complexity in Cortical Networks: An In Vitro Perturbational Study

Measuring the spatiotemporal complexity of cortical responses to direct perturbations provides a reliable index of the brain's capacity for consciousness in humans under both physiological and pathological conditions. Upon loss of consciousness, the complex pattern of causal interactions observed during wakefulness collapses into a stereotypical slow wave, suggesting that cortical bistability may play a role. Bistability is mainly expressed in the form of slow oscillations, a default pattern of activity that emerges from cortical networks in conditions of functional or anatomical disconnection. Here, we employ an in vitro model to understand the relationship between bistability and complexity in cortical circuits. We adapted the perturbational complexity index applied in humans to electrically stimulated cortical slices under different neuromodulatory conditions. At this microscale level, we demonstrate that perturbational complexity can be effectively modulated by pharmacological reduction of bistability and, albeit to a lesser extent, by enhancement of excitability, providing mechanistic insights into the macroscale measurements performed in humans.

Neural Correlates of Anesthesia in Newborn Mice and Humans

Monitoring the hypnotic component of anesthesia during surgeries is critical to prevent intraoperative awareness and reduce adverse side effects. For this purpose, electroencephalographic (EEG) methods complementing measures of autonomic functions and behavioral responses are in use in clinical practice. However, in human neonates and infants existing methods may be unreliable and the correlation between brain activity and anesthetic depth is still poorly understood. Here, we characterized the effects of different anesthetics on brain activity in neonatal mice and developed machine learning approaches to identify electrophysiological features predicting inspired or end-tidal anesthetic concentration as a proxy for anesthetic depth. We show that similar features from EEG recordings can be applied to predict anesthetic concentration in neonatal mice and humans. These results might support a novel strategy to monitor anesthetic depth in human newborns.

A fast and general method to empirically estimate the complexity of brain responses to transcranial and intracranial stimulations

BACKGROUND

The Perturbational Complexity Index (PCI) was recently introduced to assess the capacity of thalamocortical circuits to engage in complex patterns of causal interactions. While showing high accuracy in detecting consciousness in brain-injured patients, PCI depends on elaborate experimental setups and offline processing, and has restricted applicability to other types of brain signals beyond transcranial magnetic stimulation and high-density EEG (TMS/hd-EEG) recordings.

OBJECTIVE

We aim to address these limitations by introducing PCI^ST, a fast method for estimating perturbational complexity of any given brain response signal.

METHODS

PCI^ST is based on dimensionality reduction and state transitions (ST) quantification of evoked potentials. The index was validated on a large dataset of TMS/hd-EEG recordings obtained from 108 healthy subjects and 108 brain-injured patients, and tested on sparse intracranial recordings (SEEG) of 9 patients undergoing intracranial single-pulse electrical stimulation (SPES) during wakefulness and sleep.

RESULTS

When calculated on TMS/hd-EEG potentials, PCI^ST performed with the same accuracy as the original PCI, while improving on the previous method by being computed in less than a second and requiring a simpler set-up. In SPES/SEEG signals, the index was able to quantify a systematic reduction of intracranial complexity during sleep, confirming the occurrence of state-dependent changes in the effective connectivity of thalamocortical circuits, as originally assessed through TMS/hd-EEG.

CONCLUSIONS

PCI^ST represents a fundamental advancement towards the implementation of a reliable and fast clinical tool for the bedside assessment of consciousness as well as a general measure to explore the neuronal mechanisms of loss/recovery of brain complexity across scales and models.

Sensory Processing Across Conscious and Nonconscious Brain States: From Single Neurons to Distributed Networks for Inferential Representation

Neuronal activity is markedly different across brain states: it varies from desynchronized activity during wakefulness to the synchronous alternation between active and silent states characteristic of deep sleep. Surprisingly, limited attention has been paid to investigating how brain states affect sensory processing. While it was long assumed that the brain was mostly disconnected from external stimuli during sleep, an increasing number of studies indicates that sensory stimuli continue to be processed across all brain states-albeit differently. In this review article, we first discuss what constitutes a brain state. We argue that-next to global, behavioral states such as wakefulness and sleep-there is a concomitant need to distinguish bouts of oscillatory dynamics with specific global/local activity patterns and lasting for a few hundreds of milliseconds, as these can lead to the same sensory stimulus being either perceived or not. We define these short-lasting bouts as micro-states. We proceed to characterize how sensory-evoked neural responses vary between conscious and nonconscious states. We focus on two complementary aspects: neuronal ensembles and inter-areal communication. First, we review which features of ensemble activity are conducive to perception, and how these features vary across brain states. Properties such as heterogeneity, sparsity and synchronicity in neuronal ensembles will especially be considered as essential correlates of conscious processing. Second, we discuss how inter-areal communication varies across brain states and how this may affect brain operations and sensory processing. Finally, we discuss predictive coding (PC) and the concept of multi-level representations as a key framework for understanding conscious sensory processing. In this framework the brain implements conscious representations as inferences about world states across multiple representational levels. In this representational hierarchy, low-level inference may be carried out nonconsciously, whereas high levels integrate across different sensory modalities and larger spatial scales, correlating with conscious processing. This inferential framework is used to interpret several cellular and population-level findings in the context of brain states, and we briefly compare its implications to two other theories of consciousness. In conclusion, this review article, provides foundations to guide future studies aiming to uncover the mechanisms of sensory processing and perception across brain states.

Intrinsic Functional Connectivity Resembles Cortical Architecture at Various Levels of Isoflurane Anesthesia

Cortical single neuron activity and local field potential patterns change at different depths of general anesthesia. Here, we investigate the associated network level changes of functional connectivity. We recorded ongoing electrocorticographic (ECoG) activity from temporo-parieto-occipital cortex of 6 ferrets at various levels of isoflurane/nitrous oxide anesthesia and determined functional connectivity by computing amplitude envelope correlations. Through hierarchical clustering, we derived typical connectivity patterns corresponding to light, intermediate and deep anesthesia. Generally, amplitude correlation strength increased strongly with depth of anesthesia across all cortical areas and frequency bands. This was accompanied, at the deepest level, by the emergence of burst-suppression activity in the ECoG signal and a change of the spectrum of the amplitude envelope. Normalization of functional connectivity to the distribution of correlation coefficients showed that the topographical patterns remained similar across depths of anesthesia, reflecting the functional association of the underlying cortical areas. Thus, while strength and temporal properties of amplitude co-modulation vary depending on the activity of local neural circuits, their network-level interaction pattern is presumably most strongly determined by the underlying structural connectivity.

Highly Dynamic Spatiotemporal Organization of Low-Frequency Activities During Behavioral States in the Mouse Cerebral Cortex

Although low-frequency (LF < 10 Hz) activities have been considered as a hallmark of nonrapid eye movement (NREM) sleep, several studies have recently reported LF activities in the membrane potential of cortical neurons from different areas in awake mice. However, little is known about the spatiotemporal organization of LF activities across cortical areas during wakefulness and to what extent it differs during NREM sleep. We have thus investigated the dynamics of LF activities across cortical areas in awake and sleeping mice using chronic simultaneous local field potential recordings. We found that LF activities had higher amplitude in somatosensory and motor areas during quiet wakefulness and decreased in most areas during active wakefulness, resulting in a global state change that was overall correlated with motor activity. However, we also observed transient desynchronization of cortical states between areas, indicating a more local state regulation. During NREM sleep, LF activities had higher amplitude in all areas but slow-wave activity was only poorly correlated across cortical areas. Despite a maximal amplitude during NREM sleep, the coherence of LF activities between areas that are not directly connected dropped from wakefulness to NREM sleep, potentially reflecting a breakdown of long-range cortical integration associated with loss of consciousness.

Dysfunctional Autism Risk Genes Cause Circuit-Specific Connectivity Deficits With Distinct Developmental Trajectories

Autism spectrum disorders (ASD) are a set of complex neurodevelopmental disorders for which there is currently no targeted therapeutic approach. It is thought that alterations of genes regulating migration and synapse formation during development affect neural circuit formation and result in aberrant connectivity within distinct circuits that underlie abnormal behaviors. However, it is unknown whether deviant developmental trajectories are circuit-specific for a given autism risk-gene. We used MRI to probe changes in functional and structural connectivity from childhood to adulthood in Fragile-X (Fmr1-/y) and contactin-associated (CNTNAP2-/-) knockout mice. Young Fmr1-/y mice (30 days postnatal) presented with a robust hypoconnectivity phenotype in corticocortico and corticostriatal circuits in areas associated with sensory information processing, which was maintained until adulthood. Conversely, only small differences in hippocampal and striatal areas were present during early postnatal development in CNTNAP2-/- mice, while major connectivity deficits in prefrontal and limbic pathways developed between adolescence and adulthood. These findings are supported by viral tracing and electron micrograph approaches and define 2 clearly distinct connectivity endophenotypes within the autism spectrum. We conclude that the genetic background of ASD strongly influences which circuits are most affected, the nature of the phenotype, and the developmental time course of the associated changes.

Inter-subject FDG PET Brain Networks Exhibit Multi-scale Community Structure with Different Normalization Techniques

Inter-subject networks are used to model correlations between brain regions and are particularly useful for metabolic imaging techniques, like 18F-2-deoxy-2-(18F)fluoro-D-glucose (FDG) positron emission tomography (PET). Since FDG PET typically produces a single image, correlations cannot be calculated over time. Little focus has been placed on the basic properties of inter-subject networks and if they are affected by group size and image normalization. FDG PET images were acquired from rats (n = 18), normalized by whole brain, visual cortex, or cerebellar FDG uptake, and used to construct correlation matrices. Group size effects on network stability were investigated by systematically adding rats and evaluating local network connectivity (node strength and clustering coefficient). Modularity and community structure were also evaluated in the differently normalized networks to assess meso-scale network relationships. Local network properties are stable regardless of normalization region for groups of at least 10. Whole brain-normalized networks are more modular than visual cortex- or cerebellum-normalized network (p < 0.00001); however, community structure is similar at network resolutions where modularity differs most between brain and randomized networks. Hierarchical analysis reveals consistent modules at different scales and clustering of spatially-proximate brain regions. Findings suggest inter-subject FDG PET networks are stable for reasonable group sizes and exhibit multi-scale modularity.

Emergence EEG pattern classification in sevoflurane anesthesia

OBJECTIVE

Significant spectral electroencephalogram (EEG) pattern characteristics exist in individual patients during the re-establishment of consciousness after general anesthesia. However, these EEG patterns cannot be quantitatively identified using commercially available depth of anesthesia (DoA) monitors. This study proposes an effective classification method and indices to classify these patterns among patients.

APPROACH

Four types of emergence EEG patterns were identified based on the EEG data set from 52 patients undergoing sevoflurane general anesthesia from two hospitals. Then, the relative power spectrum density (RPSD) of five frequency sub-bands of clinical interest (delta, theta, alpha, beta and gamma) were selected for emergence state analysis. Finally, a genetic algorithm support vector machine (GA-SVM) was used to identify the emergence EEG patterns. The performance was reported in terms of sensitivity (SE), specificity (SP) and accuracy (AC).

MAIN RESULTS

The combination of the mean and mode of RPSD in the delta and alpha band (P (delta)/P (alpha) performed the best in the GA-SVM classification. The AC indices obtained by GA-SVM across the four patterns were 90.64  ±  7.61, 81.79  ±  5.84, 82.14  ±  7.99 and 72.86  ±  11.11 respectively. Furthermore, the emergence time of the patients with EEG emergence patterns I and III increased as the patients' age increased. However, for patients with EEG emergence pattern IV, the emergence time positively correlates with the patients' age when they are under 50, and negatively correlates with it when they are over 50.

SIGNIFICANCE

The mean and mode of P (delta)/P (alpha) is a useful index to classify the different emergence EEG patterns. In addition, these patterns may correlate with an underlying neural substrate which is related to the patients' age. Highlights ► Four emergence EEG patterns were found in γ-amino-butyric acid (GABA)-ergic anesthetic drugs. ► A genetic algorithm combined with a support vector machine (GA-SVM) was proposed to identify the emergence EEG patterns. ► The relative power spectrum density (RPSD) was used as a feature to classify the emergence EEG patterns and good accuracy was achieved. ► The statistics shows that the emergence EEG patterns are age-related and may have value in assessing postoperative brain states.

Menthol enhances nicotine-induced locomotor sensitization and in vivo functional connectivity in adolescence

Mentholated cigarettes capture a quarter of the US market, and are disproportionately smoked by adolescents. Menthol allosterically modulates nicotinic acetylcholine receptor function, but its effects on the brain and nicotine addiction are unclear. To determine if menthol is psychoactive, we assessed locomotor sensitization and brain functional connectivity. Adolescent male Sprague Dawley rats were administered nicotine (0.4 mg/kg) daily with or without menthol (0.05 mg/kg or 5.38 mg/kg) for nine days. Following each injection, distance traveled in an open field was recorded. One day after the sensitization experiment, functional connectivity was assessed in awake animals before and after drug administration using magnetic resonance imaging. Menthol (5.38 mg/kg) augmented nicotine-induced locomotor sensitization. Functional connectivity was compared in animals that had received nicotine with or without the 5.38 mg/kg dosage of menthol. Twenty-four hours into withdrawal after the last drug administration, increased functional connectivity was observed for ventral tegmental area and retrosplenial cortex with nicotine+menthol compared to nicotine-only exposure. Upon drug re-administration, the nicotine-only, but not the menthol groups, exhibited altered functional connectivity of the dorsal striatum with the amygdala. Menthol, when administered with nicotine, showed evidence of psychoactive properties by affecting brain activity and behavior compared to nicotine administration alone.

Stereotypical modulations in dynamic functional connectivity explained by changes in BOLD variance

Spontaneous activity measured in human subject under the absence of any task exhibits complex patterns of correlation that largely correspond to large-scale functional topographies obtained with a wide variety of cognitive and perceptual tasks. These "resting state networks" (RSNs) fluctuate over time, forming and dissolving on the scale of seconds to minutes. While these fluctuations, most prominently those of the default mode network, have been linked to cognitive function, it remains unclear whether they result from random noise or whether they index a nonstationary process which could be described as state switching. In this study, we use a sliding windows-approach to relate temporal dynamics of RSNs to global modulations in correlation and BOLD variance. We compare empirical data, phase-randomized surrogate data, and data simulated with a stationary model. We find that RSN time courses exhibit a large amount of coactivation in all three cases, and that the modulations in their activity are closely linked to global dynamics of the underlying BOLD signal. We find that many properties of the observed fluctuations in FC and BOLD, including their ranges and their correlations amongst each other, are explained by fluctuations around the average FC structure. However, we also report some interesting characteristics that clearly support nonstationary features in the data. In particular, we find that the brain spends more time in the troughs of modulations than can be expected from stationary dynamics.

BOLD temporal variability differentiates wakefulness from anesthesia-induced unconsciousness

Even though a number of findings, based on information content or information integration, are shown to define neural underpinnings characteristic of a conscious experience, the neurophysiological mechanism of consciousness is still poorly understood. Here, we investigated the brain activity and functional connectivity changes that occur in the isoflurane-anesthetized unconscious state in contrast to the awake state in rats (awake and/or anesthetized, n = 68 rats). We examined nine information measures previously shown to distinguish between conscious states: blood oxygen level-dependent (BOLD) variability, functional connectivity strength, modularity, weighted modularity, efficiency, clustering coefficient, small-worldness, and spatial and temporal Lempel-Ziv complexity measure. We also identified modular membership, seed-based network connectivity, and absolute and normalized power spectrums to assess the integrity of the BOLD functional networks between awake and anesthesia. fMRI BOLD variability and related absolute power were the only information measures significantly higher during the awake state compared with isoflurane anesthesia across animals, and with varying levels of anesthesia, after correcting for motion and respiration confounds. Thus, we conclude that, at least under the specific conditions examined here, global measures of information integration/sharing do not properly distinguish the anesthetized state from wakefulness, and heightened overall, global and local, BOLD variability is the most reliable determinant of conscious brain activity relative to isoflurane anesthesia. NEW & NOTEWORTHY Multiple metrics previously suggested to be able to distinguish between states of consciousness were compared, within and across rats in awake and isoflurane anesthesia-induced unconsciousness. All measures tested showed sensitivity to confounds, correcting for motion and for respiration changes due to anesthesia. Resting state local BOLD variability and the related absolute power were the only information measures that robustly differentiated wakefulness states. These results caution against the general applicability of global information measures in identifying levels of consciousness, thus challenging the popular concept that these measures reflect states of consciousness, and also pointing to local signal variability as a more reliable indicator of states of wakefulness.

Neural and metabolic basis of dynamic resting state fMRI

Resting state fMRI (rsfMRI) as a technique showed much initial promise for use in psychiatric and neurological diseases where diagnosis and treatment were difficult. To realize this promise, many groups have moved towards examining "dynamic rsfMRI," which relies on the assumption that rsfMRI measurements on short time scales remain relevant to the underlying neural and metabolic activity. Many dynamic rsfMRI studies have demonstrated differences between clinical or behavioral groups beyond what static rsfMRI measured, suggesting a neurometabolic basis. Correlative studies combining dynamic rsfMRI and other physiological measurements have supported this. However, they also indicate multiple mechanisms and, if using correlation alone, it is difficult to separate cause and effect. Hypothesis-driven studies are needed, a few of which have begun to illuminate the underlying neurometabolic mechanisms that shape observed differences in dynamic rsfMRI. While the number of potential noise sources, potential actual neurometabolic sources, and methodological considerations can seem overwhelming, dynamic rsfMRI provides a rich opportunity in systems neuroscience. Even an incrementally better understanding of the neurometabolic basis of dynamic rsfMRI would expand rsfMRI's research and clinical utility, and the studies described herein take the first steps on that path forward.

Evidence of Altered Brain Responses to Nicotine in an Animal Model of Attention Deficit/Hyperactivity Disorder

Introduction

Individuals with attention deficit/hyperactivity disorder (ADHD) are susceptible to earlier and more severe nicotine addiction. To shed light on the relationship between nicotine and ADHD, we examined nicotine's effects on functional brain networks in an animal model of ADHD.

Methods

Awake magnetic resonance imaging was used to compare functional connectivity in adolescent (post-natal day 44 ± 2) males of the spontaneously hypertensive rat (SHR) strain and two control strains, Wistar-Kyoto and Sprague-Dawley (n = 16 each). We analyzed functional connectivity immediately before and after nicotine exposure (0.4 mg/kg base) in naïve animals, using a region-of-interest approach focussing on 16 regions previously implicated in reward and addiction.

Results

Relative to the control groups, the SHR strain demonstrated increased functional connectivity between the ventral tegmental area (VTA) and retrosplenial cortex in response to nicotine, suggesting an aberrant response to nicotine. In contrast, increased VTA-substantia nigra connectivity in response to a saline injection in the SHR was absent following a nicotine injection, suggesting that nicotine normalized function in this circuit.

Conclusions

In the SHR, nicotine triggered an atypical response in one VTA circuit while normalizing activity in another. The VTA has been widely implicated in drug reward. Our data suggest that increased susceptibility to nicotine addiction in individuals with ADHD may involve altered responses to nicotine involving VTA circuits.

Implications

Nicotine addiction is more common among individuals with ADHD. We found that two circuits involving the VTA responded differently to nicotine in animals that model ADHD in comparison to two control strains. In one circuit, nicotine normalized activity that was abnormal in the ADHD animals, while in the other circuit nicotine caused an atypical brain response in the ADHD animals. The VTA has been implicated in drug reward. Our results would be consistent with an interpretation that nicotine may normalize abnormal brain activity in ADHD, and that nicotine may be more rewarding for individuals with ADHD.

Temporal Processing in the Visual Cortex of the Awake and Anesthetized Rat

The activity pattern and temporal dynamics within and between neuron ensembles are essential features of information processing and believed to be profoundly affected by anesthesia. Much of our general understanding of sensory information processing, including computational models aimed at mathematically simulating sensory information processing, rely on parameters derived from recordings conducted on animals under anesthesia. Due to the high variety of neuronal subtypes in the brain, population-based estimates of the impact of anesthesia may conceal unit- or ensemble-specific effects of the transition between states. Using chronically implanted tetrodes into primary visual cortex (V1) of rats, we conducted extracellular recordings of single units and followed the same cell ensembles in the awake and anesthetized states. We found that the transition from wakefulness to anesthesia involves unpredictable changes in temporal response characteristics. The latency of single-unit responses to visual stimulation was delayed in anesthesia, with large individual variations between units. Pair-wise correlations between units increased under anesthesia, indicating more synchronized activity. Further, the units within an ensemble show reproducible temporal activity patterns in response to visual stimuli that is changed between states, suggesting state-dependent sequences of activity. The current dataset, with recordings from the same neural ensembles across states, is well suited for validating and testing computational network models. This can lead to testable predictions, bring a deeper understanding of the experimental findings and improve models of neural information processing. Here, we exemplify such a workflow using a Brunel network model.

Shaping the Default Activity Pattern of the Cortical Network

Slow oscillations have been suggested as the default emergent activity of the cortical network. This is a low complexity state that integrates neuronal, synaptic, and connectivity properties of the cortex. Shaped by variations of physiological parameters, slow oscillations provide information about the underlying healthy or pathological network. We review how this default activity is shaped, how it acts as a powerful attractor, and how getting out of it is necessary for the brain to recover the levels of complexity associated with conscious states. We propose that slow oscillations provide a robust unifying paradigm for the study of cortical function.

Spike-Based Functional Connectivity in Cerebral Cortex and Hippocampus: Loss of Global Connectivity Is Coupled to Preservation of Local Connectivity During Non-REM Sleep

Behavioral states are commonly considered global phenomena with homogeneous neural determinants. However, recent studies indicate that behavioral states modulate spiking activity with neuron-level specificity as a function of brain area, neuronal subtype, and preceding history. Although functional connectivity also strongly depends on behavioral state at a mesoscopic level and is globally weaker in non-REM (NREM) sleep and anesthesia than wakefulness, it is unknown how neuronal communication is modulated at the cellular level. We hypothesize that, as for neuronal activity, the influence of behavioral states on neuronal coupling strongly depends on type, location, and preceding history of involved neurons. Here, we applied nonlinear, information-theoretical measures of functional connectivity to ensemble recordings with single-cell resolution to quantify neuronal communication in the neocortex and hippocampus of rats during wakefulness and sleep. Although functional connectivity (measured in terms of coordination between firing rate fluctuations) was globally stronger in wakefulness than in NREM sleep (with distinct traits for cortical and hippocampal areas), the drop observed during NREM sleep was mainly determined by a loss of inter-areal connectivity between excitatory neurons. Conversely, local (intra-area) connectivity and long-range (inter-areal) coupling between interneurons were preserved during NREM sleep. Furthermore, neuronal networks that were either modulated or not by a behavioral task remained segregated during quiet wakefulness and NREM sleep. These results show that the drop in functional connectivity during wake-sleep transitions globally holds true at the cellular level, but confine this change mainly to long-range coupling between excitatory neurons.

SIGNIFICANCE STATEMENT

Studies performed at a mesoscopic level of analysis have shown that communication between cortical areas is disrupted in non-REM sleep and anesthesia. However, the neuronal determinants of this phenomenon are not known. Here, we applied nonlinear, information-theoretical measures of functional coupling to multi-area tetrode recordings from freely moving rats to investigate whether and how brain state modulates coordination between individual neurons. We found that the previously observed drop in functional connectivity during non-REM (NREM) sleep can be explained by a decrease in coupling between excitatory neurons located in distinct brain areas. Conversely, intra-area communication and coupling between interneurons are preserved. Our results provide significant new insights into the neuron-level mechanisms responsible for the loss of consciousness occurring in NREM sleep.

How structure sculpts function: Unveiling the contribution of anatomical connectivity to the brain's spontaneous correlation structure

Intrinsic brain activity is characterized by highly organized co-activations between different regions, forming clustered spatial patterns referred to as resting-state networks. The observed co-activation patterns are sustained by the intricate fabric of millions of interconnected neurons constituting the brain's wiring diagram. However, as for other real networks, the relationship between the connectional structure and the emergent collective dynamics still evades complete understanding. Here, we show that it is possible to estimate the expected pair-wise correlations that a network tends to generate thanks to the underlying path structure. We start from the assumption that in order for two nodes to exhibit correlated activity, they must be exposed to similar input patterns from the entire network. We then acknowledge that information rarely spreads only along a unique route but rather travels along all possible paths. In real networks, the strength of local perturbations tends to decay as they propagate away from the sources, leading to a progressive attenuation of the original information content and, thus, of their influence. Accordingly, we define a novel graph measure, topological similarity, which quantifies the propensity of two nodes to dynamically correlate as a function of the resemblance of the overall influences they are expected to receive due to the underlying structure of the network. Applied to the human brain, we find that the similarity of whole-network inputs, estimated from the topology of the anatomical connectome, plays an important role in sculpting the backbone pattern of time-average correlations observed at rest.

Structural connectome topology relates to regional BOLD signal dynamics in the mouse brain

Brain dynamics are thought to unfold on a network determined by the pattern of axonal connections linking pairs of neuronal elements; the so-called connectome. Prior work has indicated that structural brain connectivity constrains pairwise correlations of brain dynamics ("functional connectivity"), but it is not known whether inter-regional axonal connectivity is related to the intrinsic dynamics of individual brain areas. Here we investigate this relationship using a weighted, directed mesoscale mouse connectome from the Allen Mouse Brain Connectivity Atlas and resting state functional MRI (rs-fMRI) time-series data measured in 184 brain regions in eighteen anesthetized mice. For each brain region, we measured degree, betweenness, and clustering coefficient from weighted and unweighted, and directed and undirected versions of the connectome. We then characterized the univariate rs-fMRI dynamics in each brain region by computing 6930 time-series properties using the time-series analysis toolbox, hctsa. After correcting for regional volume variations, strong and robust correlations between structural connectivity properties and rs-fMRI dynamics were found only when edge weights were accounted for, and were associated with variations in the autocorrelation properties of the rs-fMRI signal. The strongest relationships were found for weighted in-degree, which was positively correlated to the autocorrelation of fMRI time series at time lag τ = 34 s (partial Spearman correlation ρ=0.58), as well as a range of related measures such as relative high frequency power (f > 0.4 Hz: ρ=-0.43). Our results indicate that the topology of inter-regional axonal connections of the mouse brain is closely related to intrinsic, spontaneous dynamics such that regions with a greater aggregate strength of incoming projections display longer timescales of activity fluctuations.

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.

N-methyl-D-aspartate receptor encephalitis mediates loss of intrinsic activity measured by functional MRI

Spontaneous brain activity is required for the development and maintenance of normal brain function. Many disease processes disrupt the organization of intrinsic brain activity, but few pervasively reduce the amplitude of resting state blood oxygen level dependent (BOLD) fMRI fluctuations. We report the case of a female with anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis, longitudinally studied during the course of her illness to determine the contribution of NMDAR signaling to spontaneous brain activity. Resting state BOLD fMRI was measured at the height of her illness and 18 weeks following discharge from hospital. Conventional resting state networks were defined using established methods. Correlation and covariance matrices were calculated by extracting the BOLD time series from regions of interest and calculating either the correlation or covariance quantity. The intrinsic activity was compared between visits, and to expected activity from 45 similarly aged healthy individuals. Near the height of the illness, the patient exhibited profound loss of consciousness, high-amplitude slowing of the electroencephalogram, and a severe reduction in the amplitude of spontaneous BOLD fMRI fluctuations. The patient's neurological status and measures of intrinsic activity improved following treatment. We conclude that NMDAR-mediated signaling plays a critical role in the mechanisms that give rise to organized spontaneous brain activity. Loss of intrinsic activity is associated with profound disruptions of consciousness and cognition.

Large-scale functional connectivity networks in the rodent brain

Resting-state functional Magnetic Resonance Imaging (rsfMRI) of the human brain has revealed multiple large-scale neural networks within a hierarchical and complex structure of coordinated functional activity. These distributed neuroanatomical systems provide a sensitive window on brain function and its disruption in a variety of neuropathological conditions. The study of macroscale intrinsic connectivity networks in preclinical species, where genetic and environmental conditions can be controlled and manipulated with high specificity, offers the opportunity to elucidate the biological determinants of these alterations. While rsfMRI methods are now widely used in human connectivity research, these approaches have only relatively recently been back-translated into laboratory animals. Here we review recent progress in the study of functional connectivity in rodent species, emphasising the ability of this approach to resolve large-scale brain networks that recapitulate neuroanatomical features of known functional systems in the human brain. These include, but are not limited to, a distributed set of regions identified in rats and mice that may represent a putative evolutionary precursor of the human default mode network (DMN). The impact and control of potential experimental and methodological confounds are also critically discussed. Finally, we highlight the enormous potential and some initial application of connectivity mapping in transgenic models as a tool to investigate the neuropathological underpinnings of the large-scale connectional alterations associated with human neuropsychiatric and neurological conditions. We conclude by discussing the translational potential of these methods in basic and applied neuroscience.

The power of using functional fMRI on small rodents to study brain pharmacology and disease

Functional magnetic resonance imaging (fMRI) is an excellent tool to study the effect of pharmacological modulations on brain function in a non-invasive and longitudinal manner. We introduce several blood oxygenation level dependent (BOLD) fMRI techniques, including resting state (rsfMRI), stimulus-evoked (st-fMRI), and pharmacological MRI (phMRI). Respectively, these techniques permit the assessment of functional connectivity during rest as well as brain activation triggered by sensory stimulation and/or a pharmacological challenge. The first part of this review describes the physiological basis of BOLD fMRI and the hemodynamic response on which the MRI contrast is based. Specific emphasis goes to possible effects of anesthesia and the animal’s physiological conditions on neural activity and the hemodynamic response. The second part of this review describes applications of the aforementioned techniques in pharmacologically induced, as well as in traumatic and transgenic disease models and illustrates how multiple fMRI methods can be applied successfully to evaluate different aspects of a specific disorder. For example, fMRI techniques can be used to pinpoint the neural substrate of a disease beyond previously defined hypothesis-driven regions-of-interest. In addition, fMRI techniques allow one to dissect how specific modifications (e.g., treatment, lesion etc.) modulate the functioning of specific brain areas (st-fMRI, phMRI) and how functional connectivity (rsfMRI) between several brain regions is affected, both in acute and extended time frames. Furthermore, fMRI techniques can be used to assess/explore the efficacy of novel treatments in depth, both in fundamental research as well as in preclinical settings. In conclusion, by describing several exemplary studies, we aim to highlight the advantages of functional MRI in exploring the acute and long-term effects of pharmacological substances and/or pathology on brain functioning along with several methodological considerations.