Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Featured ArticleArticles, Cellular/Molecular

Interdependent Conductances Drive Infraslow Intrinsic Rhythmogenesis in a Subset of Accessory Olfactory Bulb Projection Neurons

Monika Gorin, Chryssanthi Tsitoura, Anat Kahan, Katja Watznauer, Daniela R. Drose, Martijn Arts, Rudolf Mathar, Simon O'Connor, Ileana L. Hanganu-Opatz, Yoram Ben-Shaul and Marc Spehr
Journal of Neuroscience 16 March 2016, 36 (11) 3127-3144; DOI: https://doi.org/10.1523/JNEUROSCI.2520-15.2016
Monika Gorin
1Department of Chemosensation, Institute for Biology II,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chryssanthi Tsitoura
1Department of Chemosensation, Institute for Biology II,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anat Kahan
3School of Medicine, Department of Medical Neurobiology, The Hebrew University of Jerusalem, Jerusalem, 91120 Israel,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katja Watznauer
1Department of Chemosensation, Institute for Biology II,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniela R. Drose
1Department of Chemosensation, Institute for Biology II,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Martijn Arts
2Institute for Theoretical Information Technology, RWTH Aachen University, D-52074 Aachen, Germany,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rudolf Mathar
2Institute for Theoretical Information Technology, RWTH Aachen University, D-52074 Aachen, Germany,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Simon O'Connor
4Biocomputation Group, University of Hertfordshire, Hatfield, Hertfordshire AL10 9AB, United Kingdom, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ileana L. Hanganu-Opatz
5Developmental Neurophysiology, Institute of Neuroanatomy, University Medical Center Hamburg-Eppendorf, D-20251 Hamburg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yoram Ben-Shaul
3School of Medicine, Department of Medical Neurobiology, The Hebrew University of Jerusalem, Jerusalem, 91120 Israel,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marc Spehr
1Department of Chemosensation, Institute for Biology II,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marc Spehr
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

The accessory olfactory system controls social and sexual behavior. However, key aspects of sensory signaling along the accessory olfactory pathway remain largely unknown. Here, we investigate patterns of spontaneous neuronal activity in mouse accessory olfactory bulb mitral cells, the direct neural link between vomeronasal sensory input and limbic output. Both in vitro and in vivo, we identify a subpopulation of mitral cells that exhibit slow stereotypical rhythmic discharge. In intrinsically rhythmogenic neurons, these periodic activity patterns are maintained in absence of fast synaptic drive. The physiological mechanism underlying mitral cell autorhythmicity involves cyclic activation of three interdependent ionic conductances: subthreshold persistent Na+ current, R-type Ca2+ current, and Ca2+-activated big conductance K+ current. Together, the interplay of these distinct conductances triggers infraslow intrinsic oscillations with remarkable periodicity, a default output state likely to affect sensory processing in limbic circuits.

SIGNIFICANCE STATEMENT We show for the first time that some rodent accessory olfactory bulb mitral cells—the direct link between vomeronasal sensory input and limbic output—are intrinsically rhythmogenic. Driven by ≥3 distinct interdependent ionic conductances, infraslow intrinsic oscillations show remarkable periodicity both in vitro and in vivo. As a novel default state, infraslow autorhythmicity is likely to affect limbic processing of pheromonal information.

  • intrinsic oscillations
  • mitral cell
  • olfactory bulb
  • pacemaker
  • rhythmogenesis
  • vomeronasal system

Introduction

In rodents, the accessory olfactory system controls conspecific chemical communication during social interactions. Behaviorally relevant chemosignals are detected by neurons in the vomeronasal organ and processed along sensory pathways that include the accessory olfactory bulb (AOB), amygdala, and hypothalamus (Stowers and Logan, 2010). While two of these three primary processing centers have been intensely studied (LeDoux, 2000; Simerly, 2002), central aspects of AOB physiology remain largely unexplored and functional analogies with neurons of the main olfactory bulb (MOB) are mostly speculative (Dulac and Wagner, 2006). Principal neurons, conventionally referred to as mitral cells, are the sole projection neurons of the AOB. These neurons extend complex, branched primary dendrites that receive excitatory synaptic input from vomeronasal sensory neurons in multiple glomeruli (Larriva-Sahd, 2008). This morphology indicates that sensory computation in AOB mitral cells is fundamentally different from the processing of “glomerulus-specific” information by their MOB counterparts (Dulac and Torello, 2003). Despite important recent insights into the organizational principles of connectivity, sensory input, and integration in the AOB (Del Punta et al., 2002; Ma and Lowe, 2004; Sugai et al., 2005; Wagner et al., 2006; Castro et al., 2007; Ben-Shaul et al., 2010; Smith and Araneda, 2010; Hovis et al., 2012; Leszkowicz et al., 2012; Shpak et al., 2012; Tolokh et al., 2013; Hammen et al., 2014), a conceptual understanding of how the biophysical properties of AOB mitral cells affect their computations is lacking.

Spontaneous activity is a major determinant of a neuron's coding capacity and information transfer function (Rieke et al., 1997). Spontaneous discharge may be sparse or dense, irregular or periodic, generating a continuum of neuronal firing patterns from Poisson-like discharge (Shadlen and Newsome, 1998) to rhythmic clock-like firing (Blankenship and Feller, 2010). Patterned discharge, such as burst firing and/or oscillatory activity, is particularly relevant for neural processing: bursts often represent units of information (Izhikevich et al., 2003), whereas oscillations provide precise temporal windows of excitability for circuit computations (Mizuseki et al., 2009). Different frequencies have been reported to encode specific brain states (Buzsáki et al., 2013). At the lower end of the time scale, oscillations extend into the slow (0.1–1 Hz) and infraslow (<0.1 Hz) range (Schroeder and Lakatos, 2009). Slow stereotypical episodic “up” and “down” states (Blethyn et al., 2006) could either result from regular recurrence of excitatory and inhibitory synaptic barrages (Crunelli and Hughes, 2010) or, alternatively, represent a network-independent intrinsic neuronal property, caused by a cyclical interplay of autonomous depolarizing and hyperpolarizing conductances (Blankenship and Feller, 2010). Such intrinsic pacemaker-like discharge has been implicated both in controlling rhythmic behaviors (Peña et al., 2004; Bucher et al., 2006; Koizumi and Smith, 2008; Tazerart et al., 2008) and in driving oscillatory circuits that play important roles in sensory perception, attention, memory formation, and decision making (Gutierrez et al., 2013). In the MOB, sensory input evokes oscillatory network activity across a wide frequency range (Kay et al., 2009). Autonomous, episodic burst firing of external tufted (ET) cells within a small spectral window (∼0.2–10 Hz; Hayar et al., 2004) has been attributed a role in setting sniff cycle-dependent glomerular synchrony (Hayar et al., 2005; Shao et al., 2009) and, consequently, distinct phase-locking of MOB principle neurons (Smear et al., 2011; Fukunaga et al., 2012). Whether aspects of sensory coding in the AOB are also affected by single cell/network rhythmicity, however, is unclear.

Here, we identify patterns of coordinated spontaneous activity among AOB mitral cells. We demonstrate that a subpopulation of these neurons reveals slow stereotypical rhythmic discharge both in vitro and in vivo. Among these neurons, a group of mitral cells is intrinsically rhythmogenic. Their autorhythmicity is driven by a precisely orchestrated ensemble of interdependent Na+, Ca2+, and K+ conductances. The cyclic (in)activation of these three distinct conductances provides a mechanistic basis for infraslow intrinsic oscillations of remarkable periodicity in a subset of AOB projection neurons.

Materials and Methods

Animals.

All animal procedures were approved by local authorities and were in compliance with European Union legislation (Directive 86/609/EEC) and recommendations by the Federation of European Laboratory Animal Science Associations. C57BL/6 mice (Charles River Laboratories) were housed in groups of both sexes [room temperature (RT); 12 h light/dark cycle; food and water available ad libitum]. For in vivo recordings, sexually naive male BALB/c mice were used (Harlan Laboratories). All in vivo experiments were performed in compliance with the Hebrew University Animal Care and Use Committee.

Chemicals and solutions.

The following solutions, labeled S1–S9, were used: S1, HEPES buffered extracellular solution containing (in mm) 145 NaCl, 5 KCl, 1 CaCl2, 1 MgCl2, 10 HEPES, pH 7.3 (adjusted with NaOH), 300 mOsm (adjusted with glucose); S2, oxygenated (95% O2, 5% CO2) artificial CSF (aCSF) containing (in mm) 124 NaCl, 26 NaHCO3, 3 KCl, 1.25 NaH2PO4, 1.3 MgSO4, 1.3 CaCl2, 10 glucose, pH 7.3, 300 mOsm (adjusted with glucose); S3, oxygenated (95% O2, 5% CO2) cutting solution containing (in mm) 220 sucrose, 26 NaHCO3, 3 KCl, 1.25 NaH2PO4, 2.6 MgSO4, 10 glucose, pH 7.3, 300 mOsm (adjusted with glucose); S4, extracellular Ba2+ solution containing (in mm) 100 NaCl, 25 tetraethylammonium chloride (TEA-Cl), 10 4-aminopyridine, 10 BaCl2, 1 MgCl2, 10 HEPES, 0.001 tetrodotoxin (TTX), pH 7.3 (adjusted with HCl), 300 mOsm; S5, standard pipette solution containing (in mm) 10.5 KCl, 125 KOH, 125 gluconic acid, 2 MgCl2, 1 EGTA, 0.3 CaCl2, 10 HEPES, 2 Mg-ATP, 1 Na-GTP (free Mg2+, 2 mm; free Ca2+, 130 nm), pH 7.1 (adjusted with KOH), 290 mOsm; S6, symmetrical chloride pipette solution containing (in mm) 143 KCl, 2 KOH, 1 EGTA, 0.3 CaCl2, 10 HEPES, 2 Mg-ATP, 1 Na-GTP (free Ca2+, 130 nm), pH 7.1 (adjusted with KOH), 290 mOsm; S7, Cs+-based pipette solution containing (in mm) 125 gluconic acid, 105 CsOH, 20 NaOH,10.5 CsCl, 2 MgCl2, 1 EGTA, 0.3 CaCl2, 10 HEPES, 2 Mg-ATP, 1 Na-GTP (free Mg2+, 2 mm; free Ca2+, 130 nm), pH 7.1 (adjusted with CsOH), 290 mOsm; S8, standard blocking solution containing 5% normal bovine serum (Dianova), 0.8% Triton X-100, and 0.05% NaN3 in Ca2+/Mg2+-free PBS (PBS−/−; 100 mm); S9, standard staining solution containing 3% bovine serum albumin (IgG-free, protease-free), 0.05% NaN3, and Alexa Fluor 488 or 633 streptavidin conjugate (1:800; Life Technologies) in PBS−/− (100 mm).

