Simulating global properties of electroencephalograms with minimal random neural networks

The human electroencephalogram (EEG) is globally characterized by a 1/f power spectrum superimposed with certain peaks, whereby the "alpha peak" in a frequency range of 8-14 Hz is the most prominent one for relaxed states of wakefulness. We present simulations of a minimal dynamical network model of leaky integrator neurons attached to the nodes of an evolving directed and weighted random graph (an Erdos-Renyi graph). We derive a model of the dendritic field potential (DFP) for the neurons leading to a simulated EEG that describes the global activity of the network. Depending on the network size, we find an oscillatory transition of the simulated EEG when the network reaches a critical connectivity. This transition, indicated by a suitably defined order parameter, is reflected by a sudden change of the network's topology when super-cycles are formed from merging isolated loops. After the oscillatory transition, the power spectra of simulated EEG time series exhibit a 1/f continuum superimposed with certain peaks. (c) 2007 Elsevier B.V. All rights reserved.

Cannabidivarin is anticonvulsant in mouse and rat in vitro and in seizure models

Summary Background and purpose: Phytocannabinoids in Cannabis sativa have diverse pharmacological targets extending beyond cannabinoid receptors and several exert notable anticonvulsant effects. For the first time, we investigated the anticonvulsant profile of the phytocannabinoid cannabidivarin (CBDV) in vitro and in in vivo seizure models. Experimental approach: The effect of CBDV (1-100μM) on epileptiform local field potentials (LFPs) induced in rat hippocampal brain slices by 4-AP application or Mg2+-free conditions was assessed by in vitro multi-electrode array recordings. Additionally, the anticonvulsant profile of CBDV (50-200 mg kg-1) in vivo was investigated in four rodent seizure models: maximal electroshock (mES) and audiogenic seizures in mice, and pentylenetetrazole (PTZ) and pilocarpine-induced seizures in rat. CBDV effects in combination with commonly-used antiepileptic drugs were investigated in rat seizures. Finally, the motor side effect profile of CBDV was investigated using static beam and gripstrength assays. Key results: CDBV significantly attenuated status epilepticus-like epileptiform LFPs induced by 4-AP and Mg2+-free conditions. CBDV had significant anticonvulsant effects in mES (≥100 mg kg-1), audiogenic (≥50 mg kg-1) and PTZ-induced seizures (≥100 mg kg-1). CBDV alone had no effect against pilocarpine-induced seizures, but significantly attenuated these seizures when administered with valproate or phenobarbital at 200 mg kg-1 CBDV. CBDV had no effect on motor function. Conclusions and Implications: These results indicate that CBDV is an effective anticonvulsant across a broad range of seizure models, does not significantly affect normal motor function and therefore merits further investigation in chronic epilepsy models to justify human trials.

Classification of cortical microcircuits based on micro-electrode-array data from slices of rat barrel cortex

The bewildering complexity of cortical microcircuits at the single cell level gives rise to surprisingly robust emergent activity patterns at the level of laminar and columnar local field potentials (LFPs) in response to targeted local stimuli. Here we report the results of our multivariate data-analytic approach based on simultaneous multi-site recordings using micro-electrode-array chips for investigation of the microcircuitary of rat somatosensory (barrel) cortex. We find high repeatability of stimulus-induced responses, and typical spatial distributions of LFP responses to stimuli in supragranular, granular, and infragranular layers, where the last form a particularly distinct class. Population spikes appear to travel with about 33 cm/s from granular to infragranular layers. Responses within barrel related columns have different profiles than those in neighbouring columns to the left or interchangeably to the right. Variations between slices occur, but can be minimized by strictly obeying controlled experimental protocols. Cluster analysis on normalized recordings indicates specific spatial distributions of time series reflecting the location of sources and sinks independent of the stimulus layer. Although the precise correspondences between single cell activity and LFPs are still far from clear, a sophisticated neuroinformatics approach in combination with multi-site LFP recordings in the standardized slice preparation is suitable for comparing normal conditions to genetically or pharmacologically altered situations based on real cortical microcircuitry.

A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

The resting-state neurovascular coupling relationship: rapid changes in spontaneous neural activity in the somatosensory cortex are associated with haemodynamic fluctuations that resemble stimulus-evoked haemodynamics

Although promise exists for patterns of resting-state blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) brain connectivity to be used as biomarkers of early brain pathology, a full understanding of the nature of the relationship between neural activity and spontaneous fMRI BOLD fluctuations is required before such data can be correctly interpreted. To investigate this issue, we combined electrophysiological recordings of rapid changes in multi-laminar local field potentials from the somatosensory cortex of anaesthetized rats with concurrent two-dimensional optical imaging spectroscopy measurements of resting-state haemodynamics that underlie fluctuations in the BOLD fMRI signal. After neural ‘events’ were identified, their time points served to indicate the start of an epoch in the accompanying haemodynamic fluctuations. Multiple epochs for both neural ‘events’ and the accompanying haemodynamic fluctuations were averaged. We found that the averaged epochs of resting-state haemodynamic fluctuations taken after neural ‘events’ closely resembled the temporal profile of stimulus-evoked cortical haemodynamics. Furthermore, we were able to demonstrate that averaged epochs of resting-state haemodynamic fluctuations resembling the temporal profile of stimulus-evoked haemodynamics could also be found after peaks in neural activity filtered into specific electroencephalographic frequency bands (theta, alpha, beta, and gamma). This technique allows investigation of resting-state neurovascular coupling using methodologies that are directly comparable to that developed for investigating stimulus-evoked neurovascular responses.

