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 gradual depth-dependent change of connectivity features of supragranular pyramidal cells in rat barrel cortex

Recent experimental evidence suggests a finer genetic, structural and functional subdivision of the layers which form a cortical column. The classical layer II/III (LII/III) of rodent neocortex integrates ascending sensory information with contextual cortical information for behavioral read-out. We systematically investigated to which extent regular-spiking supragranular pyramidal neurons, located at different depths within the cortex, show different input-output connectivity patterns. Combining glutamate-uncaging with whole-cell recordings and biocytin filling, we revealed a novel cellular organization of LII/III: (i) “Lower LII/III” pyramidal cells receive a very strong excitatory input from lemniscal LIV and much fewer inputs from paralemniscal LVa. They project to all layers of the home column, including a feedback projection to LIV whereas transcolumnar projections are relatively sparse. (ii) “Upper LII/III” pyramidal cells also receive their strongest input from LIV, but in addition, a very strong and dense excitatory input from LVa. They project extensively to LII/III as well as LVa and Vb of their home and neighboring columns, (iii) “Middle LII/III” pyramidal cell show an intermediate connectivity phenotype that stands in many ways in-between the features described for lower versus upper LII/III. “Lower LII/III” intracolumnarly segregates and transcolumnarly integrates lemniscal information whereas “upper LII/III” seems to integrate lemniscal with paralemniscal information. This suggests a finegrained functional subdivision of the supragranular compartment containing multiple circuits without any obvious cytoarchitectonic, other structural or functional correlate of a laminar border in rodent barrel cortex.

Balanced excitation and inhibition: Model based analysis of local field potentials

We have developed a model of the local field potential (LFP) based on the conservation of charge, the independence principle of ionic flows and the classical Hodgkin–Huxley (HH) type intracellular model of synaptic activity. Insights were gained through the simulation of the HH intracellular model on the nonlinear relationship between the balance of synaptic conductances and that of post-synaptic currents. The latter is dependent not only on the former, but also on the temporal lag between the excitatory and inhibitory conductances, as well as the strength of the afferent signal. The proposed LFP model provides a method for decomposing the LFP recordings near the soma of layer IV pyramidal neurons in the barrel cortex of anaesthetised rats into two highly correlated components with opposite polarity. The temporal dynamics and the proportional balance of the two components are comparable to the excitatory and inhibitory post-synaptic currents computed from the HH model. This suggests that the two components of the LFP reflect the underlying excitatory and inhibitory post-synaptic currents of the local neural population. We further used the model to decompose a sequence of evoked LFP responses under repetitive electrical stimulation (5 Hz) of the whisker pad. We found that as neural responses adapted, the excitatory and inhibitory components also adapted proportionately, while the temporal lag between the onsets of the two components increased during frequency adaptation. Our results demonstrated that the balance between neural excitation and inhibition can be investigated using extracellular recordings. Extension of the model to incorporate multiple compartments should allow more quantitative interpretations of surface Electroencephalography (EEG) recordings into components reflecting the excitatory, inhibitory and passive ionic current flows generated by local neural populations.

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 model of neurovascular coupling: Implications for blood vessel dilation and constriction

Neurovascular coupling in response to stimulation of the rat barrel cortex was investigated using concurrent multichannel electrophysiology and laser Doppler flowmetry. The data were used to build a linear dynamic model relating neural activity to blood flow. Local field potential time series were subject to current source density analysis, and the time series of a layer IV sink of the barrel cortex was used as the input to the model. The model output was the time series of the changes in regional cerebral blood flow (CBF). We show that this model can provide excellent fit of the CBF responses for stimulus durations of up to 16 s. The structure of the model consisted of two coupled components representing vascular dilation and constriction. The complex temporal characteristics of the CBF time series were reproduced by the relatively simple balance of these two components. We show that the impulse response obtained under the 16-s duration stimulation condition generalised to provide a good prediction to the data from the shorter duration stimulation conditions. Furthermore, by optimising three out of the total of nine model parameters, the variability in the data can be well accounted for over a wide range of stimulus conditions. By establishing linearity, classic system analysis methods can be used to generate and explore a range of equivalent model structures (e.g., feed-forward or feedback) to guide the experimental investigation of the control of vascular dilation and constriction following stimulation.

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.

Long duration stimuli and nonlinearities in the neural–haemodynamic coupling

