Human Synapses Show a Wide Temporal Window for Spike-Timing-Dependent Plasticity

Throughout our lifetime, activity-dependent changes in neuronal connection strength enable the brain to refine neural circuits and learn based on experience. Synapses can bi-directionally alter strength and the magnitude and sign depend on the millisecond timing of presynaptic and postsynaptic action potential firing. Recent findings on laboratory animals have shown that neurons can show a variety of temporal windows for spike-timing-dependent plasticity (STDP). It is unknown what synaptic learning rules exist in human synapses and whether similar temporal windows for STDP at synapses hold true for the human brain. Here, we directly tested in human slices cut from hippocampal tissue removed for surgical treatment of deeper brain structures in drug-resistant epilepsy patients, whether adult human synapses can change strength in response to millisecond timing of pre- and postsynaptic firing. We find that adult human hippocampal synapses can alter synapse strength in response to timed pre- and postsynaptic activity. In contrast to rodent hippocampal synapses, the sign of plasticity does not sharply switch around 0-ms timing. Instead, both positive timing intervals, in which presynaptic firing preceded the postsynaptic action potential, and negative timing intervals, in which postsynaptic firing preceded presynaptic activity down to −80 ms, increase synapse strength (tLTP). Negative timing intervals between −80 to −130 ms induce a lasting reduction of synapse strength (tLTD). Thus, similar to rodent synapses, adult human synapses can show spike-timing-dependent changes in strength. The timing rules of STDP in human hippocampus, however, seem to differ from rodent hippocampus, and suggest a less strict interpretation of Hebb's predictions.

plasticity (STDP) depending on parameters such as brain area and neuron type (Froemke et al., 2005;Wittenberg and Wang, 2006;Meredith et al., 2007;Caporale and Dan, 2008). However, it is unknown what synaptic learning rules exist in human synapses and whether similar temporal windows for STDP at synapses hold true for the human brain.
There are few studies on synaptic LTP and LTD in humans. With field potential recordings from hippocampal and neocortical tissue excised from human patients it was shown that high frequency stimulation (100 Hz) induces LTP in human synapses (Chen et al., 1996;Beck et al., 2000). Blocking NMDA receptors with APV prevents LTP induction, indicating that plasticity of human synapses shares molecular mechanisms with animal models. Low frequency stimulation (1 Hz) resulted in LTD (Chen et al., 1996), showing that the strength of human synapses can be regulated bidirectionally. Indirect evidence suggests that coincident millisecond timing of activity is likely to govern synaptic changes in humans as well (Stefan et al., 2000;Wolters et al., 2003). Pairing transcranial magnetic stimulation (TMS) of motor cortex with peripheral nerve stimulation in vivo can alter motor-evoked potentials in muscles, and precise timing determines the sign of this plasticity (Wolters

IntroductIon
One of the central questions in neuroscience is how memories are formed and stored in the human brain. From a large number of studies on laboratory animals it is known that learning and memory are most likely mediated by activity-dependent neuronal circuit modifications resulting from synaptic plasticity (Bliss and Lomo, 1973;Cooke and Bliss, 2006;Whitlock et al., 2006;Letzkus et al., 2007). The ability to shape neuronal connections in an activitydependent manner enables the brain to functionally refine neural circuits in response to sensory experience and adapt to changing environments (for review see Caporale and Dan, 2008). In line with predictions made by Hebb (1949), synapse strength can be modified depending on the millisecond timing of action potential firing and the sign of synaptic plasticity depends on the spike order of presynaptic and postsynaptic neurons (Levy and Steward, 1983;Gustafsson et al., 1987;Bell et al., 1997;Magee and Johnston, 1997;Markram et al., 1997). By varying the timing and order of preand postsynaptic spiking, it was found that critical time windows exist for synaptic modification on the order of tens of milliseconds Poo, 1998, 2001). In recent years, it has become clear that diversity exists of temporal windows for spike-timing-dependent  et al., 2003, 2005). Timed TMS of motor cortex and peripheral stimulation may induce synaptic plasticity in cortical circuits, but a direct demonstration thereof is lacking. In this study, we set out to directly test whether human synapses can alter strength in response to millisecond timing of pre-and postsynaptic firing. In addition, we explored the temporal window for STDP at these synapses. We find that in contrast to rodent hippocampus, adult human hippocampal synapses show a wide temporal window for STDP.

