LFP publications

Auditory-cued behavioral training can alter neural circuits in primary auditory cortex (A1), but the mechanisms and consequences of experience-dependent cortical plasticity are not fully understood. To address this issue, we trained adult rats to detect a 5 kHz target in order to receive a food reward. After 14 days training we identified three locations within A1: (i) the region representing the characteristic frequency (CF) 5 kHz, (ii) a nearby region with CF ∼10 kHz, and (iii) a more distant region with CF ∼20 kHz. In order to compare functional connectivity in A1 near to, vs. far from, the representation of the target frequency, we placed a 16-channel multiprobe in middle- (∼10 kHz) and high- (∼20 kHz) CF regions and obtained current-source density (CSD) profiles evoked by a range of tone stimuli (CF±1-3 octaves in quarter-octave steps). Our aim was to construct "CSD receptive fields" (CSD RFs) in order to determine the laminar and spectral profile of tone-evoked current sinks, and infer changes to thalamocortical and intracortical inputs. Behavioral training altered CSD RFs at the 10 kHz, but not 20 kHz, site relative to CSD RFs in untrained control animals. At the 10 kHz site, current sinks evoked by the target frequency were enhanced in layer 2/3, but the initial current sink in layer 4 was not altered. The results imply training-induced plasticity along intracortical pathways connecting the target representation with nearby cortical regions. Finally, we related behavioral performance (sensitivity index, d') to CSD responses in individual animals, and found a significant correlation between the development of d' over training and the amplitude of the target-evoked current sink in layer 2/3. The results suggest that plasticity along intracortical pathways is important for auditory learning.

The efficacy of spinal cord stimulators is dependent on the ability of the device to functionally activate targeted structures within the spinal cord, while avoiding activation of near-by non-targeted structures. In theory, these objectives can best be achieved by delivering electrical stimuli directly to the surface of the spinal cord. The current experiments were performed to study the influence of different stimulating electrode positions on patterns of spinal cord electrophysiological activation. A custom-designed spinal cord neurostimulator was used to investigate the effects of lead position and stimulus amplitude on cortical electrophysiological responses to spinal cord stimulation. Brain recordings were obtained from subdural grids placed in four adult sheep. We systematically varied the position of the stimulating lead relative to the spinal cord and the voltage delivered by the device at each position, and then examined how these variables influenced cortical responses. A clear relationship was observed between voltage and electrode position, and the magnitude of high gamma-band oscillations. Direct stimulation of the dorsal column contralateral to the grid required the lowest voltage to evoke brain responses to spinal cord stimulation. Given the lower voltage thresholds associated with direct stimulation of the dorsal column, and its possible impact on the therapeutic window, this intradural modality may have particular clinical advantages over standard epidural techniques now in routine use.

The central thalamus plays an important role in the regulation of arousal and allocation of attentional resources in the performance of even simple tasks. To assess the contribution of central thalamic neurons to short-term adjustments of attentional effort, we analyzed 166 microelectrode recordings obtained from two rhesus monkeys performing a visuomotor simple reaction time task with a variable foreperiod. Multiunit responses showed maintained firing rate elevations during the variable delay period of the task in ∼24% of recording sites. Simultaneously recorded local field potentials demonstrated significant decreases in power at ∼10–20 Hz and increases in power at 30–100 Hz during the delay period when compared against precue baselines. Comparison of the spectral power of local field potentials during the delay period of correct and incorrect trials showed that, during incorrect trials, similar, but reduced, shifts of spectral power occurred within the same frequency bands. Sustained performance of even simple tasks requires regulation of arousal and attention that combine in the concept of “attentional effort”. Our findings suggest that central thalamic neurons regulate task performance through brief changes in firing rates and spectral power changes during task-relevant short-term shifts of attentional effort. Increases in attentional effort may be reflected in changes within the central thalamic local populations, where correct task performance associates with more robust maintenance of firing rates during the delay period. Such ongoing fluctuations of central thalamic activity likely reflect a mix of influences, including variations in moment-to-moment levels of motivation, arousal, and availability of cognitive resources.

Detailed understanding of the haemodynamic changes that underlie non-invasive neuroimaging techniques such as blood oxygen level dependent functional magnetic resonance imaging is essential if we are to continue to extend the use of these methods for understanding brain function and dysfunction. The use of animal and in particular rodent research models has been central to these endeavours as they allow in-vivo experimental techniques that provide measurements of the haemodynamic response function at high temporal and spatial resolution. A limitation of most of this research is the use of anaesthetic agents which may disrupt or mask important features of neurovascular coupling or the haemodynamic response function. In this study we therefore measured spatiotemporal cortical haemodynamic responses to somatosensory stimulation in awake rats using optical imaging spectroscopy. Trained, restrained animals received non-noxious stimulation of the whisker pad via chronically implanted stimulating microwires whilst optical recordings were made from the contralateral somatosensory cortex through a thin cranial window. The responses we measure from un-anaesthetised animals are substantially different from those reported in previous studies which have used anaesthetised animals. These differences include biphasic response regions (initial increases in blood volume and oxygenation followed by subsequent decreases) as well as oscillations in the response time series of awake animals. These haemodynamic response features do not reflect concomitant changes in the underlying neuronal activity and therefore reflect neurovascular or cerebrovascular processes. These hitherto unreported hyperemic response dynamics may have important implications for the use of anaesthetised animal models for research into the haemodynamic response function.

Photoacoustic tomography (PAT) offers three-dimensional (3D) structural and functional imaging of living biological tissue with label-free, optical absorption contrast. These attributes lend PAT imaging to a wide variety of applications in clinical medicine and preclinical research. Despite advances in live animal imaging with PAT, there is still a need for 3D imaging at centimeter depths in real-time. We report the development of four dimensional (4D) PAT, which integrates time resolutions with 3D spatial resolution, obtained using spherical arrays of ultrasonic detectors. The 4D PAT technique generates motion pictures of imaged tissue, enabling real time tracking of dynamic physiological and pathological processes at hundred micrometer-millisecond resolutions. The 4D PAT technique is used here to image needle-based drug delivery and pharmacokinetics. We also use this technique to monitor 1) fast hemodynamic changes during inter-ictal epileptic seizures and 2) temperature variations during tumor thermal therapy.

