Synaptic and dendritic architecture of different types of hippocampal somatostatin interneurons.

GABAergic inhibitory neurons fundamentally shape the activity and plasticity of cortical circuits. A major subset of these neurons contains somatostatin (SOM); these cells play crucial roles in neuroplasticity, learning, and memory in many brain areas including the hippocampus, and are implicated in several neuropsychiatric diseases and neurodegenerative disorders. Two main types of SOM-containing cells in area CA1 of the hippocampus are oriens-lacunosum-moleculare (OLM) cells and hippocampo-septal (HS) cells. These cell types show many similarities in their soma-dendritic architecture, but they have different axonal targets, display different activity patterns in vivo, and are thought to have distinct network functions. However, a complete understanding of the functional roles of these interneurons requires a precise description of their intrinsic computational properties and their synaptic interactions. In the current study we generated, analyzed, and make available several key data sets that enable a quantitative comparison of various anatomical and physiological properties of OLM and HS cells in mouse. The data set includes detailed scanning electron microscopy (SEM)-based 3D reconstructions of OLM and HS cells along with their excitatory and inhibitory synaptic inputs. Combining this core data set with other anatomical data, patch-clamp electrophysiology, and compartmental modeling, we examined the precise morphological structure, inputs, outputs, and basic physiological properties of these cells. Our results highlight key differences between OLM and HS cells, particularly regarding the density and distribution of their synaptic inputs and mitochondria. For example, we estimated that an OLM cell receives about 8,400, whereas an HS cell about 15,600 synaptic inputs, about 16% of which are GABAergic. Our data and models provide insight into the possible basis of the different functionality of OLM and HS cell types and supply essential information for more detailed functional models of these neurons and the hippocampal network.

The medial septum controls hippocampal supra-theta oscillations.

Hippocampal theta oscillations orchestrate faster beta-to-gamma oscillations facilitating the segmentation of neural representations during navigation and episodic memory. Supra-theta rhythms of hippocampal CA1 are coordinated by local interactions as well as inputs from the entorhinal cortex (EC) and CA3 inputs. However, theta-nested gamma-band activity in the medial septum (MS) suggests that the MS may control supra-theta CA1 oscillations. To address this, we performed multi-electrode recordings of MS and CA1 activity in rodents and found that MS neuron firing showed strong phase-coupling to theta-nested supra-theta episodes and predicted changes in CA1 beta-to-gamma oscillations on a cycle-by-cycle basis. Unique coupling patterns of anatomically defined MS cell types suggested that indirect MS-to-CA1 pathways via the EC and CA3 mediate distinct CA1 gamma-band oscillations. Optogenetic activation of MS parvalbumin-expressing neurons elicited theta-nested beta-to-gamma oscillations in CA1. Thus, the MS orchestrates hippocampal network activity at multiple temporal scales to mediate memory encoding and retrieval.

Huygens synchronization of medial septal pacemaker neurons generates hippocampal theta oscillation.

Episodic learning and memory retrieval are dependent on hippocampal theta oscillation, thought to rely on the GABAergic network of the medial septum (MS). To test how this network achieves theta synchrony, we recorded MS neurons and hippocampal local field potential simultaneously in anesthetized and awake mice and rats. We show that MS pacemakers synchronize their individual rhythmicity frequencies, akin to coupled pendulum clocks as observed by Huygens. We optogenetically identified them as parvalbumin-expressing GABAergic neurons, while MS glutamatergic neurons provide tonic excitation sufficient to induce theta. In accordance, waxing and waning tonic excitation is sufficient to toggle between theta and non-theta states in a network model of single-compartment inhibitory pacemaker neurons. These results provide experimental and theoretical support to a frequency-synchronization mechanism for pacing hippocampal theta, which may serve as an inspirational prototype for synchronization processes in the central nervous system from Nematoda to Arthropoda to Chordate and Vertebrate phyla.

Hippocampal sharp wave-ripples and the associated sequence replay emerge from structured synaptic interactions in a network model of area CA3.

