- State Key Laboratory of Medical Neurobiology, School of Life Science, Human Phenome Institute, Institute of Brain Science, Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai, China
What type of principle features intrinsic inside of the fluctuated input signals could drive neurons with the maximal excitations is one of the crucial neural coding issues. In this article, we examined both experimentally and theoretically the cortical neuronal responsivity (including firing rate and spike timing reliability) to input signals with different intrinsic correlational statistics (e.g., white-type noise, showed 1/f0 power spectrum, pink noise 1/f, and brown noises 1/f2) and different frequency ranges. Our results revealed that the response sensitivity and reliability of cortical neurons is much higher in response to 1/f noise stimuli with long-term correlations than 1/f0 with short-term correlations for a broad frequency range, and also higher than 1/f2 for all frequency ranges. In addition, we found that neuronal sensitivity diverges to opposite directions for 1/f noise comparing with 1/f0 white noise as a function of cutoff frequency of input signal. As the cutoff frequency is progressively increased from 50 to 1,000 Hz, the neuronal responsiveness increased gradually for 1/f noise, while decreased exponentially for white noise. Computational simulations of a general cortical model revealed that, neuronal sensitivity and reliability to input signal statistics was majorly dominated by fast sodium inactivation, potassium activation, and membrane time constants.
Introduction
For a signal, the inherent frequency structure shown in the Fourier frequency domain characterizes its second-order statistics. The power spectrum of various natural signals typically exhibits the power law 1/fβ in the frequency domain, with β close to one (Voss and Clarke, 1978; Gilden et al., 1995; Musha and Yamamoto, 1997; De Coensel et al., 2003). Moreover, this 1/f property within a specific frequency range is widely observed in neural activities at all levels, as evidenced in recordings of the membrane potential and current (Diba, 2004; Jacobson et al., 2005; Bédard et al., 2006; Yaron-Jakoubovitch, 2008; El Boustani et al., 2009), EEG (Novikov et al., 1997; Bhattacharya and Petsche, 2001; Bédard et al., 2006; Dehghani et al., 2010; Voytek et al., 2015), MEG (Novikov et al., 1997; Dehghani et al., 2010), LFPs (Bédard and Destexhe, 2009; Bedard et al., 2017; Maex, 2018), and fMRI signals (Bullmore et al., 2001; He, 2011; Ciuciu et al., 2014). In practice, the white noise (the power law 1/f0) with generally low cutoff frequencies is widely used to detect neuronal input-and-output functions (Sakai, 1992; Fairhall et al., 2001; Cook et al., 2007; Vilela and Lindner, 2009). Besides, 1/f2 noise is also observed in neural field potentials (Freeman and Zhai, 2009; Miller et al., 2009; Milstein et al., 2009; He et al., 2010; Halnes et al., 2016) and membrane currents under special conditions (Diba, 2004).
Mammalian sensory neural systems exhibit better responses to naturalistic signals rather than white-type noise signals in a specific frequency domain (Aertsen and Johannesma, 1981; Baddeley et al., 1997; de Ruyter van Steveninck et al., 1997; Yu et al., 2005; Garcia-Lazaro et al., 2006, 2011). It has been speculated that the 1/f property might be the key in shaping the neuronal function preference to the naturalistic input (Yu et al., 2005; Garcia-Lazaro et al., 2006, 2011); this preference may extend to the atomic level of neural organization, namely the single neuronal input-output function (Gal and Marom, 2013).
In addition, white and colored noises with various cutoff frequencies enhance the detection of weak signals by neuronal systems via stochastic resonance (SR) (Nozaki and Yamamoto, 1998,?; Nozaki et al., 1999b; Hutcheon and Yarom, 2000; Jia et al., 2001; Ruszczynski et al., 2001; Chizhov and Graham, 2008; Mino and Durand, 2008; Gutkin et al., 2009; Higgs and Spain, 2009; Sekine et al., 2009; Guo and Li, 2011; Sobie et al., 2011; Duan et al., 2014; Zhao et al., 2017). However, previous studies have regarded the signals with 1/fβ statistics as mere background noise. In particular, although the signal frequency range has been confirmed to significantly affect neuronal excitability (Nozaki et al., 1999b; Higgs and Spain, 2009), the exact role of the frequency range in neuronal responsiveness to the 1/fβ statistic remains unclear. As neurons generally encode information according to the rate and/or the precise timing of spikes (Nowak et al., 1997; Reinagel and Reid, 2000; Fellous et al., 2001; Brette and Guigon, 2003; Avissar et al., 2007; Freund and Cerquera, 2012), we sought to examine the neuronal firing rate and spike-timing reliability to input statistics. In this study, we focus on the effects of the signal frequency range to identify the neuronal responsiveness to 1/fβ (β = 0, 1, and 2) noises with respect to the firing rate and spike-timing reliability. Specifically, we reveal the mechanism underlying neuronal responsivity to 1/fβ and the frequency range using theoretical experiments, as this mechanism has not been clearly determined in previous model-based studies (Nozaki and Yamamoto, 1998; Nozaki et al., 1999b; Brunel and Latham, 2003; Mino and Durand, 2008; Sekine et al., 2009; Sobie et al., 2011; Ostojic et al., 2015; Schwalger et al., 2015; Zhao et al., 2017). We conducted in vitro whole-cell patch clamp recording experiments on mouse cortical pyramidal neurons to examine the neuronal firing rate and reliability to 1/f, 1/f0, and 1/f2 noises with various cutoff frequencies (Fcut). We have also carried out a set of computational simulations of a general Hodgkin-Huxley neuronal model (Yu et al., 2012) to reproduce our experimental observations, and revealed the critical factors underlying the neuronal responsiveness to second-order statistics at the cellular level.
Methods
Signal Production
Input noise stimuli of each type of 1/fβ–white noise (1/f0), pink noise (1/f), and brown noise (1/f2) were first generated digitally by computer programming. The noise stimuli were then filtered by different low-pass filter with aimed cutoff frequency range. Each signal intensity [represented by the standard deviation (SD)] was set to the defined value. All the above were done in Matlab R2017a software (Mathworks, USA), and then the signal was loaded to Micro 1401 (CED, UK), where the signals were converted from digital to analog. In the experimental study, a unique set of noise stimuli with designed cutoff frequencies and noise intensities were used as the input signals to recorded neurons. The three signal types have different slopes in PSDs (the white noise has a slope of 0 and 1/f noise has a slope of −1, while 1/f2 noise has a slope of −2), as shown in Figure 1A. The 1/f and 1/f2 noises have more power in the low frequency components and less power in the high frequency components, while white noise has equal energy at each frequency interval. All three types of signals have equal total power within the examined frequency range.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/449601/fncel-13-00311-HTML/image_m/fncel-13-00311-g001.jpg)
Figure 1. Neuronal firing rate in response to different types of input signal. (A) The PSD plots for white noise (blue), 1/f noise (red), and 1/f2 noise (black) are shown. Noises with a 1,000 Hz cutoff frequency are shown here. (B) In vitro patch clamp recording of a pyramidal cell. The scale bar indicates 50 μm. (C) Plot of the recorded membrane potentials in the resting state vs. time. Inset: The power spectrum density (PSD) plot of the resting potentials shows a 1/f property in a log-log plot. (D,E) Representative trace of input 1/f0 white noise (D) and 1/f noise (E) with different intensities over time. The signal produced at each intensity lasted for 1 s, followed by a 3 s no-stimulus interval. The bottom panel shows the membrane potential and action potential produced by a recorded pyramidal cell in response to the input signal shown in the top panel. (F) Neuronal responsive firing rate for white noise (blue), 1/f (red), and 1/f2 noise (black) with various intensities (82.37, 164.80, 247.23, 329.67, 412.10, and 494.53 pA) and a cutoff frequency of 1000 Hz. Fcut, cutoff frequency. (G) Neuronal responsive firing rate for white noise (blue), 1/f (red), and 1/f2 noise (black) at various cutoff frequencies (50, 100, 200, 500, and 1,000 Hz) with a intensity of 494.53 pA.
Brain Slice Preparation
The experimental procedures involving animal experiments in this study were approved by Animal Ethics Committee of Fudan University School of Life Science. Using 0.7% chloral hydrate, 14–28 day-old mice were anesthetized and brain slices prepared by a protective slicing and recovery method reported previously (Ting et al., 2014). Briefly, anesthetized mice were perfused intracardially with ice-cold oxygenated (95% O2, 5% CO2) NMDG-based cutting solution containing (in mM) 93 mM NMDG, 2.5 mM KCl, 1.2 mM NaH2PO4, 30 mM NaHCO3, 20 mM HEPES, 25 mM glucose, 5 mM sodium ascorbate, 2 mM thiourea, 3 mM sodium pyruvate, 10 mM MgSO4, 2 mM CaCl2, and 12 mM NAC (pH 7.3–7.4, 300–305 mOsm). Brains were carefully removed from the skull and cut coronally at a thickness of 300 μm with a vibratome VT1000S (Leica, Germany) in chilled oxygenated (95% O2, 5% CO2) NMDG-based cutting solution. Slices were initially recovered in NMDG-based cutting solution at 32°C for 10 min. Slices were then incubated in oxygenated (95% O2, 5% CO2) HEPES-modified solution containing (in mM) 94 mM NaCl, 2.5 mM KCl, 1.2 mM NaH2PO4, 30 mM NaHCO3, 20 mM HEPES, 25 mM glucose, 5 mM sodium ascorbate, 2 mM thiourea, 3 mM sodium pyruvate, 2 mM MgSO4, 2 mM CaCl2, and 6 mM NAC (pH 7.3–7.4, 300–305 mOsm) at room temperature for 30 min. Finally, slices were incubated in oxygenated (95% O2, 5% CO2) ACSF at room temperature for at least 1 h before recording. The ACSF contained 126 mM NaCl, 2.5 mM KCl, 2 mM MgSO4, 2 mM CaCl2, 26 mM NaHCO3, 1.25 mM NaH2PO4, and 12.5 mM dextrose (pH 7.3–7.4, 300–305 mOsm).
Electrophysiological Recordings
Whole-cell slice recordings were performed on the cell body of layer 5 pyramidal neurons in the prefrontal cortex (Figure 1B). In total, 12 neurons in 12 slices of six mice were examined. Oxygenated (95% O2, 5% CO2) ASCF was used as the recording solution. Recordings were conducted on an upright infrared-differential interference contrast (IR-DIC) microscope (Zeiss Axioskop 2 FS plus) with a recording solution temperature of 36 °C. Cortical slices were suspended on a net to allow an oxygenated solution to flow over both the upper and lower surfaces at a rate of 3–4 ml/min. The membrane potential in the whole-cell recordings was corrected for Donnan liquid junction potentials of 15 mV. The temperature was regulated by a Warner Instruments Corporation two channel temperature regulator (Model TC344B).
Whole-cell recordings from the soma were achieved with the help of a Multiclamp 700B amplifier (Axon Instruments, Union City, CA) and Micro 1401 converter. Pipettes had an impedance of 5–6 MΩ and were filled with an intracellular solution that contained 140 mM K-gluconate, 3 mM KCl, 2 mM MgCl2, 2 mM Na2ATP, 10 mM HEPES, and 0.2 mM EGTA, and the pH was adjusted to 7.3 with KOH–(270 mOsm). The sample rate of the membrane potential data was 25 kHz for both experiments and computational simulations. Those recorded neurons with input resistance within 100–200 MΩ were saved for analysis in this paper while others were discarded so that the studied neurons have the similar biophysical membrane properties.
Stimulation Mode
To experimentally examine the responsive firing rate to the input signal, the current signal with each intensity (σ = 82.37, 164.80, 247.23, 329.67, 412.10, and 494.53 pA) and each cutoff frequency (Fcut = 50, 100, 200, 500, and 1,000 Hz) was injected into neurons three times (Figures 1D,E). The signals (σ = 494.53 pA; Fcut = 50, 200, 500, and 1,000 Hz) in firing rate detection were also used for spike reliability examination, in which 50 repetitions of signals were injected into neurons. Every neuron in our study received all types of input signal, with its firing rate and reliability for each signal cutoff frequency examined. In the model simulation, the neuron received input signals (σ = 1.5, 3, 4.5, 6, 7.5, 9 μA/cm2; Fcut = 50, 200, 500, and 1,000 Hz) thrice for firing rate detection. The signal (σ = 9 pA; Fcut = 50, 200, 500, and 1,000 Hz) was also injected 50 times for the spike reliability examination. In each test, to simulate the background noise from external or intrinsic fluctuations of the neuron itself, a 1/f noise with 500 Hz cutoff frequency and an intensity of 0.2169 μA/cm2 was added to the input signal. The duration of stimuli signal of each cutoff frequency and intensity was 1 s.
Firing Rate and Spike Time Reliability
For each neuron under an input signal of each cutoff frequency and intensity, the spike number within the 1s stimuli duration was the firing rate. We averaged the firing rates in the three repetitions to get each neuronal firing rate value. For spike time reliability, 50 repetitions of a 1 s stimuli were injected into neurons, and cross covariance was calculated between all pair-wise combinations of trials (within stimulus size) within each neuron, on binary spike timing trains (with “1” representing an action potential with a 2 ms time bin and “0” representing a non-spiking neuron). The black-covariance function estimates the mean-removed cross-correlation between the two sequences of random processes, thus avoiding the contribution of the mean firing rate to the spike timing reliability. The resulting cross-covariance values at zero-lag (normalized by the average autocovariance function) for each neuron were used to quantify the spike reliability for each input signal statistic (Haider et al., 2010).
Hodgkin-Huxley-Style Cortical Neuronal Model
Three major ionic voltage-dependent currents were used in our cortical model: fast Na+, INa, fast K+, IK, and a leak current, IL. The equations describing the voltage and time dependence of the Na+ and K+ conductance have been reported in previous publications (McCormick and Huguenard, 1992) and the channel kinetics were modified based on models of cortical neurons (Mainen et al., 1995; Mainen and Sejnowski, 1996; Yu et al., 2008) and experimental studies (Huguenard et al., 1989; Colbert and Pan, 2002; Yu et al., 2008; Schmidt-Hieber and Bischofberger, 2010). The following equations describe the cortical axon single compartment model:
where the Q10 effect is described by Φ on regulating the temperature dependence of the biochemical reaction rate with Q10 = 2.3 (Frankenhaeuser and Moore, 1963; Matteson and Armstrong, 1982). The relationships between temperature and INa and IK activation and inactivation are not monotonic and vary in different species (Fohlmeister et al., 2010). The reverse potentials for Na+ and K+ currents were adjusted for change in temperature according to the Nernst equation (not shown). Similar results were obtained with a variety of values for Q10. For example, the use of a Q10 of 3 yielded similar results for spike efficiency and changes in spike rate with temperature. In our cortical model, Na+ kinetics were determined based on recent experimental observations (Kole et al., 2008). The parameters used were: membrane capacitance = 0.75 μF/cm2, gNa = 1950 pS/μm2, density of gK = 40 pS/μm2, and gleak = 0.25 pS/μm2, based on recent experimental results (Mainen and Sejnowski, 1995; Colbert and Pan, 2002; Kole et al., 2008; Hu et al., 2009; Fleidervish et al., 2010; Schmidt-Hieber and Bischofberger, 2010). The reversal potentials were VL = −70 mV, VNa = 60 mV, and VK = −90 mV for leak, sodium, and potassium channels, respectively.
Statistical Analysis
To detect the differential response under various signal intensities, cutoff frequencies, or signal types, we performed paired t-test and Wilcoxon rank sum test. First, Kolmogorov-Smirnov goodness-of-fit hypothesis test was done on the paired difference between the two compared data groups. Depending on whether the normality assumption held for the data, we used a paired t-test or the Wilcoxon test for statistical significant comparison. A p < 0.05 was considered statistically significant. The data presented in the figures are reported as the mean ± standard error, and the significant level are labeled by *p < 0.05, **p < 0.01, and ***p < 0.001, respectively.
Results
Neuronal Responsiveness to Input Stimuli
In the absence of external stimuli, neuronal resting membrane potentials (~-75 to −70 mV) showed fluctuating synaptic potentials that were received from somatic and dendritic synapses (Figure 1C). A power spectrum density (PSD) analysis of resting potentials displayed a typical 1/f property in a log-log plot (see the inset in Figure 1C), suggesting that the membrane potentials observed in vitro in neurons with synaptic inputs also exhibit a similar statistical property as neurons in vivo.
First, we examined the neuronal firing rate in response to the input signal (Figures 1D,E). As shown in Figure 1F and Supplementary Figure S1, the neuronal firing rate increases non-linearly as the noise intensity increases within a firing rate of 0–30 Hz. When the noise cutoff frequency is within 50 Hz, neurons are sensitive to white noise, firing with a high frequency at high noise intensities (Figure 1G and Supplementary Figure S1A). However, as the cutoff frequency increases, neurons lose their responsiveness to the white noise (Figure 1G and Supplementary Figure S1) and produce a low firing rate, even at a very high noise intensity levels and when the cutoff frequency is 1,000 Hz (Figures 1F,G). The neuronal firing rate increases minimally at a cutoff frequency of 1/f2 noise (Figure 1G). In contrast, the firing rate increases gradually with the cutoff frequency for 1/f noise (Figure 1G). When the cutoff frequency is >200 Hz, firing rates for 1/f and 1/f2 noise are higher than for the white noise (Figures 1F,G, and Supplementary Figures S1C,D). At low intensities, 1/f noise evokes lower firing rates than the 1/f2 noise; but at high intensities, 1/f noise evokes higher firing rates than the 1/f2 noise. The switch in firing rates emerges under various cutoff frequencies (Figure 1F and Supplementary Figure S1). Notably, at σ = 494.53 pA and a 1,000 Hz cutoff frequency, 1/f noise evokes the highest firing rate (25.47 ± 2.24 Hz) in all our recordings. Overall, 1/f noise evokes relatively high firing rates under all conditions, particularly at high input intensities, and high cutoff frequencies.
We next examined spike reliability (Figure 2A) at four cutoff frequencies (Fcut = 50, 200, 500, and 1,000 Hz). The neuronal responses are shown in the post-stimulus time histogram (PSTH) and raster plot (Figure 2A). Figures 2B,C show the spike-timing reliability, which quantifies the extent of repeatability of spike timing in response to the input signal (Mainen and Sejnowski, 1995; Haider et al., 2010) (calculated from the representative data shown in Figure 1A). As the cutoff frequency increases to above 200 Hz, the reliability of the response to white noise decreases significantly (Figure 2B). For both the 1/f and 1/f2 noises, the reliability increases as the cutoff frequency increases (Figure 2B). When the cutoff frequency is 50 Hz, the reliability for the 1/f noise is greater than the 1/f2 noise and less than that for the white noise (Figure 2B and Supplementary Figure S2A). When the cutoff frequency increases to 200 Hz, the reliability for white noise remains greater than for the 1/f and 1/f2 noises (Figure 2B and Supplementary Figure S2B). Interestingly, at 500 and 1,000 Hz, the 1/f noise evokes the greatest reliability (Figures 2B,C, and Supplementary Figure S2C), suggesting that neurons respond with relatively high reliability to the 1/f property in signals, particularly when the input signal has a high cutoff frequency.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/449601/fncel-13-00311-HTML/image_m/fncel-13-00311-g002.jpg)
Figure 2. Neuronal spike reliability to different types of input signal. (A) Raster plot of the responses from 50 trials (bottom panel) and PSTH (middle panel) for the repeated input stimulus (top panel) reflecting the reliability of cortical neuron firing patterns evoked by white noise (left panel), 1/f noise (middle panel), and 1/f2 noise (right panel). (B) Plot of reliability vs. cutoff frequencies for the three signal types. Fcut, cutoff frequency. (C) Plot of reliability vs. β for the three types of 1/fβ input signals with a cutoff frequency of 500 Hz.
Our data demonstrate that neurons respond to white noise with low sensitivity and reliability, except when the noise is composed of low frequency components only. When the noise intensity is weak, 1/f2 noise more readily evokes neuronal firing, but high firing rate is difficult to evoke. For the 1/f noise, neurons respond with relatively high sensitivity and high reliability under all conditions. In particular, at high cutoff frequencies or high input intensities, the 1/f noise evokes the highest firing rate and the most reliable neuronal response among all signal types tested. These findings indicate a clear neuronal preference for responses to 1/f noise with high cutoff frequencies.
Computer Simulations of Hodgkin-Huxley Model
A cortical Hodgkin-Huxley-type neuronal model was constructed to study the key factors and dynamic mechanisms underlying neuronal responsiveness to signals with different statistics. We only considered the fast sodium and potassium channels in the model and ignored other subtypes of sodium/potassium and calcium channels to identify the common mechanism underlying the neuronal preference for 1/f statistics.
The model neuron first reproduced the similar responses to input signals with various noise intensities described in our experimental study (Figure 3A). For cutoff frequency of 1,000 Hz, the firing rates of neuronal model to white noise is significantly lower than that of 1/f and 1/f2 noise, and 1/f noise evokes the highest firing rate for the high intensities (Figures 3B,C). When the cutoff frequency increases from 50 to 1,000 Hz, the neuronal firing rate of model to input white noise decreases gradually. However, it increases slightly for 1/f noise, while it keeps almost invariant for 1/f2 noise (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/449601/fncel-13-00311-HTML/image_m/fncel-13-00311-g003.jpg)
Figure 3. Model neuronal response to different types of input signal. (A) The top panel shows white noise (left) and 1/f noise (right) with step-changed variance introduced into the neuronal model. The bottom panel shows the membrane potentials of the model neuron in response to the input. The noise intensity shown here is 3, 6, and 9 μA/cm2. (B) Histogram of the neuronal firing rate vs. input signal SD for white noise (blue), 1/f noise (red), and 1/f2 noise (black) at the cutoff frequency of 1,000 Hz. (C) Histogram of the neuronal firing rate vs. cutoff frequencies (50, 100, 200, 500, and 1,000 Hz) for white noise (blue), 1/f (red), and 1/f2 noise (black) when the noise SD was 9 μA/cm2. (D) Plot of model neuronal reliability vs. cutoff frequencies (50, 100, 200, 500, and 1,000 Hz) for white noise (blue), 1/f (red), and 1/f2 noise (black). Fcut, cutoff frequency. (E) Plot of reliability vs. β for the three types of 1/fβ input signals with a cutoff frequency of 500 Hz.
We also evaluated the spike reliability of the model neuron based on the evoked spike trains by repeatedly feeding the model neuron with each type of input signal with the different cutoff frequencies. Neuronal reliability for white noise is high only at a low cutoff frequency, while it decreases significantly with an increase in cutoff frequency (Figure 3D). On the contrary, neuronal response reliability to 1/f noise keeps at a much higher reliability level for almost the whole range of cutoff frequency (Figure 3D). In addition, the neuronal response reliability to 1/f2 noise is lower than the other two type of inputs (Figure 3E). However, its reliability increases as the cutoff frequency of 1/f2 noise increases (Figure 3D). These model simulation results are consistent with the experimental observations.
Next, we investigated the critical factors of neurons affecting the neuronal preference to 1/f signals. Considering that the responseness of neuronal model is dominated by the sodium and potassium channel kinetics, especially the channel opening velocity constant (e.g., αm, αh, αn) and channel closing velocity constant (e.g., βm βh, and βn) (Yu et al., 2012), we systematically changed the values of these parameters to study their effects on neuronal preference to inputs. We have also varied the values of membrane input resistance Rinput and the membrane capacitance Cm to understand the effect of the passive membrane property on the neuronal responseness. First we examined the effect of open (αm) and close (βm) velocity constants of sodium activation variable. Although neuronal firing rate could be significantly affected by an increase or decrease of the values of αm and βm, the changing amount of the firing rate is almost same for all the three signal types (Supplementary Figures S3A,B).
However, when the value of αh is halved, the neuronal firing rate as a function of noise intensity (Figure 4A1) or cutoff frequency (Figure 4A2) decreases dramatically for both 1/f and 1/f2 noise, but slightly for white noise. Figure 4A3 summarizes the contribution effect. When theαh value increases, the neuronal firing rate increases gradually to a saturation level for both 1/f and 1/f2 noises while keeps invariant for white noise input. This suggests that αh plays a role in the neuronal preference to 1/f signals. Next, we varied the value of βh and observed almost no effect on neuronal firing rate (please see Supplementary Figure S3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/449601/fncel-13-00311-HTML/image_m/fncel-13-00311-g004.jpg)
Figure 4. Effects of ionic channel kinetic constants on the responsive firing rate of model neurons. (A1–C1) Plots of firing rate vs. intensity for white noise (blue), 1/f (red), and 1/f2 noise (black) for a normal neuron (dotted line) compared with the results from the model neuron in the case of 0.5 αh (solid line, A1), 2 αn (solid line, B1), 2 βn (solid line, C1). Fcut = 1,000 Hz. (A2–C2) Plots of firing rate vs. cutoff frequency for white noise (blue), 1/f (red), and 1/f2 noise (black) for a normal neuron (dotted line) compared with the results from the model neuron in the case of 0.5 αh (solid line, A2), 2 αn (solid line, B2), 2 βn (solid line, C2). Noise σ = 9 μA/cm2. (A3–C3) Plots of firing rate vs. times of αh (solid line, A3), αn (solid line, B3), and βn (solid line, C3) for white noise (blue), 1/f (red), and 1/f2 noise (black) for a normal neuron (dotted line) compared with the results from the model neuron. Fcut = 1,000 Hz. Noise σ = 9 μA/cm2.
Next, we examined the effect of αn and βn for the potassium channel activation variable. When the value of αn was doubled (Figures 4B1,B2), the neuronal firing rate as a function of noise intensity (Figure 4B1) or cutoff frequency (Figure 4B2) decreased dramatically for both 1/f and 1/f2 noise, but slightly for white noise. Figure 4B3 summarizes the contribution effect, and shows that the firing rates of neuronal model decrease gradually when the αn value increased for both 1/f and 1/f2 noises while deceased slightly for the white noise input. This suggests that αn also contributes to the neuronal preference to 1/f signals.
Moreover, when βn is increased, firing rate for 1/f2 noise increases greater than for 1/f noise, much more than for white noise, at each noise intensity (Figure 4C1) as well as for cutoff frequency (Figure 4C2). As a result, in the case of large βn, neuronal firing rate for 1/f2 noise is larger than 1/f noise, while in the case of small βn, neuronal firing rate for 1/f2 noise is lower than 1/f noise (Figure 4C3). This indicated that relatively large βn value may dominate the neuronal preference to 1/f and 1/f2 noise. For the much small βn value, neuronal preference to 1/f noise decreases dramatically.
Based on the equation, the membrane time constant τc = Rinput * Cm, τc could be changed by either changing Rinput or Cm. Halved Rinput increased neuronal firing rate to all the three type of input, and the neuronal sensitivity reached the maximum to white noise in the low cutoff frequency around 200 Hz, which is higher than to 1/f and 1/f2 noise stimuli (Figures 5A1,A2). Doubled Cm decreases the firing rate and makes the neuronal firing rate for 1/f2 noise larger than for 1/f noise at various noise intensity (Figure 5B1) and each cutoff frequency, respectively (Figure 5B2). As shown in Figures 5A,B, both Cm and Rinput changed the neuronal firing rate for white noise more than for 1/f noise, while the change for 1/f2 noise is not that strong. This is contrary to the situation of changing βn. As βnand Rinput or Cm have different weighted effect on neuronal response, they may act together to support the neuronal preferential response to 1/f noise with a broad frequency range. So we concluded βn, Rinput and Cm play important roles in neuronal responsive firing rate to 1/fβ noises.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/449601/fncel-13-00311-HTML/image_m/fncel-13-00311-g005.jpg)
Figure 5. Effects of membrane time constants on the responsive firing rate of model neurons. (A1,B1) Plots of firing rate vs. intensity for white noise (blue), 1/f (red), and 1/f2 noise (black) for a normal neuron (dotted line) compared with the results from the model neuron in the case of 0.5 Rinput (solid line, A1), 2 Cm (solid line, B1). Fcut = 1,000 Hz. (A2,B2) Plots of firing rate vs. cutoff frequency for white noise (blue), 1/f (red), and 1/f2 noise (black) for a normal neuron (dotted line) compared with the results from the model neuron in the case of 0.5 Rinput (solid line, A2), 2 Cm (solid line, B2). Noise σ = 9 μA/cm2. (A3,B3) Plots of firing rate vs. times of Rinput (solid line, A3), and Cm (solid line, B3) for white noise (blue), 1/f (red), and 1/f2 noise (black) for a normal neuron (dotted line) compared with the results from the model neuron. Fcut = 1000 Hz. Noise σ = 9 μA/cm2.
We also examined changes in the spike reliability in response to noises with altered neuronal biophysical parameters. As shown in Figure 6A, an increased value ofβm results in a significant decreased spike reliability for both 1/f noise (in all frequency range) and white noise (mainly with cutoff frequency higher than 200 Hz). On the contrary, increased βm results in significantly increased spike reliability for 1/f2 noise. With doubled αh, αn, and βn (Figures 6B–D) and decreased Rinput (Figure 6E), neuronal spike reliability for 1/f2 noise decreases significantly for most of frequency range. These effects for 1/f2 noise are stronger than for 1/f noise, and much stronger for white noise. On the contrary, doubled Cm decreases the reliability significantly for each type of input, especially for the white noise (Figure 6F). Noted here, there is no clear change observed for the other rate constants (i.e., αm and βh).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/449601/fncel-13-00311-HTML/image_m/fncel-13-00311-g006.jpg)
Figure 6. Effects of ionic and membrane time constants on the spike reliability of model neurons. (A–F) Plots of reliability vs. cutoff frequency for white noise (blue), 1/f (red) and 1/f2 noise (black) for a normal neuron (dotted line, A–F) compared with the results from the model neuron in the case of 2 βm (solid line, A), 2 αh (solid line, B), 2 αn (solid line, C), 2 βn (solid line, D), 0.5 Rinput (solid line, E), and 2 Cm (solid line, F). Fcut, cutoff frequency. Rinput, membrane input resistance. Cm, membrane capacitance.
Discussion
Recordings of neuronal membranes exhibit time-dependent voltage fluctuations (Destexhe et al., 2003; El Boustani et al., 2009), which have not been extensively studied in the last few decades due to the lack of a functional understanding of the underlying noise composition. In addition, 1/fβ noise has mainly been considered as background noise in previous studies (Mainen and Sejnowski, 1995; Nozaki and Yamamoto, 1998; Nozaki et al., 1999b; Jia et al., 2001; Ruszczynski et al., 2001; Richardson et al., 2003; Chizhov and Graham, 2008; Köndgen et al., 2008; Mino and Durand, 2008; Gutkin et al., 2009; Higgs and Spain, 2009; Sekine et al., 2009; Guo and Li, 2011; Sobie et al., 2011; Duan et al., 2014; Zhao et al., 2017). Here, for the first time, we directly examine cortical pyramidal neuronal responsiveness to 1/fβ (β = 0, 1, and 2) input signals. Our results shed light on the mechanism by which the statistical structure of input signals affects the dynamics of spike generation in cortical pyramidal cells.
The Effect of the Frequency Range on the 1/fβ Noise Response
We compared neuronal responses to 1/f0, 1/f, and 1/f2 noises with various cutoff frequencies (50–1,000 Hz) and observed that cortical pyramidal cells exhibit a substantial loss of reliability and sensitivity to white noise when the cutoff frequency exceeds 200 Hz. Neuronal responsiveness to 1/f is also better than 1/f2 at various cutoff frequencies. In previous studies, white noise has been widely utilized to detect neuronal input and output functions (Sakai, 1992; Fairhall et al., 2001; Cook et al., 2007; Vilela and Lindner, 2009). Based on our results, 1/f noise should be a much better probe for determining the response properties of neurons at various input frequency ranges, instead of white noise or 1/f2 noise, and the signal cutoff frequency should be carefully established in future studies.
Our work suggests that the neuronal firing rate and reliability for 1/f noise are both enhanced with an increase in the cutoff frequency. Thus, each frequency interval has an effect on neuronal responses, and future work should examine the optimum signal frequency range for neuronal responsiveness and derive the neuronal input-output function for each frequency interval. Although a previous study has found that injecting white noise to the DC component increases neuronal spike reliability (Mainen and Sejnowski, 1995), our findings with zero mean signals imply that neuronal excitability may determine the effects of signal frequency range, and ultimately modulate neuronal response.
Preferential Neuronal Response to 1/f Noise With a Large Frequency Range
The 1/f noise is a common phenomenon in nature (Bak et al., 1987; Gilden et al., 1995; Musha and Yamamoto, 1997; Novikov et al., 1997; De Coensel et al., 2003; Bédard et al., 2006; He, 2014). In the brain, EEG and ECoG recordings reveal that when neuronal populations exhibit the 1/f characteristic, the neural network is capable of highly efficient information processing (Soma et al., 2003; Lin and Chen, 2005; Shew and Plenz, 2013; Wood et al., 2016). At the cellular level, 1/f noise was found to originate from intact network inputs (El Boustani et al., 2009), and can enhance neuronal excitability and the stochastic resonanGarciace effect (Nozaki and Yamamoto, 1998; Nozaki et al., 1999a,b; El Boustani et al., 2009). Based on neuronal sensitivity to low-frequency sine wave signals (Hutcheon et al., 1996; Hunter et al., 1998; Volgushev et al., 1998; Fellous et al., 2001; Yu et al., 2001; Brumberg, 2002; Richardson et al., 2003; Köndgen et al., 2008; Levi et al., 2015), an individual neuron uses signals with a low frequency range to process information. However, recent experiments have observed high frequency components in the recorded membrane potentials of cortical neurons in vivo (El Boustani et al., 2009; Bedard et al., 2017).
Our results, for the first time, show that neurons respond preferentially to the 1/f noise in large frequency ranges, with respect to the firing rate and spike timing reliability. Because the rate and temporal encoding are two major encoding methods for neurons (Nowak et al., 1997; Reinagel and Reid, 2000; Fellous et al., 2001; Brette and Guigon, 2003; Avissar et al., 2007; Freund and Cerquera, 2012) and neuronal activities at all levels show 1/f characteristics in vivo (Bédard et al., 2006; El Boustani et al., 2009; Freeman and Zhai, 2009; Milstein et al., 2009; Dehghani et al., 2010; He et al., 2010; Ciuciu et al., 2012; Pettersen et al., 2014; Voytek et al., 2015), our results may indicate the positive effects of the 1/f property and large frequency range on the neuronal responses for normal physiological functioning. In addition, since the 1/f type of signal substantially drives neuronal response, our findings suggest new operating rules for synaptic transmission, neuronal plasticity and other activities relying on neuronal firing. It is likely that future studies of neuronal interactions and network activities will discover additional distinct effects of 1/f probe signals on neural information process and neural computation.
Neuronal Dynamics Underlying the 1/f Preference
Given the low-pass filter property of neuronal membrane, once the signal contain more power in low frequency range, it may induce more neuronal firing, as seen from the neuronal response to 1/f compared with to white noise (Fellous et al., 2001; Brumberg, 2002; Levi et al., 2015). However, this is in contradiction with our finding that 1/f noise evokes higher firing rate than 1/f2 noise in a large frequency rage. Especially, we observed that neuronal firing rate for 1/f noise increases with the cutoff frequency of 1/f noise increases. In addition, it was reported previously that the refractory period endows neuron with high-pass filtering (Nozaki et al., 1999a), increasing the response complexity. We performed numerical simulations by using the Hodgkin-Huxley model to gain deep insights into the mechanism underlying the neuronal responsive preference for the 1/f noise. We found that βn, Rinput, and Cm determine the neuronal preference to input signal with different type of statistics. Figure 4 shows that the effect of changing βn on the neuronal firing rate to 1/f2 is larger to 1/f, much larger than to 1/f0 noise. On the contrary, the effects of changing Rinput or Cm on the neuronal firing rate is stronger to 1/f0 noise than 1/f, even stronger than to 1/f2 inputs. Here, close rate constant of potassium activation βn seems to form a high-pass filter effect, while both Rinput and Cm act as low-pass filter effect. Their appropriately combined action may result in neuronal responsive preference for 1/f noise with a broad frequency range.
In addition, we found the neuronal responsive reliability for 1/fβ signal types is majorly dominated by several key factors of ion channels and membrane time constants, respectively. Specifically, as shown in results section, the spike reliability for white noise is mainly determined by Cm. Spike reliability for 1/f and 1/f2 noise is mainly dominated by sodium activation close rate constant (αm), inactivation open rate constant (αh) and both open (αn), and close (βn) rate constants of potassium channel, as well as membrane time constant. With doubled αh, αn, and βn (Figures 6B–D) and decreased Rinput (Figure 6E), neuronal spike reliability for 1/f2 noise decreases significantly for most of frequency range. These effects for 1/f2 noise are stronger than for 1/f noise, and much stronger for white noise. These results revealed that the membrane capacitance dominates the neuronal preference to low frequency component, while the combined ion channel kinetics dominates the neuronal preference to high frequency components in the 1/fβ type noise stimuli.
A prevailing hypothesis in neuronal response is that the stochastic opening and closing of individual ion channels endows cortical neurons an inherent noise. When the correlation time of the external input matches the time scale of the inherent noise, the neuronal responsiveness is maximized. Indeed, previous studies have demonstrated the existence of the optimum time scale of input signals for neuronal spiking reliability (Galán et al., 2008; McGinley et al., 2015). The neuron-preferred 1/f signal in our finding is consistent with this hypothesis.
In addition, cortical neurons are also enriched in other subtypes of sodium, potassium, and calcium channels, which may play important roles in precisely controlling the cellular sensitivity to the input signals with different temporal correlations (Wang et al., 2003). Additionally, a balance between excitatory and inhibitory synaptic inputs may adjust neuronal excitability to different signal statistics (Brunel et al., 2001; Chance et al., 2002; Wang, 2010). The dendritic morphology of the cell may even enhance the neuronal sensitivity to some frequency components within input signals (Eyal et al., 2014; Ostojic et al., 2015). Future studies should examine how the synaptic balance and the interaction between the neuronal intrinsic dynamics and modulations from the recurrent network contribute to the neuronal responses to input signals with various higher order statics and cutoff frequencies.
In summary, as the signal cutoff frequency progressively increases from 50 to 1,000 Hz, the neuronal responsive firing rate and reliability increase for 1/f noise, but decrease for white noise. Ion channel kinetic and membrane time constants endow neurons with a preferential response for 1/f noise with high cutoff frequencies. These results suggest that the 1/f noise is important in determining the computational rules and operating principles of cortical circuits.
Author Contributions
YY supervised the research. YY, GQ, and BF designed the research. GQ and XF performed the experimental research. GQ and YY wrote the paper. All authors performed data analysis and reviewed the manuscript.
Funding
YY thanks for the support from the National Natural Science Foundation of China (81761128011, 31571070), Shanghai Science and Technology Committee support (16410722600), the program for the Professor of Special Appointment (Eastern Scholar SHH1140004) at Shanghai Institutions of Higher Learning, the Research Fund for the Doctoral Program of Higher Education of China (1322051), and Omics-based precision medicine of epilepsy entrusted by the Key Research Project of the Ministry of Science and Technology of China (Grant No. 2016YFC0904400) for their support.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2019.00311/full#supplementary-material
References
Aertsen, A. M., and Johannesma, P. I. (1981). A comparison of the spectro-temporal sensitivity of auditory neurons to tonal and natural stimuli. Biol. Cybern. 42, 145–156. doi: 10.1007/BF00336732
Avissar, M., Furman, A. C., Saunders, J. C., and Parsons, T. D. (2007). Adaptation reduces spike-count reliability, but not spike-timing precision, of auditory nerve responses. J. Neurosci. 27, 6461–6472. doi: 10.1523/JNEUROSCI.5239-06.2007
Baddeley, R., Abbott, L. F., Booth, M. C., Sengpiel, F., Freeman, T., Wakeman, E. A., et al. (1997). Responses of neurons in primary and inferior temporal visual cortices to natural scenes. Proc. Biol. Sci. 264, 1775–1783. doi: 10.1098/rspb.1997.0246
Bak, P., Tang, C., and Wiesenfeld, K. (1987). Self-organized criticality: An explanation of the 1/f noise. Phys. Rev. Lett. 59, 381–384. doi: 10.1103/PhysRevLett.59.381
Bédard, C., and Destexhe, A. (2009). Macroscopic models of local field potentials and the apparent 1/f noise in brain activity. Biophys. J. 96, 2589–2603. doi: 10.1016/j.bpj.2008.12.3951
Bedard, C., Gomes, J.-M., Bal, T., and Destexhe, A. (2017). A framework to reconcile frequency scaling measurements, from intracellular recordings, local-field potentials, up to EEG and MEG signals. J. Integr. Neurosci. 16, 3–18. doi: 10.3233/JIN-160001
Bédard, C., Kröger, H., and Destexhe, A. (2006). Does the 1/f frequency scaling of brain signals reflect self-organized critical states? Phys. Rev. Lett. 97:118102. doi: 10.1103/PhysRevLett.97.118102
Bhattacharya, J., and Petsche, H. (2001). Universality in the brain while listening to music. Proc. Biol. Sci. 268, 2423–2433. doi: 10.1098/rspb.2001.1802
Brette, R., and Guigon, E. (2003). Reliability of spike timing is a general property of spiking model neurons. Neural Comput. 15, 279–308. doi: 10.1162/089976603762552924
Brumberg, J. C. (2002). Firing pattern modulation by oscillatory input in supragranular pyramidal neurons. Neuroscience 114, 239–246. doi: 10.1016/S0306-4522(02)00252-X
Brunel, N., Chance, F. S., Fourcaud, N., and Abbott, L. F. (2001). Effects of synaptic noise and filtering on the frequency response of spiking neurons. Phys. Rev. Lett. 86, 2186–2189. doi: 10.1103/PhysRevLett.86.2186
Brunel, N., and Latham, P. E. (2003). Firing rate of the noisy quadratic integrate-and-fire neuron. Neural Comput. 15, 2281–2306. doi: 10.1162/089976603322362365
Bullmore, E., Long, C., Suckling, J., Fadili, J., Calvert, G., Zelaya, F., et al. (2001). Colored noise and computational inference in neurophysiological (fMRI) time series analysis: resampling methods in time and wavelet domains. Hum. Brain Mapp. 12, 61–78. doi: 10.1002/1097-0193(200102)12:23.0.CO;2-W.3.0.CO;2-W
Chance, F. S., Abbott, L. F., and Reyes, A. D. (2002). Gain modulation from background synaptic input. Neuron 35, 773–782. doi: 10.1016/S0896-6273(02)00820-6
Chizhov, A. V., and Graham, L. J. (2008). Efficient evaluation of neuron populations receiving colored-noise current based on a refractory density method. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 77:011910. doi: 10.1103/PhysRevE.77.011910
Ciuciu, P., Abry, P., and He, B. J. (2014). Interplay between functional connectivity and scale-free dynamics in intrinsic fMRI networks. NeuroImage 95, 248–263. doi: 10.1016/j.neuroimage.2014.03.047
Ciuciu, P., Varoquaux, G., Abry, P., Sadaghiani, S., and Kleinschmidt, A. (2012). Scale-free and multifractal time dynamics of fMRI signals during rest and task. Front. Physiol. 3:186. doi: 10.3389/fphys.2012.00186
Colbert, C. M., and Pan, E. (2002). Ion channel properties underlying axonal action potential initiation in pyramidal neurons. Nat. Neurosci. 5, 533–538. doi: 10.1038/nn857
Cook, E. P., Wilhelm, A. C., Guest, J. A., Liang, Y., Masse, N. Y., and Colbert, C. M. (2007). The neuronal transfer function: contributions from voltage- and time-dependent mechanisms. Prog. Brain Res. 165, 1–12. doi: 10.1016/S0079-6123(06)65001-2
De Coensel, B., Botteldooren, D., and De Muer, T. (2003). “1/f noise in rural and urban soundscapes. Acta Acustica United Acustica 89, 287–295. doi: 10.1134/1.1560388
de Ruyter van Steveninck, R. R., Lewen, G. D., Strong, S. P., Koberle, R., and Bialek, W. (1997). Reproducibility and variability in neural spike trains. Science 275, 1805–1808. doi: 10.1126/science.275.5307.1805
Dehghani, N., Bédard, C., Cash, S. S., Halgren, E., and Destexhe, A. (2010). Comparative power spectral analysis of simultaneous elecroencephalographic and magnetoencephalographic recordings in humans suggests non-resistive extracellular media. J. Comput. Neurosci. 29, 405–421. doi: 10.1007/s10827-010-0263-2
Destexhe, A., Rudolph, M., and Paré, D. (2003). The high-conductance state of neocortical neurons in vivo. Nat. Rev. Neurosci. 4, 739–751. doi: 10.1038/nrn1198
Diba, K. (2004). Intrinsic noise in cultured hippocampal neurons: experiment and modeling. J. Neurosci. 24, 9723–9733. doi: 10.1523/JNEUROSCI.1721-04.2004
Duan, F., Chapeau-Blondeau, F., and Abbott, D. (2014). Stochastic resonance with colored noise for neural signal detection. PLoS ONE 9:e91345. doi: 10.1371/journal.pone.0091345
El Boustani, S., Marre, O., Béhuret, S., Baudot, P., Yger, P., Bal, T., et al. (2009). Network-state modulation of power-law frequency-scaling in visual cortical neurons. PLoS Comput. Biol. 5:e1000519. doi: 10.1371/journal.pcbi.1000519
Eyal, G., Mansvelder, H. D., de Kock, C. P. J., and Segev, I. (2014). Dendrites impact the encoding capabilities of the axon. J. Neurosci. Off. J. Soc. Neurosci. 34, 8063–8071. doi: 10.1523/JNEUROSCI.5431-13.2014
Fairhall, A. L., Lewen, G. D., Bialek, W., and de Ruyter Van Steveninck, R. R. (2001). Efficiency and ambiguity in an adaptive neural code. Nature 412, 787–792. doi: 10.1038/35090500
Fellous, J. M., Houweling, A. R., Modi, R. H., Rao, R. P., Tiesinga, P. H., and Sejnowski, T. J. (2001). Frequency dependence of spike reliability in cortical pyramidal cells and interneurons. J. Neurophysiol. 85, 1782–1787. doi: 10.1152/jn.2001.85.4.1782
Fleidervish, I. A., Lasser-Ross, N., Gutnick, M. J., and Ross, W. N. (2010). Na+ imaging reveals little difference in action potential-evoked Na+ influx between axon and soma. Nat. Neurosci. 13, 852–860. doi: 10.1038/nn.2574
Fohlmeister, J. F., Cohen, E. D., and Newman, E. A. (2010). Mechanisms and distribution of ion channels in retinal ganglion cells: using temperature as an independent variable. J. Neurophysiol. 103, 1357–1374. doi: 10.1152/jn.00123.2009
Frankenhaeuser, B., and Moore, L. E. (1963). The effect of temperature on the sodium and potassium permeability changes in myelinated nerve fibres of xenopus laevis. J. Physiol. 169, 431–437. doi: 10.1113/jphysiol.1963.sp007269
Freeman, W. J., and Zhai, J. (2009). Simulated power spectral density (PSD) of background electrocorticogram (ECoG). Cogn. Neurodyn. 3, 97–103. doi: 10.1007/s11571-008-9064-y
Freund, J. A., and Cerquera, A. (2012). How spurious correlations affect a correlation-based measure of spike timing reliability. Neurocomputing 81, 97–103. doi: 10.1016/j.neucom.2011.10.014
Gal, A., and Marom, S. (2013). Entrainment of the intrinsic dynamics of single isolated neurons by natural-like input. J. Neurosci. 33, 7912–7918. doi: 10.1523/JNEUROSCI.3763-12.2013
Galán, R. F., Ermentrout, G. B., and Urban, N. N. (2008). Optimal time scale for spike-time reliability: theory, simulations, and experiments. J. Neurophysiol. 99, 277–283. doi: 10.1152/jn.00563.2007
Garcia-Lazaro, J. A., Ahmed, B., and Schnupp, J. W. (2006). Tuning to natural stimulus dynamics in primary auditory cortex. Curr. Biol. CB 16, 264–271. doi: 10.1016/j.cub.2005.12.013
Garcia-Lazaro, J. A., Ahmed, B., and Schnupp, J. W. (2011). Emergence of tuning to natural stimulus statistics along the central auditory pathway. PLoS ONE 6:e22584. doi: 10.1371/journal.pone.0022584
Gilden, D. L., Thornton, T., and Mallon, M. W. (1995). 1/f noise in human cognition. Science 267, 1837–1839. doi: 10.1126/science.7892611
Guo, D., and Li, C. (2011). Signal propagation in feedforward neuronal networks with unreliable synapses. J. Comput. Neurosci. 30, 567–587. doi: 10.1007/s10827-010-0279-7
Gutkin, B. S., Jost, J., and Tuckwell, H. C. (2009). Inhibition of rhythmic neural spiking by noise: the occurrence of a minimum in activity with increasing noise. Naturwissenschaften 96, 1091–1097. doi: 10.1007/s00114-009-0570-5
Haider, B., Krause, M. R., Duque, A., Yu, Y., Touryan, J., Mazer, J. A., et al. (2010). Synaptic and network mechanisms of sparse and reliable visual cortical activity during nonclassical receptive field stimulation. Neuron 65, 107–121. doi: 10.1016/j.neuron.2009.12.005
Halnes, G., Mäki-Marttunen, T., Keller, D., Pettersen, K. H., Andreassen, O. A., and Einevoll, G. T. (2016). Effect of ionic diffusion on extracellular potentials in neural tissue. PLoS Comput. Biol. 12:e1005193. doi: 10.1371/journal.pcbi.1005193
He, B. J. (2011). Scale-free properties of the functional magnetic resonance imaging signal during rest and task. J. Neurosci. Off. J. Soc. Neurosci. 31, 13786–13795. doi: 10.1523/JNEUROSCI.2111-11.2011
He, B. J. (2014). Scale-free brain activity: past, present, and future. Trends Cogn. Sci. 18, 480–487. doi: 10.1016/j.tics.2014.04.003
He, B. J., Zempel, J. M., Snyder, A. Z., and Raichle, M. E. (2010). The temporal structures and functional significance of scale-free brain activity. Neuron 66, 353–369. doi: 10.1016/j.neuron.2010.04.020
Higgs, M. H., and Spain, W. J. (2009). Conditional bursting enhances resonant firing in neocortical layer 2-3 pyramidal neurons. J. Neurosci. 29, 1285–1299. doi: 10.1523/JNEUROSCI.3728-08.2009
Hu, W., Tian, C., Li, T., Yang, M., Hou, H., and Shu, Y. (2009). Distinct contributions of Nav1.6 and Nav1.2 in action potential initiation and backpropagation. Nat. Neurosci. 12, 996–1002. doi: 10.1038/nn.2359
Huguenard, J. R., Hamill, O. P., and Prince, D. A. (1989). Sodium channels in dendrites of rat cortical pyramidal neurons. Proc. Natl. Acad. Sci. U. S. A. 86, 2473–2477. doi: 10.1073/pnas.86.7.2473
Hunter, J. D., Milton, J. G., Thomas, P. J., and Cowan, J. D. (1998). Resonance effect for neural spike time reliability. J. Neurophysiol. 80, 1427–1438. doi: 10.1152/jn.1998.80.3.1427
Hutcheon, B., Miura, R. M., and Puil, E. (1996). Subthreshold membrane resonance in neocortical neurons. J. Neurophysiol. 76, 683–697. doi: 10.1152/jn.1996.76.2.683
Hutcheon, B., and Yarom, Y. (2000). Resonance, oscillation and the intrinsic frequency preferences of neurons. Trends Neurosci. 23, 216–222. doi: 10.1016/S0166-2236(00)01547-2
Jacobson, G. A., Diba, K., Yaron-Jakoubovitch, A., Oz, Y., Koch, C., Segev, I., et al. (2005). Subthreshold voltage noise of rat neocortical pyramidal neurones. J. Physiol. 564, 145–160. doi: 10.1113/jphysiol.2004.080903
Jia, Y., Zheng, X., Hu, X., and Li, J. (2001). Effects of colored noise on stochastic resonance in a bistable system subject to multiplicative and additive noise. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 63:031107. doi: 10.1103/PhysRevE.63.031107
Kole, M. H., Ilschner, S. U., Kampa, B. M., Williams, S. R., Ruben, P. C., and Stuart, G. J. (2008). Action potential generation requires a high sodium channel density in the axon initial segment. Nat. Neurosci. 11, 178–186. doi: 10.1038/nn2040
Köndgen, H., Geisler, C., Fusi, S., Wang, X.-J., Lüscher, H.-R., and Giugliano, M. (2008). The dynamical response properties of neocortical neurons to temporally modulated noisy inputs in vitro. Cereb. Cortex 18, 2086–2097. doi: 10.1093/cercor/bhm235
Levi, R., Akanyeti, O., Ballo, A., and Liao, J. C. (2015). Frequency response properties of primary afferent neurons in the posterior lateral line system of larval zebrafish. J. Neurophysiol. 113, 657–668. doi: 10.1152/jn.00414.2014
Lin, M., and Chen, T. (2005). Self-organized criticality in a simple model of neurons based on small-world networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 71:016133. doi: 10.1103/PhysRevE.71.016133
Maex, R. (2018). An interneuron circuit reproducing essential spectral features of field potentials. Neural Comput. 30, 1296–1322. doi: 10.1162/NECO_a_01068
Mainen, Z. F., Joerges, J., Huguenard, J. R., and Sejnowski, T. J. (1995). A model of spike initiation in neocortical pyramidal neurons. Neuron 15, 1427–1439. doi: 10.1016/0896-6273(95)90020-9
Mainen, Z. F., and Sejnowski, T. J. (1995). Reliability of spike timing in neocortical neurons. Science 268, 1503–1506. doi: 10.1126/science.7770778
Mainen, Z. F., and Sejnowski, T. J. (1996). Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 382, 363–366. doi: 10.1038/382363a0
Matteson, D. R., and Armstrong, C. M. (1982). Evidence for a population of sleepy sodium channels in squid axon at low temperature. J. Gen. Physiol. 79, 739–758. doi: 10.1085/jgp.79.5.739
McCormick, D. A., and Huguenard, J. R. (1992). A model of the electrophysiological properties of thalamocortical relay neurons. J. Neurophysiol. 68, 1384–1400. doi: 10.1152/jn.1992.68.4.1384
McGinley, M. J., David, S. V., and McCormick, D. A. (2015). Cortical membrane potential signature of optimal states for sensory signal detection. Neuron 87, 179–192. doi: 10.1016/j.neuron.2015.05.038
Miller, K. J., Sorensen, L. B., Ojemann, J. G., and den Nijs, M. (2009). Power-law scaling in the brain surface electric potential. PLoS Comput. Biol. 5:e1000609. doi: 10.1371/journal.pcbi.1000609
Milstein, J., Mormann, F., Fried, I., and Koch, C. (2009). Neuronal shot noise and Brownian 1/f2 behavior in the local field potential. PLoS ONE 4:e4338. doi: 10.1371/journal.pone.0004338
Mino, H., and Durand, D. M. (2008). Stochastic resonance can induce oscillation in a recurrent Hodgkin-Huxley neuron model with added Gaussian noise. Conf Proc IEEE Eng Med Biol Soc. 2008, 2457–2460. doi: 10.1109/IEMBS.2008.4649697
Musha, T., and Yamamoto, M. (1997). “1/f fluctuations in biological systems,” in Proceedings of the 19th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. ‘Magnificent Milestones and Emerging Opportunities in Medical Engineering’ (Cat. No. 97CH36136), Vol. 6 (Chicago, IL: IEEE), 2692–2697. doi: 10.1109/IEMBS.1997.756890
Novikov, E., Novikov, A., Shannahoff-Khalsa, D., and Schwartz, B. (1997). Scale-similar activity in the brain. Phys. Rev. 56, R2387–R2389. doi: 10.1103/PhysRevE.56.R2387
Nowak, L. G., Sanchez-Vives, M. V., and McCormick, D. A. (1997). Influence of low and high frequency inputs on spike timing in visual cortical neurons. Cereb. Cortex 7, 487–501. doi: 10.1093/cercor/7.6.487
Nozaki, D., Collins, J. J., and Yamamoto, Y. (1999a). Mechanism of stochastic resonance enhancement in neuronal models driven by 1/f noise. Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Top. 60, 4637–4644. doi: 10.1103/PhysRevE.60.4637
Nozaki, D., Marder, D. J., Grigg, P., and Collins, J. J. (1999b). Effects of colored noise on stochastic resonance in sensory neurons. Phys. Rev. Lett. 82, 2402–2405. doi: 10.1103/PhysRevLett.82.2402
Nozaki, D., and Yamamoto, Y. (1998). Enhancement of stochastic resonance in a FitzHugh-Nagumo neuronal model driven by colored noise. Phys. Lett. A 243, 281–287. doi: 10.1016/S0375-9601(98)00247-3
Ostojic, S., Szapiro, G., Schwartz, E., Barbour, B., Brunel, N., and Hakim, V. (2015). Neuronal morphology generates high-frequency firing resonance. J. Neurosci. Off. J. Soc. Neurosci. 35, 7056–7068. doi: 10.1523/JNEUROSCI.3924-14.2015
Pettersen, K. H., Lindén, H., Tetzlaff, T., and Einevoll, G. T. (2014). Power laws from linear neuronal cable theory: power spectral densities of the soma potential, soma membrane current and single-neuron contribution to the EEG. PLoS Comput. Biol. 10:e1003928. doi: 10.1371/journal.pcbi.1003928
Reinagel, P., and Reid, R. C. (2000). Temporal coding of visual information in the thalamus. J. Neurosci. Off. J. Soc. Neurosci. 20, 5392–5400. doi: 10.1523/JNEUROSCI.20-14-05392.2000
Richardson, M. J., Brunel, N., and Hakim, V. (2003). From subthreshold to firing-rate resonance. J. Neurophysiol. 89, 2538–2554. doi: 10.1152/jn.00955.2002
Ruszczynski, P. S., Kish, L. B., and Bezrukov, S. M. (2001). Noise-assisted traffic of spikes through neuronal junctions. Chaos Woodbury N. 11, 581–586. doi: 10.1063/1.1379308
Sakai, H. M. (1992). White-noise analysis in neurophysiology. Physiol. Rev. 72, 491–505. doi: 10.1152/physrev.1992.72.2.491
Schmidt-Hieber, C., and Bischofberger, J. (2010). Fast sodium channel gating supports localized and efficient axonal action potential initiation. J. Neurosci. Off. J. Soc. Neurosci. 30, 10233–10242. doi: 10.1523/JNEUROSCI.6335-09.2010
Schwalger, T., Droste, F., and Lindner, B. (2015). Statistical structure of neural spiking under non-Poissonian or other non-white stimulation. J. Comput. Neurosci. 39, 29–51. doi: 10.1007/s10827-015-0560-x
Sekine, M., Mino, H., and Durand, D. M. (2009). Noise induced oscillations in recurrent neural networks. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2009, 1521–1524. doi: 10.1109/IEMBS.2009.5333069
Shew, W. L., and Plenz, D. (2013). The functional benefits of criticality in the cortex. Neuroscientist 19, 88–100. doi: 10.1177/1073858412445487
Sobie, C., Babul, A., and de Sousa, R. (2011). Neuron dynamics in the presence of 1/f noise. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 83:051912. doi: 10.1103/PhysRevE.83.051912
Soma, R., Nozaki, D., Kwak, S., and Yamamoto, Y. (2003). 1/f noise outperforms white noise in sensitizing baroreflex function in the human brain. Phys. Rev. Lett. 91:078101. doi: 10.1103/PhysRevLett.91.078101
Ting, J. T., Daigle, T. L., Chen, Q., and Feng, G. (2014). “Acute brain slice methods for adult and aging animals: application of targeted patch clamp analysis and optogenetics,” in Patch-Clamp Methods and Protocols, eds. M. Martina and S. Taverna (New York, NY: Springer New York), 221–242. Available online at: http://link.springer.com/10.1007/978-1-4939-1096-0_14 (accessed April 24, 2016). doi: 10.1007/978-1-4939-1096-0_14
Vilela, R. D., and Lindner, B. (2009). Are the input parameters of white noise driven integrate and fire neurons uniquely determined by rate and CV? J. Theor. Biol. 257, 90–99. doi: 10.1016/j.jtbi.2008.11.004
Volgushev, M., Chistiakova, M., and Singer, W. (1998). Modification of discharge patterns of neocortical neurons by induced oscillations of the membrane potential. Neuroscience 83, 15–25. doi: 10.1016/S0306-4522(97)00380-1
Voss, R. G., and Clarke, J. (1978). “1/f noise”in music: music from 1/f noise. J. Acoust. Soc. Am. 63, 258–263. doi: 10.1121/1.381721
Voytek, B., Kramer, M. A., Case, J., Lepage, K. Q., Tempesta, Z. R., Knight, R. T., et al. (2015). Age-related changes in 1/f neural electrophysiological noise. J. Neurosci. Off. J. Soc. Neurosci. 35, 13257–13265. doi: 10.1523/JNEUROSCI.2332-14.2015
Wang, X.-J. (2010). Neurophysiological and computational principles of cortical rhythms in cognition. Physiol. Rev. 90, 1195–1268. doi: 10.1152/physrev.00035.2008
Wang, X.-J., Liu, Y., Sanchez-Vives, M. V., and McCormick, D. A. (2003). Adaptation and temporal decorrelation by single neurons in the primary visual cortex. J. Neurophysiol. 89, 3279–3293. doi: 10.1152/jn.00242.2003
Wood, T., Osredkar, D., Puchades, M., Maes, E., Falck, M., Flatebø, T., et al. (2016). Treatment temperature and insult severity influence the neuroprotective effects of therapeutic hypothermia. Sci. Rep. 6:23430. doi: 10.1038/srep23430
Yaron-Jakoubovitch, A. (2008). A paradoxical isopotentiality: a spatially uniform noise spectrum in neocortical pyramidal cells. Front. Cell. Neurosci. 2:3. doi: 10.3389/neuro.03.003.2008
Yu, Y., Hill, A. P., and McCormick, D. A. (2012). Warm body temperature facilitates energy efficient cortical action potentials. PLoS Comput. Biol. 8:e1002456. doi: 10.1371/journal.pcbi.1002456
Yu, Y., Liu, F., Wang, J., and Wang, W. (2001). Spike timing precision for a neuronal array with periodic signal. Phys. Lett. 282, 23–30. doi: 10.1016/S0375-9601(01)00167-0
Yu, Y., Romero, R., and Lee, T. S. (2005). Preference of sensory neural coding for 1/f signals. Phys. Rev. Lett. 94:108103. doi: 10.1103/PhysRevLett.94.108103
Yu, Y., Shu, Y., and McCormick, D. A. (2008). Cortical action potential backpropagation explains spike threshold variability and rapid-onset kinetics. J. Neurosci. 28, 7260–7272. doi: 10.1523/JNEUROSCI.1613-08.2008
Keywords: 1/fβ noise, cortical neurons, patch clamp recording, long-term correlation, Hodgkin-Huxley model, response reliability
Citation: Qu G, Fan B, Fu X and Yu Y (2019) The Impact of Frequency Scale on the Response Sensitivity and Reliability of Cortical Neurons to 1/fβ Input Signals. Front. Cell. Neurosci. 13:311. doi: 10.3389/fncel.2019.00311
Received: 22 January 2019; Accepted: 25 June 2019;
Published: 11 July 2019.
Edited by:
Mario Eduardo Guido, Center for Research in Biological Chemistry Córdoba (CIQUIBIC), ArgentinaReviewed by:
Richardson N. Leão, Federal University of Rio Grande do Norte, BrazilThierry Ralph Nieus, Luigi Sacco Hospital, Italy
Copyright © 2019 Qu, Fan, Fu and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuguo Yu, yuyuguo@fudan.edu.cn
†These authors have contributed equally to this work