Search

[EN] Cognition involves using attended information, maintained in working memory (WM), to guide action. During a cognitive task, a correct response requires flexible, selective gating so that only the appropriate information flows from WM to downstream effectors that carry out the response. In this work, we used biophysically-detailed modeling to explore the hypothesis that network oscillations in prefrontal cortex (PFC), leveraging local inhibition, can independently gate responses to items in WM. The key role of local inhibition was to control the period between spike bursts in the outputs, and to produce an oscillatory response no matter whether the WM item was maintained in an asynchronous or oscillatory state. We found that the WM item that induced an oscillatory population response in the PFC output layer with the shortest period between spike bursts was most reliably propagated. The network resonant frequency (i.e., the input frequency that produces the largest response) of the output layer can be flexibly tuned by varying the excitability of deep layer principal cells. Our model suggests that experimentally-observed modulation of PFC beta-frequency (15-30 Hz) and gamma -frequency (30-80 Hz) oscillations could leverage network resonance and local inhibition to govern the flexible routing of signals in service to cognitive processes like gating outputs from working memory and the selection of rule-based actions. Importantly, we show for the first time that nonspecific changes in deep layer excitability can tune the output gate's resonant frequency, enabling the specific selection of signals encoded by populations in asynchronous or fast oscillatory states. More generally, this represents a dynamic mechanism by which adjusting network excitability can govern the propagation of asynchronous and oscillatory signals throughout neocortex., This work was supported by the U.S. Army Research Office under award number ARO W911NF-12-R-0012-02 to N. K., the U.S. Office of Naval Research under award number ONR MURI N00014-16-1-2832 to M. H. and E. M., the National Institute of Mental Health under award number NIMH R37MH087027 to E. M., and The MIT Picower Institute Faculty Innovation Fund to E. M. We would like to acknowledge Joachim Hass and Michelle McCarthy for early discussions of our modeling results, as well as Andre Bastos and Mikael Lundqvist for discussions relating our modeling work to their experiments.

Persistent activity, the maintenance of neural activation over short periods of time in cortical networks, is widely thought to underlie the cognitive function of working memory. A large body of modeling studies has reproduced this kind of activity using cell assemblies with strengthened synaptic connections. However, almost all of these studies have considered persistent activity within networks with homogeneous neurons and synapses, making it difficult to judge the validity of such model results for cortical dynamics, which is based on highly heterogeneous neurons. Here, we consider persistent activity in a detailed, strongly data-driven network model of the prefrontal cortex with heterogeneous neuron and synapse parameters. Surprisingly, persistent activity could not be reproduced in this model without incorporating further constraints. We identified three factors that prevent successful persistent activity: heterogeneity in the cell parameters of interneurons, heterogeneity in the parameters of short-term synaptic plasticity and heterogeneity in the synaptic weights. Our model predicts that persistent activity is recovered if the heterogeneity in the activity of individual interneurons is diminished, which could be achieved by a homeostatic plasticity mechanism. Such a plasticity scheme could also compensate the heterogeneities in the synaptic weights and short-term plasticity when applied to the inhibitory synapses. Cell assemblies shaped in this way may be potentially targeted by distinct inputs or become more responsive to specific tuning or spectral properties. Furthermore, the model predicts that a network that exhibits persistent activity is not able to dynamically produce intrinsic in vivo-like irregular activity at the same time, because heterogeneous synaptic connections are required for these dynamics. Thus, the background noise in such a network must either be produced by external input or constitutes an entirely different state of the network, which is brought about, e.g., by neuromodulation.Author summaryTo operate effectively in a constantly changing world, it is crucial to keep relevant information in mind for short periods of time. This ability, called working memory, is commonly assumed to rest on reverberating activity among members of cell assemblies. While effective in reproducing key results of working memory, most cell assembly models rest on major simplifications such as using the same parameters for all neurons and synapses, i.e., assuming homogeneity in these parameters. Here, we show that this homogeneity assumption is necessary for persistent activity to arise, specifically for inhibitory interneurons and synapses. Using a strongly data-driven network model of the prefrontal cortex, we show that the heterogeneities in the above parameters that are implied by in vitro studies prevent persistent activity. When homogeneity is imposed on inhibitory neurons and synapses, persistent activity is recovered. We propose that the homogeneity constraints can be implemented in the brain by means of homeostatic plasticity, a form of learning that keeps the activity of a network in a constant, homeostatic state. The model makes a number of predictions for biological networks, including a structural separation of networks responsible for generating persistent activity and spontaneous, noise-like activity.