Hippocampal place cells are activated sequentially as an animal explores its environment. These activity sequences are internally recreated ('replayed'), either in the same or reversed order, during bursts of activity (sharp wave-ripples [SWRs]) that occur in sleep and awake rest. SWR-associated replay is thought to be critical for the creation and maintenance of long-term memory. In order to identify the cellular and network mechanisms of SWRs and replay, we constructed and simulated a data-driven model of area CA3 of the hippocampus. Our results show that the chain-like structure of recurrent excitatory interactions established during learning not only determines the content of replay, but is essential for the generation of the SWRs as well. We find that bidirectional replay requires the interplay of the experimentally confirmed, temporally symmetric plasticity rule, and cellular adaptation. Our model provides a unifying framework for diverse phenomena involving hippocampal plasticity, representations, and dynamics, and suggests that the structured neural codes induced by learning may have greater influence over cortical network states than previously appreciated.

HippoUnit: A software tool for the automated testing and systematic comparison of detailed models of hippocampal neurons based on electrophysiological data.

Anatomically and biophysically detailed data-driven neuronal models have become widely used tools for understanding and predicting the behavior and function of neurons. Due to the increasing availability of experimental data from anatomical and electrophysiological measurements as well as the growing number of computational and software tools that enable accurate neuronal modeling, there are now a large number of different models of many cell types available in the literature. These models were usually built to capture a few important or interesting properties of the given neuron type, and it is often unknown how they would behave outside their original context. In addition, there is currently no simple way of quantitatively comparing different models regarding how closely they match specific experimental observations. This limits the evaluation, re-use and further development of the existing models. Further, the development of new models could also be significantly facilitated by the ability to rapidly test the behavior of model candidates against the relevant collection of experimental data. We address these problems for the representative case of the CA1 pyramidal cell of the rat hippocampus by developing an open-source Python test suite, which makes it possible to automatically and systematically test multiple properties of models by making quantitative comparisons between the models and electrophysiological data. The tests cover various aspects of somatic behavior, and signal propagation and integration in apical dendrites. To demonstrate the utility of our approach, we applied our tests to compare the behavior of several different rat hippocampal CA1 pyramidal cell models from the ModelDB database against electrophysiological data available in the literature, and evaluated how well these models match experimental observations in different domains. We also show how we employed the test suite to aid the development of models within the European Human Brain Project (HBP), and describe the integration of the tests into the validation framework developed in the HBP, with the aim of facilitating more reproducible and transparent model building in the neuroscience community.

Activity and Coupling to Hippocampal Oscillations of Median Raphe GABAergic Cells in Awake Mice.

Ascending serotonergic/glutamatergic projection from the median raphe region (MRR) to the hippocampal formation regulates both encoding and consolidation of memory and the oscillations associated with them. The firing of various types of MRR neurons exhibits rhythmic modulation coupled to hippocampal oscillatory activity. A possible intermediary between rhythm-generating forebrain regions and entrained ascending modulation may be the GABAergic circuit in the MRR, known to be targeted by a diverse array of top-down inputs. However, the activity of inhibitory MRR neurons in an awake animal is still largely unexplored. In this study, we utilized whole cell patch-clamp, single cell, and multichannel extracellular recordings of GABAergic and non-GABAergic MRR neurons in awake, head-fixed mice. First, we have demonstrated that glutamatergic and serotonergic neurons receive both transient, phasic, and sustained tonic inhibition. Then, we observed substantial heterogeneity of GABAergic firing patterns but a marked modulation of activity by brain states and fine timescale coupling of spiking to theta and ripple oscillations. We also uncovered a correlation between the preferred theta phase and the direction of activity change during ripples, suggesting the segregation of inhibitory neurons into functional groups. Finally, we could detect complementary alteration of non-GABAergic neurons' ripple-coupled activity. Our findings support the assumption that the local inhibitory circuit in the MRR may synchronize ascending serotonergic/glutamatergic modulation with hippocampal activity on a subsecond timescale.

Publisher Correction: Distinct synchronization, cortical coupling and behavioral function of two basal forebrain cholinergic neuron types.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Distinct synchronization, cortical coupling and behavioral function of two basal forebrain cholinergic neuron types.

Basal forebrain cholinergic neurons (BFCNs) modulate synaptic plasticity, cortical processing, brain states and oscillations. However, whether distinct types of BFCNs support different functions remains unclear. Therefore, we recorded BFCNs in vivo, to examine their behavioral functions, and in vitro, to study their intrinsic properties. We identified two distinct types of BFCNs that differ in their firing modes, synchronization properties and behavioral correlates. Bursting cholinergic neurons (Burst-BFCNs) fired synchronously, phase-locked to cortical theta activity and fired precisely timed bursts after reward and punishment. Regular-firing cholinergic neurons (Reg-BFCNs) were found predominantly in the posterior basal forebrain, displayed strong theta rhythmicity and responded with precise single spikes after behavioral outcomes. In an auditory detection task, synchronization of Burst-BFCNs to the auditory cortex predicted the timing of behavioral responses, whereas tone-evoked cortical coupling of Reg-BFCNs predicted correct detections. We propose that differential recruitment of two basal forebrain cholinergic neuron types generates behavior-specific cortical activation.