MaterIals and Methods huMan hIppocaMpal slIce preparatIon
All procedures on human tissue were performed with the approval of the Medical Ethical Committee of the VU University Medical Center and in accordance with Dutch license procedures and the declaration of Helsinki. Human slices were cut from hippocampal tissue that had to be removed for the surgical treatment of deeper brain structures for epilepsy with written informed consent of the patients (aged 20-66 years) prior to surgery. Anesthesia was induced with intravenous fentanyl 1-3 μg/kg and a bolus dose of propofol (2-10 mg/kg) and was maintained with remyfentanyl 250 μg/kg/ min and propofol 4-12 mg/kg. After resection, the hippocampal tissue was placed within 30 s in ice-cold artificial cerebrospinal fluid (ACSF) slicing solution which contained in (mM): 110 choline chloride, 26 NaHCO 3 , 10 d-glucose, 11.6 sodium ascorbate, 7 MgCl 2 , 3.1 sodium pyruvate, 2.5 KCl, 1.25 NaH 2 PO 4 , and 0.5 CaCl 2 −300 mOsm. (Bureau et al., 2006) and transported to the neurophysiology laboratory, which is located within 200 m distance from the operating room. The transition time between resection of the tissue and the start of preparing slices was less than 15 min.
Hippocampal slices (300 μm) were prepared in ice-cold slicing solution, and were then transferred to holding chambers in which they were stored for 30 min at 34 o C and for 30 min at room temperature before recording in ACSF which contained (in mM): NaCl 125; KCl 3; NaH 2 PO 4 1.25; MgSO 4 2; CaCl 2 2; NaHCO 3 26; glucose 10, bubbled with carbogen gas (95% O 2 /5% CO 2 ). electrophysIology Hippocampal slices were visualized using either infrared differential interference contrast (IR-DIC) microscopy or Hoffman modulation contrast. After the whole cell configuration was established, membrane potential responses to steps of current injection were used to classify each cell electrophysiologically. Cells were loaded with biocytin through the recording pipette for post hoc identification. All experiments were performed at 32-35 o C. Although the hippocampal tissue was resected from the brains of epilepsy patients, none of the neurons recorded from showed spontaneous epileptiform spiking activity. Resting membrane potentials were in line with previous reports on recordings from human neurons (-64±5 mV for Pyramidal and -65±7 mV for Non-pyramidal cells). All experiments were performed in the absence of blockers of GABAergic synaptic transmission.
Post hoc visualization and neuron identification using biocytin labeling was performed as described previously (Horikawa and Armstrong, 1988;de Kock et al., 2007). Image-stacks were created using a 20× air objective using Surveyor software (Objective Imaging Ltd., Stow cum Quy, Cambridge UK available from Chromaphor, Bottrop, Germany) and ImageJ. Pyramidal and nonpyramidal neurons were classified based on morphological and electrophysiological criteria. Pyramidal neurons typically had an input resistance below 100 MΩ; non-pyramidal neurons typically had an input resistance above 100 MΩ (Tables 2 and 3).

