Interacting rhythms enhance sensitivity of target detection in a fronto-parietal computational model of visual attention.

Even during sustained attention, enhanced processing of attended stimuli waxes and wanes rhythmically, with periods of enhanced and relatively diminished visual processing (and subsequent target detection) alternating at 4 or 8 Hz in a sustained visual attention task. These alternating attentional states occur alongside alternating dynamical states, in which lateral intraparietal cortex (LIP), the frontal eye field (FEF), and the mediodorsal pulvinar (mdPul) exhibit different activity and functional connectivity at α, ß, and γ frequencies-rhythms associated with visual processing, working memory, and motor suppression. To assess whether and how these multiple interacting rhythms contribute to periodicity in attention, we propose a detailed computational model of FEF and LIP. When driven by θ-rhythmic inputs simulating experimentally-observed mdPul activity, this model reproduced the rhythmic dynamics and behavioral consequences of observed attentional states, revealing that the frequencies and mechanisms of the observed rhythms allow for peak sensitivity in visual target detection while maintaining functional flexibility.

Gender bias in academia: A lifetime problem that needs solutions.

Despite increased awareness of the lack of gender equity in academia and a growing number of initiatives to address issues of diversity, change is slow, and inequalities remain. A major source of inequity is gender bias, which has a substantial negative impact on the careers, work-life balance, and mental health of underrepresented groups in science. Here, we argue that gender bias is not a single problem but manifests as a collection of distinct issues that impact researchers' lives. We disentangle these facets and propose concrete solutions that can be adopted by individuals, academic institutions, and society.

Differential contributions of synaptic and intrinsic inhibitory currents to speech segmentation via flexible phase-locking in neural oscillators.

Current hypotheses suggest that speech segmentation-the initial division and grouping of the speech stream into candidate phrases, syllables, and phonemes for further linguistic processing-is executed by a hierarchy of oscillators in auditory cortex. Theta (∼3-12 Hz) rhythms play a key role by phase-locking to recurring acoustic features marking syllable boundaries. Reliable synchronization to quasi-rhythmic inputs, whose variable frequency can dip below cortical theta frequencies (down to ∼1 Hz), requires "flexible" theta oscillators whose underlying neuronal mechanisms remain unknown. Using biophysical computational models, we found that the flexibility of phase-locking in neural oscillators depended on the types of hyperpolarizing currents that paced them. Simulated cortical theta oscillators flexibly phase-locked to slow inputs when these inputs caused both (i) spiking and (ii) the subsequent buildup of outward current sufficient to delay further spiking until the next input. The greatest flexibility in phase-locking arose from a synergistic interaction between intrinsic currents that was not replicated by synaptic currents at similar timescales. Flexibility in phase-locking enabled improved entrainment to speech input, optimal at mid-vocalic channels, which in turn supported syllabic-timescale segmentation through identification of vocalic nuclei. Our results suggest that synaptic and intrinsic inhibition contribute to frequency-restricted and -flexible phase-locking in neural oscillators, respectively. Their differential deployment may enable neural oscillators to play diverse roles, from reliable internal clocking to adaptive segmentation of quasi-regular sensory inputs like speech.

Interactions of multiple rhythms in a biophysical network of neurons.

Neural oscillations, including rhythms in the beta1 band (12-20 Hz), are important in various cognitive functions. Often neural networks receive rhythmic input at frequencies different from their natural frequency, but very little is known about how such input affects the network's behavior. We use a simplified, yet biophysical, model of a beta1 rhythm that occurs in the parietal cortex, in order to study its response to oscillatory inputs. We demonstrate that a cell has the ability to respond at the same time to two periodic stimuli of unrelated frequencies, firing in phase with one, but with a mean firing rate equal to that of the other. We show that this is a very general phenomenon, independent of the model used. We next show numerically that the behavior of a different cell, which is modeled as a high-dimensional dynamical system, can be described in a surprisingly simple way, owing to a reset that occurs in the state space when the cell fires. The interaction of the two cells leads to novel combinations of properties for neural dynamics, such as mode-locking to an input without phase-locking to it.

Prefrontal oscillations modulate the propagation of neuronal activity required for working memory.

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.

A biophysical model of striatal microcircuits suggests gamma and beta oscillations interleaved at delta/theta frequencies mediate periodicity in motor control.

