Flexible resonance in prefrontal networks with strong feedback inhibition

Functional Dissection of Ipsilateral and Contralateral Neural Activity Propagation Using Voltage-Sensitive Dye Imaging in Mouse Prefrontal Cortex

Prefrontal cortex (PFC) intrahemispheric activity and the interhemispheric connection have a significant impact on neuropsychiatric disorder pathology. This study aimed to generate a functional map of FC intrahemispheric and interhemispheric connections. Functional dissection of mouse PFCs was performed using the voltage-sensitive dye (VSD) imaging method with high speed (1 ms/frame), high resolution (256 × 256 pixels), and a large field of view (∼10 mm). Acute serial 350 μm slices were prepared from the bregma covering the PFC and numbered 1-5 based on their distance from the bregma (i.e., 1.70, 1.34, 0.98, 0.62, and 0.26 mm) with reference to the Mouse Brain Atlas (Paxinos and Franklin, 2008). The neural response to electrical stimulation was measured at nine sites and then averaged, and a functional map of the propagation patterns was created. Intracortical propagation was observed in slices 3-5, encompassing the anterior cingulate cortex (ACC) and corpus callosum (CC). The activity reached area 33 of the ACC. Direct white matter stimulation activated area 33 in both hemispheres. Similar findings were obtained via DiI staining of the CC. Imaging analysis revealed directional biases in neural signals traveling within the ACC, whereby the signal transmission speed and probability varied based on the signal direction. Specifically, the spread of neural signals from cg2 to cg1 was stronger than that from cingulate cortex area 1(cg1) to cingulate cortex area 2(cg2), which has implications for interhemispheric functional connections. These findings highlight the importance of understanding the PFC functional anatomy in evaluating neuromodulators like serotonin and dopamine, as well as other factors related to neuropsychiatric diseases.

Uncovering the organization of neural circuits with Generalized Phase Locking Analysis

Despite the considerable progress of in vivo neural recording techniques, inferring the biophysical mechanisms underlying large scale coordination of brain activity from neural data remains challenging. One obstacle is the difficulty to link high dimensional functional connectivity measures to mechanistic models of network activity. We address this issue by investigating spike-field coupling (SFC) measurements, which quantify the synchronization between, on the one hand, the action potentials produced by neurons, and on the other hand mesoscopic "field" signals, reflecting subthreshold activities at possibly multiple recording sites. As the number of recording sites gets large, the amount of pairwise SFC measurements becomes overwhelmingly challenging to interpret. We develop Generalized Phase Locking Analysis (GPLA) as an interpretable dimensionality reduction of this multivariate SFC. GPLA describes the dominant coupling between field activity and neural ensembles across space and frequencies. We show that GPLA features are biophysically interpretable when used in conjunction with appropriate network models, such that we can identify the influence of underlying circuit properties on these features. We demonstrate the statistical benefits and interpretability of this approach in various computational models and Utah array recordings. The results suggest that GPLA, used jointly with biophysical modeling, can help uncover the contribution of recurrent microcircuits to the spatio-temporal dynamics observed in multi-channel experimental recordings.

Emotional Memory Processing during REM Sleep with Implications for Post-Traumatic Stress Disorder

REM sleep is important for the processing of emotional memories, including fear memories. Rhythmic interactions, especially in the theta band, between the medial prefrontal cortex (mPFC) and limbic structures are thought to play an important role, but the ways in which memory processing occurs at a mechanistic and circuits level are largely unknown. To investigate how rhythmic interactions lead to fear extinction during REM sleep, we used a biophysically based model that included the infralimbic cortex (IL), a part of the mPFC with a critical role in suppressing fear memories. Theta frequency (4-12 Hz) inputs to a given cell assembly in IL, representing an emotional memory, resulted in the strengthening of connections from the IL to the amygdala and the weakening of connections from the amygdala to the IL, resulting in the suppression of the activity of fear expression cells for the associated memory. Lower frequency (4 Hz) theta inputs effected these changes over a wider range of input strengths. In contrast, inputs at other frequencies were ineffective at causing these synaptic changes and did not suppress fear memories. Under post-traumatic stress disorder (PTSD) REM sleep conditions, rhythmic activity dissipated, and 4 Hz theta inputs to IL were ineffective, but higher-frequency (10 Hz) theta inputs to IL induced changes similar to those seen with 4 Hz inputs under normal REM sleep conditions, resulting in the suppression of fear expression cells. These results suggest why PTSD patients may repeatedly experience the same emotionally charged dreams and suggest potential neuromodulatory therapies for the amelioration of PTSD symptoms.SIGNIFICANCE STATEMENT Rhythmic interactions in the theta band between the mPFC and limbic structures are thought to play an important role in processing emotional memories, including fear memories, during REM sleep. The infralimbic cortex (IL) in the mPFC is thought to play a critical role in suppressing fear memories. We show that theta inputs to the IL, unlike other frequency inputs, are effective in producing synaptic changes that suppress the activity of fear expression cells associated with a given memory. Under PTSD REM sleep conditions, lower-frequency (4 Hz) theta inputs to the IL do not suppress the activity of fear expression cells associated with the given memory but, surprisingly, 10 Hz inputs do. These results suggest potential neuromodulatory therapies for PTSD.