Amyloid ß induces interneuron-specific changes in the hippocampus of APPNL-F mice.

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by cognitive decline and amyloid-beta (Aß) depositions generated by the proteolysis of amyloid precursor protein (APP) in the brain. In APPNL-F mice, APP gene was humanized and contains two familial AD mutations, and APP-unlike other mouse models of AD-is driven by the endogenous mouse APP promoter. Similar to people without apparent cognitive dysfunction but with heavy Aß plaque load, we found no significant decline in the working memory of adult APPNL-F mice, but these mice showed decline in the expression of normal anxiety. Using immunohistochemistry and 3D block-face scanning electron microscopy, we found no changes in GABAA receptor positivity and size of somatic and dendritic synapses of hippocampal interneurons. We did not find alterations in the level of expression of perineuronal nets around parvalbumin (PV) interneurons or in the density of PV- or somatostatin-positive hippocampal interneurons. However, in contrast to other investigated cell types, PV interneuron axons were occasionally mildly dystrophic around Aß plaques, and the synapses of PV-positive axon initial segment (AIS)-targeting interneurons were significantly enlarged. Our results suggest that PV interneurons are highly resistant to amyloidosis in APPNL-F mice and amyloid-induced increase in hippocampal pyramidal cell excitability may be compensated by PV-positive AIS-targeting cells. Mechanisms that make PV neurons more resilient could therefore be exploited in the treatment of AD for mitigating Aß-related inflammatory effects on neurons.

Median raphe controls acquisition of negative experience in the mouse.

Adverse events need to be quickly evaluated and memorized, yet how these processes are coordinated is poorly understood. We discovered a large population of excitatory neurons in mouse median raphe region (MRR) expressing vesicular glutamate transporter 2 (vGluT2) that received inputs from several negative experience-related brain centers, projected to the main aversion centers, and activated the septohippocampal system pivotal for learning of adverse events. These neurons were selectively activated by aversive but not rewarding stimuli. Their stimulation induced place aversion, aggression, depression-related anhedonia, and suppression of reward-seeking behavior and memory acquisition-promoting hippocampal theta oscillations. By contrast, their suppression impaired both contextual and cued fear memory formation. These results suggest that MRR vGluT2 neurons are crucial for the acquisition of negative experiences and may play a central role in depression-related mood disorders.

Brainstem nucleus incertus controls contextual memory formation.

Hippocampal pyramidal cells encode memory engrams, which guide adaptive behavior. Selection of engram-forming cells is regulated by somatostatin-positive dendrite-targeting interneurons, which inhibit pyramidal cells that are not required for memory formation. Here, we found that γ-aminobutyric acid (GABA)-releasing neurons of the mouse nucleus incertus (NI) selectively inhibit somatostatin-positive interneurons in the hippocampus, both monosynaptically and indirectly through the inhibition of their subcortical excitatory inputs. We demonstrated that NI GABAergic neurons receive monosynaptic inputs from brain areas processing important environmental information, and their hippocampal projections are strongly activated by salient environmental inputs in vivo. Optogenetic manipulations of NI GABAergic neurons can shift hippocampal network state and bidirectionally modify the strength of contextual fear memory formation. Our results indicate that brainstem NI GABAergic cells are essential for controlling contextual memories.

The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow.

Every neuron is part of a network, exerting its function by transforming multiple spatiotemporal synaptic input patterns into a single spiking output. This function is specified by the particular shape and passive electrical properties of the neuronal membrane, and the composition and spatial distribution of ion channels across its processes. For a variety of physiological or pathological reasons, the intrinsic input/output function may change during a neuron's lifetime. This process results in high variability in the peak specific conductance of ion channels in individual neurons. The mechanisms responsible for this variability are not well understood, although there are clear indications from experiments and modeling that degeneracy and correlation among multiple channels may be involved. Here, we studied this issue in biophysical models of hippocampal CA1 pyramidal neurons and interneurons. Using a unified data-driven simulation workflow and starting from a set of experimental recordings and morphological reconstructions obtained from rats, we built and analyzed several ensembles of morphologically and biophysically accurate single cell models with intrinsic electrophysiological properties consistent with experimental findings. The results suggest that the set of conductances expressed in any given hippocampal neuron may be considered as belonging to two groups: one subset is responsible for the major characteristics of the firing behavior in each population and the other is responsible for a robust degeneracy. Analysis of the model neurons suggests several experimentally testable predictions related to the combination and relative proportion of the different conductances that should be expressed on the membrane of different types of neurons for them to fulfill their role in the hippocampus circuitry.