Motor tics, a cardinal symptom of Tourette syndrome (TS), are hypothesized to arise from abnormalities within cerebro-basal ganglia circuits. Yet noninvasive neuroimaging of TS has previously identified robust activation in the cerebellum. To date, electrophysiological properties of cerebellar activation and its role in basal ganglia-mediated tic expression remain unknown. We performed multisite, multielectrode recordings of single-unit activity and local field potentials from the cerebellum, basal ganglia, and primary motor cortex using a pharmacologic monkey model of motor tics/TS. Following microinjections of bicuculline into the sensorimotor putamen, periodic tics occurred predominantly in the orofacial region, and a sizable number of cerebellar neurons showed phasic changes in activity associated with tic episodes. Specifically, 64% of the recorded cerebellar cortex neurons exhibited increases in activity, and 85% of the dentate nucleus neurons displayed excitatory, inhibitory, or multiphasic responses. Critically, abnormal discharges of cerebellar cortex neurons and excitatory-type dentate neurons mostly preceded behavioral tic onset, indicating their central origins. Latencies of pathological activity in the cerebellum and primary motor cortex substantially overlapped, suggesting that aberrant signals may be traveling along divergent pathways to these structures from the basal ganglia. Furthermore, the occurrence of tic movement was most closely associated with local field potential spikes in the cerebellum and primary motor cortex, implying that these structures may function as a gate to release overt tic movements. These findings indicate that tic-generating networks in basal ganglia mediated tic disorders extend beyond classical cerebro-basal ganglia circuits, leading to global network dysrhythmia including cerebellar circuits.

Objective: Hierarchical processing of auditory sensory information is believed to occur in two streams: a ventral stream responsible for stimulus identity and a dorsal stream responsible for processing spatial elements of a stimulus. The objective of the current study is to examine neural coding in this processing stream in the context of understanding the possibility for an auditory cortical neural prosthesis. Approach: We examined the selectivity for species-specific primate vocalizations in the ventral auditory processing stream by applying a statistical classifier to neural data recorded from microelectrode arrays. Multi-unit activity (MUA) and local field potential (LFP) data recorded simultaneously from primary auditory complex (AI) and rostral parabelt (PBr) were decoded on a trial-by-trial basis. Main results: While decode performance in AI was well above chance, mean performance in PBr did not deviate >15% from chance level. Mean performance levels were similar for MUA and LFP decodes. Increasing the spectral and temporal resolution improved decode performance; while inter-electrode spacing could be as large as 1.14 mm without degrading decode performance. Significance: These results serve as preliminary guidance for a human auditory cortical neural prosthesis; instructing interface implementation, microstimulation patterns and anatomical placement.