Striatal oscillatory activity is associated with movement, reward, and decision-making, and observed in several interacting frequency bands. Local field potential recordings in rodent striatum show dopamine- and reward-dependent transitions between two states: a "spontaneous" state involving ß (∼15-30 Hz) and low γ (∼40-60 Hz), and a state involving θ (∼4-8 Hz) and high γ (∼60-100 Hz) in response to dopaminergic agonism and reward. The mechanisms underlying these rhythmic dynamics, their interactions, and their functional consequences are not well understood. In this paper, we propose a biophysical model of striatal microcircuits that comprehensively describes the generation and interaction of these rhythms, as well as their modulation by dopamine. Building on previous modeling and experimental work suggesting that striatal projection neurons (SPNs) are capable of generating ß oscillations, we show that networks of striatal fast-spiking interneurons (FSIs) are capable of generating δ/θ (ie, 2 to 6 Hz) and γ rhythms. Under simulated low dopaminergic tone our model FSI network produces low γ band oscillations, while under high dopaminergic tone the FSI network produces high γ band activity nested within a δ/θ oscillation. SPN networks produce ß rhythms in both conditions, but under high dopaminergic tone, this ß oscillation is interrupted by δ/θ-periodic bursts of γ-frequency FSI inhibition. Thus, in the high dopamine state, packets of FSI γ and SPN ß alternate at a δ/θ timescale. In addition to a mechanistic explanation for previously observed rhythmic interactions and transitions, our model suggests a hypothesis as to how the relationship between dopamine and rhythmicity impacts motor function. We hypothesize that high dopamine-induced periodic FSI γ-rhythmic inhibition enables switching between ß-rhythmic SPN cell assemblies representing the currently active motor program, and thus that dopamine facilitates movement in part by allowing for rapid, periodic shifts in motor program execution.

Parietal low beta rhythm provides a dynamical substrate for a working memory buffer.

Working memory (WM) is a component of the brain's memory systems vital for interpretation of sequential sensory inputs and consequent decision making. Anatomically, WM is highly distributed over the prefrontal cortex (PFC) and the parietal cortex (PC). Here we present a biophysically detailed dynamical systems model for a WM buffer situated in the PC, making use of dynamical properties believed to be unique to this area. We show that the natural beta1 rhythm (12 to 20 Hz) of the PC provides a substrate for an episodic buffer that can synergistically combine executive commands (e.g., from PFC) and multimodal information into a flexible and updatable representation of recent sensory inputs. This representation is sensitive to distractors, it allows for a readout mechanism, and it can be readily terminated by executive input. The model provides a demonstration of how information can be usefully stored in the temporal patterns of activity in a neuronal network rather than just synaptic weights between the neurons in that network.

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.

Striatal cholinergic receptor activation causes a rapid, selective and state-dependent rise in cortico-striatal ß activity.

Cortico-basal ganglia-thalamic (CBT) ß oscillations (15-30 Hz) are elevated in Parkinson's disease and correlated with movement disability. To date, no experimental paradigm outside of loss of dopamine has been able to specifically elevate ß oscillations in the CBT loop. Here, we show that activation of striatal cholinergic receptors selectively increased ß oscillations in mouse striatum and motor cortex. In individuals showing simultaneous ß increases in both striatum and M1, ß partial directed coherence (PDC) increased from striatum to M1 (but not in the reverse direction). In individuals that did not show simultaneous ß increases, ß PDC increased from M1 to striatum (but not in the reverse direction), and M1 was characterized by persistent ß-high frequency oscillation phase-amplitude coupling. Finally, the direction of ß PDC distinguished between ß sub-bands. This suggests that (1) striatal cholinergic tone exerts state-dependent and frequency-selective control over CBT ß power and coordination; (2) ongoing rhythmic dynamics can determine whether elevated ß oscillations are expressed in striatum and M1; and (3) altered striatal cholinergic tone differentially modulates distinct ß sub-bands.

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.

Low-Beta Oscillations Turn Up the Gain During Category Judgments.