Co-transmission of acetylcholine and GABA regulates hippocampal states.

The basal forebrain cholinergic system is widely assumed to control cortical functions via non-synaptic transmission of a single neurotransmitter. Yet, we find that mouse hippocampal cholinergic terminals invariably establish GABAergic synapses, and their cholinergic vesicles dock at those synapses only. We demonstrate that these synapses do not co-release but co-transmit GABA and acetylcholine via different vesicles, whose release is triggered by distinct calcium channels. This co-transmission evokes composite postsynaptic potentials, which are mutually cross-regulated by presynaptic autoreceptors. Although postsynaptic cholinergic receptor distribution cannot be investigated, their response latencies suggest a focal, intra- and/or peri-synaptic localisation, while GABAA receptors are detected intra-synaptically. The GABAergic component alone effectively suppresses hippocampal sharp wave-ripples and epileptiform activity. Therefore, the differentially regulated GABAergic and cholinergic co-transmission suggests a hitherto unrecognised level of control over cortical states. This novel model of hippocampal cholinergic neurotransmission may lead to alternative pharmacotherapies after cholinergic deinnervation seen in neurodegenerative disorders.

Hippocampal Network Dynamics during Rearing Episodes.

Animals build a model of their surroundings on the basis of information gathered during exploration. Rearing on the hindlimbs changes the vantage point of the animal, increasing the sampled area of the environment. This environmental knowledge is suggested to be integrated into a cognitive map stored by the hippocampus. Previous studies have found that damage to the hippocampus impairs rearing. Here, we characterize the operational state of the hippocampus during rearing episodes. We observe an increase of theta frequency paralleled by a sink in the dentate gyrus and a prominent theta-modulated fast gamma transient in the middle molecular layer. On the descending phase of rearing, a decrease of theta power is detected. Place cells stop firing during rearing, while a different subset of putative pyramidal cells is activated. Our results suggest that the hippocampus switches to a different operational state during rearing, possibly to update spatial representation with information from distant sources.

Regulation of Hippocampal 5-HT Release by P2X7 Receptors in Response to Optogenetic Stimulation of Median Raphe Terminals of Mice.

Serotonergic and glutamatergic neurons of median raphe region (MRR) play a pivotal role in the modulation of affective and cognitive functions. These neurons synapse both onto themselves and remote cortical areas. P2X7 receptors (P2rx7) are ligand gated ion channels expressed by central presynaptic excitatory nerve terminals and involved in the regulation of neurotransmitter release. P2rx7s are implicated in various neuropsychiatric conditions such as schizophrenia and depression. Here we investigated whether 5-HT release released from the hippocampal terminals of MRR is subject to modulation by P2rx7s. To achieve this goal, an optogenetic approach was used to selectively activate subpopulation of serotonergic terminals derived from the MRR locally, and one of its target area, the hippocampus. Optogenetic activation of neurons in the MRR with 20 Hz was correlated with freezing and enhanced locomotor activity of freely moving mice and elevated extracellular levels of 5-HT, glutamate but not GABA in vivo. Similar optical stimulation (OS) significantly increased [3H]5-HT and [3H]glutamate release in acute MRR and hippocampal slices. We examined spatial and temporal patterns of [3H]5-HT release and the interaction between the serotonin and glutamate systems. Whilst [3H]5-HT release from MRR neurons was [Ca2+]o-dependent and sensitive to TTX, CNQX and DL-AP-5, release from hippocampal terminals was not affected by the latter drugs. Hippocampal [3H]5-HT released by electrical but not OS was subject to modulation by 5- HT1B/D receptors agonist sumatriptan (1 µM), whereas the selective 5-HT1A agonist buspirone (0.1 µM) was without effect. [3H]5-HT released by electrical and optical stimulation was decreased in mice genetically deficient in P2rx7s, and after perfusion with selective P2rx7 antagonists, JNJ-47965567 (0.1 µM), and AZ-10606120 (0.1 µM). Optical and electrical stimulation elevated the extracellular level of ATP. Our results demonstrate for the first time the modulation of 5-HT release from hippocampal MRR terminals by the endogenous activation of P2rx7s. P2rx7 mediated modulation of 5-HT release could contribute to various physiological and pathophysiological phenomena, related to hippocampal serotonergic transmission.