Brain-machine interfaces (BMI) were born to control "actions from thoughts" in order to recover motor capability of patients with impaired functional connectivity between the central and peripheral nervous system. The final goal of our studies is the development of a new proof-of-concept BMI-a neuromorphic chip for brain repair-to reproduce the functional organization of a damaged part of the central nervous system. To reach this ambitious goal, we implemented a multidisciplinary "bottom-up" approach in which in vitro networks are the paradigm for the development of an in silico model to be incorporated into a neuromorphic device. In this paper we present the overall strategy and focus on the different building blocks of our studies: (i) the experimental characterization and modeling of "finite size networks" which represent the smallest and most general self-organized circuits capable of generating spontaneous collective dynamics; (ii) the induction of lesions in neuronal networks and the whole brain preparation with special attention on the impact on the functional organization of the circuits; (iii) the first production of a neuromorphic chip able to implement a real-time model of neuronal networks. A dynamical characterization of the finite size circuits with single cell resolution is provided. A neural network model based on Izhikevich neurons was able to replicate the experimental observations. Changes in the dynamics of the neuronal circuits induced by optical and ischemic lesions are presented respectively for in vitro neuronal networks and for a whole brain preparation. Finally the implementation of a neuromorphic chip reproducing the network dynamics in quasi-real time (10 ns precision) is presented.

Axons from the nucleus magnocellularis (NM) and their targets in nucleus laminaris (NL) form the circuit responsible for encoding interaural time differences (ITDs). In barn owls, NL receives bilateral inputs from NM such that axons from the ipsilateral NM enter NL dorsally, while contralateral axons enter from the ventral side. These afferents and their synapses on NL neurons generate a tone-induced local field potential, or neurophonic, that varies systematically with position in NL. From dorsal to ventral within the nucleus, the best interaural time difference (ITD) of the neurophonic shifts from contralateral space to best ITDs around 0 µs. Earlier recordings suggested that in NL, iso-delay contours ran parallel to the dorsal and ventral borders of NL (Sullivan WE, Konishi M. Proc Natl Acad Sci U S A 83:8400-8404, 1986). This axis is orthogonal to that seen in chicken NL, where a single map of ITD runs from around 0 µs ITD medially to contralateral space laterally (Köppl C, Carr CE. Biol Cyber 98:541-559, 2008). Yet the trajectories of the NM axons are similar in owl and chicken (Seidl AH, Rubel EW, Harris DM, J Neurosci 30:70-80, 2010). We therefore used clicks to measure conduction time in NL and made lesions to mark the 0 µs iso-delay contour in multiple penetrations along an isofrequency slab. Iso-delay contours were not parallel to the dorsal and ventral borders of NL; instead the 0 µs iso-delay contour shifted systematically from a dorsal position in medial NL to a ventral position in lateral NL. Could different conduction delays account for the mediolateral shift in the representation of 0 µs ITD? We measured conduction delays using the neurophonic potential and developed a simple linear model of the delay-line conduction velocity. We then raised young owls with time-delaying earplugs in one ear (Gold JI, Knudsen EI, J Neurophysiol 82:2197-2209, 1999) to examine map plasticity.

It is not known how visual cortical neurons react to several moving objects and how their firing to the motion of one object is affected by neurons firing to another moving object. Here we combine imaging of voltage sensitive dye (VSD) signals, reflecting the population membrane potential from ferret visual areas 17, 18, 19, and 21, with laminar recordings of multiunit activity, (MUA), when two bars moved toward each other in the visual field, occluded one another, and continued on in opposite directions. Two zones of peak MUA, mapping the bars' motion, moved toward each other along the area 17/18 border, which in the ferret maps the vertical meridian of the field of view. This was reflected also in the VSD signal, at both the 17/18 border as well as at the 19/21 border with a short delay. After some 125 ms at the area 19/21 border, the VSD signal increased and became elongated in the direction of motion in front of both of the moving representations. This was directly followed by the phase of the signal reversing and travelling back from the 19/21 border toward the 17/18 border, seemingly without respect for retinotopic boundaries, where it arrived at 150 ms after stimulus onset. At this point the VSD signal in front of the moving bar representations along the 17/18 border also increased and became elongated in the direction of object motion; the signal now being the linear sum of what has been observed in response to single moving bars. When the neuronal populations representing the bars were some 600 μm apart on the cortex, the dye signal and laminar MUA decreased strongly, with the MUA scaling to that of a single bar during occlusion. Despite a short rebound of the dye signal and MUA, the MUA after the occlusion was significantly depressed. The interactions between the neuronal populations mapping the bars' position, and the neurons in between these populations were, apart from 19/21 to 17/18 interaction, mainly lateral-horizontal; first excitatory and inducing firing at the site of future occlusion, then inhibitory just prior to occlusion. After occlusion the neurons that had fired already to the first bar showed delayed and prolonged inhibition in response to the second bar. Thus, the interactions that were particular to the occlusion condition in these experiments were local and inhibitory at short cortical range, and delayed and inhibitory after the occlusion when the bars moved further apart.