Recent studies have shown that the haemodynamic responses to brief (<2 secs) stimuli can be well characterised as a linear convolution of neural activity with a suitable haemodynamic impulse response. In this paper, we show that the linear convolution model cannot predict measurements of blood flow responses to stimuli of longer duration (>2 secs), regardless of the impulse response function chosen. Modifying the linear convolution scheme to a nonlinear convolution scheme was found to provide a good prediction of the observed data. Whereas several studies have found a nonlinear coupling between stimulus input and blood flow responses, the current modelling scheme uses neural activity as an input, and thus implies nonlinearity in the coupling between neural activity and blood flow responses. Neural activity was assessed by current source density analysis of depth-resolved evoked field potentials, while blood flow responses were measured using laser Doppler flowmetry. All measurements were made in rat whisker barrel cortex after electrical stimulation of the whisker pad for 1 to 16 secs at 5 Hz and 1.2 mA (individual pulse width 0.3 ms).

Long-latency reductions in gamma power predict hemodynamic changes that underlie the negative BOLD signal

Studiesthat use prolonged periods of sensory stimulation report associations between regional reductions in neural activity and negative blood oxygenation level-dependent (BOLD) signaling. However, the neural generators of the negative BOLD response remain to be characterized. Here, we use single-impulse electrical stimulation of the whisker pad in the anesthetized rat to identify components of the neural response that are related to “negative” hemodynamic changes in the brain. Laminar multiunit activity and local field potential recordings of neural activity were performed concurrently withtwo-dimensional optical imaging spectroscopy measuring hemodynamic changes. Repeated measurements over multiple stimulation trials revealed significant variations in neural responses across session and animal datasets. Within this variation, we found robust long-latency decreases (300 and 2000 ms after stimulus presentation) in gammaband power (30 – 80 Hz) in the middle-superficial cortical layers in regions surrounding the activated whisker barrel cortex. This reduction in gamma frequency activity was associated with corresponding decreases in the hemodynamic responses that drive the negative BOLD signal. These findings suggest a close relationship between BOLD responses and neural events that operate over time scales that outlast the initiating sensory stimulus, and provide important insights into the neurophysiological basis of negative neuroimaging signals.

Effects of Delta9-tetrahydrocannabivarin on [35S]GTPgammaS binding in mouse brain cerebellum and piriform cortex membranes

BACKGROUND AND PURPOSE: We have recently shown that the phytocannabinoid Delta9-tetrahydrocannabivarin (Delta9-THCV) and the CB1 receptor antagonist AM251 increase inhibitory neurotransmission in mouse cerebellum and also exhibit anticonvulsant activity in a rat piriform cortical (PC) model of epilepsy. Possible mechanisms underlying cannabinoid actions in the CNS include CB1 receptor antagonism (by displacing endocannabinergic tone) or inverse agonism at constitutively active CB1 receptors. Here, we investigate the mode of cannabinoid action in [35S]GTPgammaS binding assays. EXPERIMENTAL APPROACH: Effects of Delta9-THCV and AM251 were tested either alone or against WIN55,212-2-induced increases in [35S]GTPgammaS binding in mouse cerebellar and PC membranes. Effects on non-CB receptor expressing CHO-D2 cell membranes were also investigated. KEY RESULTS :Delta9-THCV and AM251 both acted as potent antagonists of WIN55,212-2-induced increases in [35S]GTPgammaS binding in cerebellar and PC membranes (Delta9-THCV: pA2=7.62 and 7.44 respectively; AM251: pA2=9.93 and 9.88 respectively). At micromolar concentrations, Delta9-THCV or AM251 alone caused significant decreases in [35S]GTPgammaS binding; Delta9-THCV caused larger decreases than AM251. When applied alone in CHO-D2 membranes, Delta9-THCV and AM251 also caused concentration-related decreases in G protein activity. CONCLUSIONS AND IMPLICATIONS: Delta9-THCV and AM251 act as CB1 receptors antagonists in the cerebellum and PC, with AM251 being more potent than Delta9-THCV in both brain regions. Individually, Delta9-THCV or AM251 exhibited similar potency at CB1 receptors in the cerebellum and the PC. At micromolar concentrations, Delta9-THCV and AM251 caused a non-CB receptor-mediated depression of basal [35S]GTPgammaS binding.

Sox1-deficient mice suffer from epilepsy associated with abnormal ventral forebrain development and olfactory cortex hyperexcitability