Free Ca2+ concentrations were calculated using WEBMAXC (http://web.stanford.edu/∼cpatton/webmaxcE.htm). If not stated otherwise, chemicals were purchased from Sigma-Aldrich. Alexa Fluor hydrazide was purchased from Life Technologies; ω-agatoxin IVA and SNX 482 were purchased from Tocris Bioscience; paxilline was purchased from Enzo Life Science; 2-(3-carboxypropyl)-3-amino-6-(4 methoxyphenyl)pyridazinium bromide (gabazine), D-(-)-2-amino-5-phosphonopentanoic acid (AP5), 2,3-dioxo-6-nitro-1,2,3,4-tetrahydrobenzo[f]quinoxaline-7-sulfonamide (NBQX), ω-conotoxin GVIA, and 4-ethylphenylamino-1,2-dimethyl-6-methylaminopyrimidinium chloride (ZD7288) were purchased from Abcam. Final solvent concentrations were ≤0.1%. Solutions and pharmacological agents were applied either by bath or from air pressure-driven reservoirs via an eight-in-one multibarrel “perfusion pencil” (Science Products). Changes in focal superfusion (Veitinger et al., 2011) were software-controlled and, if required, synchronized with data acquisition by transistor–transistor logic input to 12 V DC solenoid valves using a TIB 14S digital output trigger interface (HEKA Elektronik).

Slice preparation.

Mice were killed by brief exposure to a CO2 atmosphere and decapitation. The left and right olfactory bulbs were rapidly removed, separated with a razor blade, embedded in 4% low-gelling temperature agarose (VWR International), and placed in ice-cold oxygenated cutting solution (S3). Parasagittal slices (250 μm) were cut with a VT1000S vibratome (Leica Biosystems). Two slices per bulb, each including the AOB, were transferred to a submerged, oxygenated storage container and allowed to recover for ≥1 h in aCSF (S2). Slices were then stored at RT until use.

Histology and immunochemistry.

For immunohistochemical labeling of AOB sections, olfactory bulbs were fixed in 4% (w/v) paraformaldehyde (PFA) in PBS−/− [10 mm, pH 7.4 (3 h, 4°C)] and subsequently cryoprotected in PBS−/− containing 30% sucrose (≥24 h, 4°C). Samples were embedded in Tissue Freezing Medium (Leica Biosystems), sectioned at 20 μm on a Leica CM1950 cryostat (Leica Biosystems), and mounted on Superfrost Plus slides (Menzel). Next, sections were incubated in 4% PFA (15 min, 4°C), washed several times in PBS−/−, and then incubated (1 h, RT) in PBS−/−-based blocking solution containing 0.8% Triton X-100, 0.05% NaN3, and normal bovine serum. Primary antibodies were anti-Slo1 (KCa1.1, 1:500; Alomone Labs), anti-Sloβ4 (1:100; Alomone Labs), and anti-CaV2.3 (1:100; Alomone Labs). Cryosections were incubated with primary antibodies (19 h, 4°C) in a dark humidified chamber. After washing in PBS−/− (3×, 10 min), sections were incubated in PBS−/− (1 h, RT) containing Alexa Fluor-conjugated secondary antibodies (1:500; Life Technologies). Excess antibodies were removed by washing (5×, 10 min). Fluorescent images were taken using an upright fixed-stage scanning confocal microscope (Leica TCS SP5 DM6000 CFS, Leica Microsystems) equipped with a 20× 1.0 numerical aperture water-immersion objective (HCX APO L, Leica Microsystems). To control for nonspecific staining and to demonstrate antibody specificity, we performed (1) antigen preadsorption controls and (2) experiments in which the primary antibodies were omitted in parallel with each procedure. Nuclei were counterstained using 1, 5-bis{[2-(di-methylamino)ethyl]amino}-4, 8-dihydroxyanthracene-9, 10-dione (DRAQ5, 1:200; ThermoFisher Scientific). Digital images were uniformly adjusted for brightness and contrast using Adobe Photoshop CS6 (Adobe Systems).

For immunoblotting, both AOB and control tissue were homogenized in lysis buffer (100 μl, 0.1% Triton X-100, 4°C) in the presence of Complete Mini protease inhibitor mixture tablets (Roche). The homogenate was sonicated and centrifuged for 10 min (4°C, 1000 × g). The supernatant was resuspended in lysis buffer and protein concentration was determined (BioPhotometer Plus, Eppendorf). Thirty microliters of Laemmli buffer (20% glycerol, 4% SDS, 125 mm Tris-HCl, and 0.02% bromphenol blue, pH 6.8) was added and equal amounts of protein were fractionated by SDS-PAGE. Separated proteins were transferred using a Bio-Rad Criterion Blotter wet-blotting system. Membranes were washed, stained with Ponceau S to control for protein transfer, and again washed with TBST (61 mm Tris-HCl, 88 mm NaCl, 0.1% Tween 20, pH 7.5). Blocking was performed in 5% non-fat dry milk/TBST overnight. Blots were then incubated in 2.5% non-fat dry milk/TBST (1 h, 4°C) with anti-Slo1 (1:200; Alomone Labs) and anti-Sloβ4 (1:50; Alomone Labs) antibodies. Membranes were then washed (4 × 15 min) in TBST and incubated with horseradish peroxidase-conjugated goat anti-rabbit IgG (Bio-Rad; 1 h, RT, 1:5000) in 2.5% non-fat dry milk/TBST. Blots were again washed in TBST (4 × 15 min) and antibody binding was detected using 2.5 ml of Lumi-Light solution (Roche).

For post hoc visualization, 3D reconstruction, and morphometric analysis of AOB mitral cells diffusion-loaded (≥20 min) with biocytin during whole-cell patch-clamp recordings, slices were kept for another 15 min in the recording chamber to wash out excess biocytin (Marx et al., 2012). AOB slices were then fixed in 4% PFA in PBS−/− [0.1 m; pH 7.4 (4°C, ≥12 h)]. Next, slices were washed (4 × 10 min) in PBS−/− and incubated in blocking/permeabilization solution (S8; 4°C, 90 min). Subsequently, biocytin-filled neurons were stained in a dark humidified chamber (S9; 60 min, RT) with fluorophore-conjugated streptavidin. After washing in PBS−/− (4 × 10 min), slices were mounted on slides and coverslipped. Confocal fluorescent z-stack images were taken using an upright fixed-stage scanning confocal microscope (TCS SP5 DM6000CFS, Leica Microsystems) equipped with a 20×/1.0 numerical aperture water-immersion objective (HCX APO L, Leica Microsystems). Alexa Fluor 488 or 633 streptavidin conjugate were excited using the 488 nm line of an argon laser or a 633 nm HeNe laser, respectively. Rendering of 3D data and morphometric measurements were performed using Imaris 8.0 software (Bitplane).

In vivo recordings.

Electrophysiological recordings of AOB neurons were performed as previously described in detail (Ben-Shaul et al., 2010; Cichy et al., 2015). Mice were anesthetized with 100 mg/kg ketamine and 10 mg/kg xylazine and a tracheotomy was performed using a polyethylene tube (inner diameter, 0.76 mm; outer diameter, 1.22 mm). Though not used during the recordings presented here, a cuff electrode was placed on the sympathetic nerve trunk. Incisions were closed and the mouse was placed in a custom-built stereotaxic apparatus where anesthesia was maintained throughout the entire experiment (0.5–1% isoflurane in O2). A craniotomy was made immediately rostral to the rhinal sinus, the dura was removed around the penetration site, and electrophysiological probes were advanced into the AOB using an electronic micromanipulator (MP-285, Sutter Instruments). All recordings were made with 32-channel multisite electrodes (A4x8-5mm-100-200-5 177, NeuroNexus Technologies). Spontaneous neuronal activity was recorded without vomeronasal stimulation for 3–10 min. Before each recording, electrodes were dipped in fluorescent dye (DiI; Invitrogen) and targeting of the AOB external cellular layer was confirmed from electrode tracts post mortem. During the experiment, large well isolated spikes were found exclusively when electrode tracks traversed the AOB external cellular layer. By contrast, negligible neural activity was observed from electrode tracks located in the granule cell layer. Recordings from MOB units were identified by both electrode tracks and signal–response patterns, and were thus removed from the analysis. Using an RZ2 processor, PZ2 preamplifier, and two RA16CH head-stage amplifiers (Tucker-Davis Technologies), neuronal activity was sampled at 25 kHz and bandpass filtered at 0.3–5 kHz. Candidate spike events were derived from continuous electrode data using the Matlab thselect function with a minimaxi criterion. The threshold was set conservatively to avoid loss of candidate spikes. Because each shank of the probe contains eight recording sites, each spike event is defined by eight simultaneously measured waveforms (spanning 3.5 ms; detected at any site). Then, for each recording site, we calculated the first two principal components of all recorded waveforms. Each spike event was then characterized according to its projections on the first two principal components of each of the eight channels per shank. Thus, each spike waveform was defined by 16 numbers. These numbers were then fed to the automatic spike-sorting algorithm KlustaKwik (Harris et al., 2000), which assigned each event to one of several clusters. Clusters were then manually examined and adjusted using the Klusters program (Hazan et al., 2006). This latter control procedure is required since the KlustaKwik function tends to overseparate clusters. Criteria for cluster evaluation included the spike waveforms themselves, their projections on the principal component space, and the interspike interval (ISI) histogram. Spike clusters were defined as single units if (1) they had a distinct spike shape, (2) they were fully separated from both the origin (noise) and other clusters in ≥1 principal component projection, and (3) their ISI histogram demonstrated a clear trough of ≥10 ms duration around time 0. If a cluster was composed of >1 single unit, it was designated as multiunit activity and excluded from the analysis. All spike-sorting procedures (except KlustaKwik and Klusters) were implemented using custom-written Matlab code.

In vitro electrophysiology.

Olfactory bulb slices were transferred to a recording chamber (Luigs & Neumann), positioned with stainless steel anchors, and visualized using an upright fixed-stage video-microscope (DM6000FS, Leica Microsystems) equipped for infrared-optimized differential interference contrast. Slices were continuously superfused with oxygenated S2 (∼3 ml/min, gravity flow, 25°C). Bath temperature was measured and controlled using the temperature controller TC07 and the PTC Minibath Chamber IV Plus Upgrade kit (Luigs & Neumann). Neurons were visualized using a 5× (N Plan 5x/0.12) and 25× (HCX IRAPO L25x/0.95W) objective, a three-position magnification changer (0.35×, 1.25×, and 4.0×) and a cooled CCD camera (DFC360FX, Leica Microsystems). Patch pipettes (5–8 MΩ) were pulled from borosilicate glass capillaries (outer diameter, 1.50 mm; inner diameter, 0.86 mm; Science Products) on a PC-10 micropipette puller (Narishige Instruments), fire-polished (MF-830 Microforge, Narishige Instruments), and filled with pipette solution (S5–S7, depending on experimental design). Alexa Fluor 488 hydrazide (20 μm) and, in some recordings, biocytin [0.3% (w/v)] was routinely added to the pipette solution to enable on-line evaluation of cell morphology and post hoc 3D reconstruction of recorded neurons, respectively. Neither chemical showed an evident effect on mitral cell electrophysiology. An agar bridge (150 mm KCl) connected the reference electrode and bath solution. An EPC-10 USB amplifier controlled by Patchmaster 2.67 software (HEKA Elektronik) was used for data acquisition. We monitored and compensated pipette and membrane capacitance (Cmem) as well as series resistance. Only neurons exhibiting relatively low (<30 MΩ) and stable access resistances were used for analysis. Liquid junction potentials were calculated using JPCalcW software (Barry, 1994) and corrected on-line. Signals were low-pass filtered [analog three-pole and four-pole Bessel filters (−3 dB); adjusted to one-quarter to one-fifth of the sampling rate (10 kHz)]. If not stated otherwise, holding potential (Vhold) was −70 mV. All electrophysiological data were recorded at RT. Mitral cells were identified according to their location [residing in the external cellular layer between the AOB glomerular layer and the lateral olfactory tract (Larriva-Sahd, 2008)], soma size (large somata; average Cmem, 14.5 ± 0.4 pF), and dendritic morphology (multiple apical/primary dendrites that terminate as tufts in the glomerular layer). “Loose-patch” recordings were performed from intact mitral cell somata to prevent dialysis of intracellular components. Action potential-driven capacitive currents were recorded in loose-seal cell-attached configuration (seal resistance, 30–150 MΩ; pipettes filled with S1). Spikes were analyzed using Igor Pro functions (SpAcAn, written by Guillaume Dugué and Charly Rousseau) for detection and analysis of spontaneous events by a threshold detection algorithm. Passive membrane properties [i.e., input resistance (Rinput), Cmem, and membrane time constant (τmem)] were obtained immediately after membrane rupture. Treated, to a first approximation, as a “biological constant” with a value of ∼1 μF/cm2 (Gentet et al., 2000), Cmem was determined using a square-pulse (5 mV, 10 ms) routine. Rinput at the mitral cell soma was determined by measuring the steady-state voltage response to a hyperpolarizing current step of −20 pA. Linear passive voltage responses were also used to estimate τmem from monoexponential fits to the voltage responses (from onset to steady state). General action potential parameters were calculated from averaged spike waveforms. Spike amplitude was measured as the threshold-to-peak distance, spike duration was calculated as the full duration at half-maximum (FDHM), and spike-generating kinetics were measured as the time-to-peak. All electrophysiological in vitro experiments used slices from young adults of either sex. We did not observe obvious gender-dependent differences.

Modeling.

The model design and specification was set up in neuroConstruct software (Gleeson et al., 2007). neuroConstruct facilitates the building and validation of models conforming to the NeuroML specification (www.neuroml.org). Model scripts were generated at multiple stages during the development process by neuroConstruct to run simulations using the Neuron simulator software (Hines and Carnevale, 1997, 2001). Simulator software contains solvers that run all the differential equations that are spatially distributed within the morphological framework of the model at a specified time step. It also allows the output of these differential equations to be recorded at specified locations and for specified parameters.

For model development, the digitized 3D reconstruction of a representative intrinsically oscillating AOB mitral cell was chosen. Morphology parameters were saved as a .hoc file and imported into neuroConstruct. Compartments were separated into groups (soma, dendrites, axon) for specifying channel allocation.

The model's specific passive parameters were tuned to target mean experimental measurements/estimates, such as somatic Rinput, τmem, and Cmem. Based on simulated somatic current injections (70 pA, 400 ms), Rinput was calculated as follows: ΔVmem/Iinject, with Vmem defined as cell membrane potential and Iinject defined as injection current. τmem Was estimated by monoexponential fits using the Octave expfit function (http://octave.sourceforge.net/). Then, a simulated voltage-clamp test pulse (5 ms, 5 mV, 20 MΩ) was applied to the soma and Cmem was estimated by fitting an exponential to the capacitive transient (Octave expfit). Iterations were repeated until passive model parameters replicated experimental values. The passive parameters were as follows: Cmem = 0.79 μF cm−2; Rmem = 59,524 Ω μm−2; and Rinput = 560 Ωcm.

New models for outward-rectifying K+-channel current (IKdr), transient K+-channel current (IKA), persistent Na+ current (INaP), and voltage-gated Ca2+ current (ICaV) were developed from current–voltage (I–V) plots, steady-state plots, and time-course datasets for activation, inactivation, and deactivation measured in voltage-clamp experiments. These models were used in combination with previously published models for transient Na+ current (INaT; Migliore et al., 2005) and big conductance KCa (BK) current (IBK; Maex and De Schutter, 1998). Ionic reversal potentials were (in millivolts) as follows: Eleak = −74.1, ENa = 67, EK = −86.5, and ECa = 80. All equations and channel densities (Gmax) for each model are found at https://github.com/Simon-at-Ely/Channel_Kinetics/blob/master/AccessoryOlfactoryBulb/MitralCell/ipython/MultipleChannelModel/iAMT_model_specification.ipynb (any tuning adjustments have been highlighted in red or blue).

Data analysis.

All in vitro data were obtained from independent experiments performed on ≥3 d using ≥3 different animals. Individual numbers of cells/experiments (n) are denoted in the figure and/or captions. If not stated otherwise, results are presented as means ± SEM. Statistical analyses were performed using paired or unpaired t tests, one-way ANOVA with Tukey's HSD post hoc test, Mann–Whitney U tests, Wilcoxon signed-rank tests, or Kruskal–Wallis tests (as dictated by data distribution and experimental design). Tests and corresponding p values that report statistical significance (≤0.05) are individually specified in the figure legends. Data were analyzed off-line using FitMaster 2.67 (HEKA Elektronik), Imaris 8.0 (Bitplane), IGOR Pro 6.3 (WaveMetrics), Matlab (Mathworks), and Excel 2013 (15.0.4779.1001, Microsoft) software. Activation curves were fitted by the Boltzmann equation to calculate the Vmem of half-maximal activation (V1/2). Time constants (τ) were calculated by fitting individual traces to monoexponential functions I(t) = I1 [exp (−t/τ)] + I0. Linear correlation of two parameters was analyzed by calculating the Pearson correlation coefficient using IGOR Pro's linear correlation procedure. Data fitting for model development was carried out using IPython Notebook (ipython.org/notebook.html) and a least-square fitting method.

To account for the variability in single-unit firing rates obtained from in vivo recordings, unit classification required statistical analysis. Bursting activity was defined statistically as follows. First, for each unit, an expected ISI distribution assuming random independent (Poisson) spiking was derived. Poisson distribution is defined by a single parameter λ, which corresponds to that unit's observed average firing rate, calculated over the entire recording period. From this distribution, the median ISI value for each unit (i.e., intervals smaller or larger than this interval were equally likely) was derived. A burst was defined as a sequence of ≥4 consecutive spikes separated by intervals smaller than the median ISI (Fig. 1Ei–Hi, red horizontal bars). A unit was classified as “nonbursting” if <50% of all spikes occurred in bursts. By contrast, a unit was classified as “bursting” if (1) >50% of all spikes occurred in bursts, and (2) that unit's ISI distribution was significantly different (p < 0.05; two-sample Kolmogorov–Smirnov test) from the expected Poisson distribution under the null hypothesis of random firing at the same mean firing rate. For additional classification of (ir)regularity, spike time autocorrelation histograms (ACHs) were calculated for each unit across multiple time periods: 5, 10, 20, 30, 60, 90, and 120 s. Each period was divided into 200 equal bins (5 s ACH ≙ 25 ms bins, 120 s ACH ≙ 600 ms bins). Based on each unit's autocorrelograms, we determined ACH-specific trough indices. The trough index (TI) was defined as follows: TI = CVNcross, where CV is the coefficient of variation of the autocorrelogram, and Ncross is the number of times that the autocorrelogram crosses its mean value. The TI is designed to capture numerically what is intuitively perceived as an oscillatory ACH (Moore et al., 1966). Thus, the TI attains high values if an autocorrelogram exhibits large fluctuations from its mean that are concentrated in a small number of peaks, a feature that is associated with robust oscillatory activity. By contrast, flat ACHs or autocorrelograms with many deviations from the mean are associated with low TI values. The rationale for examining multiple autocorrelograms per unit was to detect oscillatory activity over a range of time scales. To identify the time scale exhibiting the strongest oscillatory activity, we compared, for each neuron separately, the set of TIs for all autocorrelograms. For comparison across time scales, we used a standardized TI, defined as TI−median(TI)STD(TI). This standardization was performed across all units for each time scale separately. Thus, potential oscillatory activity was classified by a single number, i.e., the maximal standardized TI across all time periods. One caveat to the use of the TI as a classifier of oscillatory activity is that nonstationary firing (e.g., strong transient activity) could result in large TI values. From visual inspection of ISI distributions and ACHs, we found that these cases were associated with small Ncross values. Therefore, we excluded autocorrelograms of Ncross < 9. Based on the above parameters, bursting units were categorized as oscillatory/regular if their maximal standardized TI was >0.2.

In addition to each unit's classification as either nonbursting, irregular bursting, or regular bursting (i.e., oscillating), oscillatory firing patterns were further characterized by various statistics. These included the number and duration of bursts as well as the number and frequency of spikes within bursts.

Results

Slow oscillatory bursting of AOB neurons in vivo

Brain circuit computations are not only determined by sensory input, but also by intrinsically generated spatiotemporally structured patterns of spontaneous activity (Romano et al., 2015), both at the single-neuron and network levels. As information about spontaneous AOB activity is lacking, we performed in vivo recordings from the AOB mitral cell layer of anesthetized mice (Fig. 1A–C). We continuously monitored spontaneous multiunit activity of AOB neurons and validated putative single units by feature-based clustering (Fig. 1C). Inspection of spike trains over time revealed large firing-pattern heterogeneity in AOB neurons, ranging from apparently random (Fig. 1Ci) to more regular firing (Fig. 1Ciii).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Firing patterns of AOB neurons in vivo. A, Schematic drawing illustrating the multisite recording electrode arrangement in the mouse AOB. B, Digital photomicrograph reconstructing the location of a DiI-labeled 32-site recording electrode from a single fluorescent tract (red) in a 30 μm sagittal section of the olfactory bulb (nuclear staining; blue; DAPI). Dashed white line delimits the AOB. The electrode was accurately targeted to the AOB mitral cell layer. Ci, Ciii, Example extracellular recordings of spontaneous activity (10 min duration; 0.3–5 kHz bandpass filtered) in the AOB of anesthetized mice. Red boxes illustrate 20 s time windows shown at an enlarged scale below. Cii, Civ, Off-line spike waveform sorting isolated single-unit activity from two presumptive mitral cells with irregular and bursting discharge, respectively, clustered from the multiunit activity recordings displayed in Ci and Ciii. Ensembles of overlaid spike waveforms (black) and average spike shapes (red) are shown. D–H, Raster plots of spike times (Di, Ei, Fi, Gi, Hi) recorded from five individual representative AOB neurons in vivo. The corresponding ISI distribution (Dii, Eii, Fii, Gii, Hii) and autocorrelation (Diii, Eiii, Fiii, Giii, Hiii) histograms enable the classification of neurons as either irregular nonbursting (D; corresponding to Ci, Cii), irregular bursting (E), or oscillating (F–H; F corresponds to Ciii, Civ). Dii, Eii, Fii, Gii, Hii, ISI histograms superimposed with the probability (red curve) and cumulative probability (gray curve) functions, corresponding to random Poisson spiking with the same mean rates. The median ISI interval (ISI0.5) is indicated by the light red vertical line. Note that regularity [or the lack thereof (D, E)] is mirrored by the shape of the ACH (F–H; 200 ms bin width). Insets, Autocorrelograms (3 s duration; 50 ms bin width) show narrow initial peaks close to zero (E–G), indicative of bursting units. Red asterisks (*) and boxes mark the time windows of either spike times (Di, Ei, Fi, Gi, Hi) or autocorrelograms (Diii, Eiii, Fiii, Giii, Hiii) displayed at an expanded time scale; red horizontal bars in Ei, Fi, Gi, and Hi denote epochs detected as bursts (see Materials and Methods). I, Scatter dot plots depicting the distribution of IBIs, burst durations, and within-burst firing rates for all bursts recorded from 78 oscillating units in vivo. Average values are shown as means ± SD (red): IBIs, 22.1 ± 22.8 s (median, 14.0 s); durations, 2.9 ± 2.7 s (median, 2.0 s); firing rates, 5.1 ± 4.0 Hz (median, 4.0 Hz). J, Histogram plots of various descriptors of firing behavior. Diagrams depict the percentage of a unit's spikes that occurred within a burst (Ji), the deviation of a unit's ISI distribution from random firing (Jii), and the relation of both parameters to a unit's discharge regularity as indicated by its TI (Jiii; 3D scatter plot). In Jiii, units classified as oscillating are shown in red. Note the skewed distribution of burstiness (Ji; single peak fitting of a Gaussian function reveals a population of units with ≥50% spikes in bursts; red bars), the overrepresentation of both large and small p values (Jii; red bar denotes p < 0.05), and the cluster of red dots within the 3D parameter space (Jiii).

Based on these initial findings, we next quantified and classified the observed patterns (Fig. 1D–H). Consideration of multiple measures, including deviation from Poissonian spiking, burst-firing parameters, autocorrelation analysis, and each unit's TI (see Materials and Methods), revealed two different discharge categories: 339 of 466 (72.8%) units were classified as irregular nonbursting (Fig. 1D), whereas 127 of 466 (27.2%) units were classified as bursting (Fig. 1E–H). Among the bursting subpopulation of AOB neurons, 49 (10.5%) displayed irregular burst firing characterized by an exponential decay of the ISI histogram and the lack of pronounced autocorrelogram “side” peaks. These units were classified as irregular bursting neurons (Fig. 1E). By contrast, 78 (16.7%) bursting neurons showed highly rhythmic oscillatory discharge with prominent ISI and ACH peaks, corresponding to large TI values (see Materials and Methods). These neurons were classified as oscillating (Fig. 1F–H). Within this subpopulation of oscillating AOB neurons, individual discharge patterns were heterogeneous and, accordingly, interburst intervals (IBIs), burst durations, and firing rates were highly variable (Fig. 1I). IBIs ranged from 1.4 to 112.2 s and bursts lasted for ≤15.0 s. The non-normal distribution of these values indicates the absence of a single characteristic rhythm governing AOB unit oscillations. In general, the activity of oscillating neurons was markedly increased compared with previously observed low baseline rates (1–2 Hz) of irregularly firing neurons (Luo et al., 2003; Ben-Shaul et al., 2010). To identify whether the characteristics of oscillating units define a distinct subpopulation of AOB neurons or they instead represent one extreme of a continuous and normally distributed dataset, we plotted histograms of burstiness (Fig. 1Ji), deviation from Poissonian spiking (Fig. 1Jii), and regularity (Fig. 1Jiii) for the entire sample population (n = 466). Our analyses indicate that oscillating AOB neurons indeed form a distinct group. Together, these results demonstrate that rhythmic patterns of slow oscillatory activity characterize a subpopulation of AOB neurons in vivo.

A subpopulation of AOB mitral cells displays slow spontaneous oscillations

AOB neuron oscillatory discharge may result from interactions within an intact network and potential top-down modulation or, alternatively, from intrinsic rhythmogenic properties of individual neurons. To resolve the mechanistic basis of mitral cell rhythmicity and avoid the experimental drawbacks that complicate the biophysical interpretation of in vivo whole-cell patch-clamp data (Maier et al., 2011), we turned to an in vitro model. We recorded spontaneous activity from individual AOB mitral cells in sagittal sections of the mouse olfactory bulb (Fig. 2A). When continuously monitoring mitral cell Vmem for prolonged periods of time under control conditions (0 pA current injection; S2, S5; see Materials and Methods), the vast majority of neurons generated spontaneous discharge [308 of 324 cells; resting Vmem (Vrest) = −73.2 ± 0.6 mV; firing rate, 2.9 ± 0.4 Hz; means ± SEM]. Similar to our in vivo findings, AOB mitral cells displayed one of two distinct activity patterns (Fig. 2B): either irregular firing with no apparent periodicity (irregular discharge; 105 of 308 cells; Vrest = −71.9 ± 0.7 mV; Fig. 2Bi) or “phasic” firing patterns with alternating periods of activity and silence (oscillatory discharge; 203 of 308 cells; Vrest = −74.1 ± 0.8 mV; Fig. 2Bii). These slow oscillations of recurring up and down states (−63.9 ± 2.5 vs −74.5 ± 3.7 mV, mean ± SD) usually remained stable throughout the recording (≤60 min). Typically, bursts of action potentials were superimposed on the slow depolarizing envelope (burst duration, 6.0 ± 5.3 s; IBI, 11.3 ± 7.6 s; within-burst firing rate, 4.2 ± 2.3 Hz, mean ± SD). Essentially identical patterns of spontaneous activity were observed when we recorded mitral cell activity in “loose-seal” cell-attached configuration (n = 39) to prevent dialysis of cytosolic components and maintain unperturbed Vrest (Fig. 2Biii). As a measure of regularity and bursting behavior (Moore et al., 1966), we inspected the ISI distribution and the ACHs (Fig. 2C–F). Figure 2C exemplifies an irregular nonbursting neuron, whereas oscillatory mitral cells are depicted in Figure 2D–F. While discharge periodicity for individual cells was strong, patterns were heterogeneous across the population, revealing no predominant rhythm. We then asked whether these two patterns of spontaneous activity correlate with distinct morphological phenotypes (Fig. 2G). Volume-rendered 3D reconstructions of individual biocytin-filled mitral cells (Fig. 2Gi,Gii) revealed large polymorphic somata and several branched primary dendrites that terminate as multiple tufts within the homonymous AOB half (Wagner et al., 2006). However, morphometric analysis documented no obvious morphological differences between irregularly firing and oscillating neurons (Fig. 2Giii). Together, these data demonstrate that, both under in vivo and in vitro conditions, infraslow oscillatory discharge represents the default activity pattern of a substantial population of AOB mitral cells.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Patterns of spontaneous activity in AOB mitral cells in vitro. A, Photomicrographs illustrating the experimental approach. Top, Left, Parasagittal section (250 μm) of the mouse rostral brain, including the AOB (dotted line). Bottom, Left, Infrared differential interference contrast image showing the layered structure of the AOB as well as typical experimental configurations of patch electrode and perfusion pencil. Granule cell layer (GCL), glomerular layer (GL), and mitral cell layer (MCL) are indicated. Top, Right, Infrared differential interference contrast photomicrograph showing a region of the MCL at higher magnification. Dotted lines indicate the patch pipette targeting a mitral cell (MC) soma. Bottom, Right, Maximum projection of a confocal z-stack illustrates four MCs that were loaded with Alexa Fluor 488 during patch-clamp recordings (see Materials and Methods). B, Representative whole-cell current-clamp recordings (S2/S5; no pharmacological perturbation of network activity) illustrating the two distinct types of spontaneous discharge. Mitral cells either fire irregularly (Bi), or exhibit oscillatory discharge patterns (Bii). Rhythmicity (or the lack thereof) is also evident from either one or two peaks in the corresponding Vmem all-points histograms (right). Prominent oscillations are also observed in extracellular loose-seal recordings (Biii). C, Original representative trace (Ci), ISI distribution (Cii), and autocorrelation (Ciii) histogram of an irregularly firing MC. In the ISI histogram (Cii), the gray curve depicts the cumulative probability function and the light red vertical line illustrates the interval corresponding to a cumulative probability of 0.5 (ISI0.5). The red exponential curve represents the “theoretical” stochastic ISI distribution of a Poissonian spike train with a uniform interval probability of ISI0.5. Both the uniform shape of the ACH [200 ms bin width (Ciii)] and the lack of a peak near zero indicate irregular nonbursting firing. Inset, Autocorrelogram (3 s duration, 50 ms bin width) confirming the lack of a narrow initial peak. Red asterisk (*) and box denote the time window displayed at an enlarged time scale. D–F, Representative original recordings (Di, Ei, Fi), ISI histograms (Dii, Eii, Fii) and autocorrelograms (Diii, Eiii, Fiii) of MC oscillatory discharge reveal distinct types of autorhythmicity. Cumulative ISI probability curve, ISI0.5, and Poisson-type ISI curve are derived as in Cii. Regularity of bursting and its diversity is mirrored by the shape of the autocorrelograms (Diii, Eiii, Fiii) with a narrow initial peak close to zero and several “side” peaks (vs the peak-less autocorrelogram indicative of irregular firing in C). Red asterisks (*) and dashed boxes correspond to the initial 3 s time segments of the autocorrelograms that are displayed at an enlarged scale (50 ms bin width; insets). G, Representative confocal z-stack images (1.14 μm intervals; maximum intensity projection) and volume-rendered 3D reconstructions of an irregularly firing (Gi) and an oscillating (Gii) MC, respectively. Reconstructed somata are depicted in cyan; dendrites are shown in green. Morphometric analysis (Giii) included soma diameter [20.6 ± 0.6 μm (irregular) vs 18.7 ± 1.0 μm (oscillatory)], numbers of branching points [5.9 ± 0.6 (irregular) vs 5.8 ± 0.7 (oscillatory)], and dendrites [3.9 ± 0.3 (irregular) vs 3.7 ± 0.3 (oscillatory)], as well as total surface [8503 ± 847 μm2 (irregular) vs 7828 ± 1028 μm2 (oscillatory)]. Bar graphs show means ± SEM. Numbers of MCs are indicated above bars. Morphological parameters are independent of discharge type (p > 0.05; unpaired t test).

A group of AOB mitral cells are intrinsically rhythmogenic

Neural rhythmogenesis is either a consequence of network activity and thus results from regular recurrence of excitatory and inhibitory synaptic barrages (Crunelli and Hughes, 2010) or, alternatively, oscillations are generated intrinsically by the repetitive, rhythmic discharge of pacemaker-like neurons (Blankenship and Feller, 2010). To distinguish between these mechanisms, we tested whether spontaneous oscillations in AOB mitral cells depend on synaptic drive. Pharmacological inhibition of GABAergic synaptic transmission did not qualitatively affect rhythmic discharge (Fig. 3A). Furthermore, when both GABAergic and glutamatergic fast synaptic transmission were blocked (Fig. 3B,C), stable oscillations persisted in a substantial fraction of AOB mitral cells (43 of 185 cells; 23.2%). While selective block of inhibitory input altered some oscillation parameters (Fig. 3D), complete isolation from fast synaptic transmission did not change oscillatory patterns in those AOB neurons that proved resistant to pharmacological treatment (Fig. 3E). Thus, we next examined whether these neurons shared another hallmark of pacemakers, i.e., a positive causal correlation between oscillation frequency and “baseline” Vmem (Crunelli and Hughes, 2010). Similar to intrinsic oscillators described in other circuits (Hayar et al., 2004; Blethyn et al., 2006; Tazerart et al., 2008), oscillation frequency changed as a function of depolarizing or hyperpolarizing current injection in synaptically isolated oscillating AOB mitral cells (Fig. 3F–H). Hyperpolarization increased, whereas depolarization reduced IBIs (Fig. 3F,Hi). Moreover, each neuron (nine of nine cells) exhibited a characteristic Vmem threshold below which the pattern of periodically recurring up and down states switched to a stable resting state (Fig. 3F,G). Notably, within the voltage range from rest to firing threshold, neither positive nor negative current injections changed a given neuron's general discharge type (i.e., from oscillatory to irregular or vice versa; n = 10; Fig. 3Fii). Together, these results suggest that the mouse AOB contains a group of intrinsically rhythmogenic mitral cells that generate slow Vmem oscillations independent of fast synaptic input.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Rhythmic discharge of AOB mitral cells. A, B, Extracellular loose-seal recordings from two oscillating mitral cells under control conditions as well as during inhibition of either GABAergic [A; gabazine (10 μm)] or glutamatergic [B; AP5 (100 μm) + NBQX (10 μm)] fast synaptic transmission, respectively. Horizontal red bars indicate drug incubation. Bottom traces, Expanded view of two consecutive bursts during either pharmacological treatment. Ci, Original whole-cell recordings (Vhold = −65 mV; S2/S6; ECl = ∼0 mV) of spontaneous postsynaptic currents (spPSCs). Traces exemplify spPSCs as downward deflections of varying amplitudes under control conditions (I; red arrowheads indicate relatively large events). Moreover, recordings show the effects of inhibiting ionotropic glutamate receptors [II; AP5 (100 μm) + NBQX (10 μm)], GABAA receptors [III; gabazine (10 μm)], and fast synaptic transmission in general (IV; gabazine + AP5 + NBQX; drug preincubation, ≥10 s). Cii, A histogram of spPSC amplitude distribution illustrates that both AP5/NBQX (II) and gabazine (III) reduce spPSC amplitude and frequency. In combination, the three drugs block essentially all remaining spPSCs (IV). Di, Ei, Representative Vmem recordings from two mitral cells under control conditions as well as during either inhibition of GABAA receptors [D; gabazine (10 μm)] or a complete block of fast synaptic transmission [E; gabazine (10 μm) + AP5 (100 μm) + NBQX (10 μm)]. Note that oscillatory discharge persists when this mitral cell is isolated from both excitatory and inhibitory synaptic input. Horizontal red bars indicate drug incubation. Dii, Eii, Bar charts comparing oscillatory discharge parameters from rhythmogenic AOB mitral cells under control conditions (ctr; black) versus inhibition of synaptic transmission (red; gabazine, n = 18; synaptic isolation (iso) by gabazine + AP5 + NBQX, n = 12). Both normalized IBIs and burst durations as well as within-burst firing rates and membrane up and down states (Vd/Vu) are plotted as means ± SEM. Block of GABAergic fast synaptic input prolonged IBIs (2.4-fold; ±0.6) and hyperpolarized the down state (−74.9 ± 1.1 vs −76.6 ± 0.8 mV). Asterisks (*) denote statistical significance, p < 0.01 (Wilcoxon signed-rank test). By contrast, burst durations (1.3-fold; ± 0.3), within-burst firing frequencies (4.0 ± 0.8 vs 5.4 ± 0.9 Hz), and up-state potentials (−64.3 ± 0.7 vs −63.4 ± 0.9 mV) were essentially unaltered. Moreover, no parameter was changed in absence of synaptic drive (p > 0.05; Wilcoxon signed-rank test): IBIs, 1.3-fold; ± 0.2 s; burst durations, 1.1-fold; ± 0.1 s; firing frequencies, 4.3 ± 0.8 vs 4.3 ± 1.0 Hz; Vd, −75.1 ± 1.1 vs −74.1 ± 1.1 mV; Vu; −64.3 ± 0.7 vs −63.9 ± 0.9 mV. F, Original whole-cell current-clamp recordings from two representative oscillating AOB mitral cells during continuous depolarizing or hyperpolarizing current injections of variable amplitude [−15 to 15 pA (Fi)/50 pA (Fii)]. Gi, Expanded view of the current-clamp recordings shown in Fi (dashed red rectangle) illustrating the “switch” from a stable Vmem state (top trace) to slow periodic bursting (bottom trace). Red arrowheads indicate the (sub)threshold (Vt) and down-state (Vd) Vmem values, respectively, for this particular mitral cell. Gii, Quantification of Vt versus Vd values. Data points correspond to baseline Vmem measurements in absence (Vt) and presence (Vd) of oscillating discharge, respectively. Data from individual mitral cells are connected by black lines (n = 10). Compared with Vd (−77.4 ± 1.8 mV), average Vt (−80.3 ± 1.3 mV) is significantly more hyperpolarized (means ± SEM; red), indicating that the mechanisms underlying oscillatory discharge operate in a Vmem range more positive than −80 mV. Asterisk (*) denotes statistical significance; p < 0.0001 (paired-sample t test). Hi, Hii, Negative (Hi) and positive (Hii) correlation of current injection amplitude with IBIs and burst durations, respectively (n = 7). Data are normalized to individual maxima.

Pattern variability among intrinsically oscillating mitral cells

Since pacemaker-like neurons can exert profound effects on coding and computation in sensory systems (Izhikevich et al., 2003; Mizuseki et al., 2009), we focused on the subpopulation of intrinsically oscillating AOB mitral cells (iAMCs). Accordingly, we performed subsequent experiments under tonic synaptic inhibition (gabazine, AP5, NBQX). First, we asked whether iAMCs and irregularly discharging neurons differ in their passive membrane properties and/or their spike-generation machinery. Neither Cmem nor τmem differed significantly (Fig. 4Ai,Aii), suggesting that basic biophysical properties do not distinguish iAMCs from the “general” AOB mitral cell population. Rinput, however, was slightly increased in iAMCs (Fig. 4Aiii). Comparison of mean instantaneous spike frequencies as a function of stationary current input (f–I curve; 5–100 pA; Fig. 4Aiv) revealed indistinguishable curves for both neuron types with response saturation at amplitudes >80 pA and maximum average firing frequencies of ∼25 Hz. Moreover, averaged spike waveforms from iAMCs and irregularly firing neurons were similar in spike amplitude, duration, and kinetics (Fig. 4Av). Together with the results shown in Figure 2G, these data demonstrate that the unique spontaneous activity of iAMCs does not correspond to a readily distinguishable morphological or biophysical phenotype.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Discharge variability among the iAMC population. Ai–Aiii, Bar charts comparing passive membrane properties of irregularly firing (black) versus intrinsically oscillating (red) AOB mitral cells. Analysis of Cmem (15.3 ± 0.6 vs 13.9 ± 0.6 pF) and time constant (τmem; 47.6 ± 2.2 vs 46.4 ± 2.4 ms) reveals no significant (n.s.) differences. By contrast, Rinput is significantly increased in iAMCs (560.2 ± 30 vs 424 ± 25 MΩ). Asterisk (*) denotes statistical significance, p = 0.001 (Mann–Whitney U test). Data are means ± SEM. Numbers of experiments are indicated above bars. Aiv, f–I curve depicting average instantaneous discharge frequencies evoked by stationary current injection (1 s duration; 5 pA intervals; range, 5–100 pA) in oscillating (red; n = 24) and irregularly spiking mitral cells (black; n = 30). Maximum frequencies are 24.8 ± 1.6 Hz (iAMCs) and 24.5 ± 1.0 Hz (irregular), respectively. Individual data points are means ± SEM. Curves are monoexponential fits. Av, Averaged action potential waveform observed in oscillating (red) versus irregularly spiking (black) mitral cells. Mean traces reveal an average amplitude of 84.6 mV (red) versus 93.1 mV (black) and a FDHM of 2.2 ms each. Bi–Biii, Original Vmem recordings (Bi and Bii) and firing rate histograms (Biii) illustrating the heterogeneity of iAMC discharge. Bi, Bii, Bottom traces, Expanded view of the recording period delimited by the dashed rectangles. Red arrowheads indicate Vd and Vu, respectively. Biii, Spontaneous discharge frequency histogram of the two iAMCs shown above [gray; IBI, 1.5 s (Bi); black; IBI, 49.4 s (Bii)] and a neuron exhibiting an intermediate average IBI of 10.0 s (red). C, Quantification of oscillatory discharge parameters. Ci, Current-clamp recording of two consecutive bursts schematically exemplifies the discharge properties subject to analysis (horizontal red lines). Cii, Scatter dot plots of individual measurements (open black circles; n = 148). Average values are shown ± SD (red): IBIs, 10.7 ± 8.1 s (median, 9.3 s); burst durations, 5.4 ± 4.1 s (median, 3.9 s); within-burst firing rates, 4.6 ± 3.1 Hz (median, 3.7 Hz); Vd, −76.3 ± 3.8 mV (median, −76.4 mV); Vu, −64.6 ± 3.8 mV (median, −64.7 mV). Asterisk (*) denotes statistical significance, p < 0.0001 (Wilcoxon signed-rank test). D–F, Representative whole-cell (S2; S5) current-clamp recordings (Di, Ei, Fi), ISI histograms (Dii, Eii, Fii), and autocorrelograms (Diii, Eiii, Fiii) of three intrinsically oscillating neurons. Recordings were performed in isolation from both excitatory and inhibitory fast synaptic drive [gabazine (10 μm) + AP5 (100 μm) + NBQX (10 μm)]. Cumulative ISI probability curves, ISI0.5 values, and Poisson-type ISI curves are derived as described in Figure 2. Regularity is evident in the shape of the autocorrelograms (Diii, Eiii, Fiii). Phenotypic heterogeneity among the iAMC population becomes apparent from the wide distribution of IBIs, which range from ∼2 s (D) to >45 s (F). Diii, Eiii, Fiii, Insets, ACHs (3 s duration; 50 ms bin width); red asterisks (*) and dashed boxes denote “expanded” views. Gi, Gii, Distribution of IBIs (Gi) and burst durations (Gii) among the iAMC population [n = 148; 1 s (Gi) and 0.5 s (Gii) bin width, respectively]. While data points preferentially scatter on the lower end of the scale, their distribution is skewed to the right and is thus not well fit by Gaussian equations (red lines). Hi, Hii, Pairwise iAMC correlation analyses of burst durations with IBIs (Hi) and up-state voltage (Hii), respectively (n = 148). Linear regression indicates no obvious correlation of either parameter pair [Pearson's correlation coefficient, r = 0.37 (Hi) and r = −0.26 (Hii)].

Next, we analyzed discharge variability among the iAMC population. Similar to in vivo observations (Fig. 1I), measures including IBIs, burst durations, and firing rates within a burst varied considerably across neurons and were not normally distributed (Fig. 4B,C). Silent periods between bursts lasted on average 10.7 ± 8.1 s (±SD) with individual IBI values ranging from 1.5 to 49.3 s (Fig. 4B–F). Moreover, we observed bursts as brief as 0.3 s, but also prolonged firing periods of ≤22.2 s (average burst duration, 5.4 ± 3.9 s; mean ± SD; Fig. 4B–F). Increasing the bath solution's temperature from 25°C (RT) to 37°C decreased both IBIs and burst duration (n = 3; data not shown). Distribution analysis of both IBIs and burst durations, however, revealed no obviously segregated subpopulations among iAMCs (Fig. 4G). The vast majority of iAMCs (99.4%) fire >5 spikes per burst, translating into an average within-burst firing rate of 4.6 ± 3.0 Hz (mean ± SD; range, 0.9–18.6 Hz). Notably, the average voltage change during transitions from down to up states is 11.5 mV (−76.1 ± 3.9 to −64.6 ± 4.0 mV; means ± SD), indicating that any current(s) providing the excitatory drive to iAMCs must operate within this voltage regime. To examine whether oscillation parameters are interdependent, we carried out pairwise correlation analyses of IBIs, burst durations, and up-state voltage (Vu), respectively (Fig. 4H). Our results, however, indicate no obvious correlation between any parameter pair, suggesting that the antagonistic mechanisms underlying the recurring transition and maintenance of depolarization and hyperpolarization are phenotypically independent.

Together, these findings show that iAMCs are a physiologically distinct subgroup of AOB mitral cells that share pacemaker-like discharge properties, spanning a wide range of mechanistically independent characteristics within the oscillation parameter space.

INaP promotes rhythmogenesis

The cyclical interplay of depolarizing and hyperpolarizing pacemaker currents drives autorhythmicity and burst generation in a variety of neurons (Grillner, 2006; Blankenship and Feller, 2010). Conceptually, these currents provide either the transitional excitatory drive from down to up state, the translation into a regular firing pattern during the up-state plateau, or the hyperpolarization that terminates the burst (Colwell, 2011). The hyperpolarization-activated current (Ih), the low-voltage-activated T-type Ca2+ current (IT) and INaP are prototypical depolarizing pacemaker currents and, thus, major determinants of autorhythmicity (Crill, 1996; Perez-Reyes et al., 1998; Chan et al., 2004). Ih, in particular, is a predominant driving force of rhythmic oscillatory activity (Maccaferri and McBain, 1996; Liu and Shipley, 2008). Therefore, we tested whether Ih is involved in oscillatory iAMC discharge. The bradycardic agent ZD7288, a specific but isoform-independent blocker of hyperpolarization-activated cyclic nucleotide-gated channels (BoSmith et al., 1993), did not alter iAMC rhythmic activity (Fig. 5A,B). Ih expression typically manifests as a rebound depolarization (sag) at hyperpolarized potentials (Cichy et al., 2015). In current-clamp recordings, we thus measured the hyperpolarization-evoked sag potential in iAMCs under control conditions (Fig. 5C). Plotting the sag potential amplitude as a function of peak hyperpolarization (Fig. 5D) revealed that pronounced voltage sags were only detected upon membrane hyperpolarization to values substantially more negative than the average down-state potential (Vd). Thus, these results indicate that Ih serves a minor, if any, role in iAMC rhythmogenesis.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Ih and IT serve minor roles in iAMC autorhythmicity. A, Representative whole-cell current-clamp recording of oscillatory iAMC discharge before (left) and during (right) ZD7288 treatment (40 μm; horizontal blue bar; preincubation, ≥3 min). Drug efficacy was confirmed in parallel recordings from vomeronasal sensory neurons (Cichy et al., 2015). B, Bar charts comparing iAMC discharge (n = 12) under synaptically isolated conditions (iso; black) versus ZD7288-dependent hyperpolarization-activated cyclic nucleotide-gated channel inhibition (ZD; blue). Bi, Bii, Burst durations (Bi) and firing rates (Bii) are plotted as means ± SEM. Both parameters are essentially unaffected by ZD7288 treatment (p > 0.05; paired sample t test): burst durations (3.8 ± 1.1 s vs 5.5 ± 2.0 s); burst firing frequencies (6.5 ± 1.5 vs 6.2 ± 1.3 Hz). C, Representative whole-cell current-clamp recordings in response to hyperpolarizing current injections (−50 to −500 pA; −50 pA increments; 300 ms duration). Exemplary original traces of Vmem hyperpolarizations evoked by −50, −300, −400, and −500 pA (red lines in pulse protocol; inset) demonstrate the voltage dependence of the sag potential (Vsag) amplitude. No sag was observed in presence of ZD7288 (data not shown). D, Scatter plot depicting the relationship between Vsag and peak hyperpolarization. For quantification, Vsag is determined as the voltage difference between the peak hyperpolarization and the steady-state Vmem (inset). Data points correspond to intrinsically oscillating (gray) and irregularly firing (black) mitral cells (n = 12). Note that considerable hyperpolarization (≤−130 mV) is required to elicit rebound depolarization. Ei, Representative recordings of I–V relationships as determined by fast ascending voltage ramps (inset: −100 to + 80 mV, 100 ms duration; Vmem range shown, −60 to 65 mV) in absence (black) and presence (gray) of Cd2+ (200 μm). The red trace depicts the difference between I–V plots under control conditions and during Cd2+ incubation and thus corresponds to the Cd2+-sensitive Ca2+ current. Eiii, Average I–V relationship of Cd2+-sensitive Ca2+ currents (red trace; n = 18). Gray “shadow” represents SEM. Peak currents are detected at ∼0 mV (−563 ± 66 pA). Fi, Representative recordings in response to step depolarization (from −100 to −30 mV; inset) in absence (black) and presence (gray) of mibefradil (10 μm). Negligible mibefradil-sensitive current is isolated by digital subtraction (violet trace). Fii, Quantification of fractional mibefradil-sensitive current density. When normalized to total peak current, mibefradil treatment blocked 11.2 ± 4.8% (mean ± SEM; n = 4) of inward current. Gi, Representative whole-cell current-clamp recording of oscillatory iAMC discharge before (left) and during (right) mibefradil treatment (10 μm; horizontal violet bar). Gii, Bar diagrams comparing iAMC discharge parameters (n = 4) under synaptically isolated control conditions (iso; black) versus mibefradil treatment (mib; violet). Burst duration, IBI, firing rate, and Vmem values are shown before (black) and during mibefradil incubation (violet). Data are plotted as means ± SD. All parameters are essentially unaffected by mibefradil (p > 0.05; paired sample t test): burst durations, 8.9 ± 2.7 versus 8.5 ± 3.3 s; IBIs, 14.6 ± 4.9 versus 13.7 ± 2.2 s; burst firing frequencies, 3.0 ± 1.3 versus 2.5 ± 0.5 Hz; Vmem: Vd = −80.5 ± 7.7 versus −83.2 ± 8.8 mV; Vu = −71.4 ± 7.9 versus −75.2 ± 10.0 mV.

Another common subthreshold current that is frequently involved in Vmem bistability and burst generation is IT (Bean, 2007; Crunelli and Hughes, 2010). Especially in neurons characterized by relatively hyperpolarized Vrest, these low-threshold currents frequently accelerate the down-state-to-up-state transition (Llinás and Yarom, 1986) by providing additional subthreshold depolarizing drive. Thus, we asked whether iAMCs express low-threshold T-type Ca2+ channels. In fast voltage ramp recordings, Cd2+-sensitive Ca2+ currents first manifested at suprathreshold Vmem (>−40 mV) and I–V curves did not display a characteristic “dent” at subthreshold voltages (Fig. 5E). Moreover, Ca2+ currents were essentially insensitive to mibefradil (Fig. 5F), a drug that preferentially inhibits T-type Ca2+ channels (Martin et al., 2000). Consistent with this observation, rhythmic discharge patterns of oscillating iAMCs remained essentially unchanged upon exposure to mibefradil (Fig. 5G). Thus, our data strongly suggest that iAMCs do not display substantial IT.

As a regenerative depolarizing current between Vrest and spike threshold, TTX-sensitive INaP affects both firing pattern and frequency (Crill, 1996; Tazerart et al., 2008). Thus, we next examined whether TTX, in addition to disrupting action potential discharge, would also affect subthreshold Vmem oscillations. Indeed, incubation with TTX abolished Vmem fluctuations (Fig. 6Ai). In presence of the toxin, iAMCs exhibited stable Vmem values that were statistically indistinguishable from Vd (Fig. 6Aii). Since selective pharmacological block of INaP has proven difficult (Ramirez et al., 2011; Richter et al., 2014), we evoked and isolated TTX-sensitive steady-state currents by slow depolarizing voltage ramps (Carter et al., 2012). A TTX-sensitive current was first evident at ∼−75 mV and increased steeply with voltage (Fig. 6Aiii). As known for INaP (Bean, 2007), maximum steady-state currents were only a small fraction of the maximum INaT (2.2%; n = 7; data not shown). At subthreshold voltages, however, such currents of only tens of picoamperes may prove functionally significant. If so, average Vu and Vd values must be within the voltage-operating range of INaP, which is indeed the case (Fig. 6Aiv).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

INaP is a major determinant of iAMC autorhythmicity. Ai, Representative whole-cell current-clamp recording under isolation from fast synaptic input (gabazine + AP5 + NBQX) in absence and presence of TTX (300 nm; horizontal red bar). Aii, Bar diagram quantifying average Vmem values (n = 13) before (black; Vd = −77.2 ± 0.9 mV; Vu = −65.8 ± 0.8 mV) and after TTX incubation (red; −73.5 ± 0.8 mV). Data are means ± SEM. Asterisks (*) denote statistical significance, p < 0.0001 (Kruskal–Wallis test); n.s., not significantly different. Aiii, Isolation of INaP: representative I–V relationships of average whole-cell currents (3 consecutive recordings) evoked by a slow (10 mV/300 ms) ascending voltage ramp from −93 to −53 mV (inset) in absence (black trace) and presence (gray trace) of TTX (300 nm). When data recorded during drug treatment were digitally subtracted from control recordings and plotted as a function of the command voltage, the TTX-sensitive negative slope current becomes apparent (red trace). Note the TTX-insensitive linear leak current. Aiv, At different [Ca2+]ex, TTX-sensitive currents (revealed by off-line subtraction) were converted to conductances and normalized quasi-steady-state activation curves were fit by the Boltzmann equation (0.1 mm Ca2+) or extrapolated from monoexponential functions (solid lines on top of shaded traces). Dashed horizontal red line denotes a 5% conductance threshold. The gray rectangle indicates the “oscillatory” Vmem range delimited by average Vd and Vu values (Fig. 4Cii). Av, Bar chart depicting Vmem values (n = 7; means ± SEM) that correspond to 5% activation of persistent Na+ conductance. Data are plotted as a function of [Ca2+]ex: −80.2 ± 0.4 mV (0.1 mm; orange); −73.4 ± 2.4 mV (0.5 mm; green); −70.3 ± 2.2 mV (1 mm; black); −65.2 ± 2.2 mV (5 mm; gray). Asterisks (*) denote statistical significance, p < 0.05 (one-way ANOVA with Tukey's HSD post hoc test). Bi, Bii, Original whole-cell current-clamp recordings from two representative iAMCs at different [Ca2+]ex levels (preincubation, ≥3 min). Bii, Note that oscillations recover after restoration of physiological [Ca2+]ex. Biii, Bar diagram illustrating IBIs as a function of [Ca2+]ex. Data (n = 5; means ± SEM) are normalized to control conditions (1 mm Ca2+): 0.44 ± 0.15 (0.1 mm; orange); 0.88 ± 0.14 (0.5 mm; green). Asterisk (*) denotes statistical significance, p < 0.05 (one-way ANOVA with Tukey's HSD post hoc test).

A shift in voltage dependence caused by changes in extracellular Ca2+ concentration ([Ca2+]ex) is another hallmark of INaP (Su et al., 2001). Thus, we tested whether activation of TTX-sensitive steady-state currents was altered by changes in [Ca2+]ex (Fig. 6Aiv,Av). Likely as a result of surface charge screening (Hille, 2001), increasing [Ca2+]ex indeed shifted INaP activation to more depolarized values. When converted to conductances and fit by a Boltzmann function (Fig. 6Aiv), INaP activation (calculated as 5% conductance) varied as a function of [Ca2+]ex, ranging from −65.2 ± 2.2 mV (5 mm Ca2+) to −80.2 ± 0.4 mV (0.1 mm Ca2+; Fig. 6Av). Accordingly, if INaP is a major determinant of iAMC autorhythmicity, changes in [Ca2+]ex should alter oscillation patterns. Indeed, IBIs progressively decreased with lower [Ca2+]ex levels (Fig. 6Bi–Biii). By contrast, elevated concentrations (5 mm Ca2+) resulted in cessation of oscillations and stable Vmem values (Fig. 6Bi), an effect that was reversible upon restoration of physiological conditions (1 mm Ca2+; Fig. 6Bii). Together, we conclude that TTX-sensitive INaP is the major excitatory element that drives Vd–Vu transitions in iAMCs.

R-type Ca2+ channels determine iAMC oscillation patterns

During the up state, distinct ionic conductances translate the subthreshold excitatory drive into regular trains of action potentials (Colwell, 2011). Thus, we next investigated the mechanistic basis of burst firing during the depolarizing envelope of iAMCs. While the predominant role of INaT in AOB mitral cell spike generation is well established (Castro et al., 2007; Smith and Araneda, 2010; Shpak et al., 2012), voltage-gated Ca2+ (CaV) channels additionally shape various neuronal discharge parameters (Bean, 2007). Using a pharmacological approach, we first characterized the repertoire of high voltage-activated Ca2+ currents expressed in iAMCs. L-type, P-type/Q-type, and N-type Ca2+ currents were isolated using selective dihydropyridine or peptide toxin CaV-channel antagonists (McCleskey et al., 1987; Bossert and Vater, 1989; Adams et al., 1993). Notably, iAMC autorhythmicity persisted essentially unaltered in the presence of either CaV-channel blocker (data not shown). By contrast, addition of SNX-482, a selective R-type Ca2+ channel inhibitor (Bourinet et al., 2001), dramatically altered iAMC oscillations (Fig. 7A,B). While within-burst firing rates remained unchanged, inhibition of R-type Ca2+ currents (IR) significantly decreased IBIs and prolonged burst duration (Fig. 7Aii). These results indicate that, while iAMCs express at least one member of each high voltage-activated Ca2+ channel subfamily, only R-type/CaV2.3 channels play a significant role in orchestrating autonomous iAMC oscillations.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Oscillatory discharge involves R-type CaV channels. A, Inhibition of R-type/CaV2.3 channels alters iAMC oscillations. Ai, Representative whole-cell current-clamp recording of oscillatory iAMC discharge before (left) and during (right) SNX-482 treatment (100 nm; preincubation, ≥3 min). Aii, Bar charts comparing iAMC discharge parameters (n = 9) under fast synaptic isolation (iso; black) versus additional inhibition of R-type CaV channels (SNX-482; violet). IBIs, burst durations, and firing rates are plotted as means ± SEM. Block of R-type Ca2+ currents shortened IBIs (20.5 ± 3.7 vs 11.9 ± 2.0 s) and increased burst duration (4.2 ± 0.7 vs 6.4 ± 1.1 s). By contrast, within-burst firing frequencies (4.5 ± 0.9 vs 4.5 ± 0.9 Hz) were essentially unaffected. Asterisks (*) denote statistical significance, p < 0.01 (Wilcoxon signed-rank test). Bi, Bii, Both representative traces from voltage step recordings (Bi; whole-cell; S4/S7) and the corresponding I–V curve (Bii) reveal a voltage-dependent slowly inactivating inward current sensitive to SNX-482 (digital off-line subtraction). Steady-state current densities (means ± SEM; n = 4) are plotted as a function of membrane depolarization. Step protocol as indicated. Maximum currents (−19.3 ± 4.3 pA/pF) are detected at −30 mV. Ci, Spike waveforms recorded under synaptic isolation (iso; left; n = 131) and in additional presence of SNX-482 (right; n = 493). Recording duration, 180 s. Average action potentials are shown in red; gray shadows indicate SD. Inset, Schematic spike illustrating action potential analysis parameters: rise time, peak Vmem, FDHM. Cii, Bar charts comparing characteristic spike parameters in absence (iso; black) and presence (SNX; violet) of SNX-482. Data are means ± SEM (n = 9). While rise time is unchanged (1.38 ± 0.1 vs 1.42 ± 0.1 ms), SNX-482 treatment results in reduced peak Vmem values (32.3 ± 0.6 vs 28.8 ± 1.4 mV) and prolonged repolarization (FDHM: 2.1 ± 0.2 vs 2.4 ± 0.1 ms; integral: 212 ± 16 vs 241 ± 14 mV*ms). Asterisks (*) denote statistical significance, p < 0.01 (Wilcoxon signed-rank test).

In iAMCs, Vu is below the activation threshold of SNX-482-sensitive IR (>−50 mV; Fig. 7B). Therefore, action potentials must provide the depolarization required for IR activation during the up state. Accordingly, we next examined whether action potential amplitude or shape changed upon IR inhibition (Fig. 7C). Rise time analysis revealed that SNX-482-sensitive currents contribute little to the action potential rising phase (Fig. 7Cii). By contrast, SNX-482 significantly reduced the amplitude and prolonged the FDHM of average iAMC action potentials (Fig. 7Cii). These data indicate that CaV2.3 channels activate near the action potential peak and that IR is largest during the falling phase when the driving force on Ca2+ increases.

Big conductance K+ channel-dependent negative feedback regulates burst duration

A prolonged action potential FDHM as a result of blocking a Ca2+ inward current (Fig. 7Cii) is somewhat counterintuitive. Apparently, the net effect of blocking Ca2+ entry is to inhibit a net outward current. Ca2+-activated K+ (KCa) channels are likely mediators of such effects (Bean, 2007). Therefore, we next asked whether Ca2+-activated K+ conductances are involved in iAMC oscillations.

Underlying the afterhyperpolarization that follows single spikes or bursts, for example in the hippocampus (Alger and Nicoll, 1980), small conductance KCa (SK) channels are major determinants of firing rate. Thus, we first investigated whether the selective SK channel antagonist apamin (Blatz and Magleby, 1986) affects rhythmic iAMC discharge. Surprisingly, SK channel block did not affect iAMC oscillations (Fig. 8A), consistent with the lack of a substantial apamin-sensitive KCa current in voltage-clamp experiments (Fig. 8B).

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Functional interaction of CaV and BK channels determines burst properties. Ai, Representative current-clamp recording of oscillatory discharge before (top) and during (bottom) apamin treatment (300 nm; preincubation, ≥3 min). Aii, Bar charts comparing iAMC discharge parameters (n = 10) under control conditions (black) versus inhibition of SK channels (apamin; brown). IBIs, burst durations, within-burst firing rates, and Vmem state are plotted as means ± SEM. Apamin treatment shows essentially no effect on either parameter (p > 0.05; paired sample t test): IBIs, 10.8 ± 0.9 versus 9.9 ± 1.1 s; burst durations, 7.0 ± 1.3 versus 6.6 ± 1.4 s; within-burst firing frequencies, 4.6 ± 1.3 versus 4.4 ± 1.1 Hz; Vd, −81.1 ± 1.1 versus −80.0 ± 0.9 mV; Vu, −68.7 ± 0.6 versus −67.8 ± 0.9 mV. B, Representative recording of Ca2+-dependent whole-cell currents from AOB mitral cells in response to a “pseudo” tail current pulse protocol (inset) designed to isolate SK channel-mediated currents: “brief” step depolarization (+10 mV, 100 ms) triggers Ca2+ influx via CaV channels; recordings are performed upon repolarization to −50 mV (Hallworth et al., 2003; Stocker, 2004). Using digital subtraction (brown trace), essentially no apamin-sensitive current is isolated (n = 6). Ci, Representative Vmem recording illustrating oscillatory activity before and after, as well as irregular discharge during wash-in of relatively low TEA concentrations (1 mm; horizontal bar indicates drug incubation). Cii, Ciii, Expanded view of periodic bursting before (Cii) and irregular activity during (Ciii) TEA treatment, respectively. Rhythmicity (as well as the lack thereof) is also evident from the corresponding Vmem all-points histograms. Red arrowheads indicate Vmem states (Vd; Vu) under synaptic isolation. Di, Using a prerecorded prototypical burst as the voltage-clamp command waveform (template; top), individual current components are isolated owing to their specific pharmacological profile. Currents sensitive to TEA (1 mm; middle) and Cd2+ (200 μm; bottom) are obtained by successive digital subtraction. Red arrowheads indicate the first and a randomly chosen intermediate (29th) action potential (see below). Dii, Expanded view and overlay of TEA-sensitive and Cd2+-sensitive currents evoked by the first (light colors) and 29th (dark colors) action potential, respectively. Note that, while the initial action potential shows the largest amplitude in the prerecorded burst template (Di), both peak current and charge transfer (trace integrals) evoked by the first spike are substantially smaller compared with later events. When TEA-sensitive (Diii) and Cd2+-sensitive (Div) charge transfer is normalized to each cell's maximum and plotted as a function of spike count (event number), the charge transfer function increased monotonically. Black dots represent results from individual measurements; red dots indicate means [n = 7 (Diii) and 6 (Div)]. Curves were fit by monoexponential functions (red lines). Dv, Scatter plot depicting the relationship between TEA-sensitive and Cd2+-sensitive absolute charge transfer during a burst. Linear regression indicates strong correlation of both parameters (Pearson's correlation coefficient, r = 0.977; n = 6). Ei, Representative Vmem recording of oscillatory discharge before (top) and during (bottom) paxilline treatment (5 μm; preincubation, ≥3 min). Eii, Bar charts comparing iAMC discharge parameters (n = 8) under conditions of synaptic isolation (iso; black) versus additional inhibition of BK channels (paxilline; green). IBIs, burst durations, within-burst firing rates, and Vmem states are plotted as means ± SEM. Paxilline treatment shortened IBIs (13.9 ± 2.0 vs 8.7 ± 0.8 s) and increased burst durations (5.5 ± 1.1 vs 10.6 ± 1.9 s). Firing frequencies (3.4 ± 0.6 vs 3.9 ± 0.5 Hz) as well as Vd and Vu (−75.6 ± 0.8 vs −75.1 ± 1.4 mV; −65.2 ± 0.9 vs −62.0 ± 1.1 mV) were essentially unaffected. Asterisks (*) denote statistical significance, p < 0.05 (paired sample t test).

Other candidates for mediating a Ca2+-activated K+ efflux during discharge are big conductance KCa (BK) channels. BK channels, which are cooperatively activated by depolarization and increased cytoplasmic Ca2+ (Fakler and Adelman, 2008), are inhibited by both TEA and the selective organic blocker paxilline (Brenner et al., 2005). In iAMCs, relatively low TEA concentrations (1 mm) reversibly switched oscillatory discharge to irregular tonic firing (Fig. 8C), suggesting that coupling of Ca2+ entry to activation and “build up” of IBK might cause burst termination and, consequently, the following transition to Vd (Crunelli et al., 2012). To quantify the relative contributions of IBK and ICaV to pacemaking, we then used an iAMC's own representative burst firing pattern as a voltage command and recorded both TEA-sensitive and Cd2+-sensitive currents (Fig. 8Di). When single spike-dependent currents were plotted as a function of event number (Fig. 8Dii–Div), both putative IBK and ICaV monotonically increased and gradually reached saturation after ∼15–20 events. Thus, while the recorded template showed spike-amplitude adaptation, contributions of BK and CaV channels were largest in later stages of up-state bursting. Moreover, TEA-sensitive and Cd2+-sensitive charge transfer was strongly correlated (Fig. 8Dv). Together, these data suggest that, during prolonged firing, a gradual increase in cytoplasmic Ca2+ is coupled to progressive activation of BK channels, which ultimately causes burst termination. If so, selective inhibition of IBK should increase up-state duration and reduce IBIs. Low micromolar concentrations of paxilline indeed caused corresponding changes in iAMC rhythmogenesis (Fig. 8E), supporting the notion that BK channels play a major role in setting iAMC oscillation frequency.

A “triangle” of interdependent distinct conductances is sufficient to drive autorhythmicity

The above results indicate that mitral cell autorhythmicity emerges from coordinated reciprocal interaction of INaP, IR, and IBK. While the molecular correlate of INaP is debated and subthreshold persistent currents might even originate from the same channels that carry suprathreshold transient currents (Carter et al., 2012), the CaV2.3 channel and the BK channel α-1 subunit Sloα1 have been shown to mediate IR and IBK, respectively. Indeed, antibodies raised against CaV2.3 and Sloα1 label a large population of cells in the AOB mitral cell layer (Fig. 9A,B). Notably, the pharmacological profile of IBK in AOB mitral cells—sensitivity to TEA and paxilline, but resistance to iberiotoxin and charybdotoxin (data not shown)—suggests coexpression of an auxiliary β4 subunit (Sloβ4) that confers type-II BK channel properties (Brenner et al., 2005), such as relatively slow gating kinetics. Immunolabeling against Sloβ4 confirmed its specific expression in AOB mitral cells (Fig. 9C). Moreover, both Sloα1 and Sloβ4 were expressed in neurons identified as iAMCs in prior electrophysiological recordings (Fig. 9D,E).

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Expression of BK and CaV2.3 channels in AOB mitral cells. A–C, Confocal laser-scanning microscopy images of the mouse AOB (parasagittal cryosections). Immunofluorescence labeling (green) of CaV2.3 (A), the BK channel subunit α-1 (Sloα1; B), and the BK channel subunit β-4 (Sloβ4; C). Nuclei are counterstained with DRAQ5 (red). Dashed lines indicate both the lateral olfactory tract (LOT) and the glomerular layer (GL) in Aiii, Biii, and Ciii. Representative single optical section images show an overview of the AOB (Ai, Bi, Ci) and a high-magnification view (Aii, Bii, Cii) of the mitral cell layer (MCL). Note the robust labeling of cells in the MCL as well as the absence of immunofluorescence following peptide preadsorption (Aiii, Biii, Ciii; no DRAQ5 counterstaining). D, E, Post hoc immunolabeling of previously recorded intrinsically rhythmic neurons. Confocal images of AOB cryosections containing biocytin-filled iAMC somata and proximal dendrites (visualized by Alexa Fluor 488/633 streptavidin conjugate). Sections were counterstained against the BK channel α-1 (D) and β-4 (E) subunits, respectively. Merged images (Diii, Eiii) reveal BK expression in these representative iAMCs. Div, Eiv, Western blot analysis of total olfactory bulb (right lanes) and posterior brain (left lanes) protein extracts. Immunoblots were probed with anti-Sloα1 (D) and anti-Sloβ4 (E) antibodies, respectively. Expected band size is indicated above both blots. d, Dorsal; FC, frontal cortex; v, ventral.

Are INaP, IR, and IBK unique properties of rhythmogenic AOB mitral cells or are these conductances broadly expressed in both iAMCs and nonrhythmogenic neurons? To address this question, we compared electrophysiological profiles of INaP, IR, and IBK in both mitral cell populations (Fig. 10A–G). All three currents are found in both iAMCs and irregularly firing mitral cells. However, with the exception of SNX-482-sensitive IR (Fig. 10D,E), properties of rhythmogenic currents markedly differed between both neuronal populations. The INaP activation threshold was significantly lower in iAMCs than in nonrhythmogenic neurons (Fig. 10A,Bi). In addition to the resulting left shift in the I–V relationship (Fig. 10C), maximum INaP amplitudes—calculated from sigmoidal fits to the individual I–V curves—were significantly increased in iAMCs (Fig. 10Bii). TEA-sensitive K+ currents also differed between oscillating and irregularly firing mitral cells (Fig. 10F,G). Larger IBK amplitudes are recorded from iAMCs and steady-state activation curves in oscillating neurons are considerably shifted toward lower Vmem values (Fig. 10G). These results show that, while not exclusively expressed in the intrinsically oscillating mitral cell population, INaP, IR, and IBK confer a distinct electrophysiological phenotype on iAMCs.

Figure 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 10.

A distinct combination of currents is sufficient to generate iAMC autorhythmicity. A, I–V relationship of TTX-sensitive whole-cell currents (300 nm; mean of 3 consecutive recordings; INaP isolated by digital subtraction) evoked by slow (10 mV/300 ms) ascending voltage ramps from −90 to −50 mV (inset) in a representative iAMC under synaptic isolation [gabazine (10 μm) + AP5 (100 μm) + NBQX (10 μm)]. Bi, Bii, Bar diagrams quantifying either (Bi) the average activation threshold of INaP in iAMCs (oscill.; red; Vthreshold = −80.1 ± 1.2 mV) versus nonrhythmogenic neurons (irreg.; black; Vthreshold = −74.6 ± 0.7 mV) or (Bii) the maximum INaP amplitude in oscillating (red; −164.0 ± 36.8 pA) versus irregular neurons (black; −87.9 ± 9.7 pA). Data are means ± SEM. Asterisks (*) denote statistical significance, p < 0.0001 (Bi) or p < 0.05 (Bii; unpaired sample t test). Numbers of experiments are denoted above bars. C, Average I–V relationships of TTX-sensitive INaP recorded from oscillating (red; n = 6–17) versus irregular mitral cells (black; n = 16–23). Data points depict means ± SEM. Curves are sigmoidal Boltzmann-type equation fits assuming threshold Vmem values shown in Bi. Vmem values corresponding to half-maximal activation are −60.4 ± 0.8 mV versus −60.6 ± 14.5 mV in oscillating and irregular neurons, respectively. D, SNX-482 (100 nm) sensitive IR recorded from an irregularly firing mitral cell under synaptic isolation in response to a fast ascending voltage ramp (−100 to + 80 mV; 100 ms duration) and plotted as a function of command voltage. E, Bar charts comparing IR characteristics in oscillating (violet) versus irregular (black) mitral cells. Ei–Eiv, SNX-sensitive peak current (Ei; 10.3 ± 3.1 vs 9.9 ± 1.6 pA/pF), activation threshold (Eii; −38.8 ± 3.1 vs −35.1 ± 4.2 mV), Vmem inducing maximum currents (Eiii; −20.6 ± 6.9 vs −17.1 ± 4.3 mV), and R-type channel conductance (Eiv; 0.73 ± 0.2 vs 0.64 ± 0.5 nS) are not significantly different between both mitral cell populations (unpaired sample t test). Data are means ± SEM. Numbers of experiments are denoted above bars. F, K+ current sensitive to TEA (1 mm) recorded from a representative iAMC (whole-cell voltage clamp). Inset depicts pulse protocol. G, Normalized IBK activation curve recorded from oscillating (green, n = 3) versus irregular (black, n = 3) mitral cells. Steady-state current amplitudes are measured upon stepping back to 0 mV (+ in F) and plotted as a function of prepulse Vmem. Solid lines represent sigmoidal fits to data points (mean ± SD; V1/2 = −2.0 ± 0.8 vs 19.0 ± 0.7 mV). V1/2 shift as indicated. Inset depicts average steady-state current densities of TEA-sensitive IBK (* in F) as a function of command voltage (n = 3; means ± SD). Hi, Volume-rendered 3D reconstruction of a representative iAMC used for modeling spontaneous discharge. Hii, Hiii, Model-based Vmem simulations (top) qualitatively reproduce experimental recordings (bottom) with respect to both IBIs and burst durations. Note that the different Vmem simulation outcomes result mainly from model variation in Kdr channel density and deactivation time constant.

We next asked whether the interplay of INaP, IR, and IBK is sufficient to drive autorhythmicity in a model AOB mitral cell. Based on the volume-rendered 3D reconstruction of a representative iAMC (Fig. 10Hi), Vmem simulations were run using the Neuron simulation environment (Hines and Carnevale, 1997). In addition to INaP, IR, and IBK, the model neuron also expressed IKdr and IKA (based on our own voltage-clamp experiments; see Materials and Methods), as well as INaT (Migliore et al., 2005). Strikingly, model-based simulations qualitatively reproduced the experimental data, with relatively small changes in Kdr channel density accounting for much of the experimentally observed heterogeneity in both IBIs and burst duration (Fig. 10Hii–Hiii). Together these data demonstrate that iAMCs express a “triangle” of conductances with distinct biophysical properties—INaP, IR, and IBK —that are sufficient to generate autorhythmicity in mouse AOB mitral cells.

Discussion

The accessory olfactory system is central to social information processing in various mammalian species. However, many basic physiological principles underlying sensory processing in the AOB remain poorly understood (Dulac and Wagner, 2006). Here, we investigate spontaneous activity in AOB mitral cells, the sole output/projection neurons of the AOB. We report that a group of AOB mitral cells displays slow stereotypical rhythmic discharge both in vitro and in vivo. Among these neurons, a subpopulation is intrinsically rhythmogenic. In these pacemaker-like cells, reciprocal interactions between depolarizing and hyperpolarizing conductances generate stable patterns of Vmem oscillations within a broad spectral window.

The prevalence of slow oscillatory discharge was larger in AOB slices than observed in vivo. While we recorded rhythmic bursting in 16.7% of all units in vivo, 69.7% of AOB mitral cells displayed periodic discharge in vitro. Several factors could account for this apparent discrepancy. First, AOB slices are isolated from both peripheral sensory input and top-down modulation. When connectivity is intact, either factor could add substantial “noise” to a given mitral cell's output. Second, experimental in vitro conditions might favor oscillatory discharge. While we performed loose-seal extracellular recordings in parallel to whole-cell patch-clamp measurements to prevent dialysis of intracellular components, extracellular in vitro conditions might not exactly mirror endogenous ionic concentrations. For example, our data indicate that small differences in [Ca2+]ex could exert profound effects on rhythmogenesis. Regardless of the exact prevalence, slow/infraslow mitral cell rhythmicity will likely have considerable physiological impact on sensory processing along the accessory olfactory pathway.

At the lower end of the “bandwidth” scale (Buzsáki, 2006), infraslow oscillations with a periodicity of tens of seconds have been described in a variety of neurons and circuits (Hughes et al., 2002; Blethyn et al., 2006; Crunelli and Hughes, 2010; Buzsáki et al., 2013). In iAMCs, rhythmogenesis appears to be driven by cyclical interplay of ≥3 voltage-dependent and/or Ca2+-dependent conductances. Low threshold INaP functions as the major excitatory element that drives iAMC transition from Vd to Vu. R-type CaV channels play a significant role in oscillation maintenance and shape, while the resulting increase in cytoplasmic Ca2+ is coupled to progressive activation of BK channels which, in concert with slow voltage-dependent INaP inactivation (Jasinski et al., 2013), ultimately causes burst termination. Here, we observed that INaP becomes evident in iAMCs at voltages ≥−80 mV. This finding supports recent reports of INaP activation at more negative voltages than previously appreciated (Huang and Trussell, 2008; Carter et al., 2012). The INaP pacemaking “engine” can therefore be engaged or disengaged by small shifts in net current at the down-state Vmem range (Yamada-Hanff and Bean, 2013). In nonrhythmogenic neurons, the INaP activation threshold is shifted to ∼−75 mV, a Vmem value that is more positive than the down state in most iAMCs (median Vd = −76.4 mV). An interplay of R-type CaV and BK channels also supports relatively prolonged bursting. IR, which contributes to bursting in hippocampal (Metz et al., 2005) and thalamic reticular neurons (Zaman et al., 2011), is similar to IT but activates above threshold (Foehring et al., 2000). Importantly, among high-voltage-activated Ca2+ channels, only R-type channels do not associate with BK channels in Ca2+ nanodomains (Fakler and Adelman, 2008), indicating the requirement of a slower, more global increase in cytoplasmic Ca2+ to activate substantial IBK. Together with the pronounced left shift in the Vmem dependence of IBK activation that we observed in iAMCs, the specific IR–IBK combination expressed in oscillating neurons appears ideally suited to shape autorhythmicity.

Our results do not exclude a functional role of other conductances or electrogenic pumps (Krey et al., 2010; Zylbertal et al., 2015). While we could not demonstrate a significant contribution of Ih, IT, or L-type, P-type/Q-type, and N-type Ca2+ currents or SK channel-mediated K+ currents, the activity-dependent cytoplasmic build-up of Na+ and/or Ca2+ during the burst could additionally “recruit” other currents. Potential candidates include Na+-activated K+ currents (Hage and Salkoff, 2012) and Ca2+-activated nonselective cation currents (Peña et al., 2004; Blethyn et al., 2006; Shpak et al., 2012). The Na+-activated K+ current, in particular, could prove important by mitigating INaP-dependent depolarization, effectively decelerating the Vd-to-Vu transition. Indeed, the biophysical foundation of pacemaker-like activity (e.g., in central pattern generators) is often redundant, such that identical firing patterns can be achieved by multiple current combinations (Marder and Goaillard, 2006). Notwithstanding, when added to a model mitral cell that otherwise expresses “standard” Hodgkin–Huxley-type voltage-gated Na+ and K+ channels, the combination of INaP, IR, and IBK is sufficient to drive autorhythmicity.

To what extent intrinsically rhythmogenic mitral cells shape information processing and sensory coding in the AOB remains so far unexplored. Future experiments will thus have to address the impact of iAMC oscillatory discharge on the intact AOB network. The wide spectrum of discharge patterns, as well as the absence of a dominant oscillation frequency among iAMCs, argues against a distinct prevalent AOB rhythm. Instead, parallel pacemaker-like activity of phenotypically different iAMCs might generate synchronous discharge in several AOB microcircuits. By temporally linking AOB neurons into functional assemblies, synchronous firing could facilitate synaptic plasticity and input selection (Buzsáki and Draguhn, 2004). Rhythmic cycles between high and low postsynaptic excitability states would add a temporal dimension to a given circuit's sensory coding space (Schroeder and Lakatos, 2009). Similar to intrinsic theta frequency oscillations that entrain to the sniffing cycle in ET cells of the MOB (Hayar et al., 2004, 2005), operation of the vomeronasal pump (Meredith and O'Connell, 1979) could entrain iAMC oscillations. In the hamster vomeronasal organ, peristaltic vasoconstriction cycles of 0.2–0.5 Hz were recorded in vivo (Meredith, 1994). In addition to input entrainment, orchestration of AOB neurons into functional ensembles could ensure communication reliability and selectivity (Izhikevich et al., 2003) by controlling signal flow among anatomically connected networks. Notably, downstream processing modules for vomeronasal stimuli include several nuclei that mediate pulsatile neuroendocrine release by synchronized slow rhythmic bursting of, for example, gonadotropin-releasing hormone neurons (Chu et al., 2012) or vasopressin magnocellular neurosecretory cells (Brown, 2004).

