Interface dynamics in planar neural field models

Neural field models describe the coarse-grained activity of populations of interacting neurons. Because of the laminar structure of real cortical tissue they are often studied in two spatial dimensions, where they are well known to generate rich patterns of spatiotemporal activity. Such patterns have been interpreted in a variety of contexts ranging from the understanding of visual hallucinations to the generation of electroencephalographic signals. Typical patterns include localized solutions in the form of traveling spots, as well as intricate labyrinthine structures. These patterns are naturally defined by the interface between low and high states of neural activity. Here we derive the equations of motion for such interfaces and show, for a Heaviside firing rate, that the normal velocity of an interface is given in terms of a non-local Biot-Savart type interaction over the boundaries of the high activity regions. This exact, but dimensionally reduced, system of equations is solved numerically and shown to be in excellent agreement with the full nonlinear integral equation defining the neural field. We develop a linear stability analysis for the interface dynamics that allows us to understand the mechanisms of pattern formation that arise from instabilities of spots, rings, stripes and fronts. We further show how to analyze neural field models with linear adaptation currents, and determine the conditions for the dynamic instability of spots that can give rise to breathers and traveling waves.

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.

Realistic mean field forward predictions for the integration of co-registered EEG/fMRI

Brain activity can be measured non-invasively with functional imaging techniques. Each pixel in such an image represents a neural mass of about 105 to 107 neurons. Mean field models (MFMs) approximate their activity by averaging out neural variability while retaining salient underlying features, like neurotransmitter kinetics. However, MFMs incorporating the regional variability, realistic geometry and connectivity of cortex have so far appeared intractable. This lack of biological realism has led to a focus on gross temporal features of the EEG. We address these impediments and showcase a "proof of principle" forward prediction of co-registered EEG/fMRI for a full-size human cortex in a realistic head model with anatomical connectivity, see figure 1. MFMs usually assume homogeneous neural masses, isotropic long-range connectivity and simplistic signal expression to allow rapid computation with partial differential equations. But these approximations are insufficient in particular for the high spatial resolution obtained with fMRI, since different cortical areas vary in their architectonic and dynamical properties, have complex connectivity, and can contribute non-trivially to the measured signal. Our code instead supports the local variation of model parameters and freely chosen connectivity for many thousand triangulation nodes spanning a cortical surface extracted from structural MRI. This allows the introduction of realistic anatomical and physiological parameters for cortical areas and their connectivity, including both intra- and inter-area connections. Proper cortical folding and conduction through a realistic head model is then added to obtain accurate signal expression for a comparison to experimental data. To showcase the synergy of these computational developments, we predict simultaneously EEG and fMRI BOLD responses by adding an established model for neurovascular coupling and convolving "Balloon-Windkessel" hemodynamics. We also incorporate regional connectivity extracted from the CoCoMac database [1]. Importantly, these extensions can be easily adapted according to future insights and data. Furthermore, while our own simulation is based on one specific MFM [2], the computational framework is general and can be applied to models favored by the user. Finally, we provide a brief outlook on improving the integration of multi-modal imaging data through iterative fits of a single underlying MFM in this realistic simulation framework.

Cortical thickness changes in the non-lesioned hemisphere associated with non-paretic arm immobilization in modified CI therapy

Recent evidence suggests that immobilization of the upper limb for 2–3 weeks induces changes in cortical thickness as well as motor performance. In constraint induced (CI) therapy, one of the most effective interventions for hemiplegia, the non-paretic arm is constrained to enforce the use of the paretic arm in the home setting. With the present study we aimed to explore whether non-paretic arm immobilization in CI therapy induces structural changes in the non-lesioned hemisphere, and how these changes are related to treatment benefit. 31 patients with chronic hemiparesis participated in CI therapy with (N = 14) and without (N = 17) constraint. Motor ability scores were acquired before and after treatment. Diffusion tensor imaging (DTI) data was obtained prior to treatment. Cortical thickness was measured with the Freesurfer software. In both groups cortical thickness in the contralesional primary somatosensory cortex increased and motor function improved with the intervention. However the cortical thickness change was not associated with the magnitude of motor function improvement. Moreover, the treatment effect and the cortical thickness change were not significantly different between the constraint and the non-constraint groups. There was no correlation between fractional anisotropy changes in the non-lesioned hemisphere and treatment outcome. CI therapy induced cortical thickness changes in contralesional sensorimotor regions, but this effect does not appear to be driven by the immobilization of the non-paretic arm, as indicated by the absence of differences between the constraint and the non-constraint groups. Our data does not suggest that the arm immobilization used in CI therapy is associated with noticeable cortical thinning.