The role of early auditory processing may be to extract some elementary features from an acoustic mixture in order to organize the auditory scene. To accomplish this task, the central auditory system may rely on the fact that sensory objects are often composed of spectral edges, i.e., regions where the stimulus energy changes abruptly over frequency. The processing of acoustic stimuli may benefit from a mechanism enhancing the internal representation of spectral edges. While the visual system is thought to rely heavily on this mechanism (enhancing spatial edges), it is still unclear whether a related process plays a significant role in audition. We investigated the cortical representation of spectral edges, using acoustic stimuli composed of multi-tone pips whose time-averaged spectral envelope contained suppressed or enhanced regions. Importantly, the stimuli were designed such that neural responses properties could be assessed as a function of stimulus frequency during stimulus presentation. Our results suggest that the representation of acoustic spectral edges is enhanced in the auditory cortex, and that this enhancement is sensitive to the characteristics of the spectral contrast profile, such as depth, sharpness and width. Spectral edges are maximally enhanced for sharp contrast and large depth. Cortical activity was also suppressed at frequencies within the suppressed region. To note, the suppression of firing was larger at frequencies nearby the lower edge of the suppressed region than at the upper edge. Overall, the present study gives critical insights into the processing of spectral contrasts in the auditory system.

Cerebrocortical injuries such as stroke are a major source of disability. Maladaptive consequences can result from post-injury local reorganization of cortical circuits. For example, epilepsy is a common sequela of cortical stroke, but the mechanisms responsible for seizures following cortical injuries remain unknown. In addition to local reorganization, long-range, extra-cortical connections might be critical for seizure maintenance. In rats, we found that the thalamus, a structure that is remote from, but connected to, the injured cortex, was required to maintain cortical seizures. Thalamocortical neurons connected to the injured epileptic cortex underwent changes in HCN channel expression and became hyperexcitable. Targeting these neurons with a closed-loop optogenetic strategy revealed that reducing their activity in real-time was sufficient to immediately interrupt electrographic and behavioral seizures. This approach is of therapeutic interest for intractable epilepsy, as it spares cortical function between seizures, in contrast with existing treatments, such as surgical lesioning or drugs.

The mismatch negativity (MMN) is a preattentive component of the auditory event-related potential that is elicited by a change in a repetitive acoustic pattern. While MMN has been extensively used in human electrophysiological studies of auditory processing, the neural mechanisms and brain regions underlying its generation remain unclear. We investigate possible homologs of the MMN in macaque primary auditory cortex (A1) using a frequency oddball paradigm in which rare “deviant” tones are randomly interspersed among frequent “standard” tones. Standards and deviants had frequencies equal to the best frequency (BF) of the recorded neural population or to a frequency that evoked a response half the amplitude of the BF response. Early and later field potentials, current source density components, multiunit activity, and induced high-gamma band responses were larger when elicited by deviants than by standards of the same frequency. Laminar analysis indicated that differences between deviant and standard responses were more prominent in later activity, thus suggesting cortical amplification of initial responses driven by thalamocortical inputs. However, unlike the human MMN, larger deviant responses were characterized by the enhancement of “obligatory” responses rather than the introduction of new components. Furthermore, a control condition wherein deviants were interspersed among many tones of variable frequency replicated the larger responses to deviants under the oddball condition. Results suggest that differential responses under the oddball condition in macaque A1 reflect stimulus-specific adaptation rather than deviance detection per se. We conclude that neural mechanisms of deviance detection likely reside in cortical areas outside of A1.