In summary, we identify a subpopulation of AOB mitral cells that exhibit slow/infraslow spontaneous rhythmicity both in vivo and in vitro. Some of these projection neurons are intrinsically rhythmogenic. The mechanistic basis of mitral cell autorhythmicity is the coordinated interplay of a regenerative subthreshold Na+ current, a high voltage-activated R-type Ca2+ current, and a burst-terminating Ca2+-activated K+ current that is cooperatively activated by depolarization and cytoplasmic Ca2+ accumulation.

Footnotes

  • This work was supported by grants from the Volkswagen Foundation (I/83533; to M.S.), the Deutsche Forschungsgemeinschaft (SP724/8-1, to M.S.; MA1184/20-2, to R.M.), DFG priority program SPP 1392: Integrative Analysis of Olfaction, Marie Skłodowska-Curie Actions (PCIG10-GA-2011-303785, to Y.B.-S.), the Lady Davis Fellowship Trust (to A.K.), and the German-Israeli Foundation for Scientific Research and Development (1-1193-153.13/2012, to M.S. and Y.B.-S.). M.S. is a Lichtenberg Professor of the Volkswagen Foundation. We thank both Corinna H. Engelhardt and Susanne Lipartowski (RWTH Aachen University) for excellent technical assistance.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Marc Spehr, RWTH Aachen University, Institute for Biology II/Dept. Chemosensation, Worringerweg 3, D-52074 Aachen, Germany. m.spehr{at}sensorik.rwth-aachen.de