Median raphe region stimulation alone generates remote, but not recent fear memory traces.

The median raphe region (MRR) is believed to control the fear circuitry indirectly, by influencing the encoding and retrieval of fear memories by amygdala, hippocampus and prefrontal cortex. Here we show that in addition to this established role, MRR stimulation may alone elicit the emergence of remote but not recent fear memories. We substituted electric shocks with optic stimulation of MRR in C57BL/6N male mice in an optogenetic conditioning paradigm and found that stimulations produced agitation, but not fear, during the conditioning trial. Contextual fear, reflected by freezing was not present the next day, but appeared after a 7 days incubation. The optogenetic silencing of MRR during electric shocks ameliorated conditioned fear also seven, but not one day after conditioning. The optogenetic stimulation patterns (50Hz theta burst and 20Hz) used in our tests elicited serotonin release in vitro and lead to activation primarily in the periaqueductal gray examined by c-Fos immunohistochemistry. Earlier studies demonstrated that fear can be induced acutely by stimulation of several subcortical centers, which, however, do not generate persistent fear memories. Here we show that the MRR also elicits fear, but this develops slowly over time, likely by plastic changes induced by the area and its connections. These findings assign a specific role to the MRR in fear learning. Particularly, we suggest that this area is responsible for the durable sensitization of fear circuits towards aversive contexts, and by this, it contributes to the persistence of fear memories. This suggests the existence a bottom-up control of fear circuits by the MRR, which complements the top-down control exerted by the medial prefrontal cortex.

Cellular architecture and transmitter phenotypes of neurons of the mouse median raphe region.

The median raphe region (MRR, which consist of MR and paramedian raphe regions) plays a crucial role in regulating cortical as well as subcortical network activity and behavior, while its malfunctioning may lead to disorders, such as schizophrenia, major depression, or anxiety. Mouse MRR neurons are classically identified on the basis of their serotonin (5-HT), vesicular glutamate transporter type 3 (VGLUT3), and gamma-aminobutyric acid (GABA) contents; however, the exact cellular composition of MRR regarding transmitter phenotypes is still unknown. Using an unbiased stereological method, we found that in the MR, 8.5 % of the neurons were 5-HT, 26 % were VGLUT3, and 12.8 % were 5-HT and VGLUT3 positive; whereas 37.2 % of the neurons were GABAergic, and 14.4 % were triple negative. In the whole MRR, 2.1 % of the neurons were 5-HT, 7 % were VGLUT3, and 3.6 % were 5-HT and VGLUT3 positive; whereas 61 % of the neurons were GABAergic. Surprisingly, 25.4 % of the neurons were triple negative and were only positive for the neuronal marker NeuN. PET-1/ePET-Cre transgenic mouse lines are widely used to specifically manipulate only 5-HT containing neurons. Interestingly, however, using the ePET-Cre transgenic mice, we found that far more VGLUT3 positive cells expressed ePET than 5-HT positive cells, and about 38 % of the ePET cells contained only VGLUT3, while more than 30 % of 5-HT cells were ePET negative. These data should facilitate the reinterpretation of PET-1/ePET related data in the literature and the identification of the functional role of a putatively new type of triple-negative neuron in the MRR.

The Effects of Realistic Synaptic Distribution and 3D Geometry on Signal Integration and Extracellular Field Generation of Hippocampal Pyramidal Cells and Inhibitory Neurons.