Synchrony between local field potential (LFP) rhythms is thought to boost the signal of attended sensory inputs. Other cognitive functions could benefit from such gain control. One is categorization where decisions can be difficult if categories differ in subtle ways. Monkeys were trained to flexibly categorize smoothly varying morphed stimuli, using orthogonal boundaries to carve up the same stimulus space in 2 different ways. We found evidence for category-specific patterns of low-beta (16-20 Hz) synchrony in the lateral prefrontal cortex (PFC). This synchrony was stronger when a given category scheme was relevant. We also observed an overall increase in low-beta LFP synchrony for stimuli that were near the category boundary and thus more difficult to categorize. Beta category selectivity was evident in partial field-field coherence measurements, which measure local synchrony, but the boundary enhancement was not. Thus, it seemed that category selectivity relied on local interactions while boundary enhancement was a more global effect. The results suggest that beta synchrony helps form category ensembles and may reflect recruitment of additional cortical resources for categorizing challenging stimuli, thus serving as a form of gain control.

Mesoscale-duration activated states gate spiking in response to fast rises in membrane voltage in the awake brain.

Seconds-scale network states, affecting many neurons within a network, modulate neural activity by complementing fast integration of neuron-specific inputs that arrive in the milliseconds before spiking. Nonrhythmic subthreshold dynamics at intermediate timescales, however, are less well characterized. We found, using automated whole cell patch clamping in vivo, that spikes recorded in CA1 and barrel cortex in awake mice are often preceded not only by monotonic voltage rises lasting milliseconds but also by more gradual (lasting tens to hundreds of milliseconds) depolarizations. The latter exert a gating function on spiking, in a fashion that depends on the gradual rise duration: the probability of spiking was higher for longer gradual rises, even when controlled for the amplitude of the gradual rises. Barrel cortex double-autopatch recordings show that gradual rises are shared across some, but not all, neurons. The gradual rises may represent a new kind of state, intermediate both in timescale and in proportion of neurons participating, which gates a neuron's ability to respond to subsequent inputs.NEW & NOTEWORTHY We analyzed subthreshold activity preceding spikes in hippocampus and barrel cortex of awake mice. Aperiodic voltage ramps extending over tens to hundreds of milliseconds consistently precede and facilitate spikes, in a manner dependent on both their amplitude and their duration. These voltage ramps represent a "mesoscale" activated state that gates spike production in vivo.

Aberrant Network Activity in Schizophrenia.

Brain dynamic changes associated with schizophrenia are largely equivocal, with interpretation complicated by many factors, such as the presence of therapeutic agents and the complex nature of the syndrome itself. Evidence for a brain-wide change in individual network oscillations, shared by all patients, is largely equivocal, but stronger for lower (delta) than for higher (gamma) bands. However, region-specific changes in rhythms across multiple, interdependent, nested frequencies may correlate better with pathology. Changes in synaptic excitation and inhibition in schizophrenia disrupt delta rhythm-mediated cortico-cortical communication, while enhancing thalamocortical communication in this frequency band. The contrasting relationships between delta and higher frequencies in thalamus and cortex generate frequency mismatches in inter-regional connectivity, leading to a disruption in temporal communication between higher-order brain regions associated with mental time travel.

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.

Frontal beta-theta network during REM sleep.

We lack detailed knowledge about the spatio-temporal physiological signatures of REM sleep, especially in humans. By analyzing intracranial electrode data from humans, we demonstrate for the first time that there are prominent beta (15-35 Hz) and theta (4-8 Hz) oscillations in both the anterior cingulate cortex (ACC) and the DLPFC during REM sleep. We further show that these theta and beta activities in the ACC and the DLPFC, two relatively distant but reciprocally connected regions, are coherent. These findings suggest that, counter to current prevailing thought, the DLPFC is active during REM sleep and likely interacting with other areas. Since the DLPFC and the ACC are implicated in memory and emotional regulation, and the ACC has motor areas and is thought to be important for error detection, the dialogue between these two areas could play a role in the regulation of emotions and in procedural motor and emotional memory consolidation.

GABAB receptor-mediated, layer-specific synaptic plasticity reorganizes gamma-frequency neocortical response to stimulation.

Repeated presentations of sensory stimuli generate transient gamma-frequency (30-80 Hz) responses in neocortex that show plasticity in a task-dependent manner. Complex relationships between individual neuronal outputs and the mean, local field potential (population activity) accompany these changes, but little is known about the underlying mechanisms responsible. Here we show that transient stimulation of input layer 4 sufficient to generate gamma oscillations induced two different, lamina-specific plastic processes that correlated with lamina-specific changes in responses to further, repeated stimulation: Unit rates and recruitment showed overall enhancement in supragranular layers and suppression in infragranular layers associated with excitatory or inhibitory synaptic potentiation onto principal cells, respectively. Both synaptic processes were critically dependent on activation of GABAB receptors and, together, appeared to temporally segregate the cortical representation. These data suggest that adaptation to repetitive sensory input dramatically alters the spatiotemporal properties of the neocortical response in a manner that may both refine and minimize cortical output simultaneously.