Early and late stimulus-evoked cortical hemodynamic responses provide insight into the neurogenic nature of neurovascular coupling

Understanding neurovascular coupling is a prerequisite for the interpretation of results obtained from modern neuroimaging techniques. This study investigated the hemodynamic and neural responses in rat somatosensory cortex elicited by 16 seconds electrical whisker stimuli. Hemodynamics were measured by optical imaging spectroscopy and neural activity by multichannel electrophysiology. Previous studies have suggested that the whisker-evoked hemodynamic response contains two mechanisms, a transient ‘backwards’ dilation of the middle cerebral artery, followed by an increase in blood volume localized to the site of neural activity. To distinguish between the mechanisms responsible for these aspects of the response, we presented whisker stimuli during normocapnia (‘control’), and during a high level of hypercapnia. Hypercapnia was used to ‘predilate’ arteries and thus possibly ‘inhibit’ aspects of the response related to the ‘early’ mechanism. Indeed, hemodynamic data suggested that the transient stimulus-evoked response was absent under hypercapnia. However, evoked neural responses were also altered during hypercapnia and convolution of the neural responses from both the normocapnic and hypercapnic conditions with a canonical impulse response function, suggested that neurovascular coupling was similar in both conditions. Although data did not clearly dissociate early and late vascular responses, they suggest that the neurovascular coupling relationship is neurogenic in origin.

Temporal coupling between stimulus-evoked neural activity and hemodynamic responses from individual cortical columns

Using previously published data from the whisker barrel cortex of anesthetized rodents (Berwick et al 2008 J. Neurophysiol. 99 787–98) we investigated whether highly spatially localized stimulus-evoked cortical hemodynamics responses displayed a linear time-invariant (LTI) relationship with neural activity. Presentation of stimuli to individual whiskers of 2 s and 16 s durations produced hemodynamics and neural activity spatially localized to individual cortical columns. Two-dimensional optical imaging spectroscopy (2D-OIS) measured hemoglobin responses, while multi-laminar electrophysiology recorded neural activity. Hemoglobin responses to 2 s stimuli were deconvolved with underlying evoked neural activity to estimate impulse response functions which were then convolved with neural activity evoked by 16 s stimuli to generate predictions of hemodynamic responses. An LTI system more adequately described the temporal neuro-hemodynamics coupling relationship for these spatially localized sensory stimuli than in previous studies that activated the entire whisker cortex. An inability to predict the magnitude of an initial 'peak' in the total and oxy- hemoglobin responses was alleviated when excluding responses influenced by overlying arterial components. However, this did not improve estimation of the hemodynamic responses return to baseline post-stimulus cessation.

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.

Phenolic acid intake, delivered via moderate champagne wine consumption, improves spatial working memory via the modulation of hippocampal and cortical protein expression/activation.

Aims: While much data exist for the effects of flavonoid-rich foods on spatial memory in rodents, there are no such data for foods/beverages predominantly containing hydroxycinnamates and phenolic acids. To address this, we investigated the effects of moderate Champagne wine intake, which is rich in these components, on spatial memory and related mechanisms relative to the alcohol- and energy-matched controls. Results: In contrast to the isocaloric and alcohol-matched controls, supplementation with Champagne wine (1.78 ml/kg BW, alcohol 12.5% vol.) for 6 weeks led to an improvement in spatial working memory in aged rodents. Targeted protein arrays indicated that these behavioral effects were paralleled by the differential expression of a number of hippocampal and cortical proteins (relative to the isocaloric control group), including those involved in signal transduction, neuroplasticity, apoptosis, and cell cycle regulation. Western immunoblotting confirmed the differential modulation of brain-derived neurotrophic factor, cAMP response-element-binding protein (CREB), p38, dystrophin, 2',3'-cyclic-nucleotide 3'-phosphodiesterase, mammalian target of rapamycin (mTOR), and Bcl-xL in response to Champagne supplementation compared to the control drink, and the modulation of mTOR, Bcl-xL, and CREB in response to alcohol supplementation. Innovation: Our data suggest that smaller phenolics such as gallic acid, protocatechuic acid, tyrosol, caftaric acid, and caffeic acid, in addition to flavonoids, are capable of exerting improvements in spatial memory via the modulation in hippocampal signaling and protein expression. Conclusion: Changes in spatial working memory induced by the Champagne supplementation are linked to the effects of absorbed phenolics on cytoskeletal proteins, neurotrophin expression, and the effects of alcohol on the regulation of apoptotic events in the hippocampus and cortex. Antioxid. Redox Signal. 00, 000-000.

Isolating single primary rat hippocampal neurons & astrocytes on ultra-thin patterned parylene-C/silicon dioxide substrates