References

  1. ↵
    1. Adams ME,
    2. Myers RA,
    3. Imperial JS,
    4. Olivera BM
    (1993) Toxityping rat brain calcium channels with omega-toxins from spider and cone snail venoms. Biochemistry 32:12566–12570, doi:10.1021/bi00210a003, pmid:8251474.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Alger BE,
    2. Nicoll RA
    (1980) Epileptiform burst afterhyperolarization: calcium-dependent potassium potential in hippocampal CA1 pyramidal cells. Science 210:1122–1124, doi:10.1126/science.7444438, pmid:7444438.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Barry PH
    (1994) JPCalc, a software package for calculating liquid junction potential corrections in patch-clamp, intracellular, epithelial and bilayer measurements and for correcting junction potential measurements. J Neurosci Methods 51:107–116, doi:10.1016/0165-0270(94)90031-0, pmid:8189746.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bean BP
    (2007) The action potential in mammalian central neurons. Nat Rev Neurosci 8:451–465, pmid:17514198.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Ben-Shaul Y,
    2. Katz LC,
    3. Mooney R,
    4. Dulac C
    (2010) In vivo vomeronasal stimulation reveals sensory encoding of conspecific and allospecific cues by the mouse accessory olfactory bulb. Proc Natl Acad Sci U S A 107:5172–5177, doi:10.1073/pnas.0915147107, pmid:20194746.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Blankenship AG,
    2. Feller MB
    (2010) Mechanisms underlying spontaneous patterned activity in developing neural circuits. Nat Rev Neurosci 11:18–29, pmid:19953103.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Blatz AL,
    2. Magleby KL
    (1986) Single apamin-blocked Ca-activated K+ channels of small conductance in cultured rat skeletal muscle. Nature 323:718–720, doi:10.1038/323718a0, pmid:2430185.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Blethyn KL,
    2. Hughes SW,
    3. Tóth TI,
    4. Cope DW,
    5. Crunelli V
    (2006) Neuronal basis of the slow (<1 Hz) oscillation in neurons of the nucleus reticularis thalami in vitro. J Neurosci 26:2474–2486, doi:10.1523/JNEUROSCI.3607-05.2006, pmid:16510726.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. BoSmith RE,
    2. Briggs I,
    3. Sturgess NC
    (1993) Inhibitory actions of ZENECA ZD7288 on whole-cell hyperpolarization activated inward current (If) in guinea-pig dissociated sinoatrial node cells. Br J Pharmacol 110:343–349, doi:10.1111/j.1476-5381.1993.tb13815.x, pmid:7693281.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Bossert F,
    2. Vater W
    (1989) 1,4-Dihydropyridines-A basis for developing new drugs. Med Res Rev 9:291–324, doi:10.1002/med.2610090304, pmid:2666803.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Bourinet E,
    2. Stotz SC,
    3. Spaetgens RL,
    4. Dayanithi G,
    5. Lemos J,
    6. Nargeot J,
    7. Zamponi GW
    (2001) Interaction of SNX482 with domains III and IV inhibits activation gating of alpha(1E) (Ca(V)2.3) calcium channels. Biophys J 81:79–88, doi:10.1016/S0006-3495(01)75681-0, pmid:11423396.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Brenner R,
    2. Chen QH,
    3. Vilaythong A,
    4. Toney GM,
    5. Noebels JL,
    6. Aldrich RW
    (2005) BK channel b4 subunit reduces dentate gyrus excitability and protects against temporal lobe seizures. Nat Neurosci 8:1752–1759, doi:10.1038/nn1573, pmid:16261134.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Brown CH
    (2004) Rhythmogenesis in vasopressin cells. J Neuroendocrinol 16:727–739, doi:10.1111/j.1365-2826.2004.01227.x, pmid:15344911.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Bucher D,
    2. Taylor AL,
    3. Marder E
    (2006) Central pattern generating neurons simultaneously express fast and slow rhythmic activities in the stomatogastric ganglion. J Neurophysiol 95:3617–3632, doi:10.1152/jn.00004.2006, pmid:16495367.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Buzsáki G
    (2006) Rhythms of the brain (Oxford UP, New York).
  16. ↵
    1. Buzsáki G,
    2. Draguhn A
    (2004) Neuronal oscillations in cortical networks. Science 304:1926–1929, doi:10.1126/science.1099745, pmid:15218136.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Buzsáki G,
    2. Logothetis N,
    3. Singer W
    (2013) Scaling brain size, keeping timing: evolutionary preservation of brain rhythms. Neuron 80:751–764, doi:10.1016/j.neuron.2013.10.002, pmid:24183025.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Carter BC,
    2. Giessel AJ,
    3. Sabatini BL,
    4. Bean BP
    (2012) Transient sodium current at subthreshold voltages: activation by EPSP waveforms. Neuron 75:1081–1093, doi:10.1016/j.neuron.2012.08.033, pmid:22998875.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Castro JB,
    2. Hovis KR,
    3. Urban NN
    (2007) Recurrent dendrodendritic inhibition of accessory olfactory bulb mitral cells requires activation of group I metabotropic glutamate receptors. J Neurosci 27:5664–5671, doi:10.1523/JNEUROSCI.0613-07.2007, pmid:17522311.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Chan CS,
    2. Shigemoto R,
    3. Mercer JN,
    4. Surmeier DJ
    (2004) HCN2 and HCN1 channels govern the regularity of autonomous pacemaking and synaptic resetting in globus pallidus neurons. J Neurosci 24:9921–9932, doi:10.1523/JNEUROSCI.2162-04.2004, pmid:15525777.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Chu Z,
    2. Tomaiuolo M,
    3. Bertram R,
    4. Moenter SM
    (2012) Two types of burst firing in gonadotrophin-releasing hormone neurones. J Neuroendocrinol 24:1065–1077, doi:10.1111/j.1365-2826.2012.02313.x, pmid:22435872.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Cichy A,
    2. Ackels T,
    3. Tsitoura C,
    4. Kahan A,
    5. Gronloh N,
    6. Söchtig M,
    7. Engelhardt CH,
    8. Ben-Shaul Y,
    9. Müller F,
    10. Spehr J,
    11. Spehr M
    (2015) Extracellular pH regulates excitability of vomeronasal sensory neurons. J Neurosci 35:4025–4039, doi:10.1523/JNEUROSCI.2593-14.2015, pmid:25740530.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Colwell CS
    (2011) Linking neural activity and molecular oscillations in the SCN. Nat Rev Neurosci 12:553–569, doi:10.1038/nrn3086, pmid:21886186.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Crill WE
    (1996) Persistent sodium current in mammalian central neurons. Annu Rev Physiol 58:349–362, doi:10.1146/annurev.ph.58.030196.002025, pmid:8815799.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Crunelli V,
    2. Hughes SW
    (2010) The slow (<1 Hz) rhythm of non-REM sleep: a dialogue between three cardinal oscillators. Nat Neurosci 13:9–17, doi:10.1038/nn.2445, pmid:19966841.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Crunelli V,
    2. Lörincz ML,
    3. Errington AC,
    4. Hughes SW
    (2012) Activity of cortical and thalamic neurons during the slow (<1 Hz) rhythm in the mouse in vivo. Pflugers Arch 463:73–88, doi:10.1007/s00424-011-1011-9, pmid:21892727.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Del Punta K,
    2. Puche A,
    3. Adams NC,
    4. Rodriguez I,
    5. Mombaerts P
    (2002) A divergent pattern of sensory axonal projections is rendered convergent by second-order neurons in the accessory olfactory bulb. Neuron 35:1057–1066, doi:10.1016/S0896-6273(02)00904-2, pmid:12354396.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Dulac C,
    2. Torello AT
    (2003) Molecular detection of pheromone signals in mammals: from genes to behaviour. Nat Rev Neurosci 4:551–562, doi:10.1038/nrn1140, pmid:12838330.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Dulac C,
    2. Wagner S
    (2006) Genetic analysis of brain circuits underlying pheromone signaling. Annu Rev Genet 40:449–467, doi:10.1146/annurev.genet.39.073003.093937, pmid:16953793.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Fakler B,
    2. Adelman JP
    (2008) Control of K(Ca) channels by calcium nano/microdomains. Neuron 59:873–881, doi:10.1016/j.neuron.2008.09.001, pmid:18817728.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Foehring RC,
    2. Mermelstein PG,
    3. Song WJ,
    4. Ulrich S,
    5. Surmeier DJ
    (2000) Unique properties of R-type calcium currents in neocortical and neostriatal neurons. J Neurophysiol 84:2225–2236, pmid:11067968.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Fukunaga I,
    2. Berning M,
    3. Kollo M,
    4. Schmaltz A,
    5. Schaefer AT
    (2012) Two distinct channels of olfactory bulb output. Neuron 75:320–329, doi:10.1016/j.neuron.2012.05.017, pmid:22841316.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Gentet LJ,
    2. Stuart GJ,
    3. Clements JD
    (2000) Direct measurement of specific membrane capacitance in neurons. Biophys J 79:314–320, doi:10.1016/S0006-3495(00)76293-X, pmid:10866957.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Gleeson P,
    2. Steuber V,
    3. Silver RA
    (2007) neuroConstruct: a tool for modeling networks of neurons in 3D space. Neuron 54:219–235, doi:10.1016/j.neuron.2007.03.025, pmid:17442244.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Grillner S
    (2006) Biological pattern generation: the cellular and computational logic of networks in motion. Neuron 52:751–766, doi:10.1016/j.neuron.2006.11.008, pmid:17145498.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Gutierrez GJ,
    2. O'Leary T,
    3. Marder E
    (2013) Multiple mechanisms switch an electrically coupled, synaptically inhibited neuron between competing rhythmic oscillators. Neuron 77:845–858, doi:10.1016/j.neuron.2013.01.016, pmid:23473315.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Hage TA,
    2. Salkoff L
    (2012) Sodium-activated potassium channels are functionally coupled to persistent sodium currents. J Neurosci 32:2714–2721, doi:10.1523/JNEUROSCI.5088-11.2012, pmid:22357855.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Hallworth NE,
    2. Wilson CJ,
    3. Bevan MD
    (2003) Apamin-sensitive small conductance calcium-activated potassium channels, through their selective coupling to voltage-gated calcium channels, are critical determinants of the precision, pace, and pattern of action potential generation in rat subthalamic nucleus neurons in vitro. J Neurosci 23:7525–7542, pmid:12930791.
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Hammen GF,
    2. Turaga D,
    3. Holy TE,
    4. Meeks JP
    (2014) Functional organization of glomerular maps in the mouse accessory olfactory bulb. Nat Neurosci 17:953–961, doi:10.1038/nn.3738, pmid:24880215.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Harris KD,
    2. Henze DA,
    3. Csicsvari J,
    4. Hirase H,
    5. Buzsáki G
    (2000) Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements. J Neurophysiol 84:401–414, pmid:10899214.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Hayar A,
    2. Karnup S,
    3. Shipley MT,
    4. Ennis M
    (2004) Olfactory bulb glomeruli: external tufted cells intrinsically burst at theta frequency and are entrained by patterned olfactory input. J Neurosci 24:1190–1199, doi:10.1523/JNEUROSCI.4714-03.2004, pmid:14762137.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Hayar A,
    2. Shipley MT,
    3. Ennis M
    (2005) Olfactory bulb external tufted cells are synchronized by multiple intraglomerular mechanisms. J Neurosci 25:8197–8208, doi:10.1523/JNEUROSCI.2374-05.2005, pmid:16148227.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Hazan L,
    2. Zugaro M,
    3. Buzsáki G
    (2006) Klusters, NeuroScope, NDManager: a free software suite for neurophysiological data processing and visualization. J Neurosci Methods 155:207–216, doi:10.1016/j.jneumeth.2006.01.017, pmid:16580733.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Hille B
    (2001) Ion channels of excitable membranes (Sinauer Associates, Sunderland, MA), 3rd ed.
  45. ↵
    1. Hines ML,
    2. Carnevale NT
    (1997) The NEURON simulation environment. Neural Comput 9:1179–1209, doi:10.1162/neco.1997.9.6.1179, pmid:9248061.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Hines ML,
    2. Carnevale NT
    (2001) NEURON: a tool for neuroscientists. Neurosci 7:123–135, pmid:11496923.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Hovis KR,
    2. Ramnath R,
    3. Dahlen JE,
    4. Romanova AL,
    5. LaRocca G,
    6. Bier ME,
    7. Urban NN
    (2012) Activity regulates functional connectivity from the vomeronasal organ to the accessory olfactory bulb. J Neurosci 32:7907–7916, doi:10.1523/JNEUROSCI.2399-11.2012, pmid:22674266.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Huang H,
    2. Trussell LO
    (2008) Control of presynaptic function by a persistent Na+ current. Neuron 60:975–979, doi:10.1016/j.neuron.2008.10.052, pmid:19109905.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Hughes SW,
    2. Cope DW,
    3. Blethyn KL,
    4. Crunelli V
    (2002) Cellular mechanisms of the slow (<1 Hz) oscillation in thalamocortical neurons in vitro. Neuron 33:947–958, doi:10.1016/S0896-6273(02)00623-2, pmid:11906700.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Izhikevich EM,
    2. Desai NS,
    3. Walcott EC,
    4. Hoppensteadt FC
    (2003) Bursts as a unit of neural information: selective communication via resonance. Trends Neurosci 26:161–167, doi:10.1016/S0166-2236(03)00034-1, pmid:12591219.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Jasinski PE,
    2. Molkov YI,
    3. Shevtsova NA,
    4. Smith JC,
    5. Rybak IA
    (2013) Sodium and calcium mechanisms of rhythmic bursting in excitatory neural networks of the pre-Bötzinger complex: a computational modelling study. Eur J Neurosci 37:212–230, doi:10.1111/ejn.12042, pmid:23121313.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Kay LM,
    2. Beshel J,
    3. Brea J,
    4. Martin C,
    5. Rojas-Líbano D,
    6. Kopell N
    (2009) Olfactory oscillations: the what, how and what for. Trends Neurosci 32:207–214, doi:10.1016/j.tins.2008.11.008, pmid:19243843.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Koizumi H,
    2. Smith JC
    (2008) Persistent Na+ and K+-dominated leak currents contribute to respiratory rhythm generation in the pre-Bötzinger complex in vitro. J Neurosci 28:1773–1785, doi:10.1523/JNEUROSCI.3916-07.2008, pmid:18272697.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Krey RA,
    2. Goodreau AM,
    3. Arnold TB,
    4. Del Negro CA
    (2010) Outward currents contributing to inspiratory burst termination in preBötzinger complex neurons of neonatal mice studied in vitro. Front Neural Circuits 4:124, doi:10.3389/fncir.2010.00124, pmid:21151816.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Larriva-Sahd J
    (2008) The accessory olfactory bulb in the adult rat: a cytological study of its cell types, neuropil, neuronal modules, and interactions with the main olfactory system. J Comp Neurol 510:309–350, doi:10.1002/cne.21790, pmid:18634021.
    OpenUrlCrossRefPubMed
  56. ↵
    1. LeDoux JE
    (2000) Emotion circuits in the brain. Annu Rev Neurosci 23:155–184, doi:10.1146/annurev.neuro.23.1.155, pmid:10845062.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Leszkowicz E,
    2. Khan S,
    3. Ng S,
    4. Ved N,
    5. Swallow DL,
    6. Brennan PA
    (2012) Noradrenaline-induced enhancement of oscillatory local field potentials in the mouse accessory olfactory bulb does not depend on disinhibition of mitral cells. Eur J Neurosci 35:1433–1445, doi:10.1111/j.1460-9568.2012.08070.x, pmid:22487171.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Liu S,
    2. Shipley MT
    (2008) Multiple conductances cooperatively regulate spontaneous bursting in mouse olfactory bulb external tufted cells. J Neurosci 28:1625–1639, doi:10.1523/JNEUROSCI.3906-07.2008, pmid:18272683.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Llinás R,
    2. Yarom Y
    (1986) Oscillatory properties of guinea-pig inferior olivary neurones and their pharmacological modulation: an in vitro study. J Physiol 376:163–182, doi:10.1113/jphysiol.1986.sp016147, pmid:3795074.
    OpenUrlCrossRefPubMed
  60. ↵
    1. Luo M,
    2. Fee MS,
    3. Katz LC
    (2003) Encoding pheromonal signals in the accessory olfactory bulb of behaving mice. Science 299:1196–1201, doi:10.1126/science.1082133, pmid:12595684.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Ma J,
    2. Lowe G
    (2004) Action potential backpropagation and multiglomerular signaling in the rat vomeronasal system. J Neurosci 24:9341–9352, doi:10.1523/JNEUROSCI.1782-04.2004, pmid:15496670.
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Maccaferri G,
    2. McBain CJ
    (1996) The hyperpolarization-activated current (Ih) and its contribution to pacemaker activity in rat CA1 hippocampal stratum oriens-alveus interneurones. J Physiol 497:119–130, doi:10.1113/jphysiol.1996.sp021754, pmid:8951716.
    OpenUrlCrossRefPubMed
  63. ↵
    1. Maex R,
    2. De Schutter E
    (1998) Synchronization of Golgi and granule cell firing in a detailed network model of the cerebellar granule cell layer. J Neurophysiol 80:2521–2537, pmid:9819260.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Maier N,
    2. Tejero-Cantero A,
    3. Dorrn AL,
    4. Winterer J,
    5. Beed PS,
    6. Morris G,
    7. Kempter R,
    8. Poulet JF,
    9. Leibold C,
    10. Schmitz D
    (2011) Coherent phasic excitation during hippocampal ripples. Neuron 72:137–152, doi:10.1016/j.neuron.2011.08.016, pmid:21982375.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Marder E,
    2. Goaillard JM
    (2006) Variability, compensation and homeostasis in neuron and network function. Nat Rev Neurosci 7:563–574, doi:10.1038/nrn1949, pmid:16791145.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Martin RL,
    2. Lee JH,
    3. Cribbs LL,
    4. Perez-Reyes E,
    5. Hanck DA
    (2000) Mibefradil block of cloned T-type calcium channels. J Pharmacol Exp Ther 295:302–308, pmid:10991994.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Marx M,
    2. Günter RH,
    3. Hucko W,
    4. Radnikow G,
    5. Feldmeyer D
    (2012) Improved biocytin labeling and neuronal 3D reconstruction. Nat Protoc 7:394–407, doi:10.1038/nprot.2011.449, pmid:22301777.
    OpenUrlCrossRefPubMed
  68. ↵
    1. McCleskey EW,
    2. Fox AP,
    3. Feldman DH,
    4. Cruz LJ,
    5. Olivera BM,
    6. Tsien RW,
    7. Yoshikami D
    (1987) Omega-conotoxin: direct and persistent blockade of specific types of calcium channels in neurons but not muscle. Proc Natl Acad Sci U S A 84:4327–4331, doi:10.1073/pnas.84.12.4327, pmid:2438698.
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Meredith M
    (1994) Chronic recording of vomeronasal pump activation in awake behaving hamsters. Physiol Behav 56:345–354, doi:10.1016/0031-9384(94)90205-4, pmid:7938248.
    OpenUrlCrossRefPubMed
  70. ↵
    1. Meredith M,
    2. O'Connell RJ
    (1979) Efferent control of stimulus access to the hamster vomeronasal organ. J Physiol 286:301–316, doi:10.1113/jphysiol.1979.sp012620, pmid:439028.
    OpenUrlCrossRefPubMed
  71. ↵
    1. Metz AE,
    2. Jarsky T,
    3. Martina M,
    4. Spruston N
    (2005) R-type calcium channels contribute to afterdepolarization and bursting in hippocampal CA1 pyramidal neurons. J Neurosci 25:5763–5773, doi:10.1523/JNEUROSCI.0624-05.2005, pmid:15958743.
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Migliore M,
    2. Hines ML,
    3. Shepherd GM
    (2005) The role of distal dendritic gap junctions in synchronization of mitral cell axonal output. J Comput Neurosci 18:151–161, doi:10.1007/s10827-005-6556-1, pmid:15714267.
    OpenUrlCrossRefPubMed
  73. ↵
    1. Mizuseki K,
    2. Sirota A,
    3. Pastalkova E,
    4. Buzsáki G
    (2009) Theta oscillations provide temporal windows for local circuit computation in the entorhinal-hippocampal loop. Neuron 64:267–280, doi:10.1016/j.neuron.2009.08.037, pmid:19874793.
    OpenUrlCrossRefPubMed
  74. ↵
    1. Moore GP,
    2. Perkel DH,
    3. Segundo JP
    (1966) Statistical analysis and functional interpretation of neuronal spike data. Annu Rev Physiol 28:493–522, doi:10.1146/annurev.ph.28.030166.002425, pmid:5323104.
    OpenUrlCrossRefPubMed
  75. ↵
    1. Peña F,
    2. Parkis MA,
    3. Tryba AK,
    4. Ramirez JM
    (2004) Differential contribution of pacemaker properties to the generation of respiratory rhythms during normoxia and hypoxia. Neuron 43:105–117, doi:10.1016/j.neuron.2004.06.023, pmid:15233921.
    OpenUrlCrossRefPubMed
  76. ↵
    1. Perez-Reyes E,
    2. Cribbs LL,
    3. Daud A,
    4. Lacerda AE,
    5. Barclay J,
    6. Williamson MP,
    7. Fox M,
    8. Rees M,
    9. Lee JH
    (1998) Molecular characterization of a neuronal low-voltage-activated T-type calcium channel. Nature 391:896–900, doi:10.1038/36110, pmid:9495342.
    OpenUrlCrossRefPubMed
  77. ↵
    1. Ramirez JM,
    2. Koch H,
    3. Garcia AJ 3rd.,
    4. Doi A,
    5. Zanella S
    (2011) The role of spiking and bursting pacemakers in the neuronal control of breathing. J Biol Phys 37:241–261, doi:10.1007/s10867-011-9214-z, pmid:22654176.
    OpenUrlCrossRefPubMed
  78. ↵
    1. Richter JM,
    2. Schaefer M,
    3. Hill K
    (2014) Riluzole activates TRPC5 channels independently of PLC activity. Br J Pharmacol 171:158–170, doi:10.1111/bph.12436, pmid:24117252.
    OpenUrlCrossRefPubMed
  79. ↵
    1. Rieke F,
    2. de Ruyter van Steveninck R,
    3. Bialek W
    (1997) in Spikes: exploring the neural code, eds Sejnowski TJ, Poggio TA (MIT, Cambridge, MA).
  80. ↵
    1. Romano SA,
    2. Pietri T,
    3. Pérez-Schuster V,
    4. Jouary A,
    5. Haudrechy M,
    6. Sumbre G
    (2015) Spontaneous neuronal network dynamics reveal circuit's functional adaptations for behavior. Neuron 85:1070–1085, doi:10.1016/j.neuron.2015.01.027, pmid:25704948.
    OpenUrlCrossRefPubMed
  81. ↵
    1. Schroeder CE,
    2. Lakatos P
    (2009) Low-frequency neuronal oscillations as instruments of sensory selection. Trends Neurosci 32:9–18, doi:10.1016/j.tins.2008.09.012, pmid:19012975.
    OpenUrlCrossRefPubMed
  82. ↵
    1. Shadlen MN,
    2. Newsome WT
    (1998) The variable discharge of cortical neurons: implications for connectivity, computation, and information coding. J Neurosci 18:3870–3896, pmid:9570816.
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Shao Z,
    2. Puche AC,
    3. Kiyokage E,
    4. Szabo G,
    5. Shipley MT
    (2009) Two GABAergic intraglomerular circuits differentially regulate tonic and phasic presynaptic inhibition of olfactory nerve terminals. J Neurophysiol 101:1988–2001, doi:10.1152/jn.91116.2008, pmid:19225171.
    OpenUrlAbstract/FREE Full Text
  84. ↵
    1. Shpak G,
    2. Zylbertal A,
    3. Yarom Y,
    4. Wagner S
    (2012) Calcium-activated sustained firing responses distinguish accessory from main olfactory bulb mitral cells. J Neurosci 32:6251–6262, doi:10.1523/JNEUROSCI.4397-11.2012, pmid:22553031.
    OpenUrlAbstract/FREE Full Text
  85. ↵
    1. Simerly RB
    (2002) Wired for reproduction: organization and development of sexually dimorphic circuits in the mammalian forebrain. Annu Rev Neurosci 25:507–536, doi:10.1146/annurev.neuro.25.112701.142745, pmid:12052919.
    OpenUrlCrossRefPubMed
  86. ↵
    1. Smear M,
    2. Shusterman R,
    3. O'Connor R,
    4. Bozza T,
    5. Rinberg D
    (2011) Perception of sniff phase in mouse olfaction. Nature 479:397–400, doi:10.1038/nature10521, pmid:21993623.
    OpenUrlCrossRefPubMed
  87. ↵
    1. Smith RS,
    2. Araneda RC
    (2010) Cholinergic modulation of neuronal excitability in the accessory olfactory bulb. J Neurophysiol 104:2963–2974, doi:10.1152/jn.00446.2010, pmid:20861438.
    OpenUrlAbstract/FREE Full Text
  88. ↵
    1. Stocker M
    (2004) Ca(2+)-activated K+ channels: molecular determinants and function of the SK family. Nat Rev Neurosci 5:758–770, doi:10.1038/nrn1516, pmid:15378036.
    OpenUrlCrossRefPubMed
  89. ↵
    1. Stowers L,
    2. Logan DW
    (2010) Sexual dimorphism in olfactory signaling. Curr Opin Neurobiol 20:770–775, doi:10.1016/j.conb.2010.08.015, pmid:20833534.
    OpenUrlCrossRefPubMed
  90. ↵
    1. Su H,
    2. Alroy G,
    3. Kirson ED,
    4. Yaari Y
    (2001) Extracellular calcium modulates persistent sodium current-dependent burst-firing in hippocampal pyramidal neurons. J Neurosci 21:4173–4182, pmid:11404402.
    OpenUrlAbstract/FREE Full Text
  91. ↵
    1. Sugai T,
    2. Miyazawa T,
    3. Yoshimura H,
    4. Onoda N
    (2005) Developmental changes in oscillatory and slow responses of the rat accessory olfactory bulb. Neuroscience 134:605–616, doi:10.1016/j.neuroscience.2005.05.010, pmid:16019155.
    OpenUrlCrossRefPubMed
  92. ↵
    1. Tazerart S,
    2. Vinay L,
    3. Brocard F
    (2008) The persistent sodium current generates pacemaker activities in the central pattern generator for locomotion and regulates the locomotor rhythm. J Neurosci 28:8577–8589, doi:10.1523/JNEUROSCI.1437-08.2008, pmid:18716217.
    OpenUrlAbstract/FREE Full Text
  93. ↵
    1. Tolokh II,
    2. Fu X,
    3. Holy TE
    (2013) Reliable sex and strain discrimination in the mouse vomeronasal organ and accessory olfactory bulb. J Neurosci 33:13903–13913, doi:10.1523/JNEUROSCI.0037-13.2013, pmid:23966710.
    OpenUrlAbstract/FREE Full Text
  94. ↵
    1. Veitinger S,
    2. Veitinger T,
    3. Cainarca S,
    4. Fluegge D,
    5. Engelhardt CH,
    6. Lohmer S,
    7. Hatt H,
    8. Corazza S,
    9. Spehr J,
    10. Neuhaus EM,
    11. Spehr M
    (2011) Purinergic signalling mobilizes mitochondrial Ca2+ in mouse Sertoli cells. J Physiol 589:5033–5055, doi:10.1113/jphysiol.2011.216309, pmid:21859825.
    OpenUrlCrossRefPubMed
  95. ↵
    1. Wagner S,
    2. Gresser AL,
    3. Torello AT,
    4. Dulac C
    (2006) A multireceptor genetic approach uncovers an ordered integration of VNO sensory inputs in the accessory olfactory bulb. Neuron 50:697–709, doi:10.1016/j.neuron.2006.04.033, pmid:16731509.
    OpenUrlCrossRefPubMed
  96. ↵
    1. Yamada-Hanff J,
    2. Bean BP
    (2013) Persistent sodium current drives conditional pacemaking in CA1 pyramidal neurons under muscarinic stimulation. J Neurosci 33:15011–15021, doi:10.1523/JNEUROSCI.0577-13.2013, pmid:24048831.
    OpenUrlAbstract/FREE Full Text
  97. ↵
    1. Zaman T,
    2. Lee K,
    3. Park C,
    4. Paydar A,
    5. Choi JH,
    6. Cheong E,
    7. Lee CJ,
    8. Shin HS
    (2011) Cav2.3 channels are critical for oscillatory burst discharges in the reticular thalamus and absence epilepsy. Neuron 70:95–108, doi:10.1016/j.neuron.2011.02.042, pmid:21482359.
    OpenUrlCrossRefPubMed
  98. ↵
    1. Zylbertal A,
    2. Kahan A,
    3. Ben-Shaul Y,
    4. Yarom Y,
    5. Wagner S
    (2015) Prolonged intracellular Na+ dynamics govern electrical activity in accessory olfactory bulb mitral cells. PLoS Biol 13:e1002319, doi:10.1371/journal.pbio.1002319, pmid:26674618.
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 36 (11)
Journal of Neuroscience
Vol. 36, Issue 11
16 Mar 2016
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Interdependent Conductances Drive Infraslow Intrinsic Rhythmogenesis in a Subset of Accessory Olfactory Bulb Projection Neurons
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Interdependent Conductances Drive Infraslow Intrinsic Rhythmogenesis in a Subset of Accessory Olfactory Bulb Projection Neurons
Monika Gorin, Chryssanthi Tsitoura, Anat Kahan, Katja Watznauer, Daniela R. Drose, Martijn Arts, Rudolf Mathar, Simon O'Connor, Ileana L. Hanganu-Opatz, Yoram Ben-Shaul, Marc Spehr
Journal of Neuroscience 16 March 2016, 36 (11) 3127-3144; DOI: 10.1523/JNEUROSCI.2520-15.2016

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Interdependent Conductances Drive Infraslow Intrinsic Rhythmogenesis in a Subset of Accessory Olfactory Bulb Projection Neurons
Monika Gorin, Chryssanthi Tsitoura, Anat Kahan, Katja Watznauer, Daniela R. Drose, Martijn Arts, Rudolf Mathar, Simon O'Connor, Ileana L. Hanganu-Opatz, Yoram Ben-Shaul, Marc Spehr
Journal of Neuroscience 16 March 2016, 36 (11) 3127-3144; DOI: 10.1523/JNEUROSCI.2520-15.2016
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • intrinsic oscillations
  • mitral cell
  • olfactory bulb
  • pacemaker
  • rhythmogenesis
  • vomeronasal system

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Articles

Cellular/Molecular

  • KCNQ2/3 Gain-of-Function Variants and Cell Excitability: Differential Effects in CA1 versus L2/3 Pyramidal Neurons
  • Tonic activation of NR2D-containing NMDARs exacerbates dopaminergic neuronal loss in MPTP-injected Parkinsonian mice
  • Axon Initial Segment GABA Inhibits Action Potential Generation throughout Periadolescent Development
Show more Cellular/Molecular
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.