spIke-tIMIng-dependent plastIcIty
Excitatory postsynaptic potential (EPSPs) were evoked every 7 s using an extracellular stimulation electrode positioned at approximately 100-150 μm along the cell's somatodendritic axis (Figures 1A,B). The slope of the initial 2 ms of the EPSP was analyzed to ensure that the data reflected only the monosynaptic component of each experiment (Froemke et al., 2005;Couey et al., 2007). Synaptic gain was measured as the percent change in EPSP slope when comparing the average in the period 20-30 min postconditioning to the average baseline EPSP slope. During the induction protocol spike-timings were measured from the onset of the evoked EPSP to the peak of the postsynaptic AP. Mean baseline EPSP slopes were averaged from at least 30 sweeps (with amplitudes in the range of 6-7.5 mV). During the conditioning period pre-postsynaptic stimulus pairing was repeated 40-50 times, with a 7 s (0.14 Hz) interval between each pairing. Two distinct conditioning protocols were used, two spikes with 10-ms interval (tested in 17 cells) and single spike (tested in two cells). During experiments, cell input resistance was monitored throughout by applying a −10 to −100 pA, 200-500 ms hyperpolarizing pulse at the end of each sweep, additional determination of input resistance was done by extracting an I-V curve from each cell, in the beginning of the experiment, after pairing and at the but it has not been tested directly whether adult human synapses can change strength as a result of the precise timing of pre-and postsynaptic spiking activity. To test this, we made whole-cell recordings from pyramidal neurons and non-pyramidal cells of adult human hippocampus (20-66 years of age) and stimulated glutamatergic inputs by extracellular stimulation (Figure 1). During whole-cell recordings, neurons were labeled with biocytin for post hoc morphological identification. Hippocampal pyramidal neurons and non-pyramidal cells had distinct characteristic morphologies ( Figure 1). As described for pyramidal neurons and non-pyramidal cells in human association cortex (Foehring et al., 1991), hippocampal pyramidal neurons, and non-pyramidal cells also showed distinct basic electrophysiological properties (Tables 1-3; Figure 1). Passive and active cell properties measured in human hippocampal cells were in some aspects comparable to values obtained from rodent cells (Staff et al., 2000;Mercer et al., 2007;Routh et al., 2009). For example pyramidal cell input resistance in rat is 65.6 ± 4.4 MΩ, in mouse (C57BL/6) is 65.4 ± 1.7 MΩ, while in human pyramidals it tends to be slightly lower: 49.1 ± 8.7 MΩ. Input resistance of wide arbor basket cells in rat CA2 is 111.8 ± 36.7 MΩ, while in human, non-pyramidal cells we found a slightly higher average input resistance of 141.7 ± 6MΩ. Pyramidal cell AP amplitude agree well with reported values for rat pyramidal cells (humans: 113.8 ± 8.0 mV; rat: 112.0 ± 9.0 mV). Human cells however, have slower rate of rise compared with CA1 rat pyramidal cells. Rat: 381 ± 18 mV/ms and human: 231.3 ± 38.4 mV/ms, but similar end of the recording. Experiments were not included in the analysis if the cell input resistance varied by more than ± 30% during the experiment, which was the case for 13 out of 32 recordings. The t-test and the Mann-Whitney U test were used to assess significance on changes in slope for each experiment. Data are given as mean ± SEM, with p < 0.05 as minimum for statistical significance. Experiments lasted 20-70 min and where interrupted if the conditions of the patch degenerated or if the cell displayed repeated spiking after potentiation or failures after depression. When recording lasted less than 30 min (two of 19 recordings included in analysis), the average of baseline was compared to the average of all points in the post-pairing period.
The amplitude of the ADP component of the voltage response was measured from its maxima. The time constant of the decay of the ADP response back to the resting potential, τ decay , was estimated from the fitting to a single exponential, that is, where "A" is a dummy scaling factor.
Reversing the order of presynaptic and postsynaptic action potential firing such that the postsynaptic neuron fires before the presynaptic stimulation induces synaptic depression in rodent hippocampal synapses (Bi and Poo, 1998;Nishiyama et al., 2000;Wittenberg and Wang, 2006). The temporal window for spiketiming synaptic depression (tLTD) in these synapses is around 20 ms. To test whether adult human hippocampal synapses show tLTD in response to negative timing intervals, postsynaptic action potentials were induced 10-80 ms before the presynaptic stimulus during pairing (Figure 3). In contrast to rodent hippocampal synapses, these intervals induced a robust increase in EPSP amplitude and slope (Figures 3C,D). At −35 ms interval, the EPSP amplitude increased from 5.22 ± 1.13 to 9.44 ± 1.13 mV after pairing (an increase of 80%) and the slope from 1.57 ± 0.31 mV/ms at baseline half-width, rat: 0.93 ± 0.03 ms and human 1.1 ± 0.2 ms, which may indicate possible differences in voltage-gated sodium channel dynamics.
In rodent hippocampus, glutamatergic inputs to pyramidal neurons potentiate when the temporal order of action potential firing is such that spiking of presynaptic fibers is followed by postsynaptic firing (Bi and Poo, 1998;Meredith et al., 2003). The temporal window in which the postsynaptic neuron must spike is ∼20 ms after the presynaptic stimulus for spike-timing-dependent potentiation (tLTP) to occur (Bi and Poo, 1998). To test whether human hippocampal synapses increase strength in response to spike-timing within this positive timing window of 20 ms, we paired presynaptic stimulation with postsynaptic firing with a positive interval of 10 ms (Figure 2). Since it was reported that in rodent hippocampus the effectiveness of tLTP induction diminishes with age (Meredith et al., 2003), we induced two postsynaptic action potentials (10-ms interval) with each presynaptic stimulus to optimize tLTP induction. After recording 4-5 min of baseline EPSPs, repeated pairing of EPSPs with postsynaptic action potentials (Figure 2A; 40-50 times at 0.14 Hz) resulted in a lasting increase of both EPSP amplitude (from 5.47 ± 1.17 mV at baseline to 9.39 ± 1.90 mV after pairing, an increase of 72%) and slope (from 2.21 ± 0.39 mV/ms at baseline to 3.25 ± 0.79 mV/ms after pairing, an increase of 65%, t-test to 2.85 ± 0.34 mV/ms after pairing (an increase of 81%, t-test p < 0.0001 Figures 3C,D). On average, in pyramidal neurons in which post-before presynaptic timing was applied, synapse strength increased by 58% (n = 4). In non-pyramidal cells, a negative time interval of −10 ms induced a 116% change in EPSP slope and a negative time interval of −20 ms induced a 10% increase in slope (Figure 6). These results indicate that in contrast to rodent hippocampal synapses, adult human hippocampal synapses show tLTP at negative timing intervals up to −80 ms. The temporal window for tLTP induction covers intervals between −80 and +10 ms. Action potentials in human pyramidal neurons showed a prominent after-depolarization (ADP; Figures 3A and 4A; Table 1), possibly reflecting dendritic action potential propagation (Larkum et al., 2001). Dendritic action potential propagation is crucial for tLTP induction Couey et al., 2007;Fuenzalida et al., 2010). In the recordings with negative timing interval pairings between −10 and −80 ms, EPSPs coincided with the downward slope of the ADP (Figure 3A), which may have contributed to induction of tLTP. To test whether tLTD would be induced when the EPSP would occur after the membrane potential had returned to baseline, we increased the negative timing interval (Figure 4). When postsynaptic firing was followed by an EPSP after 130 ms, the membrane potential had returned to baseline at the time of presynaptic stimulation ( Figure 4A). At this interval, the EPSP amplitude and slope indeed showed a sustained reduction (Figures 4C,D). Both in pyramidal neurons and non-pyramidal cells negative timing intervals between −80 and −130 ms induced tLTD. On average, the EPSPs slope was reduced by −55% (n = 3) in pyramidal neurons. In one non-pyramidal cell, negative timing of -110 ms also resulted in a reduction of EPSP slope of −31% (Figure 6, Table 2). These findings show that in contrast to rodent hippocampal synapses, the sign of plasticity does not sharply switch around 0 ms timing in human hippocampal synapses. It switches around −80 ms.
During maturation of the rodent hippocampus, the effectiveness of postsynaptic spikes to induce tLTP diminishes. Pairing EPSPs with single action potentials fails to induce tLTP at ages beyond 20 days, while pairing EPSPs with a pair of action potentials suffices to induce tLTP at these ages (Meredith et al., 2003). To test whether single spikes paired with EPSPs would induce tLTP in adult human synapses, we applied a timed presynaptic stimulus with a single postsynaptic spike (Figure 5). When the action potential was followed by an EPSP after 75 ms in a pyramidal neuron, the EPSP slope showed a small but significant increase from 1.34 ± 0.27 mV/ms to 1.48 ± 0.36 mV/ms after pairing (an increase of 10%, t-test p < 0.05, Figure 5D). In a non-pyramidal neuron, when the EPSP was followed by a single postsynaptic what timing rules determine the sign of plasticity. Our main findings are: (1) Adult human hippocampal synapses can alter synapse strength in response to pairing EPSPs with postsynaptic burst activity and possibly also with single postsynaptic action potentials. (2) In contrast to rodent hippocampal synapses, the sign of plasticity does not sharply switch around 0 ms timing. Instead, both positive timing intervals, in which presynaptic firing preceded the postsynaptic action potential up to 20 ms, and negative timing intervals, in which postsynaptic firing preceded presynaptic activity down to −80 ms, induce tLTP. (3) Negative timing intervals between −80 to −130 ms induce tLTD.
In rodent hippocampus, the rules for induction of STDP change over development (Meredith et al., 2003). A developmental shift occurs in the effectiveness of a single-spike pairing protocol at inducing tLTP. In young hippocampus, repeated pairing of presynaptic activity with a single postsynaptic action potential suffices to induce tLTP, whereas in older rodents a postsynaptic burst of action potentials is necessary to induce synaptic strengthening (Meredith et al., 2003). Blocking GABAergic inhibition with bicuculline rescued the effectiveness of single-spike pairing in inducing tLTP in older animals, which suggests that rules for STDP are affected by a developmental maturation of GABAergic inhibition in the rodent hippocampus (Meredith et al., 2003). In pyramidal neurons, somatic action action potential after 5 ms, the slope increased by 69%, t-test p < 0.0001. These data may suggest that excitatory synapses in human hippocampus can change strength in response to single postsynaptic action potential pairing. However, given the low number of observations on this induction protocol a firm conclusion on this issue awaits further testing.
In rodents and other species, the size and shape of the time windows in which positive and negative synaptic weight changes occur vary for different brain regions (Bi and Poo, 2001;Caporale and Dan, 2008). Excitatory synapses in the human hippocampus showed a wide temporal window for STDP (Figure 6). Pairing intervals between −80 and +10 ms induced robust increases in synaptic weight, showing that tLTP was induced at positive and negative timing intervals. Negative timing intervals between −130 and −80 ms induced substantial tLTD. Our findings show that both for excitatory synapses on pyramidal neurons as well as on non-pyramidal neurons, synapse strength could be altered bi-directionally by spike-timing.