Network resonance can be generated independently at distinct levels of neuronal organization

Resonance is defined as maximal response of a system to periodic inputs in a limited frequency band. Resonance may serve to optimize inter-neuronal communication, and has been observed at multiple levels of neuronal organization. However, it is unknown how neuronal resonance observed at the network level is generated and how network resonance depends on the properties of the network building blocks. Here, we first develop a metric for quantifying spike timing resonance in the presence of background noise, extending the notion of spiking resonance for in vivo experiments. Using conductance-based models, we find that network resonance can be inherited from resonances at other levels of organization, or be intrinsically generated by combining mechanisms across distinct levels. Resonance of membrane potential fluctuations, postsynaptic potentials, and single neuron spiking can each be generated independently of resonance at any other level and be propagated to the network level. At all levels of organization, interactions between processes that give rise to low- and high-pass filters generate the observed resonance. Intrinsic network resonance can be generated by the combination of filters belonging to different levels of organization. Inhibition-induced network resonance can emerge by inheritance from resonance of membrane potential fluctuations, and be sharpened by presynaptic high-pass filtering. Our results demonstrate a multiplicity of qualitatively different mechanisms that can generate resonance in neuronal systems, and provide analysis tools and a conceptual framework for the mechanistic investigation of network resonance in terms of circuit components, across levels of neuronal organization.

How one part of the brain responds to periodic input from another part depends on resonant circuit properties. Resonance is a basic property of physical systems, and has been experimentally observed at various levels of neuronal organization both in vitro and in vivo. Yet how resonance is generated in neuronal networks is largely unknown. In particular, whether resonance can be generated directly at the level of a network of spiking neurons remains to be determined. Using detailed biophysical modeling, we develop a conceptual framework according to which resonance at a given level of organization is generated by the interplay of low- and high-pass filters, implemented at either the same or across levels of neuronal organization. We tease apart representative, biophysically-plausible generative mechanisms of resonance at four different levels of organization: membrane potential fluctuations, single neuron spiking, synaptic transmission, and neuronal networks. We identify conditions under which resonance at one level can be inherited to another level of organization, provide conclusive evidence that resonance at each level can be generated without resonance at any other level, and describe a number of representative routes to network resonance. The proposed framework facilitates the investigation of resonance in neuronal systems.

Layer-specific pyramidal neuron properties underlie diverse anterior cingulate cortical motor and limbic networks

The laminar cellular and circuit mechanisms by which the anterior cingulate cortex (ACC) exerts flexible control of motor and affective information for goal-directed behavior have not been elucidated. Using multimodal tract-tracing, in vitro patch-clamp recording and computational approaches in rhesus monkeys (M. mulatta), we provide evidence that specialized motor and affective network dynamics can be conferred by layer-specific biophysical and structural properties of ACC pyramidal neurons targeting two key downstream structures -the dorsal premotor cortex (PMd) and the amygdala (AMY). AMY-targeting neurons exhibited significant laminar differences, with L5 more excitable (higher input resistance and action potential firing rates) than L3 neurons. Between-pathway differences were found within L5, with AMY-targeting neurons exhibiting greater excitability, apical dendritic complexity, spine densities, and diversity of inhibitory inputs than PMd-targeting neurons. Simulations using a pyramidal-interneuron network model predict that these layer- and pathway-specific single-cell differences contribute to distinct network oscillatory dynamics. L5 AMY-targeting networks are more tuned to slow oscillations well-suited for affective and contextual processing timescales, while PMd-targeting networks showed strong beta/gamma synchrony implicated in rapid sensorimotor processing. These findings are fundamental to our broad understanding of how layer-specific cellular and circuit properties can drive diverse laminar activity found in flexible behavior.

Interneuron-specific gamma synchronization indexes cue uncertainty and prediction errors in lateral prefrontal and anterior cingulate cortex