We report here the patterning of primary rat neurons and astrocytes from the postnatal hippocampus on ultra-thin parylene-C deposited on a silicon dioxide substrate, following observations of neuronal, astrocytic and nuclear coverage on strips of different lengths, widths and thicknesses. Neuronal and glial growth was characterized ‘on’, ‘adjacent to’ and ‘away from’ the parylene strips. In addition, the article reports how the same material combination can be used to isolate single cells along thin tracks of parylene-C. This is demonstrated with a series of high magnification images of the experimental observations for varying parylene strip widths and thicknesses. Thus, the findings demonstrate the possibility to culture cells on ultra-thin layers of parylene-C and localize single cells on thin strips. Such work is of interest and significance to the Neuroengineering and Multi-Electrode Array (MEA) communities, as it provides an alternative insulating material in the fabrication of embedded micro-electrodes, which can be used to facilitate single cell stimulation and recording in capacitive coupling mode.

First human hNT neurons patterned on parylene-C/silicon dioxide substrates: Combining an accessible cell line and robust patterning technology for the study of the pathological adult human brain

In this communication, we describe a new method which has enabled the first patterning of human neurons (derived from the human teratocarcinoma cell line (hNT)) on parylene-C/silicon dioxide substrates. We reveal the details of the nanofabrication processes, cell differentiation and culturing protocols necessary to successfully pattern hNT neurons which are each key aspects of this new method. The benefits in patterning human neurons on silicon chip using an accessible cell line and robust patterning technology are of widespread value. Thus, using a combined technology such as this will facilitate the detailed study of the pathological human brain at both the single cell and network level.

Guided growth of neurons and glia using microfabricated patterns of parylene-C on a SiO2 background

This paper describes a simple technique for the patterning of glia and neurons. The integration of neuronal patterning to Multi-Electrode Arrays (MEAs), planar patch clamp and silicon based ‘lab on a chip’ technologies necessitates the development of a microfabrication-compatible method, which will be reliable and easy to implement. In this study a highly consistent, straightforward and cost effective cell patterning scheme has been developed. It is based on two common ingredients: the polymer parylene-C and horse serum. Parylene-C is deposited and photo-lithographically patterned on silicon oxide (SiO2) surfaces. Subsequently, the patterns are activated via immersion in horse serum. Compared to non-activated controls, cells on the treated samples exhibited a significantly higher conformity to underlying parylene stripes. The immersion time of the patterns was reduced from 24 to 3 h without compromising the technique. X-ray photoelectron spectroscopy (XPS) analysis of parylene and SiO2 surfaces before and after immersion in horse serum and gel based eluant analysis suggests that the quantity and conformation of proteins on the parylene and SiO2 substrates might be responsible for inducing glial and neuronal patterning.

An epigenetic signature of developmental potential in neural stem cells and early neurons.

A cardinal property of neural stem cells (NSCs) is their ability to adopt multiple fates upon differentiation. The epigenome is widely seen as a read-out of cellular potential and a manifestation of this can be seen in embryonic stem cells (ESCs), where promoters of many lineage-specific regulators are marked by a bivalent epigenetic signature comprising trimethylation of both lysine 4 and lysine 27 of histone H3 (H3K4me3 and H3K27me3, respectively). Bivalency has subsequently emerged as a powerful epigenetic indicator of stem cell potential. Here, we have interrogated the epigenome during differentiation of ESC-derived NSCs to immature GABAergic interneurons. We show that developmental transitions are accompanied by loss of bivalency at many promoters in line with their increasing developmental restriction from pluripotent ESC through multipotent NSC to committed GABAergic interneuron. At the NSC stage, the promoters of genes encoding many transcriptional regulators required for differentiation of multiple neuronal subtypes and neural crest appear to be bivalent, consistent with the broad developmental potential of NSCs. Upon differentiation to GABAergic neurons, all non-GABAergic promoters resolve to H3K27me3 monovalency, whereas GABAergic promoters resolve to H3K4me3 monovalency or retain bivalency. Importantly, many of these epigenetic changes occur before any corresponding changes in gene expression. Intriguingly, another group of gene promoters gain bivalency as NSCs differentiate toward neurons, the majority of which are associated with functions connected with maturation and establishment and maintenance of connectivity. These data show that bivalency provides a dynamic epigenetic signature of developmental potential in both NSCs and in early neurons. Stem Cells 2013;31:1868-1880.

Fibroblast growth factor 2 maintains the neurogenic capacity of embryonic neural progenitor cells in vitro but changes their neuronal subtype specification