dIscussIon
In this study, we directly tested in whole-cell recordings from human hippocampal neurons whether human synapses can alter strength in response to millisecond timing of pre-and postsynaptic activity and The size and shape of the temporal STDP windows in which positive and negative synapse strength changes can vary for different brain regions (Caporale and Dan, 2008). In rodent hippocampus, the window for synaptic modification is restricted to about 40 ms (Bi and Poo, 1998;Debanne et al., 1998;Nishiyama et al., 2000;Wittenberg and Wang, 2006) and a sharp switch of the direction of synaptic change exists at 0-ms timing interval. In adult human hippocampus we did not observe a sharp change of sign of plasticity; positive as well as negative timing intervals induce tLTP. Increases in synapse strength at negative timing intervals have been observed at excitatory synapses onto GABAergic Purkinje-like neurons in electric fish (Bell et al., 1997). Negative timing intervals up to −50 ms resulted in tLTP. Increases in synapse strength in response to negative timing intervals have also been observed in excitatory synapses at distal dendritic locations in neocortex. In neocortex, the shape of the temporal STDP windows depends on dendritic location of synapses (Froemke et al., 2005). In layer 5 pyramidal neurons, proximal and distal synapses exhibit a progressive distance-dependent shift in the timing requirements of the induction of tLTP and tLTD . Distal synapses potentiate when the EPSP arrives after the onset of an AP, in contrast to the timing requirements of proximal synapses at the same dendrites. Most likely during pairing at −10 ms, distal potentials back-propagate deep into the dendritic tree and activate voltage-gated calcium channels in proximal and distal parts of dendrites, inducing substantial amounts of calcium influx in dendrites and dendritic spines that can trigger synaptic strength changes (Yuste and Denk, 1995;Magee and Johnston, 1997;Stuart et al., 1997;Koester and Sakmann, 1998). Dendritic back-propagation of action potentials is under GABAergic inhibitory control (Tsubokawa and Ross, 1996). Coincident activation of GABAergic inputs reduces dendritic action potential amplitude and dendritic calcium signals associated with the action potential, which may explain the loss of effectiveness of singlespike pairing in inducing tLTP in older rodent hippocampus. In adult human hippocampus of 20-to 66-years old, repeated pairing of presynaptic activity with a single postsynaptic action potential did induce tLTP. This may suggest that differences may exist between human and rodent adult hippocampal dendrites either in the effectiveness of action potential propagation or in the effectiveness of postsynaptic calcium to trigger the molecular machinery for synapse strengthening. Despite the absence of epileptiform activity, brain slices used in our study were removed from brains of epileptic patients, where GABA might have a depolarizing effect (Dzhala and Staley, 2003;Rheims et al., 2008). It is valid to consider the influence that GABAergic input, potentially recruited via extracellular stimulation, might have on the time window for STDP induction. important differences also exist. In contrast to the human synapses, for both synapses on Purkinje-like neurons and rat neocortical distal synapses there was a sharp switch of the sign of plasticity at 0-ms timing interval. Uncovering the mechanisms underlying the tLTP window in human synapses requires further experimental testing.
In rodent hippocampus, the capacity for synaptic depression in synapses declines with age (Dudek and Bear, 1993;Bear and Abraham, 1996). In mouse somatosensory cortex, tLTD induced by a negative timing order in which the postsynaptic neuron fires before the EPSP disappears with age (Banerjee et al., 2009). At ages up to postnatal day 25, negative timing intervals of −10 to −15 ms elicited robust tLTD. At ages beyond 25 days, these timing intervals did not change synaptic strength at all. In adult human hippocampus we find that negative timing intervals between −130 and −80 ms elicit robust tLTD. What the mechanisms are can explain the observation that the timing window lies around −100 ms is currently not known. However, at these intervals, the ADP that follows the action potential has subsided back to baseline by the time the EPSP occurs, which may indicate that at this time the EPSP does not coincide with dendritic calcium dynamics. Our findings do show that human synapses can change bi-directionally depending on spike-timing.
EPSPs coincide with dendritic calcium dynamics induced by the backpropagating action potential Cornelisse et al., 2007). To prevent variability in the amount of synaptic gain due to spatial dependence of the stimulation site (Sjöström and Häusser, 2006), distance between the recording and stimulation electrodes was carefully controlled in the present study, as described in the Section "Materials and methods". Here, we observed that in human hippocampal synapses a wide window for tLTP exists. Possibly timing delays due to propagation of action potentials in extensive dendritic trees could play a role in determining the timing windows. Alternatively, species differences in synapse dynamics may explain wide STDP temporal window in human synapses.
Whether similar mechanisms as for rat neocortical distal synapses underlie tLTP induction at negative timing intervals in human hippocampal synapses is not clear. Human hippocampal pyramidal neurons did show a prominent ADP, which may result from dendritic action potential propagation. At negative timing intervals between 0 and −80 ms, the EPSP coincided with the falling flank of the ADP, possibly indicating that the EPSP coincided with calcium dynamics induced by the dendritic action potential. Similar mechanisms as in distal neocortical layer 5 synapses may extend the tLTP window to negative timing intervals in human synapses. However, mechanistic basis for a wide temporal window in human hippocampus is at this point not clear, but dendritic calcium dynamics induced by dendritic action potential backpropagation may be involved. If this is the case, and EPSPs coinciding with dendritic calcium dynamics induced by dendritic firing can induce tLTP, then a less strict interpretation of Hebb's theory still applies to human synapses. A wider temporal window for strengthening of synapses in the human brain may allow for the association of larger variety of events with less emphasis on the temporal order.
In conclusion, we find that the temporal window for STDP in adult human hippocampal synapses differs from rodent hippocampus, but that the core principles of spike-timing-dependent plasticity apply also to human synapses.  (Hebb, 1949)? He proposed that neurons that fire together, also wire together. With the finding that the sign of plasticity depends on the order of presynaptic and postsynaptic firing (Bell et al., 1997;Magee and Johnston, 1997;Markram et al., 1997) it has been emphasized that temporal specificity is a central feature of Hebb's postulate (Bi and Poo, 2001). Indeed, in rodent hippocampus with a temporal window of 20 ms for tLTP (Bi and Poo, 1998;Debanne et al., 1998;Nishiyama et al., 2000;Wittenberg and Wang, 2006), neurons that fire together will only wire together if pre-and postsynaptic spike timing is tightly correlated. In human synapses it seems that the timing that is necessary to wire together is less strict. The neurons we recorded from showed tLTP in a wide temporal window of about 100 ms. This may be an underestimation since we did not test intervals outside −130 and +10 ms in this study. Different timing windows have been found in different brain areas (Caporale and Dan, 2008). Asymmetric anti-Hebbian STDP, or the depression of an EPSP that occurs if a presynaptic spike follows an increase in the probability of a postsynaptic spike during pairing and reverses into potentiation if the pairing order is reversed (Roberts and Bell, 2002;Zilberter et al., 2009), has been observed in a cell-specific manner in the dorsal cochlear nucleus of the rodent brainstem (Tzounopoulos et al., 2004). The