Resting tremor classification and detection in Parkinson's disease patients

Parkinson is a neurodegenerative disease, in which tremor is the main symptom. This paper investigates the use of different classification methods to identify tremors experienced by Parkinsonian patients.Some previous research has focussed tremor analysis on external body signals (e.g., electromyography, accelerometer signals, etc.). Our advantage is that we have access to sub-cortical data, which facilitates the applicability of the obtained results into real medical devices since we are dealing with brain signals directly. Local field potentials (LFP) were recorded in the subthalamic nucleus of 7 Parkinsonian patients through the implanted electrodes of a deep brain stimulation (DBS) device prior to its internalization. Measured LFP signals were preprocessed by means of splinting, down sampling, filtering, normalization and rec-tification. Then, feature extraction was conducted through a multi-level decomposition via a wavelettrans form. Finally, artificial intelligence techniques were applied to feature selection, clustering of tremor types, and tremor detection.The key contribution of this paper is to present initial results which indicate, to a high degree of certainty, that there appear to be two distinct subgroups of patients within the group-1 of patients according to the Consensus Statement of the Movement Disorder Society on Tremor. Such results may well lead to different resultant treatments for the patients involved, depending on how their tremor has been classified. Moreover, we propose a new approach for demand driven stimulation, in which tremor detection is also based on the subtype of tremor the patient has. Applying this knowledge to the tremor detection problem, it can be concluded that the results improve when patient clustering is applied prior to detection.

Potential contribution of natural enemies to patterns of local adaptation in plants

Genetic differentiation among plant populations and adaptation to local environmental conditions are well documented. However, few studies have examined the potential contribution of plant antagonists, such as insect herbivores and pathogens, to the pattern of local adaptation. Here, a reciprocal transplant experiment was set up at three sites across Europe using two common plant species, Holcus lanatus and Plantago lanceolata. The amount of damage by the main above-ground plant antagonists was measured: a rust fungus infecting Holcus and a specialist beetle feeding on Plantago, both in low-density monoculture plots and in competition with interspecific neighbours. Strong genetic differentiation among provenances in the amount of damage by antagonists in both species was found. Local provenances of Holcus had significantly higher amounts of rust infection than foreign provenances, whereas local provenances of Plantago were significantly less damaged by the specialist beetle than the foreign provenances. The presence of surrounding vegetation affected the amount of damage but had little influence on the ranking of plant provenances. The opposite pattern of population differentiation in resistance to local antagonists in the two species suggests that it will be difficult to predict the consequences of plant translocations for interactions with organisms of higher trophic levels.

Changing management in Scottish birch woodlands: a potential threat to local invertebrate biodiversity

The silvicultural management of Scottish birch woodlands for timber production is replacing traditional low intensity management practices, such as domesticated livestock grazing. These new management practices involve thinning of existing woodlands to prescribed densities to maximize biomass and timber quality. Although presently infrequent, the wide scale adoption of this practice could affect invertebrate community diversity. The impact of these changes in management on Staphylinidae and Carabidae (Coleoptera) in 19 woodlands in Aberdeenshire, north-east Scotland was investigated. Grazing and logging practices were important determinants of beetle community structure. Woodland area had no effect on any measure of beetle community structure, although isolation did influence the abundance of one carabid species. Changes towards timber production forestry will influence the structure of invertebrate communities, although the scale at which this occurs will determine its effect.

The potential of grass field margin management for enhancing beetle diversity in intensive livestock farms

1. Declining populations of UK grassland flora and fauna have been attributed to intensification of agricultural management practices, including changes in cutting, fertilizer, grazing and drainage regimes. We aimed to develop field margin management practices that could reverse declines in intensively managed grassland biodiversity that would have application in the UK and Europe. Here we focus on one aspect of grassland biodiversity, the beetles. 2. In four intensively managed livestock farms in south-west England, 10-m wide field margins in existing grasslands were managed to create seven treatments of increasing sward architectural complexity. This was achieved through combinations of inorganic (NPK) fertilizer, cattle grazing, and timing and height of cutting. To examine the potential influence of complexity on faunal diversity, beetles were identified to species level from suction samples taken between 2003 and 2005, and their assemblage structure was related to margin management, floral assemblages and sward architecture. 3. Beetle abundance, and species richness and evenness were influenced by margin management treatment and its interaction with year. Correlations with sward architecture and the percentage cover of dominant forbs and grasses were also found. Functional groups of the beetles showed different responses to the management treatments. In particular, higher proportional abundances of seed/flower-feeding guilds were found in treatments not receiving NPK fertilizer. 4. The assemblage structure was shown to respond to margin management treatments, sward architecture and the percentage cover of dominant forbs and grasses. The most extensively managed treatments were characterized by distinct successional trajectories from the control treatment. 5. Synthesis and applications. This study provides management options suitable for use within agri-environment schemes intended to improve faunal diversity associated with intensively managed lowland grasslands. Field margins receiving either no management or a single July silage cut were shown to support greater abundances and species richness of beetles, although subtler modifications of conventional management may also be beneficial, for example the absence of NPK fertilizer while maintaining grazing and silage cutting systems.