Heterogeneous neural amplifier integration for scalable extracellular microelectrodes.

We here demonstrate multi-chip heterogeneous integration of microfabricated extracellular recording electrodes with neural amplifiers, highlighting a path to scaling electrode channel counts without the need for more complex monolithic integration. We characterize the noise and impedance performance of the heterogeneously integrated neural recording electrodes, and analyze the design parameters that enable the low-voltage neural input signals to co-exist with the high-frequency and high-voltage digital outputs on the same silicon substrate. This heterogeneous integration approach can enable future scaling efforts for microfabricated neural probes, and provides a design path for modular, fast, and independent scaling innovations in recording electrodes and neural amplifiers.

Potential Mechanisms Underlying Intercortical Signal Regulation via Cholinergic Neuromodulators.

The dynamical behavior of the cortex is extremely complex, with different areas and even different layers of a cortical column displaying different temporal patterns. A major open question is how the signals from different layers and different brain regions are coordinated in a flexible manner to support function. Here, we considered interactions between primary auditory cortex and adjacent association cortex. Using a biophysically based model, we show how top-down signals in the beta and gamma regimes can interact with a bottom-up gamma rhythm to provide regulation of signals between the cortical areas and among layers. The flow of signals depends on cholinergic modulation: with only glutamatergic drive, we show that top-down gamma rhythms may block sensory signals. In the presence of cholinergic drive, top-down beta rhythms can lift this blockade and allow signals to flow reciprocally between primary sensory and parietal cortex. SIGNIFICANCE STATEMENT: Flexible coordination of multiple cortical areas is critical for complex cognitive functions, but how this is accomplished is not understood. Using computational models, we studied the interactions between primary auditory cortex (A1) and association cortex (Par2). Our model is capable of replicating interaction patterns observed in vitro and the simulations predict that the coordination between top-down gamma and beta rhythms is central to the gating process regulating bottom-up sensory signaling projected from A1 to Par2 and that cholinergic modulation allows this coordination to occur.

A direct-to-drive neural data acquisition system.

Driven by the increasing channel count of neural probes, there is much effort being directed to creating increasingly scalable electrophysiology data acquisition (DAQ) systems. However, all such systems still rely on personal computers for data storage, and thus are limited by the bandwidth and cost of the computers, especially as the scale of recording increases. Here we present a novel architecture in which a digital processor receives data from an analog-to-digital converter, and writes that data directly to hard drives, without the need for a personal computer to serve as an intermediary in the DAQ process. This minimalist architecture may support exceptionally high data throughput, without incurring costs to support unnecessary hardware and overhead associated with personal computers, thus facilitating scaling of electrophysiological recording in the future.

Thalamic mechanisms underlying alpha-delta sleep with implications for fibromyalgia.

Alpha-delta sleep is the abnormal intrusion of alpha activity (8- to 13-Hz oscillations) into the delta activity (1- to 4-Hz oscillations) that defines slow-wave sleep. Alpha-delta sleep is especially prevalent in fibromyalgia patients, and there is evidence suggesting that the irregularities in the sleep of these patients may cause the muscle and tissue pain that characterizes the disorder. We constructed a biophysically realistic mathematical model of alpha-delta sleep. Imaging studies in fibromyalgia patients suggesting altered levels of activity in the thalamus motivated a thalamic model as the source of alpha activity. Since sodium oxybate helps to alleviate the symptoms of fibromyalgia and reduces the amount of alpha-delta sleep in fibromyalgia patients, we examined how changes in the molecular targets of sodium oxybate affected alpha-delta activity in our circuit. Our model shows how alterations in GABAB currents and two thalamic currents, Ih (a hyperpolarization-activated current) and a potassium leak current, transform a circuit that normally produces delta oscillations into one that produces alpha-delta activity. Our findings suggest that drugs that reduce Ih conductances and/or increase potassium conductances, without necessarily increasing GABAB conductances, might be sufficient to restore delta sleep. Furthermore, they suggest that delta sleep might be restored by drugs that preferentially target these currents in the thalamus; such drugs might have fewer side effects than drugs that act systemically.