Mutations in several classes of embryonically-expressed transcription factor genes are associated with behavioral disorders and epilepsies. However, there is little known about how such genetic and neurodevelopmental defects lead to brain dysfunction. Here we present the characterization of an epilepsy syndrome caused by the absence of the transcription factor SOX1 in mice. In vivo electroencephalographic recordings from SOX1 mutants established a correlation between behavioral changes and cortical output that was consistent with a seizure origin in the limbic forebrain. In vitro intracellular recordings from three major forebrain regions, neocortex, hippocampus and olfactory (piriform) cortex (OC) showed that only the OC exhibits abnormal enhanced synaptic excitability and spontaneous epileptiform discharges. Furthermore, the hyperexcitability of the OC neurons was present in mutants prior to the onset of seizures but was completely absent from both the hippocampus and neocortex of the same animals. The local inhibitory GABAergic neurotransmission remained normal in the OC of SOX1-deficient brains, but there was a severe developmental deficit of OC postsynaptic target neurons, mainly GABAergic projection neurons within the olfactory tubercle and the nucleus accumbens shell. Our data show that SOX1 is essential for ventral telencephalic development and suggest that the neurodevelopmental defect disrupts local neuronal circuits leading to epilepsy in the SOX1-deficient mice

Cholinergic modulation of intrinsic fibre-evoked excitatory transmission contains a nicotinic component in immature but not adult rat piriform cortex, in vitro

The piriform cortex (PC) is highly prone to epileptogenesis, particularly in immature animals, where decreased muscarinic modulation of PC intrinsic fibre excitatory neurotransmission is implicated as a likely cause. However, whether higher levels of acetylcholine (ACh) release occur in immature vs. adult PC remains unclear. We investigated this using in vitro extracellular electrophysiological recording techniques. Intrinsic fibre-evoked extracellular field potentials (EFPs) were recorded from layers II to III in PC brain slices prepared from immature (P14-18) and adult (P>40) rats. Adult and immature PC EFPs were suppressed by eserine (1muM) or neostigmine (1muM) application, with a greater suppression in immature ( approximately 40%) than adult ( approximately 30%) slices. Subsequent application of atropine (1muM) reversed EFP suppression, producing supranormal ( approximately 12%) recovery in adult slices, suggesting that suppression was solely muscarinic ACh receptor-mediated and that some 'basal' cholinergic 'tone' was present. Conversely, atropine only partially reversed anticholinesterase effects in immature slices, suggesting the presence of additional non-muscarinic modulation. Accordingly, nicotine (50muM) caused immature field suppression ( approximately 30%) that was further enhanced by neostigmine, whereas it had no effect on adult EFPs. Unlike atropine, nicotinic antagonists, mecamylamine and methyllycaconitine, induced immature supranormal field recovery ( approximately 20%) following anticholinesterase-induced suppression (with no effect on adult slices), confirming that basal cholinergic 'tone' was also present. We suggest that nicotinic inhibitory cholinergic modulation occurs in the immature rat PC intrinsic excitatory fibre system, possibly to complement the existing, weak muscarinic modulation, and could be another important developmentally regulated system governing immature PC susceptibility towards epileptogenesis.

Developmental changes in presynaptic muscarinic modulation of excitatory and inhibitory neurotransmission in rat piriform cortex in vitro: relevance to epileptiform bursting susceptibility

Suppression of depolarizing postsynaptic potentials and isolated GABA-A receptor-mediated fast inhibitory postsynaptic potentials by the muscarinic acetylcholine receptor agonist, oxotremorine-M (10 microM), was investigated in adult and immature (P14-P30) rat piriform cortical (PC) slices using intracellular recording. Depolarizing postsynaptic potentials evoked by layers II-III stimulation underwent concentration-dependent inhibition in oxotremorine-M that was most likely presynaptic and M2 muscarinic acetylcholine receptor-mediated in immature, but M1-mediated in adult (P40-P80) slices; percentage inhibition was smaller in immature than in adult piriform cortex. In contrast, compared with adults, layer Ia-evoked depolarizing postsynaptic potentials in immature piriform cortex slices in oxotremorine-M, showed a prolonged multiphasic depolarization with superimposed fast transients and spikes, and an increased 'all-or-nothing' character. Isolated N-methyl-d-aspartate receptor-mediated layer Ia depolarizing postsynaptic potentials (although significantly larger in immature slices) were however, unaffected by oxotremorine-M, but blocked by dl-2-amino-5-phosphonovaleric acid. Fast inhibitory postsynaptic potentials evoked by layer Ib or layers II-III-fiber stimulation in immature slices were significantly smaller than in adults, despite similar estimated mean reversal potentials ( approximately -69 and -70 mV respectively). In oxotremorine-M, only layer Ib-fast inhibitory postsynaptic potentials were suppressed; suppression was again most likely presynaptic M2-mediated in immature slices, but M1-mediated in adults. The degree of fast inhibitory postsynaptic potential suppression was however, greater in immature than in adult piriform cortex. Our results demonstrate some important physiological and pharmacological differences between excitatory and inhibitory synaptic systems in adult and immature piriform cortex that could contribute toward the increased susceptibility of this region to muscarinic agonist-induced epileptiform activity in immature brain slices.

