Explore the vast range of titles available.

To investigate the impact of experimental interventions on living biological tissue, ex vivo rodent brain slices are often used as a more controllable alternative to a live animal model. However, for meaningful results, the biological sample must be known to be healthy and viable. One of the gold-standard approaches to identifying tissue viability status is to measure the rate of tissue oxygen consumption under specific controlled conditions. Here, we work with thin (400 μm) slices of mouse cortical brain tissue which are sustained by a steady flow of oxygenated artificial cerebralspinal fluid (aCSF) at room temperature. To quantify tissue oxygen consumption (Q), we measure oxygen partial pressure (pO2) as a function of probe depth. The curvature of the obtained parabolic (or parabola-like) pO2 profiles can be used to extract Q, providing one knows the Krogh coefficient Kt, for the tissue. The oxygen trends are well described by a Fick’s law diffusion–consumption model developed by Ivanova and Simeonov, and expressed in terms of ratio (Q/K), being the rate of oxygen consumption in tissue divided by the Krogh coefficient (oxygen diffusivity × oxygen solubility) for tissue. If the fluid immediately adjacent to the tissue can be assumed to be stationary (i.e., nonflowing), one may invoke conservation of oxygen flux K·(∂P/∂x) across the interface to deduce (Kt/Kf), the ratio of Krogh coefficients for tissue and fluid. Using published interpolation formulas for the effect of salt content and temperature on oxygen diffusivity and solubility for pure water, we estimate Kf, the Krogh coefficient for aCSF, and hence deduce the Kt coefficient for tissue. We distinguish experimental uncertainty from natural biological variability by using pairs of repeated profiles at the same tissue location. We report a dimensionless Krogh ratio (Kt/Kf)=0.562±0.088 (mean ± SD), corresponding to a Krogh coefficient Kt=(1.29±0.21)×10−14 mol/(m·s·Pa) for mouse cortical tissue at room temperature, but acknowledge the experimental limitation of being unable to verify that the fluid boundary layer is truly stationary. We compare our results with those reported in the literature, and comment on the challenges and ambiguities caused by the extensive use of ‘biologically convenient’ non-SI units for tissue Krogh coefficient.

Seasonal phase disturbances in satellite Interferometric Synthetic Aperture Radar (InSAR) measurements have been reported in other studies to suggest sub-centimetre land surface terrain motion. These have been interpreted in various ways because they correlate with multiple other (sub-)seasonal signatures of, e.g., clay swelling/shrinkage and groundwater level. Recent microwave radar studies mention the occurrence of phase disturbances in different soil types and soil moisture. This study further explored this topic by modeling phase disturbances caused by both soil and vegetation surface characteristics and aimed to interpret what their possible effects on InSAR-interpreted terrain motion is. Our models, based on fundamental microwave reflection and transmission theory, found phase disturbances caused by seasonal variation of soil and vegetation that have the same magnitude as interpreted seasonal land movement in earlier InSAR studies. We showed that small, temporal differences in soil moisture and vegetation can lead to relatively large phase disturbances in InSAR measurements. These disturbances are a result of waves having to comply with boundary conditions at the interface between media with different dielectric properties. The findings of this study explain the seasonal variations found in other InSAR studies and will therefore bring new insights and alternative explanations to help improve interpretation of InSAR-derived seasonal terrain motion.

Growth of critical fluctuations prior to catastrophic state transition is generally regarded as a universal phenomenon, providing a valuable early warning signal in dynamical systems. Using an ecological fisheries model of three populations (juvenile prey J, adult prey A and predator P), a recent study has reported silent early warning signals obtained from P and A populations prior to saddle-node (SN) bifurcation, and thus concluded that early warning signals are not universal. By performing a full eigenvalue analysis of the same system we demonstrate that while J and P populations undergo SN bifurcation, A does not jump to a new state, so it is not expected to carry early warning signs. In contrast with the previous study, we capture a significant increase in the noise-induced fluctuations in the P population, but only on close approach to the bifurcation point; it is not clear why the P variance initially shows a decaying trend. Here we resolve this puzzle using observability measures from control theory. By computing the observability coefficient for the system from the recordings of each population considered one at a time, we are able to quantify their ability to describe changing internal dynamics. We demonstrate that precursor fluctuations are best observed using only the J variable, and also P variable if close to transition. Using observability analysis we are able to describe why a poorly observable variable (P) has poor forecasting capabilities although a full eigenvalue analysis shows that this variable undergoes a bifurcation. We conclude that observability analysis provides complementary information to identify the variables carrying early-warning signs about impending state transition.