[EN] 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, We thank T. Womelsdorf for helpful suggestions on an earlier version of the manuscript. We also thank the two reviewers for the constructive comments that enhanced the quality of the manuscript. In particular, their question regarding the resonant properties of SPNs under distinct mean input helped us to uncover how the resonance of D2 SPNs shifts in frequency space (Fig. 3E). Our research was supported by the Army Research Office (ARO) Grant W911NF-12-R-0012-02 (to N.K.). Additionally, S.A. and N.K. were supported by NSF Grant DMS-1042134, and M.M.M. was supported by the Collaborative Research in Computational Neuroscience (CRCNS) NIH Grant CRCNS 1R01N5081716

[EN] 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., This material is based upon research supported by the U. S. Army Research Office under award number ARO W911NF-12-R-0012-02 to N. K., the U. S. Office of Naval Research under award number ONR MURI N00014-16-1-2832 to M. H., and the National Science Foundation under award number NSF DMS-1042134 (Cognitive Rhythms Collaborative: A Discovery Network) to N. K. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Classical accounts of biased competition (BC) 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 medium spiny neurons (MSNs) receive balanced cortical input, yet their activity determines the choice between GO and NO-GO pathways in the basal ganglia. We present a corticostriatal model identifying three candidate mechanisms that rely on physiological asymmetries to effect rate- and time-coded BC in the presence of balanced inputs. First, tonic input strength determines which MSN phenotype exhibit higher mean firing rate (FR). Second, low strength oscillatory inputs induce higher FR in D2 MSNs but higher coherence between D1 MSNs. Third, high strength inputs oscillating at distinct frequencies preferentially activate D1 or D2 MSN populations. Of these candidate mechanisms, only the latter accommodates observed rhythmic activity supporting rule-based decision making in prefrontal cortex.

[EN] 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., This material is based upon research supported by the U.S. Army Research Office under award number ARO W911NF-12-R-0012-02, the U.S. Office of Naval Research under award number ONR MURI N00014-16-1-2832, and the National Science Foundation under award number NSF DMS-1042134 (Cognitive Rhythms Collaborative: A Discovery Network)

The anesthetic propofol elicits many different spectral properties on the EEG, including alpha oscillations (8–12 Hz), Slow Wave Oscillations (SWO, 0.1–1.5 Hz), and dose-dependent phase-amplitude coupling (PAC) between alpha and SWO. Propofol is known to increase GABAA inhibition and decrease H-current strength, but how it generates these rhythms and their interactions is still unknown. To investigate both generation of the alpha rhythm and its PAC to SWO, we simulate a Hodgkin-Huxley network model of a hyperpolarized thalamus and corticothalamic inputs. We find, for the first time, that the model thalamic network is capable of independently generating the sustained alpha seen in propofol, which may then be relayed to cortex and expressed on the EEG. This dose-dependent sustained alpha critically relies on propofol GABAA potentiation to alter the intrinsic spindling mechanisms of the thalamus. Furthermore, the H-current conductance and background excitation of these thalamic cells must be within specific ranges to exhibit any intrinsic oscillations, including sustained alpha. We also find that, under corticothalamic SWO UP and DOWN states, thalamocortical output can exhibit maximum alpha power at either the peak or trough of this SWO; this implies the thalamus may be the source of propofol-induced PAC. Hyperpolarization level is the main determinant of whether the thalamus exhibits trough-max PAC, which is associated with lower propofol dose, or peak-max PAC, associated with higher dose. These findings suggest: the thalamus generates a novel rhythm under GABAA potentiation such as under propofol, its hyperpolarization may determine whether a patient experiences trough-max or peak-max PAC, and the thalamus is a critical component of propofol-induced cortical spectral phenomena. Changes to the thalamus may be a critical part of how propofol accomplishes its effects, including unconsciousness., Author summary Anesthetics make patients lose consciousness, but how they affect brain dynamics is still unknown. Changes in EEG brainwaves are some of the few noninvasive signals we can use to learn about this. By analyzing such data, we can develop more targeted anesthetics, expand our knowledge of sleep circuits, and better understand how diseases impact these systems. The anesthetic propofol is known, among other effects, to increase synaptic inhibition, but it is unclear how these changes induce EEG alpha (8–12 Hz) oscillations and their interaction with slow wave (0.1–1.5 Hz) oscillations; these signals have been correlated with the state of propofol-infused consciousness. We simulated a network of thalamic cells to understand the mechanisms generating these signals. Propofol-potentiated inhibition produced a novel, sustained alpha rhythm in our network. Changes to the tonic level of depolarization enabled the alpha oscillations to occur at different phases in the slow wave oscillation, as seen clinically with increasing propofol dose. The thalamus may be critical to propofol-induced alpha oscillations and their coupling to slow wave oscillations. By understanding the mechanisms generating alpha, we may be able to design experiments to dissociate alpha from slow waves and determine their independent effects on levels of consciousness.

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.