Inhibitory interneurons are believed to realize critical gating functions in cortical circuits, but it has been difficult to ascertain the content of gated information for well-characterized interneurons in primate cortex. Here, we address this question by characterizing putative interneurons in primate prefrontal and anterior cingulate cortex while monkeys engaged in attention demanding reversal learning. We find that subclasses of narrow spiking neurons have a relative suppressive effect on the local circuit indicating they are inhibitory interneurons. One of these interneuron subclasses showed prominent firing rate modulations and (35–45 Hz) gamma synchronous spiking during periods of uncertainty in both, lateral prefrontal cortex (LPFC) and anterior cingulate cortex (ACC). In LPFC, this interneuron subclass activated when the uncertainty of attention cues was resolved during flexible learning, whereas in ACC it fired and gamma-synchronized when outcomes were uncertain and prediction errors were high during learning. Computational modeling of this interneuron-specific gamma band activity in simple circuit motifs suggests it could reflect a soft winner-take-all gating of information having high degree of uncertainty. Together, these findings elucidate an electrophysiologically characterized interneuron subclass in the primate, that forms gamma synchronous networks in two different areas when resolving uncertainty during adaptive goal-directed behavior.

Cortical gamma-band resonance preferentially transmits coherent input

Synchronization has been implicated in neuronal communication, but causal evidence remains indirect. We use optogenetics to generate depolarizing currents in pyramidal neurons of the cat visual cortex, emulating excitatory synaptic inputs under precise temporal control, while measuring spike output. The cortex transforms constant excitation into strong gamma-band synchronization, revealing the well-known cortical resonance. Increasing excitation with ramps increases the strength and frequency of synchronization. Slow, symmetric excitation profiles reveal hysteresis of power and frequency. White-noise input sequences enable causal analysis of network transmission, establishing that the cortical gamma-band resonance preferentially transmits coherent input components. Models composed of recurrently coupled excitatory and inhibitory units uncover a crucial role of feedback inhibition and suggest that hysteresis can arise through spike-frequency adaptation. The presented approach provides a powerful means to investigate the resonance properties of local circuits and probe how these properties transform input and shape transmission.

Rhythmic synchronization has been implicated in neuronal communication, yet causal evidence has remained scarce. Lewis et al. optogenetically stimulate the visual cortex to emulate synaptic input while recording spike output. Cortex resonates at the gamma band (30–90 Hz) and preferentially transmits input that is coherent to the ongoing gamma-band rhythm.

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.

Oscillatory activity in auditory cortex is believed to play an important role in auditory and speech processing. One suggested function of these rhythms is to divide the speech stream into candidate phonemes, syllables, words, and phrases, to be matched with learned linguistic templates. This requires brain rhythms to flexibly synchronize with regular acoustic features of the speech stream. How neuronal circuits implement this task remains unknown. In this study, we explored the contribution of inhibitory currents to flexible phase-locking in neuronal theta oscillators, believed to perform initial syllabic segmentation. We found that a combination of specific intrinsic inhibitory currents at multiple timescales, present in a large class of cortical neurons, enabled exceptionally flexible phase-locking, which could be used to precisely segment speech by identifying vowels at mid-syllable. This suggests that the cells exhibiting these currents are a key component in the brain’s auditory and speech processing architecture.

Layer and rhythm specificity for predictive routing

An established theoretical model, predictive coding, states that the brain is constantly building models (signifying changing predictions) of the environment. The brain does this by forming predictions and signaling sensory inputs which deviate from predictions (“prediction errors”). Various hypotheses exist about how predictive coding could be implemented in the brain. We recorded neural spiking and oscillations with laminar resolution in a network of cortical areas as monkeys performed a working memory task with changing stimulus predictability. Predictability modulated the patterns of feedforward/feedback flow, cortical layers, and oscillations used to process a visual stimulus. These data support the theory of predictive coding but suggest an alternate model for its neural implementation: predictive routing.

In predictive coding, experience generates predictions that attenuate the feeding forward of predicted stimuli while passing forward unpredicted “errors.” Different models have suggested distinct cortical layers, and rhythms implement predictive coding. We recorded spikes and local field potentials from laminar electrodes in five cortical areas (visual area 4 [V4], lateral intraparietal [LIP], posterior parietal area 7A, frontal eye field [FEF], and prefrontal cortex [PFC]) while monkeys performed a task that modulated visual stimulus predictability. During predictable blocks, there was enhanced alpha (8 to 14 Hz) or beta (15 to 30 Hz) power in all areas during stimulus processing and prestimulus beta (15 to 30 Hz) functional connectivity in deep layers of PFC to the other areas. Unpredictable stimuli were associated with increases in spiking and in gamma-band (40 to 90 Hz) power/connectivity that fed forward up the cortical hierarchy via superficial-layer cortex. Power and spiking modulation by predictability was stimulus specific. Alpha/beta power in LIP, FEF, and PFC inhibited spiking in deep layers of V4. Area 7A uniquely showed increases in high-beta (∼22 to 28 Hz) power/connectivity to unpredictable stimuli. These results motivate a conceptual model, predictive routing. It suggests that predictive coding may be implemented via lower-frequency alpha/beta rhythms that “prepare” pathways processing-predicted inputs by inhibiting feedforward gamma rhythms and associated spiking.

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 neural circuit model for a contextual association task inspired by recommender systems