Electrical recordings of brain activity during the transition from wake to anesthetic coma show temporal and spectral alterations that are correlated with gross changes in the underlying brain state. Entry into anesthetic unconsciousness is signposted by the emergence of large, slow oscillations of electrical activity (≲1Hz ) similar to the slow waves observed in natural sleep. Here we present a two-dimensional mean-field model of the cortex in which slow spatiotemporal oscillations arise spontaneously through a Turing (spatial) symmetry-breaking bifurcation that is modulated by a Hopf (temporal) instability. In our model, populations of neurons are densely interlinked by chemical synapses, and by interneuronal gap junctions represented as an inhibitory diffusive coupling. To demonstrate cortical behavior over a wide range of distinct brain states, we explore model dynamics in the vicinity of a general-anesthetic-induced transition from “wake” to “coma.” In this region, the system is poised at a codimension-2 point where competing Turing and Hopf instabilities coexist. We model anesthesia as a moderate reduction in inhibitory diffusion, paired with an increase in inhibitory postsynaptic response, producing a coma state that is characterized by emergent low-frequency oscillations whose dynamics is chaotic in time and space. The effect of long-range axonal white-matter connectivity is probed with the inclusion of a single idealized point-to-point connection. We find that the additional excitation from the long-range connection can provoke seizurelike bursts of cortical activity when inhibitory diffusion is weak, but has little impact on an active cortex. Our proposed dynamic mechanism for the origin of anesthetic slow waves complements—and contrasts with—conventional explanations that require cyclic modulation of ion-channel conductances. We postulate that a similar bifurcation mechanism might underpin the slow waves of natural sleep and comment on the possible consequences of chaotic dynamics for memory processing and learning.

[Display omitted] We present a detailed analysis of the Hindriks and van Putten thalamocortical mean-field model for propofol anesthesia [NeuroImage 60(23), 2012]. The Hindriks and van Putten (HvP) model predicts increases in delta and alpha power for moderate (up to 130%) prolongation of GABAA inhibitory response, corresponding to light anesthetic sedation. Our analysis reveals that, for deeper anesthetic effect, the model exhibits an unexpected abrupt jump in cortical activity from a low-firing state to an extremely high-firing stable state (∼250 spikes/s), and remains locked there even at GABAA prolongations as high as 300% which would be expected to induce full comatose suppression of all firing activity. We demonstrate that this unphysiological behavior can be completely suppressed with appropriate tuning of the parameters controlling the sigmoidal functions that map soma voltage to firing rate for the excitatory and inhibitory neural populations, coupled with elimination of the putative population-dependent anesthetic efficacies introduced in the HvP model. The modifications reported here constrain the anesthetized brain activity into a biologically plausible range in which the cortex now has access to a moderate-firing state (“awake”) and a low-firing (“anesthetized”) state such that the brain can transition from “awake” to “anesthetized” states at a critical level of drug concentration. The modified HvP model predicts a drug–effect hysteresis in which the drug concentration required for induction is larger than that at emergence. In addition, the revised model shows a decrease in the intensity and frequency of alpha-band fluctuations, transitioning to delta-band dominance, with deepening anesthesia. These predicted drug concentration-dependent changes in EEG dynamics are consistent with clinical reports.

Background Investigation of the nonlinear pattern dynamics of a reaction-diffusion system almost always requires numerical solution of the system’s set of defining differential equations. Traditionally, this would be done by selecting an appropriate differential equation solver from a library of such solvers, then writing computer codes (in a programming language such as C or Matlab ) to access the selected solver and display the integrated results as a function of space and time. This “code-based” approach is flexible and powerful, but requires a certain level of programming sophistication. A modern alternative is to use a graphical programming interface such as Simulink to construct a data-flow diagram by assembling and linking appropriate code blocks drawn from a library. The result is a visual representation of the inter-relationships between the state variables whose output can be made completely equivalent to the code-based solution. Results As a tutorial introduction, we first demonstrate application of the Simulink data-flow technique to the classical van der Pol nonlinear oscillator, and compare Matlab and Simulink coding approaches to solving the van der Pol ordinary differential equations. We then show how to introduce space (in one and two dimensions) by solving numerically the partial differential equations for two different reaction-diffusion systems: the well-known Brusselator chemical reactor, and a continuum model for a two-dimensional sheet of human cortex whose neurons are linked by both chemical and electrical (diffusive) synapses. We compare the relative performances of the Matlab and Simulink implementations. Conclusions The pattern simulations by Simulink are in good agreement with theoretical predictions. Compared with traditional coding approaches, the Simulink block-diagram paradigm reduces the time and programming burden required to implement a solution for reaction-diffusion systems of equations. Construction of the block-diagram does not require high-level programming skills, and the graphical interface lends itself to easy modification and use by non-experts.

