The balanced and introspective brain

Method for quantifying arousal and consciousness in healthy states and severe brain injury via EEG-based measures of corticothalamic physiology

BACKGROUND

Characterization of normal arousal states has been achieved by fitting predictions of corticothalamic neural field theory (NFT) to electroencephalographic (EEG) spectra to yield relevant physiological parameters.

NEW METHOD

A prior fitting method is extended to distinguish conscious and unconscious states in healthy and brain injured subjects by identifying additional parameters and clusters in parameter space.

RESULTS

Fits of NFT predictions to EEG spectra are used to estimate neurophysiological parameters in healthy and brain injured subjects. Spectra are used from healthy subjects in wake and sleep and from patients with unresponsive wakefulness syndrome, in a minimally conscious state (MCS), and emerged from MCS. Subjects cluster into three groups in parameter space: conscious healthy (wake and REM), sleep, and brain injured. These are distinguished by the difference X-Y between corticocortical (X) and corticothalamic (Y) feedbacks, and by mean neural response rates α and β to incoming spikes. X-Y tracks consciousness in healthy individuals, with smaller values in wake/REM than sleep, but cannot distinguish between brain injuries. Parameters α and β differentiate deep sleep from wake/REM and brain injury.

COMPARISON WITH EXISTING METHODS

Other methods typically rely on laborious clinical assessment, manual EEG scoring, or evaluation of measures like Φ from integrated information theory, for which no efficient method exists. In contrast, the present method can be automated on a personal computer.

CONCLUSION

The method provides a means to quantify consciousness and arousal in healthy and brain injured subjects, but does not distinguish subtypes of brain injury.

Discrete spectral eigenmode-resonance network of brain dynamics and connectivity

The problem of finding a compact natural representation of brain dynamics and connectivity is addressed using an expansion in terms of physical spatial eigenmodes and their frequency resonances. It is demonstrated that this discrete expansion via the system transfer function enables linear and nonlinear dynamics to be analyzed in compact form in terms of natural dynamic "atoms," each of which is a frequency resonance of an eigenmode. Because these modal resonances are determined by the system dynamics, not the investigator, they are privileged over widely used phenomenological patterns, and obviate the need for artificial discretizations and thresholding in coordinate space. It is shown that modal resonances participate as nodes of a discrete spectral network, are noninteracting in the linear regime, but are linked nonlinearly by wave-wave coalescence and decay processes. The modal resonance formulation is shown to be capable of speeding numerical calculations of strongly nonlinear interactions. Recent work in brain dynamics, especially based on neural field theory (NFT) approaches, allows eigenmodes and their resonances to be estimated from data without assuming a specific brain model. This means that dynamic equations can be inferred using system identification methods from control theory, rather than being assumed, and resonances can be interpreted as control-systems data filters. The results link brain activity and connectivity with control-systems functions such as prediction and attention via gain control and can also be linked to specific NFT predictions if desired, thereby providing a convenient bridge between physiologically based theories and experiment. Amplitudes of modes and resonances can also be tracked to provide a more direct and temporally localized representation of the dynamics than correlations and covariances, which are widely used in the field. By synthesizing many different lines of research, this work provides a way to link quantitative electrophysiological and imaging measurements, connectivity, brain dynamics, and function. This underlines the need to move between coordinate and spectral representations as required. Moreover, standard theoretical-physics approaches and mathematical methods can be used in place of ad hoc statistical measures such as those based on graph theory of artificially discretized and decimated networks, which are highly prone to selection effects and artifacts.

Neural field theory of neural avalanche exponents

The power-law exponents of observed size and lifetime distributions of near-critical neural avalanches are calculated from neural field theory using diagrammatic methods. This brings neural avalanches within the ambit of neural field theory, which has also previously explained near-critical 1/f spectra and many other observed features of neural activity. This strengthens the case for near-criticality of the brain and opens the way for these other phenomena to be interrelated with avalanches and their dynamics.

Relationships between lognormal distributions of neural properties, activity, criticality, and connectivity

Relationships between convergence of inputs onto neurons, divergence of outputs from them, synaptic strengths, nonlinear firing response properties, and randomness of axonal ranges are systematically explored by interrelating means and variances of synaptic strengths, firing rates, and soma voltages. When self-consistency is imposed, it is found that broad distributions of synaptic strength are a necessary concomitant of the known massive convergence of inputs to individual neurons, and observed widths of lognormal distributions of synaptic strength and firing rate are explained provided the brain is in a near-critical state, consistent with independent observations. The strongest individual synapses are shown to have an effect on soma voltage comparable to the effect of all others combined, which supports suggestions that they may have a key role in neural communication. Remarkably, inclusion of moderate randomness in characteristic axonal ranges is shown to account for the observed [Formula: see text]-fold variability in two-point connectivity at a given separation and [Formula: see text]-fold overall when the known mean exponential fall-off is included, consistent with observed near-lognormal distributions. Inferred axonal deviations from straight-line paths are also consistent with independent estimates.