In vivo and in vitro multichannel field and somatic intracellular recordings are frequently used to study mechanisms of network pattern generation. When interpreting these data, neurons are often implicitly considered as electrotonically compact cylinders with a homogeneous distribution of excitatory and inhibitory inputs. However, the actual distributions of dendritic length, diameter, and the densities of excitatory and inhibitory input are non-uniform and cell type-specific. We first review quantitative data on the dendritic structure and synaptic input and output distribution of pyramidal cells (PCs) and interneurons in the hippocampal CA1 area. Second, using multicompartmental passive models of four different types of neurons, we quantitatively explore the effect of differences in dendritic structure and synaptic distribution on the errors and biases of voltage clamp measurements of inhibitory and excitatory postsynaptic currents. Finally, using the 3-dimensional distribution of dendrites and synaptic inputs we calculate how different inhibitory and excitatory inputs contribute to the generation of local field potential in the hippocampus. We analyze these effects at different realistic background activity levels as synaptic bombardment influences neuronal conductance and thus the propagation of signals in the dendritic tree. We conclude that, since dendrites are electrotonically long and entangled in 3D, somatic intracellular and field potential recordings miss the majority of dendritic events in some cell types, and thus overemphasize the importance of perisomatic inhibitory inputs and belittle the importance of complex dendritic processing. Modeling results also suggest that PCs and inhibitory neurons probably use different input integration strategies. In PCs, second- and higher-order thin dendrites are relatively well-isolated from each other, which may support branch-specific local processing as suggested by studies of active dendritic integration. In the electrotonically compact parvalbumin- and cholecystokinincontaining interneurons, synaptic events are visible in the whole dendritic arbor, and thus the entire dendritic tree may form a single integrative element. Calretinin-containing interneurons were found to be electrotonically extended, which suggests the possibility of complex dendritic processing in this cell type. Our results also highlight the need for the integration of methods that allow the measurement of dendritic processes into studies of synaptic interactions and dynamics in neural networks.

Divergent in vivo activity of non-serotonergic and serotonergic VGluT3-neurones in the median raphe region.

KEY POINTS: The median raphe is a key subcortical modulatory centre involved in several brain functions, such as regulation of the sleep-wake cycle, emotions and memory storage. A large proportion of median raphe neurones are glutamatergic and implement a radically different mode of communication compared to serotonergic cells, although their in vivo activity is unknown. We provide the first description of the in vivo, brain state-dependent firing properties of median raphe glutamatergic neurones identified by immunopositivity for the vesicular glutamate transporter type 3 (VGluT3) and serotonin (5-HT). Glutamatergic populations (VGluT3+/5-HT- and VGluT3+/5-HT+) were compared with the purely serotonergic (VGluT3-/5-HT+ and VGluT3-/5-HT-) neurones. VGluT3+/5-HT+ neurones fired similar to VGluT3-/5-HT+ cells, whereas they significantly diverged from the VGluT3+/5-HT- population. Activity of the latter subgroup resembled the spiking of VGluT3-/5-HT- cells, except for their diverging response to sensory stimulation. The VGluT3+ population of the median raphe may broadcast rapidly varying signals on top of a state-dependent, tonic modulation. ABSTRACT: Subcortical modulation is crucial for information processing in the cerebral cortex. Besides the canonical neuromodulators, glutamate has recently been identified as a key cotransmitter of numerous monoaminergic projections. In the median raphe, a pure glutamatergic neurone population projecting to limbic areas was also discovered with a possibly novel, yet undetermined function. In the present study, we report the first functional description of the vesicular glutamate transporter type 3 (VGluT3)-expressing median raphe neurones. Because there is no appropriate genetic marker for the separation of serotonergic (5-HT+) and non-serotonergic (5-HT-) VGluT3+ neurones, we utilized immunohistochemistry after recording and juxtacellular labelling in anaesthetized rats. VGluT3+/5-HT- neurones fired faster, more variably and were permanently activated during sensory stimulation, as opposed to the transient response of the slow firing VGluT3-/5-HT+ subgroup. VGluT3+/5-HT- cells were also more active during hippocampal theta. In addition, the VGluT3-/5-HT- population, comprising putative GABAergic cells, resembled the firing of VGluT3+/5-HT- neurones but without any significant reaction to the sensory stimulus. Interestingly, the VGluT3+/5-HT+ group, spiking slower than the VGluT3+/5-HT- population, exhibited a mixed response (i.e. the initial transient activation was followed by a sustained elevation of firing). Phase coupling to hippocampal and prefrontal slow oscillations was found in VGluT3+/5-HT- neurones, also differentiating them from the VGluT3+/5-HT+ subpopulation. Taken together, glutamatergic neurones in the median raphe may implement multiple, highly divergent forms of modulation in parallel: a slow, tonic mode interrupted by sensory-evoked rapid transients, as well as a fast one capable of conveying complex patterns influenced by sensory inputs.