Further characterization of muscarinic agonist-induced epileptiform bursting activity in immature rat piriform cortex, in vitro

The characteristics of muscarinic acetylcholine receptor agonist-induced epileptiform bursting seen in immature rat piriform cortex slices in vitro were further investigated using intracellular recording, with particular focus on its postnatal age-dependence (P+14-P+30), pharmacology, site(s) of origin and the likely contribution of the muscarinic acetylcholine receptor agonist-induced post-stimulus slow afterdepolarization and gap junction functionality toward its generation. The muscarinic agonist, oxotremorine-M (10 microM), induced rhythmic bursting only in immature piriform cortex slices; however, paroxysmal depolarizing shift amplitude, burst duration and burst incidence were inversely related to postnatal age. No significant age-dependent changes in neuronal membrane properties or postsynaptic muscarinic responsiveness accounted for this decline. Burst incidence was higher when recorded in anterior and posterior regions of the immature piriform cortex. In adult and immature neurones, oxotremorine-M effects were abolished by M1-, but not M2-muscarinic acetylcholine receptor-selective antagonists. Rostrocaudal lesions, between piriform cortex layers I and II, or layer III and endopiriform nucleus in adult or immature slices did not influence oxotremorine-M effects; however, the slow afterdepolarization in adult (but not immature) lesioned slices was abolished. Gap junction blockers (carbenoxolone or octanol) disrupted muscarinic bursting and diminished the slow afterdepolarization in immature slices, suggesting that gap junction connectivity was important for bursting. Our data show that neural networks within layers II-III function as primary oscillatory circuits for burst initiation in immature rat piriform cortex during persistent muscarinic receptor activation. Furthermore, we propose that muscarinic slow afterdepolarization induction and gap junction communication could contribute towards the increased epileptiform susceptibility of this brain area.

A novel component of cannabis extract potentiates excitatory synaptic transmission in rat olfactory cortex in vitro

Cannabis is a potential treatment for epilepsy, although the few human studies supporting this use have proved inconclusive. Previously, we showed that a standardized cannabis extract (SCE), isolated Delta(9)-tetrahydrocannabinol (Delta(9)-THC), and even Delta(9)-THC-free SCE inhibited muscarinic agonist-induced epileptiform bursting in rat olfactory cortical brain slices, acting via CB1 receptors. The present work demonstrates that although Delta(9)-THC (1microM) significantly depressed evoked depolarizing postsynaptic potentials (PSPs) in rat olfactory cortex neurones, both SCE and Delta(9)-THC-free SCE significantly potentiated evoked PSPs (all results were fully reversed by the CB1 receptor antagonist SR141716A, 1microM); interestingly, the potentiation by Delta(9)-THC-free SCE was greater than that produced by SCE. On comparing the effects of Delta(9)-THC-free SCE upon evoked PSPs and artificial PSPs (aPSPs; evoked electrotonically following brief intracellular current injection), PSPs were enhanced, whereas aPSPs were unaffected, suggesting that the effect was not due to changes in background input resistance. Similar recordings made using CB1 receptor-deficient knockout mice (CB1(-/-)) and wild-type littermate controls revealed cannabinoid or extract-induced changes in membrane resistance, cell excitability and synaptic transmission in wild-type mice that were similar to those seen in rat neurones, but no effect on these properties were seen in CB1(-/-) cells. It appears that the unknown extract constituent(s) effects over-rode the suppressive effects of Delta(9)-THC on excitatory neurotransmitter release, which may explain some patients' preference for herbal cannabis rather than isolated Delta(9)-THC (due to attenuation of some of the central Delta(9)-THC side effects) and possibly account for the rare incidence of seizures in some individuals taking cannabis recreationally

Inverse amygdala and medial prefrontal cortex responses to surprised faces

Here we show inverse fMRI activation patterns in amygdala and medial prefrontal cortex (mPFC) depending upon whether subjects interpreted surprised facial expressions positively or negatively. More negative interpretations of surprised faces were associated with greater signal changes in the right ventral amygdala, while more positive interpretations were associated with greater signal changes in the ventral mPFC. Accordingly, signal change within these two areas was inversely correlated. Thus, individual differences in the judgment of surprised faces are related to a systematic inverse relationship between amygdala and mPFC activity, a circuitry that the animal literature suggests is critical to the assessment of stimuli that predict potential positive vs negative outcomes.