Importance of self-connections for brain connectivity and spectral connectomics

Spectral analysis and neural field theory are used to investigate the role of local connections in brain connectivity matrices (CMs) that quantify connectivity between pairs of discretized brain regions. This work investigates how the common procedure of omitting such self-connections (i.e., the diagonal elements of CMs) in published studies of brain connectivity affects the properties of functional CMs (fCMs) and the mutually consistent effective CMs (eCMs) that correspond to them. It is shown that retention of self-connections in the fCM calculated from two-point activity covariances is essential for the fCM to be a true covariance matrix, to enable correct inference of the direct total eCMs from the fCM, and to ensure their compatibility with it; the deCM and teCM represent the strengths of direct connections and all connections between points, respectively. When self-connections are retained, inferred eCMs are found to have net inhibitory self-connections that represent the local inhibition needed to balance excitation via white matter fibers at longer ranges. This inference of spatially unresolved connectivity exemplifies the power of spectral connectivity methods, which also enable transformation of CMs to compact diagonal forms that allow accurate approximation of the fCM and total eCM in terms of just a few modes, rather than the full [Formula: see text] CM entries for connections between N brain regions. It is found that omission of fCM self-connections affects both local and long-range connections in eCMs, so they cannot be omitted even when studying the large-scale. Moreover, retention of local connections enables inference of subgrid short-range inhibitory connectivity. The results are verified and illustrated using the NKI-Rockland dataset from the University of Southern California Multimodal Connectivity Database. Deletion of self-connections is common in the field; this does not affect case-control studies but the present results imply that such fCMs must have self-connections restored before eCMs can be inferred from them.

Neural Field Theory of Evoked Response Potentials With Attentional Gain Dynamics

A generalized neural field model of large-scale activity in the corticothalamic system is used to predict standard evoked potentials. This model embodies local feedbacks that modulate the gains of neural activity as part of the response to incoming stimuli and thus enables both activity changes and effective connectivity changes to be calculated as parts of a generalized evoked response, and their relative contributions to be determined. The results show that incorporation of gain modulations enables a compact and physically justifiable description of the differences in gain between background-EEG and standard-ERP conditions, with the latter able to be initiated from the background state, rather than requiring distinct parameters as in earlier work. In particular, top-down gains are found to be reduced during an ERP, consistent with recent theoretical suggestions that the role of internal models is diminished in favor of external inputs when the latter change suddenly. The static-gain and modulated-gain system transfer functions are analyzed via control theory in terms of system resonances that were recently shown to implement data filtering whose gain adjustments can be interpreted as attention. These filters are shown to govern early and late features in standard evoked responses and their gain parameters are shown to be dynamically adjusted in a way that implements a form of attention. The results show that dynamically modulated resonant filters responsible for the low-frequency oscillations in an evoked potential response have different parameters than those responsible for low-frequency resting EEG responses, while both responses share similar mid- and high-frequency resonant filters. These results provide a biophysical mechanism by which oscillatory activity in the theta, alpha, and beta frequency ranges of an evoked response are modulated as reflections of attention; notably theta is enhanced and alpha suppressed during the latter parts of the ERP. Furthermore, the model enables the part of the ERP response induced by gain modulations to be estimated and interpreted in terms of attention.

Dynamical phase transition in spike neuronal firing patterns of hippocampal cells

There is increasing evidence that the brain resides in a state of criticality. The purpose of the present work is to characterize the dynamics of individual hippocampal CA1 pyramidal cells and to investigate how it is influenced by changes in Kv7.2/7.3 (M-channel) ion channel modulation, which is known to be key in determining the neuronal excitability. We show that the resting activity of CA1 neurons exhibit random dynamics with low information content, while changes in M-channel modulation move the neuronal activity near a phase transition to richer non-trivial dynamics. We interpret these results as the basis upon which the state of self-organized criticality is built.

Neural field theory of effects of brain modifications and lesions on functional connectivity: Acute effects, short-term homeostasis, and long-term plasticity

Neural field theory is used to predict the functional connectivity effects of lesions or other modifications to effective connectivity. Widespread initial changes are predicted after localized or diffuse changes to white or gray matter, consistent with observations, and enabling lesion severity indexes to be defined. It is shown how short-term homeostasis and longer-term plasticity can reduce perturbations while maintaining brain criticality under conditions where some connections remain fixed because of damage in the lesion core. The extent to which such effects can compensate for initial connectivity changes is then explored, showing that the strongest corrective changes are concentrated toward the edges of the perturbation if it is localized and its core is fixed. The results are applicable to inferring underlying connectivity changes and to interpreting and monitoring functional connectivity modifications after lesions, injury, surgery, drugs, or brain stimulation.

Physical brain connectomics