The amygdala, important in forming and storing memories of aversive events, is believed to play a major role in debilitating tinnitus and hyperacusis. To explore this hypothesis, we recorded from the lateral amygdala (LA) and auditory cortex (AC) before and after treating rats with a dose of salicylate that induces tinnitus and hyperacusis-like behavior. Salicylate unexpectedly increased the amplitude of the local field potential (LFP) in the LA making it hyperactive to sounds≥60 dB SPL. Frequency receptive fields (FRFs) of multiunit (MU) clusters in the LA were also dramatically altered by salicylate. Neuronal activity at frequencies below 10 kHz and above 20 kHz was depressed at low intensities, but was greatly enhanced for stimuli between 10 and 20 kHz (frequencies near the pitch of the salicylate-induced tinnitus in the rat). These frequency-dependent changes caused the FRF of many LA neurons to migrate towards 10-20 kHz thereby amplifying activity from this region. To determine if salicylate-induced changes restricted to the LA would remotely affect neural activity in the AC, we used a micropipette to infuse salicylate (20 μl, 2.8 mM) into the amygdala. Local delivery of salicylate to the amygdala significantly increased the amplitude of the LFP recorded in the AC and selectively enhanced the neuronal activity of AC neurons at the mid-frequencies (10-20 kHz), frequencies associated with the tinnitus pitch. Taken together, these results indicate that systemic salicylate treatment can induce hyperactivity and tonotopic shift in the amygdala and infusion of salicylate into the amygdala can profoundly enhance sound-evoked activity in AC, changes likely to increase the perception and emotional salience of tinnitus and loud sounds. This article is part of a Special Issue entitled: Tinnitus Neuroscience.

Oscillatory neural activity within the gamma band (25–90 Hz) is generally thought to be able to provide a timing signal for harmonizing neural computations across different brain regions. Using time–frequency analyses of the dynamics of gamma-band activity in the local field potentials recorded from monkey primary visual cortex, we found identical temporal characteristics of gamma activity in both awake and anesthetized brain states, including large variability of peak frequency, brief oscillatory epochs (<100 ms on average), and stochastic statistics of the incidence and duration of oscillatory events. These findings indicate that gamma-band activity is temporally unstructured and is inherently a stochastic signal generated by neural networks. This idea was corroborated further by our neural-network simulations. Our results suggest that gamma-band activity is too random to serve as a clock signal for synchronizing neuronal responses in awake as in anesthetized monkeys. Instead, gamma-band activity is more likely to be filtered neuronal network noise. Its mean frequency changes with global state and is reduced under anesthesia.

While functional imaging is widely used in studies of the brain, how well the hemodynamic signal represents the underlying neural activity is still unclear. And there is a debate on whether hemodynamic signal is more tightly related to synaptic activity or action potentials. This study intends to address these questions by examining neurovascular coupling driven by pyramidal cells in the motor cortex of rats. Pyramidal cells in the motor cortex of rats were selectively transduced with the light sensitive cation channel channelrhodopsin-2 (ChR2). Electrophysiological recordings and optical intrinsic signal imaging were performed simultaneously and synchronously to capture the neural activity and hemodynamics induced by optical stimulation of ChR2-expressing pyramidal cells. Our results indicate that both synaptic activity (local field potential, LFP) and action potentials (multi-unit activity, MUA) are tightly related to hemodynamic signals. While LFPs in γ band are better in predicting hemodynamic signals elicited by short stimuli, MUA has better predictions to hemodynamic signals elicited by long stimuli. Our results also indicate that strong nonlinearity exists in neurovascular coupling.

Topographically organized maps of the sensory receptor epithelia are regarded as cornerstones of cortical organization as well as valuable readouts of diverse biological processes ranging from evolution to neural plasticity. However, maps are most often derived from multiunit activity recorded in the thalamic input layers of anesthetized animals using near-threshold stimuli. Less distinct topography has been described by studies that deviated from the formula above, which brings into question the generality of the principle. Here, we explicitly compared the strength of tonotopic organization at various depths within core and belt regions of the auditory cortex using electrophysiological measurements ranging from single units to delta-band local field potentials (LFP) in the awake and anesthetized mouse. Unit recordings in the middle cortical layers revealed a precise tonotopic organization in core, but not belt, regions of auditory cortex that was similarly robust in awake and anesthetized conditions. In core fields, tonotopy was degraded outside the middle layers or when LFP signals were substituted for unit activity, due to an increasing proportion of recording sites with irregular tuning for pure tones. However, restricting our analysis to clearly defined receptive fields revealed an equivalent tonotopic organization in all layers of the cortical column and for LFP activity ranging from gamma to theta bands. Thus, core fields represent a transition between topographically organized simple receptive field arrangements that extend throughout all layers of the cortical column and the emergence of nontonotopic representations outside the input layers that are further elaborated in the belt fields.