Many in vitro systems used to examine multipotential neural progenitor cells (NPCs) rely on mitogens including fibroblast growth factor 2 (FGF2) for their continued expansion. However, FGF2 has also been shown to alter the expression of transcription factors (TFs) that determine cell fate. Here, we report that NPCs from the embryonic telencephalon grown without FGF2 retain many of their in vivo characteristics, making them a good model for investigating molecular mechanisms involved in cell fate specification and differentiation. However, exposure of cortical NPCs to FGF2 results in a profound change in the types of neurons generated, switching them from a glutamatergic to a GABAergic phenotype. This change closely correlates with the dramatic upregulation of TFs more characteristic of ventral telencephalic NPCs. In addition, exposure of cortical NPCs to FGF2 maintains their neurogenic potential in vitro, and NPCs spontaneously undergo differentiation following FGF2 withdrawal. These results highlight the importance of TFs in determining the types of neurons generated by NPCs in vitro. In addition, they show that FGF2, as well as acting as a mitogen, changes the developmental capabilities of NPCs. These findings have implications for the cell fate specification of in vitro-expanded NPCs and their ability to generate specific cell types for therapeutic applications. Disclosure of potential conflicts of interest is found at the end of this article.

Expression of the guanine nucleotide exchange factor, mr-gef, is regulated during the differentiation of specific subsets of telencephalic neurons

We have performed a screen combining subtractive hybridization with PCR to isolate genes that are regulated when neuroepithelial (NE) cells differentiate into neurons. From this screen, we have isolated a number of known genes that have not previously been associated with neurogenesis, together with several novel genes. Here we report that one of these genes, encoding a guanine nucleotide exchange factor (GEF), is regulated during the differentiation of distinct neuronal populations. We have cloned both rat and mouse GEF genes and shown that they are orthologs of the human gene, MR-GEF, which encodes a GEF that specifically activates the small GTPase, Rap1. We have therefore named the rat gene rat mr-gef (rmr-gef) and the mouse gene mouse mr-gef (mmr-gef). Here, we will collectively refer to these two rodent genes as mr-gef. Expression studies show that mr-gef is expressed by young neurons of the developing rodent CNS but not by progenitor cells in the ventricular zone (VZ). The expression pattern of mr-gef during early telencephalic neurogenesis is strikingly similar to that of GABA and the LIM homeobox gene Lhx6, a transcription factor expressed by GABAergic interneurons generated in the ventral telencephalon, some of which migrate into the cortex during development. These observations suggest that mr-gef encodes a protein that is part of a signaling pathway involved in telencephalic neurogenesis; particularly in the development of GABAergic interneurons.

Mechanisms of protease-activated receptor 2-evoked hyperexcitability of nociceptive neurons innervating the mouse colon

Agonists of protease-activated receptor 2 (PAR(2)) evoke hyperexcitability of dorsal root ganglia (DRG) neurons by unknown mechanisms. We examined the cellular mechanisms underlying PAR(2)-evoked hyperexcitability of mouse colonic DRG neurons to determine their potential role in pain syndromes such as visceral hyperalgesia. Colonic DRG neurons were identified by injecting Fast Blue and DiI retrograde tracers into the mouse colon. Using immunofluorescence, we found that DiI-labelled neurons contained PAR(2) immunoreactivity, confirming the presence of receptors on colonic neurons. Whole-cell current-clamp recordings of acutely dissociated neurons demonstrated that PAR(2) activation with a brief application (3 min) of PAR(2) agonists, SLIGRL-NH(2) and trypsin, evoked sustained depolarizations (up to 60 min) which were associated with increased input resistance and a marked reduction in rheobase (50% at 30 min). In voltage clamp, SLIGRL-NH(2) markedly suppressed delayed rectifier I(K) currents (55% at 10 min), but had no effect on the transient I(A) current or TTX-resistant Na(+) currents. In whole-cell current-clamp recordings, the sustained excitability evoked by PAR(2) activation was blocked by the PKC inhibitor, calphostin, and the ERK(1/2) inhibitor PD98059. Studies of ERK(1/2) phosphorylation using confocal microscopy demonstrated that SLIGRL-NH(2) increased levels of immunoreactive pERK(1/2) in DRG neurons, particularly in proximity to the plasma membrane. Thus, activation of PAR(2) receptors on colonic nociceptive neurons causes sustained hyperexcitability that is related, at least in part, to suppression of delayed rectifier I(K) currents. Both PKC and ERK(1/2) mediate the PAR(2)-induced hyperexcitability. These studies describe a novel mechanism of sensitization of colonic nociceptive neurons that may be implicated in conditions of visceral hyperalgesia such as irritable bowel syndrome.