Brain connectivity and structure-function relationships are analyzed from a physical perspective in place of common graph-theoretic and statistical approaches that overwhelmingly ignore the brain's physical structure and geometry. Field theory is used to define connectivity tensors in terms of bare and dressed propagators, and discretized representations are implemented that respect the physical nature and dimensionality of the quantities involved, retain the correct continuum limit, and enable diagrammatic analysis. Eigenfunction analysis is used to simultaneously characterize and probe patterns of brain connectivity and activity, in place of statistical or phenomenological patterns. Physically based measures that characterize the connectivity are then developed in coordinate and spectral domains; some of which generalize or rectify graph-theoretic measures to implement correct dimensionality and continuum limits, and some replace graph-theoretic quantities. Traditional graph-based measures are shown to be highly prone to artifacts introduced by discretization and threshold, often because essential physical constraints have not been imposed, dimensionality has not been included, and/or distinctions between scalar, vector, and tensor quantities have not been considered. The results can replace them in ways that converge correctly and measure properties of brain structure, rather than of its discretization, and thus potentially enable physical interpretation of the many phenomenological results in the literature. Geometric effects are shown to dominate in determining many brain properties and care must be taken not to interpret geometric differences as differences in intrinsic neural connectivity. The results demonstrate the need to use systematic physical methods to analyze the brain and the potential of such methods to obtain new insights from data, make new predictions for experimental test, and go beyond phenomenological classification to dynamics and mechanisms.

Physiology-based ERPs in normal and abnormal states

Evoked response potentials (ERPs) and other transients are modeled as impulse responses using physiology-based neural field theory (NFT) of the corticothalamic system of neural activity in the human brain that incorporates synaptic and dendritic dynamics, firing response, axonal propagation, and corticocortical and corticothalamic pathways. The properties of model-predicted ERPs are explored throughout the stability zone of the corticothalamic system, and predicted time series and wavelet spectra are also analyzed. This provides a unified treatment of predicted ERPs for both normal and abnormal states within the brain's stability zone, including likely parameters to represent abnormal states of reduced arousal.

A New Analysis of Resting State Connectivity and Graph Theory Reveals Distinctive Short-Term Modulations due to Whisker Stimulation in Rats

Resting state (RS) connectivity has been increasingly studied in healthy and diseased brains in humans and animals. This paper presents a new method to analyze RS data from fMRI that combines multiple seed correlation analysis with graph-theory (MSRA). We characterize and evaluate this new method in relation to two other graph-theoretical methods and ICA. The graph-theoretical methods calculate cross-correlations of regional average time-courses, one using seed regions of the same size (SRCC) and the other using whole brain structure regions (RCCA). We evaluated the reproducibility, power, and capacity of these methods to characterize short-term RS modulation to unilateral physiological whisker stimulation in rats. Graph-theoretical networks found with the MSRA approach were highly reproducible, and their communities showed large overlaps with ICA components. Additionally, MSRA was the only one of all tested methods that had the power to detect significant RS modulations induced by whisker stimulation that are controlled by family-wise error rate (FWE). Compared to the reduced resting state network connectivity during task performance, these modulations implied decreased connectivity strength in the bilateral sensorimotor and entorhinal cortex. Additionally, the contralateral ventromedial thalamus (part of the barrel field related lemniscal pathway) and the hypothalamus showed reduced connectivity. Enhanced connectivity was observed in the amygdala, especially the contralateral basolateral amygdala (involved in emotional learning processes). In conclusion, MSRA is a powerful analytical approach that can reliably detect tiny modulations of RS connectivity. It shows a great promise as a method for studying RS dynamics in healthy and pathological conditions.

Theory of corticothalamic brain activity in a spherical geometry: Spectra, coherence, and correlation

Corticothalamic neural field theory is applied to a spherical geometry to better model neural activity in the human brain and is also compared with planar approximations. The frequency power spectrum, correlation, and coherence functions are computed analytically and numerically. The effects of cortical boundary conditions and resulting modal aspects of spherical corticothalamic dynamics are explored, showing that the results of spherical and finite planar geometries converge to those for the infinite planar geometry in the limit of large brain size. Estimates are made of the point at which modal series can be truncated and it is found that for physiologically plausible parameters only the lowest few spatial eigenmodes are needed for an accurate representation of macroscopic brain activity. A difference between the geometries is that there is a low-frequency 1/f spectrum in the infinite planar geometry, whereas in the spherical geometry it is 1/f^{2}. Another difference is that the alpha peak in the spherical geometry is sharper and stronger than in the planar geometry. Cortical modal effects can lead to a double alpha peak structure in the power spectrum, although the main determinant of the alpha peak is corticothalamic feedback. In the spherical geometry, the cross spectrum between two points is found to only depend on their relative distance apart. At small spatial separations the low-frequency cross spectrum is stronger than for an infinite planar geometry and the alpha peak is sharper and stronger due to the partitioning of the energy into discrete modes. In the spherical geometry, the coherence function between points decays monotonically as their separation increases at a fixed frequency, but persists further at resonant frequencies. The correlation between two points is found to be positive, regardless of the time lag and spatial separation, but decays monotonically as the separation increases at fixed time lag. At fixed distance the correlation has peaks at multiples of the period of the dominant frequency of system activity.