Subjects detected rarely-occurring shifts between two simple tone-patterns, in a paradigm that dissociated the effects of rarity from those of pitch, habituation, and attention. Whole-head magnetoencephalography suggested that rare attended pattern-shifts evoked activity first in the superior temporal plane (STP, peak ~100 ms), then superior temporal sulcus (STS, peak ~130 ms), then posteroventral prefrontal (pvpF, peak ~230 ms), and anterior temporal cortices (aT, peak ~370 ms). Activity was more prominent in the right hemisphere. After subtracting the effects of non-shift tones (balanced for pitch and habituation status), weak but consistent differential effects of pattern-shifts began in aT at 90–130 ms, spread to STS and STP at ~130 ms, then pvpF, and finally returned to aT. Cingulate activity resembled prefrontal. Responses to pattern shifts were greatly attenuated when the same stimuli were ignored, suggesting that the initial superior temporal activity reflected an attention-related Mismatch Negativity. The prefrontal activity at ~230 ms corresponded in latency and task correlates with simultaneously recorded event-related potential components N2b and P3a; the subsequent temporal activity corresponded to the P3b. These results were confirmed in sensors specific for frontal or temporal cortex, and thus are independent of the inverse method used. Overall, these results suggest that auditory working memory for temporal patterns begins with detection of the pattern change by an interaction of anterior and superior temporal structures, followed by identification of the event and its consequences led by posteroventral prefrontal and cingulate cortices, and finally, definitive encoding of the event in anterior temporal areas.

Prior neuroimaging evidence indicates that decision conflict activates medial and lateral prefrontal and parietal cortices. Theoretical accounts of cognitive control highlight anterior cingulate cortex (ACC) as a central node in this network. However, a better understanding of the relative primacy and functional contributions of these areas to decision conflict requires insight into the neural dynamics of successive processing stages including conflict detection, response selection and execution. Moderate alcohol intoxication impairs cognitive control as it interferes with the ability to inhibit dominant, prepotent responses when they are no longer correct. To examine the effects of moderate intoxication on successive processing stages during cognitive control, spatio-temporal changes in total event-related theta power were measured during Stroop-induced conflict. Healthy social drinkers served as their own controls by participating in both alcohol (0.6 g/kg ethanol for men, 0.55 g/kg women) and placebo conditions in a counterbalanced design. Anatomically-constrained magnetoencephalography (aMEG) approach was applied to complex power spectra for theta (4-7 Hz) frequencies. The principal generator of event-related theta power to conflict was estimated to ACC, with contributions from fronto-parietal areas. The ACC was uniquely sensitive to conflict during both early conflict detection, and later response selection and execution stages. Alcohol attenuated theta power to conflict across successive processing stages, suggesting that alcohol-induced deficits in cognitive control may result from theta suppression in the executive network. Slower RTs were associated with attenuated theta power estimated to ACC, indicating that alcohol impairs motor preparation and execution subserved by the ACC. In addition to their relevance for the currently prevailing accounts of cognitive control, our results suggest that alcohol-induced impairment of top-down strategic processing underlies poor self-control and inability to refrain from drinking.

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.