Properties of singular vectors using convective available potential energy as final time norm

We study the feasibility of using the singular vector technique to create initial condition perturbations for short-range ensemble prediction systems (SREPS) focussing on predictability of severe local storms and in particular deep convection. For this a new final time semi-norm based on the convective available potential energy (CAPE) is introduced. We compare singular vectors using the CAPE-norm with SVs using the more common total energy (TE) norm for a 2-week summer period in 2007, which includes a case of mesoscale extreme rainfall in the south west of Finland. The CAPE singular vectors perturb the CAPE field by increasing the specific humidity and temperature of the parcel and increase the lapse rate above the parcel in the lower troposphere consistent with physical considerations. The CAPE-SVs are situated in the lower troposphere. This in contrast to TE-SVs with short optimization times which predominantly remain in the high troposphere. By examining the time evolution of the CAPE singular values we observe that the convective event in the south west of Finland is clearly associated with high CAPE singular values.

Exploring the market potential of 'local' in food systems

Local food initiatives create a niche market in many developed countries where consumer choice is being met with an expanding offering in both conventional as well as complementary retail outlets. Supermarkets in conjunction with the food service sector currently dominate food sales and consumption, and are likely to do so for the foreseeable future. However, the local food sector offers an opportunity for implementing niche marketing strategies for many businesses. Local food activities tend to be relatively independent activities and a clearer definition for “local” food would assist in consolidating this important component of the food system. Related to this, consumers would benefit from the establishment of some form of assurance system for the ‘localness’ of food. In the UK, with its well established local food market, farmers’ markets, farm shops and box schemes are currently having the largest impact in terms of total sales. Hence further research is required to confirm that support for similar business ventures in Australia would be a viable strategy for strengthening its local food systems.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Up-gradient eddy fluxes of potential vorticity near the subtropical jet

The role of eddy fluxes in the general circulation is often approached by treating eddies as (macro)turbulence. In this approach, the eddies act to diffuse certain quasiconservative quantities, such as potential vorticity (PV), along isentropic surfaces in the free atmosphere. The eddy fluxes are determined primarily by the eddy diffusivities and are necessarily down-gradient of the basic state PV field. Support for the (macro)turbulence approach stems from the fact that the eddy fluxes of PV in the free atmosphere are generally down-gradient in the long-term mean. Here we call attention to a pronounced and significant region of upgradient eddy PV fluxes on the poleward flank of the jet core in both hemispheres. The region of up-gradient (i.e., notionally “antidiffusive”) eddy PV fluxes is most pronounced during the winter and spring seasons and partially contradicts the turbulence approach described above. Analyses of the PV variance (potential enstrophy) budget suggest that the up-gradient PV fluxes represent local wave decay and are maintained by poleward fluxes of PV variance. Finite-amplitude effects thus represent leading order contributions to the PV variance budget, whereas dissipation is only of secondary importance locally. The appearance of up-gradient PV fluxes in the long-term mean is associated with the poleward shift of the jet—and thus the region of wave decay relative to wave growth—following wave-breaking events.

Field-perpendicular and field-aligned plasma flows observed by EISCAT during a prolonged period of northward IMF

The effect of a prolonged period of strongly northward Interplanetary Magnetic Field (IMF) on the high-latitude F-region is studied using data from the EISCAT Common Programme Zero mode of operation on 11–12 August 1982. The analysis of the raw autocorrelation functions is kept to the directly derived parameters Ne, Te, Ti and velocity, and limits are defined for the errors introduced by assumptions about ion composition and by changes in the transmitted power and system constant. Simple data-cleaning criteria are employed to eliminate problems due to coherent signals and large background noise levels. The observed variations in plasma densities, temperatures and velocities are interpreted in terms of supporting data from ISEE-3 and local riometers and magnetometers. Both field-aligned and field-perpendicular plasma flows at Tromsø showed effects of the northward IMF: convection was slow and irregular and field-aligned flow profiles were characteristic of steady-state polar wind outflow with flux of order 1012 m−2 s−1. This period followed a strongly southward IMF which had triggered a substorm. The substorm gave enhanced convection, with a swing to equatorward flow and large (5 × 1012 m−2 s−1), steady-state field-aligned fluxes, leading to the possibility of O+ escape into the magnetosphere. The apparent influence of the IMF over both field-perpendicular and field-aligned flows is explained in terms of the cross-cap potential difference and the location of the auroral oval.