Behavioral data shows that humans and animals have the capacity to learn rules of associations applied to specific examples, and generalize these rules to a broad variety of contexts. This article focuses on neural circuit mechanisms to perform a context-dependent association task that requires linking sensory stimuli to behavioral responses and generalizing to multiple other symmetrical contexts. The model uses neural gating units that regulate the pattern of physiological connectivity within the circuit. These neural gating units can be used in a learning framework that performs low-rank matrix factorization analogous to recommender systems, allowing generalization with high accuracy to a wide range of additional symmetrical contexts. The neural gating units are trained with a biologically inspired framework involving traces of Hebbian modification that are updated based on the correct behavioral output of the network. This modeling demonstrates potential neural mechanisms for learning context-dependent association rules and for the change in selectivity of neurophysiological responses in the hippocampus. The proposed computational model is evaluated using simulations of the learning process and the application of the model to new stimuli. Further, human subject behavioral experiments were performed and the results validate the key observation of a low-rank synaptic matrix structure linking stimuli to responses.

State-Dependent Entrainment of Prefrontal Cortex Local Field Potential Activity Following Patterned Stimulation of the Cerebellar Vermis

The cerebellum is involved in sensorimotor, cognitive, and emotional functions through cerebello-cerebral connectivity. Cerebellar neurostimulation thus likely affects cortical circuits, as has been shown in studies using cerebellar stimulation to treat neurological disorders through modulation of frontal EEG oscillations. Here we studied the effects of different frequencies of cerebellar stimulation on oscillations and coherence in the cerebellum and prefrontal cortex in the urethane-anesthetized rat. Local field potentials were recorded in the right lateral cerebellum (Crus I/II) and bilaterally in the prefrontal cortex (frontal association area, FrA) in adult male Sprague-Dawley rats. Stimulation was delivered to the cerebellar vermis (lobule VII) using single pulses (0.2 Hz for 60 s), or repeated pulses at 1 Hz (30 s), 5 Hz (10 s), 25 Hz (2 s), and 50 Hz (1 s). Effects of stimulation were influenced by the initial state of EEG activity which varies over time during urethane-anesthesia; 1 Hz stimulation was more effective when delivered during the slow-wave state (Stage 1), while stimulation with single-pulse, 25, and 50 Hz showed stronger effects during the activated state (Stage 2). Single-pulses resulted in increases in oscillatory power in the delta and theta bands for the cerebellum, and in frequencies up to 80 Hz in cortical sites. 1 Hz stimulation induced a decrease in 0–30 Hz activity and increased activity in the 30–200 Hz range, in the right FrA. 5 Hz stimulation reduced power in high frequencies in Stage 1 and induced mixed effects during Stage 2.25 Hz stimulation increased cortical power at low frequencies during Stage 2, and increased power in higher frequency bands during Stage 1. Stimulation at 50 Hz increased delta-band power in all recording sites, with the strongest and most rapid effects in the cerebellum. 25 and 50 Hz stimulation also induced state-dependent effects on cerebello-cortical and cortico-cortical coherence at high frequencies. Cerebellar stimulation can therefore entrain field potential activity in the FrA and drive synchronization of cerebello-cortical and cortico-cortical networks in a frequency-dependent manner. These effects highlight the role of the cerebellar vermis in modulating large-scale synchronization of neural networks in non-motor frontal cortex.

The Role of Hierarchical Dynamical Functions in Coding for Episodic Memory and Cognition

Behavioral research in human verbal memory function led to the initial definition of episodic memory and semantic memory. A complete model of the neural mechanisms of episodic memory must include the capacity to encode and mentally reconstruct everything that humans can recall from their experience. This article proposes new model features necessary to address the complexity of episodic memory encoding and recall in the context of broader cognition and the functional properties of neurons that could contribute to this broader scope of memory. Many episodic memory models represent individual snapshots of the world with a sequence of vectors, but a full model must represent complex functions encoding and retrieving the relations between multiple stimulus features across space and time on multiple hierarchical scales. Episodic memory involves not only the space and time of an agent experiencing events within an episode but also features shown in neurophysiological data such as coding of speed, direction, boundaries, and objects. Episodic memory includes not only a spatio-temporal trajectory of a single agent but also segments of spatio-temporal trajectories for other agents and objects encountered in the environment consistent with data on encoding the position and angle of sensory features of objects and boundaries. We will discuss potential interactions of episodic memory circuits in the hippocampus and entorhinal cortex with distributed neocortical circuits that must represent all features of human cognition.

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.