Activation of nicotinic acetylcholine receptors (nAChRs) by systemic nicotine enhances sensory-cognitive function and sensory-evoked cortical responses. Although nAChRs mediate fast neurotransmission at many synapses in the nervous system, nicotinic regulation of cortical processing is neuromodulatory. To explore potential mechanisms of nicotinic neuromodulation, we examined whether intracellular signal transduction involving mitogen-activated protein kinase (MAPK) contributes to regulation of tone-evoked responses in primary auditory cortex (A1) in the mouse. Systemic nicotine enhanced characteristic frequency (CF) tone-evoked current-source density (CSD) profiles in A1, including the shortest-latency (presumed thalamocortical) current sink in layer 4 and longer-latency (presumed intracortical) sinks in layers 2–4, by increasing response amplitudes and decreasing response latencies. Microinjection of the MAPK kinase (MEK) inhibitor U0126 into the thalamus, targeting the auditory thalamocortical pathway, blocked the effect of nicotine on the initial (thalamocortical) CSD component but did not block enhancement of longer-latency (intracortical) responses. Conversely, microinjection of U0126 into supragranular layers of A1 blocked nicotine's effect on intracortical, but not thalamocortical, CSD components. Simultaneously with enhancement of CF-evoked responses, responses to spectrally distant (nonCF) stimuli were reduced, implying nicotinic “sharpening” of frequency receptive fields, an effect also blocked by MEK inhibition. Consistent with these physiological results, acoustic stimulation with nicotine produced immunolabel for activated MAPK in A1, primarily in layer 2/3 cell bodies. Immunolabel was blocked by intracortical microinjection of the nAChR antagonist dihydro-β-erythroidine, but not methyllycaconitine, implicating α4β2*, but not α7, nAChRs. Thus activation of MAPK in functionally distinct forebrain circuits—thalamocortical, local intracortical, and long-range intracortical—underlies nicotinic neuromodulation of A1.

The aim of this paper is twofold: first, to introduce a neural field model motivated by a well-known neural mass model; second, to show how one can estimate model parameters pertaining to spatial (anatomical) properties of neuronal sources based on EEG or LFP spectra using Bayesian inference. Specifically, we consider neural field models of cortical activity as generative models in the context of dynamic causal modeling (DCM). This paper considers the simplest case of a single cortical source modeled by the spatiotemporal dynamics of hidden neuronal states on a bounded cortical surface or manifold. We build this model using multiple layers, corresponding to cortical lamina in the real cortical manifold. These layers correspond to the populations considered in classical (Jansen and Rit) neural mass models. This allows us to formulate a neural field model that can be reduced to a neural mass model using appropriate constraints on its spatial parameters. In turn, this enables one to compare and contrast the predicted responses from equivalent neural field and mass models respectively. We pursue this using empirical LFP data from a single electrode to show that the parameters controlling the spatial dynamics of cortical activity can be recovered, using DCM, even in the absence of explicit spatial information in observed data.

Memories are thought to be encoded as a distributed representation in the neocortex. The medial prefrontal cortex (mPFC) has been shown to support the expression of memories that initially depend on the hippocampus (HPC), yet the mechanisms by which the HPC and mPFC access the distributed representations in the neocortex are unknown. By measuring phase synchronization of local field potential (LFP) oscillations, we found that learning initiated changes in neuronal communication of the HPC and mPFC with the lateral entorhinal cortex (LEC), an area that is connected with many other neocortical regions. LFPs were recorded simultaneously from the three brain regions while rats formed an association between an auditory stimulus (CS) and eyelid stimulation (US) in a trace eyeblink conditioning paradigm, as well as during retention 1 month following learning. Over the course of learning, theta oscillations in the LEC and mPFC became strongly synchronized following presentation of the CS on trials in which rats exhibited a conditioned response (CR), and this strengthened synchronization was also observed during remote retention. In contrast, CS-evoked theta synchronization between the LEC and HPC decreased with learning. Our results suggest that communication between the LEC and mPFC are strengthened with learning whereas the communication between the LEC and HPC are concomitantly weakened, suggesting that enhanced LEC-mPFC communication may be a neuronal correlate for theoretically proposed neocortical reorganization accompanying encoding and consolidation of a memory.