During slow-wave sleep, general anesthesia, and generalized seizures, there is an absence of consciousness. These states are characterized by low-frequency large-amplitude traveling waves in scalp electroencephalogram. Therefore the oscillatory state might be an indication of failure to form coherent neuronal assemblies necessary for consciousness. A generalized seizure event is a pathological brain state that is the clearest manifestation of waves of synchronized neuronal activity. Since gap junctions provide a direct electrical connection between adjoining neurons, thus enhancing synchronous behavior, reducing gap-junction conductance should suppress seizures; however there is no clear experimental evidence for this. Here we report theoretical predictions for a physiologically-based cortical model that describes the general anesthetic phase transition from consciousness to coma, and includes both chemical synaptic and direct electrotonic synapses. The model dynamics exhibits both Hopf (temporal) and Turing (spatial) instabilities; the Hopf instability corresponds to the slow (≲8 Hz) oscillatory states similar to those seen in slow-wave sleep, general anesthesia, and seizures. We argue that a delicately balanced interplay between Hopf and Turing modes provides a canonical mechanism for the default non-cognitive rest state of the brain. We show that the Turing mode, set by gap-junction diffusion, is generally protective against entering oscillatory modes; and that weakening the Turing mode by reducing gap conduction can release an uncontrolled Hopf oscillation and hence an increased propensity for seizure and simultaneously an increased sensitivity to GABAergic anesthesia.

We argue that spatial patterns of cortical activation observed with EEG, MEG and fMRI might arise from spontaneous self-organisation of interacting populations of excitatory and inhibitory neurons. We examine the dynamical behavior of a mean-field cortical model that includes chemical and electrical (gap-junction) synapses, focusing on two limiting cases: the “slow-soma” limit with slow voltage feedback from soma to dendrite, and the “fast-soma” limit in which the feedback action of soma voltage onto dendrite reversal potentials is instantaneous. For slow soma-dendrite feedback, we find a low-frequency (∼1 Hz) dynamic Hopf instability, and a stationary Turing instability that catalyzes formation of patterned distributions of cortical firing-rate activity with pattern wavelength ∼2 cm. Turing instability can only be triggered when gap-junction diffusion between inhibitory neurons is strong, but patterning is destroyed if the tonic level of subcortical excitation is raised sufficiently. Interaction between the Hopf and Turing instabilities may describe the non-cognitive background or “default” state of the brain, as observed by BOLD imaging. In the fast-soma limit, the model predicts a high-frequency Hopf (∼35 Hz) instability, and a traveling-wave gamma-band instability that manifests as a 2-D standing-wave pattern oscillating in place at ∼30 Hz. Small levels of inhibitory diffusion enhance and broaden the definition of the gamma antinodal regions by suppressing higher-frequency spatial modes, but gamma emergence is not contingent on the presence of inhibitory gap junctions; higher levels of diffusion suppress gamma activity. Fast-soma instabilities are enhanced by increased subcortical stimulation. Prompt soma-dendrite feedback may be an essential component of the genesis and large-scale cortical synchrony of gamma activity observed at the point of cognition.

Electrical recordings of brain activity show that entry into anesthetic unconsciousness is signposted by the emergence of large, slow oscillations of electrical activity (~1 Hz) that appear very similar to the slow waves observed in natural sleep. Anesthetic effect is modeled as a moderate reduction in inhibitory diffusion, paired with an increase in inhibitory postsynaptic potential. SEE PDF] Conclusion A spontaneous, spatiotemporally chaotic state—generated by nonlinear Turing–Hopf interaction—is the underlying mechanism for the slow oscillation observed in general anesthesia.