Biased competition in the absence of input bias revealed through corticostriatal computation.

Classical accounts of biased competition require an input bias to resolve the competition between neuronal ensembles driving downstream processing. However, flexible and reliable selection of behaviorally relevant ensembles can occur with unbiased stimulation: striatal D1 and D2 spiny projection neurons (SPNs) receive balanced cortical input, yet their activity determines the choice between GO and NO-GO pathways in the basal ganglia. We here present a corticostriatal model identifying three mechanisms that rely on physiological asymmetries to effect rate- and time-coded biased competition in the presence of balanced inputs. First, tonic input strength determines which one of the two SPN phenotypes exhibits a higher mean firing rate. Second, low-strength oscillatory inputs induce higher firing rate in D2 SPNs but higher coherence between D1 SPNs. Third, high-strength inputs oscillating at distinct frequencies can preferentially activate D1 or D2 SPN populations. Of these mechanisms, only the latter accommodates observed rhythmic activity supporting rule-based decision making in prefrontal cortex.

Flexible resonance in prefrontal networks with strong feedback inhibition.

Oscillations are ubiquitous features of brain dynamics that undergo task-related changes in synchrony, power, and frequency. The impact of those changes on target networks is poorly understood. In this work, we used a biophysically detailed model of prefrontal cortex (PFC) to explore the effects of varying the spike rate, synchrony, and waveform of strong oscillatory inputs on the behavior of cortical networks driven by them. Interacting populations of excitatory and inhibitory neurons with strong feedback inhibition are inhibition-based network oscillators that exhibit resonance (i.e., larger responses to preferred input frequencies). We quantified network responses in terms of mean firing rates and the population frequency of network oscillation; and characterized their behavior in terms of the natural response to asynchronous input and the resonant response to oscillatory inputs. We show that strong feedback inhibition causes the PFC to generate internal (natural) oscillations in the beta/gamma frequency range (>15 Hz) and to maximize principal cell spiking in response to external oscillations at slightly higher frequencies. Importantly, we found that the fastest oscillation frequency that can be relayed by the network maximizes local inhibition and is equal to a frequency even higher than that which maximizes the firing rate of excitatory cells; we call this phenomenon population frequency resonance. This form of resonance is shown to determine the optimal driving frequency for suppressing responses to asynchronous activity. Lastly, we demonstrate that the natural and resonant frequencies can be tuned by changes in neuronal excitability, the duration of feedback inhibition, and dynamic properties of the input. Our results predict that PFC networks are tuned for generating and selectively responding to beta- and gamma-rhythmic signals due to the natural and resonant properties of inhibition-based oscillators. They also suggest strategies for optimizing transcranial stimulation and using oscillatory networks in neuromorphic engineering.

Multimodal imaging of repetition priming: Using fMRI, MEG, and intracranial EEG to reveal spatiotemporal profiles of word processing.

Repetition priming is a core feature of memory processing whose anatomical correlates remain poorly understood. In this study, we use advanced multimodal imaging (functional magnetic resonance imaging (fMRI) and magnetoencephalography; MEG) to investigate the spatiotemporal profile of repetition priming. We use intracranial electroencephalography (iEEG) to validate our fMRI/MEG measurements. Twelve controls completed a semantic judgment task with fMRI and MEG that included words presented once (new, 'N') and words that repeated (old, 'O'). Six patients with epilepsy completed the same task during iEEG recordings. Blood-oxygen level dependent (BOLD) responses for N vs. O words were examined across the cortical surface and within regions of interest. MEG waveforms for N vs. O words were estimated using a noise-normalized minimum norm solution, and used to interpret the timecourse of fMRI. Spatial concordance was observed between fMRI and MEG repetition effects from 350 to 450 ms within bilateral occipitotemporal and medial temporal, left prefrontal, and left posterior temporal cortex. Additionally, MEG revealed widespread sources within left temporoparietal regions, whereas fMRI revealed bilateral reductions in occipitotemporal and left superior frontal, and increases in inferior parietal, precuneus, and dorsolateral prefrontal activity. BOLD suppression in left posterior temporal, left inferior prefrontal, and right occipitotemporal cortex correlated with MEG repetition-related reductions. IEEG responses from all three regions supported the timecourse of MEG and localization of fMRI. Furthermore, iEEG decreases to repeated words were associated with decreased gamma power in several regions, providing evidence that gamma oscillations are tightly coupled to cognitive phenomena and reflect regional activations seen in the BOLD signal.

DynaSim: A MATLAB Toolbox for Neural Modeling and Simulation.

DynaSim is an open-source MATLAB/GNU Octave toolbox for rapid prototyping of neural models and batch simulation management. It is designed to speed up and simplify the process of generating, sharing, and exploring network models of neurons with one or more compartments. Models can be specified by equations directly (similar to XPP or the Brian simulator) or by lists of predefined or custom model components. The higher-level specification supports arbitrarily complex population models and networks of interconnected populations. DynaSim also includes a large set of features that simplify exploring model dynamics over parameter spaces, running simulations in parallel using both multicore processors and high-performance computer clusters, and analyzing and plotting large numbers of simulated data sets in parallel. It also includes a graphical user interface (DynaSim GUI) that supports full functionality without requiring user programming. The software has been implemented in MATLAB to enable advanced neural modeling using MATLAB, given its popularity and a growing interest in modeling neural systems. The design of DynaSim incorporates a novel schema for model specification to facilitate future interoperability with other specifications (e.g., NeuroML, SBML), simulators (e.g., NEURON, Brian, NEST), and web-based applications (e.g., Geppetto) outside MATLAB. DynaSim is freely available at http://dynasimtoolbox.org. This tool promises to reduce barriers for investigating dynamics in large neural models, facilitate collaborative modeling, and complement other tools being developed in the neuroinformatics community.

Hetereogeneity in Neuronal Intrinsic Properties: A Possible Mechanism for Hub-Like Properties of the Rat Anterior Cingulate Cortex during Network Activity.

The anterior cingulate cortex (ACC) is vital for a range of brain functions requiring cognitive control and has highly divergent inputs and outputs, thus manifesting as a hub in connectomic analyses. Studies show diverse functional interactions within the ACC are associated with network oscillations in the ß (20-30 Hz) and γ (30-80 Hz) frequency range. Oscillations permit dynamic routing of information within cortex, a function that depends on bandpass filter-like behavior to selectively respond to specific inputs. However, a putative hub region such as ACC needs to be able to combine inputs from multiple sources rather than select a single input at the expense of others. To address this potential functional dichotomy, we modeled local ACC network dynamics in the rat in vitro. Modal peak oscillation frequencies in the ß- and γ-frequency band corresponded to GABAAergic synaptic kinetics as seen in other regions; however, the intrinsic properties of ACC principal neurons were highly diverse. Computational modeling predicted that this neuronal response diversity broadened the bandwidth for filtering rhythmic inputs and supported combination-rather than selection-of different frequencies within the canonical γ and ß electroencephalograph bands. These findings suggest that oscillating neuronal populations can support either response selection (routing) or combination, depending on the interplay between the kinetics of synaptic inhibition and the degree of heterogeneity of principal cell intrinsic conductances.