Abstract
Spike conduction velocity characteristically differs between myelinated and unmyelinated axons. Here we test whether spikes of myelinated and unmyelinated paths differ in other respects by measuring rat retinal ganglion cell (RGC) spike duration in the intraretinal, unmyelinated nerve fiber layer and the extraretinal, myelinated optic nerve and optic chiasm. We find that rapid spike firing and illumination broaden spikes in intraretinal axons but not in extraretinal axons. RGC axons thus initiate spikes intraretinally and normalize spike duration extraretinally. Additionally, we analyze spikes that were recorded in a previous study of rhesus macaque retinogeniculate transmission and find that rapid spike firing does not broaden spikes in optic tract. The spike normalization we find reduces the number of spike properties that can change during RGC light responses. However, this is not because identical spikes fire in all axons. Instead, our recordings show that different subtypes of RGC generate axonal spikes of different durations and that the differences resemble spike duration increases that alter neurotransmitter release from other neurons. Moreover, previous studies have shown that RGC spikes of shorter duration can fire at higher maximum frequencies. These properties should facilitate signal transfer by different mechanisms at RGC synapses onto subcortical target neurons.
Significance Statement
Myelination increases spike conduction velocity in axons. To test whether unmyelinated and myelinated segments of axons differ only in propagation speed, we compared spikes in unmyelinated intraretinal and myelinated extraretinal portions of retinal ganglion cell (RGC) axons. Short interspike intervals and illumination increased spike duration intraretinally but not extraretinally. Consequently, light can change only whether and how often extraretinal spikes fire. Notably, this does not eliminate axonal spike duration differences between subtypes of RGC. This would allow previously described properties of narrower and broader spikes to facilitate signal transfer by different mechanisms at RGC synapses onto subcortical target neurons.
Introduction
Axons are generally described as either myelinated or unmyelinated. Retinal ganglion cells (RGCs) differ markedly from this dichotomy in that their axons are formed by long unmyelinated proximal portions and, in nearly all species examined to date, myelinated distal portions (Bruesch and Arey, 1942; Hildebrand et al., 1985; Stone et al., 1995). The proximal portions are intraretinal, arise from RGC somata or primary dendrites, join fascicles in the innermost sublayer of the retina (the nerve fiber layer), converge at the optic disk, and range in length up to 3.5 mm in rats and ∼20 mm in humans (Ramón y Cajal, 1972; Fukuda, 1977; Curcio and Allen, 1990). The distal extraretinal portions exit the eye, project to numerous subcortical brain areas via the optic nerve, optic chiasm, and optic tract (Berson, 2008; Morin and Studholme, 2014), and are ∼20 mm long in rats and 100 mm long in humans (Sumitomo et al., 1969; Matheson, 1970; Lang and Reiter, 1985; Bernstein et al., 2016).
Consistent with the ability of myelinated axons to conduct spikes faster than unmyelinated axons, the conduction velocities of extraretinal RGC axons are significantly faster than those of intraretinal RGC axons (Granit, 1955; Stone and Freeman, 1971; Stanford, 1987). However, it is not known whether spikes of intraretinal and extraretinal axon segments differ only in conduction velocity. This is of interest because intraretinally recorded RGC spikes change in shape during repeated spike firing (O’Brien et al., 2002; Weick and Demb, 2011) and changes in spike shape alter synaptic efficacy in various preparations (Kusano et al., 1967; Sabatini and Regehr, 1997; Geiger and Jonas, 2000). Moreover, RGCs spike over a broad range of frequencies during changes in illumination (Sakmann and Creutzfeldt, 1969), and increases in instantaneous spike frequency increase the efficacy of signal transfer from at least some RGCs to thalamic neurons (Singer et al., 1972; Mastronarde, 1987; Usrey et al., 1998; Sincich et al., 2007; Weyand, 2007).
Here we examine spikes recorded from the nerve fiber layer, optic nerve, optic chiasm, and optic tract, and test whether the spikes of intraretinal and extraretinal RGC axons differ in three respects known from previous somatic and optic nerve recordings (Hartline, 1938; Cleland and Levick, 1974; Fukuda and Stone, 1974; Hidaka and Ishida, 1998; O’Brien et al., 2002; Kim and Rieke, 2003; Weick and Demb, 2011). Specifically, we test whether spike duration varies with interspike interval, whether the onset of light and/or dark alters spike duration, and whether the axons of subtypes of RGCs generate spikes of different duration. We find that short interspike intervals and the onset of illumination increase spike duration in intraretinal axons but not in extraretinal axons. These results identify a previously unrecognized function of the extraretinal portion of RGC axons, namely to limit the effect of light to altering whether and how often spikes fire extraretinally. In addition, we find that differences among cell types in spike duration of intraretinal axons persist in extraretinal axons. These differences show that the loss of activity- and light-induced spike broadening in extraretinal axons does not result from the generation of identical spikes in the myelinated portion of all RGC axons. This identifies a second function of the extraretinal portion of RGC axons, namely to help establish spike duration phenotypes of different cell subtypes. Notably, the spike duration differences we find between the extraretinal axons of different RGCs resemble spike duration changes that alter the amount of neurotransmitter released per spike in other systems (Sabatini and Regehr, 1997; Geiger and Jonas, 2000).
Materials and Methods
The experiments reported here were performed on ex vivo retinae, in vivo optic nerve, in vivo optic chiasm, and in vivo optic tract.
The ex vivo retina recordings were performed at the University of California, Davis. All animal procedures were performed using methods approved by the University of California, Davis, Animal Use and Care Administrative Advisory Committee in accordance with National Institutes of Health guidelines. Long–Evans rats of both sexes [n = 45; postnatal day 60 (P60) to P120; weight, 150–250 g; Envigo Bioproducts] were killed by a lethal dose of sodium pentobarbital (75 mg/kg, i.p.; Western Medical Supply). After enucleating and hemisecting the eyes, the retinae were isolated and placed in recording chambers (Partida et al., 2018). We found no consistent difference in RGC spike properties between male and female rats, and therefore pooled the results we obtained from them.
The optic nerve and optic chiasm spikes were recorded in vivo at the University of South Florida, using methods approved by the Institutional Animal Care and Use Committee of the University of South Florida in accordance with National Institutes of Health guidelines. Brown Norway rats (n = 25; male; weight, 300–400 g; Envigo Bioproducts) were anesthetized by intraperitoneal injection of ketamine hydrochloride (75 mg/kg) and xylazine (7.5 mg/kg), rested on a heating pad (37°C), and maintained under anesthesia for the remainder of the experiment via intravenous infusion of ketamine (30 mg/kg/h), xylazine (1.5 mg/kg/h), and 0.9% physiological saline. Eye movements were blocked by adding gallamine triethiodide (40 mg/kg/h) to the infusate. A tracheotomy was performed to provide mechanical ventilation during paralysis. Body temperature and heart rate were monitored with a rectal thermometer and EKG electrodes, respectively, and data were collected while vital signs remained at normal levels. Animals were killed with Euthasol (50 mg/kg; Virbac AH) at the end of experiments.
The optic tract spikes were recorded in vivo at the University of California, San Francisco, using methods approved by the University of California, San Francisco, Institutional Animal Care and Use Committee in accordance with National Institutes of Health guidelines. These spikes were recorded in rhesus monkeys (n = 2; unknown sex) during a previously published study of retinogeniculate signal transmission (Sincich et al., 2007). None of the observations presented here duplicate previously published analyses, figures, or graphs.
Intracellular recording
Spikes that were recorded intracellularly from dissociated RGC somata during a previously published study (Hayashida et al., 2009) are shown to help explain how spike duration was measured. As described in detail in that study, these spikes were elicited and measured in perforated-patch whole-cell current-clamp mode with a discontinuous single-electrode voltage-clamp amplifier (SEC-05LX, npi electronic). The conclusions of the present study are drawn entirely from the extracellularly recorded spikes described below, and none of the intracellularly recorded membrane voltage (Vm) traces and measurements shown here duplicate any previously published analyses, figures, or graphs.
Multielectrode array recording and spike identification
Spikes were recorded from isolated retinae by use of multielectrode arrays (MEAs) and classified as somatic or axonal based on the shape of their mean waveforms (Terzuolo and Araki, 1961; Meister et al., 1994). Each retina was placed ganglion cell layer down on an MEA (catalog #60MEA200/30iR-ITO, MultiChannel Systems), held in place with a slice harp (catalog #SHD-41/10, Warner Instruments), and superfused at 35–37°C with buffered Ames medium (A1372-25, United States Biologicals). The individual electrodes were 30 μm in diameter and arranged in an 8 × 8 rectilinear grid with 200 μm interelectrode spacing. Analog data were acquired simultaneously from each of the recording electrodes for 60–180 min, bandpass filtered (1 Hz to 3 kHz, 60 dB/decade), and sampled at a rate of 25 kHz/channel. Data that will be termed “spontaneous spikes” were collected in the absence of light stimuli, under steady ambient illumination of 150 lux. “Electrically stimulated spikes” were activated by biphasic current pulses delivered via one of the array electrodes (Sekirnjak et al., 2006). These spikes propagated as rapidly as in other RGCs (Zeck et al., 2011; Greschner et al., 2014; Li et al., 2015). Moreover, we excluded spikes firing at electrodes surrounding single stimulating electrodes on multiple sides (Greschner et al., 2014). The electrically stimulated spikes reported here are likely to have been recorded distal to the axon initial segments (AISs) of RGC axons because the stimulating and recording electrodes were separated by at least 400 μm (i.e., by distances greater than those between RGC somata and AISs; Boiko et al., 2003). Because we could record these spikes in a low Ca2+/high Mg2+ superfusate (see below) that blocked light-evoked responses, these spikes are likely to have been elicited by direct stimulation rather than by presynaptic inputs (Masland and Ames, 1976). “Light responses” were elicited by a ViSaGe visual stimulus generator (Cambridge Research Systems). Stimuli were presented using an M115HD projector (Dell). The following two types of visual stimuli were used to characterize responses: a full-field, spatially uniform stimulus that stepped between two radiance levels and a pseudorandom binary white-noise stimulus (Reid et al., 1997) that consisted of a 16 × 16 grid of squares that were white or black one-half of the time, as determined by an m-sequence of length 215–1. The spectral radiance was maximal at 450, 515, and 630 nm; measured 170, 60, and 129 μW/sr/m2/nm, respectively, during full-field light onset; and was <0.5 μW/sr/m2 between 380 and 760 nm during full-field light offset (model 670 SpectraScan, Photograph Research). RGC light responses were identified as On or Off by reverse correlation analysis of the response to m-sequence stimulation (Reid et al., 1997), and as transient or sustained based on whether the response diminished or not during steps of full-field illumination (Cleland et al., 1971). In experiments that tested whether spikes could be activated electrically after blocking synaptically mediated light responses, the superfusate was switched from Ames medium to Ca2+- and Mg2+-free Ames medium (A1372-25A, United States Biologicals) that was supplemented with 0.1 mm Ca2+ and 3.4 mm Mg2+ (or, in a few experiments, 0.5 mm Ca2+ and 3.0 mm Mg2+).
Offline Sorter (Plexon; RRID:SCR_000012) was used to sort the MEA-recorded waveforms based on their first three principal components. The waveforms were evaluated without post hoc filtering. Clusters were identified by using the expectation-maximization cluster algorithm in Offline Sorter and then manually edited for clustering errors. None of the spikes attributed to single spiking units in this study occurred simultaneously with other spikes (i.e., none showed amplitude multiplication attributable to coincident firing by more than one unit) or during a time period assumed to be the absolute refractory period (1 ms).
In vivo optic chiasm, optic nerve, and optic tract recording
Optic chiasm and optic nerve spikes were recorded in vivo by protocols described previously (Freeman et al., 2008) and modified as noted here. The head of each anesthetized animal was supported in a stereotaxic apparatus, and a craniotomy was performed at bregma. A tungsten-in-glass microelectrode was inserted in a metal guide needle, which doubled as the reference electrode, and the recording and reference electrode leads were connected to a high-input impedance amplifier (Xcell-3x4, FHC). For optic chiasm recordings, the microelectrode was blindly lowered into the brain at or near bregma until activity was detected in response to visual stimuli (Freeman et al., 2008). For optic nerve recordings, brain tissue was removed and the microelectrode was lowered under visual guidance into the optic nerve near the chiasm (Partida et al., 2018). Signals were amplified differentially (10,000×), bandpass filtered (second-order Butterworth filter, 10 Hz to 1 kHz), and digitally sampled at 100 kHz. (High-pass filtering at 10 Hz did not alter spike duration and was used to reduce slow drift in our recording baselines.) Light responses of individual RGCs having well isolated spikes were classified as On, Off, sustained, and transient using a battery of tests, as described previously (Freeman et al., 2008; Heine and Passaglia, 2011). Data were then collected for constant full-field illumination (30 cd/m2) or, if maintained discharge was low, for alternating step increases to 60 cd/m2 and decreases to 0 cd/m2 with a repeat period of 1 s. Stimuli were displayed on a video monitor measuring 30.2 × 40.4 cm (Trinitron GDM-C520, Sony), running at 100 Hz, and viewed at a distance of 16.5 cm from the cornea. The spectral radiance was maximal at ∼455, 520, and 630 nm, and measured at 767, 955, and 1177 μW/sr/m2/nm, respectively, at 60 cd/m2 (model STS-VIS, Ocean Optics).
The optic tract spikes included here were recorded in vivo by protocols described in detail previously (Sincich et al., 2007). The spikes were recorded extracellularly by use of Epoxylite-coated tungsten microelectrodes (5–8 MΩ; Frederick Haer). The potentials were analog filtered between 0.3 and 3 kHz, differentially amplified 1000× (model 1800, A-M Systems), and acquired digitally at 25 kHz (Power1401, Cambridge Electronic Design). The shapes of individual optic tract spikes were similar to those in previous studies (McIlwain, 1964). Moreover, previous studies have shown that spike durations as brief as those we show for optic tract are not significantly distorted by the high-pass filtering used here (Vigneswaran et al., 2011).
Somatic and axonal spike durations
We examine here the duration of extracellularly recorded spikes. To characterize the spike duration of single cells and to compare these with spikes of other cells, spike properties were typically calculated from mean waveforms. In the absence of a unique definition, and because the numbers of spikes we were able to record varied from cell to cell, we calculated spike duration in several ways.
For illustrative purposes, we first measured the duration of intracellularly recorded somatic spikes as the time elapsing between two “corner” points (Fig. 1A, c1, c2). To locate these points, we calculated the second-order time derivative (d2Vm/dt2) of the mean of the intracellular spike trace (to track how rapidly the membrane voltage changed) and then set criterion levels along the resulting waveform by comparison with studies of voltage-gated Ca2+ current flow during action potentials (Sabatini and Regehr, 1999). Because of noise in the d2Vm/dt2, the rising phase (baseline to first peak) and the falling phase (second peak to baseline) were each separately fit with a Gaussian function centered on the peak. Only the left half or the right half of the Gaussian, respectively, was used in this fit. c1 was defined as the first point of the rising phase (left half) Gaussian fit to exceed 5% of the height of the peak of the fit. c2 was defined as the last point of the falling phase (right half) Gaussian fit to exceed 50% of the peak of the fit (Fig. 1A). Similar values were obtained for individual cells by fitting the rising and falling phases of the average of all spikes recorded and by averaging the durations based on fits to the rising and falling phases of each spike recorded.
Because the c1–c2 time is not available from extracellular recordings, and because the extracellular voltage is approximately the first-order time derivative (dV/dt) of the intracellular voltage, we calculated spike widths using four other time points in the dVm/dt of the mean intracellular spike waveforms. In temporal order, we refer to these points as the inflection (inf), peak 1 (p1), zero crossing (zc), and peak 2 (p2). The timing of p1 is defined as the maximum of the dVm/dt waveform. The timing of inf is defined as the maximum of the second derivative of the dVm/dt waveform in the 1 ms preceding p1. The timing of zc is an interpolated point defined as the first zero crossing following p1 (as dVm/dt returned from p1 and crossed the baseline to begin the downward deflection). The timing of p2 is defined as the minimum point of the dVm/dt waveform within the 2 ms following p1. Figure 1B shows these points calculated from the dVm/dt of a mean intracellular waveform (and then projected onto that mean waveform). Plots of the times that elapse between these points (Fig. 1C) show that the rising phase of each spike (inf–zc), the width between the moments at which the rates of depolarization and repolarization are maximum (p1–p2), and the longest width we extracted from the dVm/dt traces (inf–p2) all increase with the duration defined as c1–c2. The r2 of these lines is 0.913 for inf–p2, 0.903 for p1–p2, and 0.840 for inf–zc. This suggests that the time courses of the extracellular waveform measured by the time points defined here are linearly related to the time course of the intracellularly recorded spikes.
Extracellularly recorded spike durations were calculated from a mean waveform obtained by averaging all of the spikes collected from a given cell under a condition of interest and then finding time points comparable to those used to measure intracellularly recorded spike widths (inf, p1, zc, and p2; Fig. 1D–F). Inf was the time of the minimum of the second derivative of the spike waveform in the 1 ms preceding p1; p1 and p2 were the minimum and maximum voltages reached during the spike waveform (i.e., the voltages at the peak of the inward and outward deflections); and zc was the first zero crossing following p1. When constructing histograms of spike widths, if the spike numbers and frequencies were not sufficient to create a reliable mean waveform for a given time bin, we calculated the spike duration from a fitted version of each individual spike. In these cases, each spike was modeled by a sum-of-three-Gaussians function, with each Gaussian corresponding to a different region of the spike. These Gaussians were of the following form: where a is the y-offset, b is the height, c is the x-location of the peak, and d is a sharpness factor. For extracellularly recorded axonal spikes, the first Gaussian was a positive-going Gaussian fitting the peak at ∼0.4 ms, the second Gaussian was a negative-going Gaussian fitting the peak at ∼0.75 ms, and the third Gaussian was a positive-going Gaussian fitting the peak at ∼1.4 ms. The fit was accepted if there was a zero crossing in the fit between the second and third Gaussians, and the value (mean squared error of the fit)/(absolute value of the height of the second peak raised to the 1.25 power) was <0.5 (Fig. 1G–I). The exponent (1.25) in the denominator was hand chosen to make the fit criterion approximately equal by eye across the range of observed spike amplitudes. Without this exponent, attempts at choosing a criterion were either too permissive for low-amplitude spikes or too restrictive for high-amplitude spikes. By these criteria, ∼70% of the fits were accepted for duration calculations. Rejected fits (e.g., Fig. 1G) differed from accepted fits (Fig. 1H,I) most notably in the shape and amplitude of the initial outward and large inward peaks of the recorded waveforms. All fitting was done with the “fmincon” function in MATLAB (R2014A, MathWorks; RRID:SCR_001622), which allows constraints to be placed on individual parameters of the fit.
Experimental design and statistical analysis
We compared the durations of spikes generated by single units and of spikes generated by different units (see below). One comparison contrasted the duration of spikes as a function of the time between each spike and the spike that preceded it (i.e., the interspike interval or “quiet period”). We will refer to the first spike after a long interspike interval as a cardinal (“C”) spike and to spikes following shorter interspike intervals as follower (“F”) spikes (Partida et al., 2018). To determine the minimum interspike interval that we used to define C spikes in this study (50 ms), we compared C spikes that were preceded by quiet periods of different lengths. For each cell, we found all C spikes that fell within bins of quiet periods ranging from 50 to 150 ms, 150 to 300 ms, 300 to 500 ms, 500 to 1000 ms, and 1000 to 2000 ms. Cells with at least five spikes in each of these five bins were included for further analysis (nerve fiber layer, n = 56; optic chiasm, n = 26). For each of these cells, a mean spike was created from the individual spikes in each bin, and the duration of each mean spike was calculated as described above. We then used interspike intervals that allowed C spike durations to recover from the effects we are reporting here. Having defined C spikes, we compared C and F spikes at interspike intervals found during light responses (Heine and Passaglia, 2011; Baden et al., 2016).
To compare the duration of C and F spikes, we plotted the duration of F spikes as a function of C spike durations. We also plotted histograms of the percentage increase in F1 spike duration (i.e., the duration of the first F spike), calculated by normalizing the Euclidean distance of the C/F spike width pairs (inf–p2) from the unity line by the mean of the C spike durations for individual somata and axons. These normalized C–F changes in width were histogrammed from −40% to 40% in bins of 1% (for optic chiasm data) or 5% (all other data).
We compared spike durations between RGCs that responded differently to illumination (see Results; Cleland and Levick, 1974; Fukuda and Stone, 1974). Because we found cell-to-cell differences in the number of spikes generated, we digitally averaged the cardinal spikes of a selected cell and calculated the duration of the averaged waveform. We then repeated this calculation for each of the other cells of the same type, yielding one duration per spiking unit. The mean of the durations of these units, and the digital average of the averaged waveforms, were taken as representative of the cell type. This enabled us to compare the spike durations of different cell types without biasing the mean waveform by large numbers of spikes from a few cells with uncharacteristically narrow or broad spikes.
Likewise, each cell generated different numbers of spikes in the dark and light, and these durations could vary from cell to cell. To test whether spike duration differed in the dark and light, we recorded spikes in response to two changes in illumination (dark–light, light–dark) and compared the duration of equal numbers of spikes in dark versus light prior to these transitions and also in dark versus light after these transitions. To do so, we (1) calculated the mean of the durations of all spikes from a given unit during a given time bin in the light (to compare spikes before the transitions) and in the dark (to compare spikes after the transitions); (2) identified the condition that produced fewer spikes (e.g., dark for an On RGC); (3) equalized the number of spikes to be compared by randomly selecting the same number of spikes under the other condition (e.g., light for that On cell); (4) normalized the duration of each of the raw spikes from step 2 and each of the selected spikes from step 3 by the mean from step 1 for the time bin under consideration; (5) repeated steps 1–4 for all units recorded from; (6) pooled the normalized durations from these units; (7) calculated the mean of these normalized durations; and (8) repeated steps 1–7 for each time bin. We then plotted the mean of the normalized durations as a function of time before and after the dark–light and light–dark transitions.
Statistical analyses were performed using MATLAB (R2014A). Unless stated otherwise, measurements are presented as the mean ± SEM and sample size. For sufficiently large sample sizes, we used two-tailed Student’s t tests to assess whether differences are statistically significant and report the mean, SEM, value range, sample size, and p value. For the data presented in scatter plots, we used Wilcoxon signed-rank tests to assess whether differences are statistically significant, plot the mean and SEM, and state the sample size and p values. The sample size (i.e., the number of spiking somata or axons) are stated in the text and figure legends.
Results
This study compares the duration of spikes that are recorded from rat nerve fiber layer, optic nerve, and optic chiasm. To test whether and how these spikes differ, we record spikes (1) as they fire spontaneously, (2) during responses to electrical stimulation, and (3) before and after the onset and cessation of illumination. We then measure spike duration as the time that elapses between pairs of time points (defined in the Materials and Methods as p1, p2, inf, and zc) during each spike waveform (Fig. 1) and compare spikes that fire at different frequencies and in functionally different cell types. When describing multiple spikes, we refer to the first spike after a long interspike interval as a cardinal (C) spike and to spikes after shorter interspike intervals thereafter as follower (F) spikes (Partida et al., 2018). Additionally, to test whether spikes broaden in other species, we compare C and F spikes that were recorded in rhesus macaque optic tract (Sincich et al., 2007). The comparisons below are based on C spikes that were each preceded by a ≥50 ms spike-free period, and on the first F (F1) spike that fired within 30 ms after each C spike, unless stated otherwise.
Activity alters intraretinal spike duration but not extraretinal spike duration
Voltage-gated Na+ current declines in amplitude as RGCs are depolarized at frequencies >10 Hz (Hidaka and Ishida, 1998; Kim and Rieke, 2003; Hayashida and Ishida, 2004). The types of K+ current that gate at the voltages traversed by RGC spikes also inactivate during sustained (Lipton and Tauck, 1987) and repeated (Geiger and Jonas, 2000; Faber and Sah, 2003) depolarizations. These activity-induced losses of current would be expected to reduce spike maximum rate-of-rise and amplitude (Fleidervish et al., 1996) and increase spike duration (Jackson et al., 1991; Giese et al., 1998). We therefore tested whether a spike preceded recently by another spike (short interspike interval) differed in shape from spikes preceded by long spike-free periods (long interspike intervals).
Because our MEA recordings captured somatic and axonal spikes, we first checked whether repeated spike firing altered somatic spikes as reported previously (O’Brien et al., 2002; Weick and Demb, 2011). Repetitive spiking broadened spikes we recorded from somata (n = 166) in situ on multielectrode arrays (Fig. 2A–D). Plotting F1 spike duration versus C spike duration shows that most F1 spikes were broader than the C spikes that preceded them. Replotting these data as a histogram of the activity-induced changes in F1 spike durations shows that a large majority of the F1 spike durations are longer than the C spike durations (Fig. 2C, inset). The differences in duration were statistically significant regardless of the time points they were measured from (Fig. 2D; p = 5.68e-4, 0.0101, and 8.34e-4 for inf–zc, p1–p2, and inf–p2, respectively). These effects were rapid in onset and reversible, in that they were seen at interspike intervals as short as 3 ms and because C spikes presented a consistent (control) shape after spike-free periods of ≥30 ms (see below).
Do short interspike intervals also broaden spikes in axons? We examined this first in the nerve fiber layer, and thereafter in the optic nerve, optic chiasm, and optic tract. Averaged across all of the intraretinal axons we recorded from (n = 201), short interspike intervals increased the width of spontaneously firing F spikes (Fig. 2E–H). The increase was statistically significant (Fig. 2H) regardless of whether spike width was measured as inf–zc (p = 4.14e-7), p1–p2 (p = 4.61e-5), or inf–p2 (p = 4.72e-6). Aligning the mean waveforms of C and F1 spikes at their p1 peaks showed that activity increased the time elapsing between p1 and p2 (Fig. 2F). The F1 spike width was broader than the C spike width in most of the nerve fiber layer axons we recorded from (Fig. 2G, inset therein). Similar effects were produced by electrically activating C and F1 spikes with an MEA electrode and varying the interval between the pulses that elicited these spikes (Fig. 3). The p1–p2 width increased (Fig. 3A) for nearly all of the nerve fiber layer axons we recorded from (e.g., Fig. 3C, 3 ms panel) and the increases in spike width were statistically significant at time intervals between 3 and 30 ms [Fig. 3B; intervals for p values are in parentheses: 7.24e-5 (3 ms), 2.28e-5 (5 ms), 8.61e-5 (8 ms), 8.40e-5 (10 ms), 2.17e-4 (15 ms), 0.0015 (30 ms), 0.3989 (50 ms), 0.5017 (100 ms), 0.0268 (500 ms)]. Figures 2 and 3 thus show that short interspike intervals broaden the spikes of RGC somata and intraretinal axons.
Do short interspike intervals also broaden spikes in extraretinal axons? We examined this by recording from optic chiasm (and, to the extent we found possible, optic nerve) and by analyzing spikes recorded from optic tract. We detected no effect of activity on optic chiasm F spike widths (Fig. 2I–L; p = 0.8851, 0.4463, and 0.4405 for inf–zc, p1–p2, and inf–p2, respectively). Scatter plots showed no cases where F1 and corresponding C spike widths differed (Fig. 2K; n = 38). Although we were unable to record for periods long enough to characterize optic nerve spikes in detail, we obtained comparable results from a limited number of optic nerve axons (n = 3; not illustrated). C and F1 axonal spikes from two different optic nerves showed no difference in duration at interspike intervals <50 ms. C and F1 spikes from one other optic nerve increased in duration at interspike intervals <30 ms. These results show that repetitive spike firing broadens spikes in intraretinal axons but not in extraretinal fibers (with one exception in our datasets).
Given that C spikes broadened intraretinal F1 spikes, we tested two potential reasons why C spikes might not have broadened extraretinal F1 spikes. One possibility is that the spike-free times used to define C spikes were long enough for extraretinal C spikes to return to control values, but were too short for intraretinal C spikes to return to control values. Our data indicate otherwise, as we found no difference between C spike inf–p2 widths measured after quiet periods between 50 and 2000 ms, for axonal spikes recorded either intraretinally or extraretinally (Fig. 4A). As determined by one-way ANOVA, there were no significant differences between these C spike widths in the nerve fiber layer (F(4,275) = 0.38, p = 0.8195) or in the optic chiasm (F(4,125) = 0.02, p = 0.9995). A second possibility is that the extraretinal C and F1 spikes in our recordings were separated by longer times than the intraretinal spikes. Again, our data indicate otherwise, as plots of inf–p2 showed that the broadening of nerve fiber layer F spikes increased as the interspike interval decreased from 35 to ≤5 ms and that optic chiasm F spikes did not broaden over the same range of interspike intervals (Fig. 4B). These results thus show that extraretinal spike widths are not altered by interspike intervals that broaden spikes in intraretinal axons.
The results in Figures 2–4 raise two questions. First, does spike broadening occur distal to the optic chiasm? Second, do short interspike intervals broaden spikes in RGC axons of other species? We addressed both questions by analyzing spikes recorded from rhesus macaque optic tract fibers (n = 23; Sincich et al., 2007). Figure 5 plots averages of the C spike waveforms and of the F1 spike waveforms recorded from an individual fiber (Fig. 5A), C and F1 spike widths averaged across all of the fibers recorded from (Fig. 5B), F1 spike width as a function of C spike width for each of the fibers recorded from (Fig. 5C), the Euclidean distance of the paired C/F values from the unity line (Fig. 5C, inset), and the difference between C and F1 spike width as a function of interspike interval (Fig. 5D). These panels in Figure 5 show that C spikes did not induce F spike broadening in optic tract (p = 0.1177, 0.6647, and 0.4264 for inf–zc, p1–p2, and inf–p2, respectively) and that, in this respect, the results we obtained from macaque optic tract closely resembled those we obtained in rat optic chiasm.
Light alters intraretinal spike duration but not extraretinal spike duration
The results in Figure 4A indicate that, given sufficiently long quiet periods, axonal C spike widths are constant. However, these spikes were recorded under constant illumination, and RGCs normally respond to changes in illumination. Because alternating steps of light and dark modulate spike firing frequency in all functional types of RGCs identified to date (Baden et al., 2016), we next tested whether the onset of light and/or dark alter axonal C spike shape. To assess this intraretinally, spikes were recorded on MEAs. Light onset increased spike firing frequency (i.e., elicited On responses; Hartline, 1938) at some electrodes (Fig. 6A). None of these light responses were sluggish in onset and offset as in intrinsically photosensitive RGCs (Berson et al., 2002). When we sorted and averaged the spike waveforms recorded during these responses (i.e., during 1.0 s before and after light offset, and during 1.0 s before and after light onset), the resulting waveforms showed that (1) after light onset, spike width was greater during illumination than during each preceding dark period (Fig. 6C); and (2) spike width did not noticeably change after light offset (i.e., it did not differ between light periods and subsequent dark periods; Fig. 6B).
Because light increased both spike frequency and spike width in Figure 6A,C, and D, we next asked whether illumination broadens spikes simply by increasing spike frequency. We found that the effects of light and short interspike intervals on spike width differ in multiple ways. One is that while short spike intervals increase F spike width (Figs. 2–4), light onset increases the width of all spikes, including C spikes (Fig. 6C,D). A second difference is that light can increase spike width even if it decreases spike frequency. For example, although spike frequency is less during light than during dark during the Off light response in Figure 7A, the spike width is greater during light (Fig. 7C). Notably, spikes from three axons were extracted by principal component analysis (PCA) in this recording and light onset increased the mean spike width in all three (e.g., Fig. 7C,E). Moreover, in RGCs that generated either On or Off light responses, short interspike intervals increased F spike width both in the light and the dark. These results (Figs. 6, 7) show that light and short interspike intervals can modulate intraretinal RGC spike duration independently and suggest, in turn, differences in underlying mechanisms.
To compare the spike widths of different axons regardless of differences in the numbers of spikes they generated, we averaged all of the spikes recorded from each axon in the light and in the dark. Across all of the axons we recorded from, the mean of these averages during the first 1.0 s after light onset was ∼15% greater than during the first 1.0 s after light offset, and the differences of these means were statistically significant (Fig. 6D, asterisks; p = 0.0025, 2.8e-08, and 2.1e-07 for inf–zc, p1–p2, and inf–p2, respectively). Viewed axon by axon, the C spike widths in the light were greater than those in the dark for most of the axons we recorded from and the largest changes in C spike width we found during illumination were increases (Fig. 6E, and inset). Plotting the averages of the C spike widths as a function of time before and after the change in illumination shows that C spike width increased markedly at light onset, remained high during the first second thereafter, and then declined during the subsequent 1–1.5 s to values that did not differ by statistically significant amounts from spike durations in the dark (Fig. 6F, blue lines). This plot also shows that C spike width did not noticeably change after light offset (i.e., it did not differ between the light periods and subsequent dark periods; Fig. 6F, red lines).
Does illumination broaden spikes in all axons? Although light significantly increased the mean spike width, some axons we recorded from showed either no change in width or a decrease in width (Fig. 6E). This might reflect systematic differences between certain subtypes of RGC or variability among axons regardless of type (e.g., in the rate of return to control values). Figure 6G–J summarizes the spike widths we measured in axons classified into four types (On, Off, sustained, and transient) for comparison with previous studies (Fukuda and Stone, 1974; Wong et al., 2012). These plots show that light onset increased spike width by statistically significant amounts in each cell class and also illustrates differences with the effects of short interspike intervals in addition to those noted above. First, the increases during illumination were a few fold larger than the effect of C spikes on F1 spikes (compare Figs. 4B, 6F–J). Second, spike width increased gradually within the first 500 ms after light onset (Fig. 6G–J) but increased maximally within 2–5 ms after C spikes (Fig. 4B). This became apparent when we pooled the duration data in 250 ms bins (Fig. 6G–J) but not in 500 ms bins (Fig. 6F). Altogether, the results in Figures 6 and 7 show multiple differences between the effects of light and interspike interval on RGC spike duration.
Having found that light broadens spikes in intraretinal axons, we next asked whether light broadens spikes in extraretinal axons. Averaged across all of the optic chiasm axons we recorded from (n = 19), none of the mean C spike widths measured in darkness differed by significant amounts from those measured in light (Fig. 8D; p = 0.9536, 0.8184, and 0.9369 for inf–zc, p1–p2, and inf–p2, respectively). None of the mean C spike widths differed between dark and light by >100 μs (Fig. 8E), and the inf–p2 widths in the dark and light showed no significant differences as a function of time before and after the changes in illumination (Fig. 8F). Thus, in marked contrast to the effects we found in nerve fiber layer axons (Figs. 6, 7), we found no differences between C spike durations in the dark and light in the optic chiasm (Fig. 8). Together, Figures 4-8 show that light and short interspike intervals broaden the spikes of intraretinal RGC axons but that neither light nor short interspike intervals broaden spikes of extraretinal RGC axons.
Cell subtype-specific spike duration differences are preserved in intraretinal and extraretinal axons
Spikes of different duration have been recorded intraretinally from some subtypes of RGC. For example, RGCs termed “transient” and “sustained” generate brief and longer bursts of spikes following the onset of illumination (Cleland et al., 1971), and electrical stimulation elicits somatic spikes of shorter duration in transient cells than in sustained cells (Fukuda and Stone, 1974). Likewise, light onset and offset elicit increases in spiking in RGCs termed “On” and “Off,” respectively, and morphologic counterparts of these cells have been found to generate individual somatic spikes of shorter duration in Off cells than in On cells (Wong et al., 2012). We therefore compared axonal C spikes of these types of RGCs, and did so in intraretinal and extraretinal axons.
The C spikes were broader in On nerve fiber layer axons than in Off axons (Fig. 9A,E). The difference between the mean widths of these On and Off units was ∼10% and was statistically significant (p = 0.0359, 0.0124, and 0.0328 for inf–zc, p1–p2, and inf–p2, respectively; n = 85 On units and n = 80 Off units). Likewise, the C spikes of On optic chiasm axons (n = 19) were ∼15% broader than those of Off optic chiasm axons (n = 19; Fig. 9C,G). This difference was statistically significant when calculated from p1–p2 (p = 0.0174). The same trend was seen when the spike widths were calculated as inf–zc or inf–p2. However, the differences did not reach statistical significance (p = 0.1375 and 0.0538, respectively).
The intraretinal axonal C spikes of sustained cells (n = 110) were ∼19% broader than those of transient cells (n = 44; Fig. 9B,F). These differences were statistically significant (p = 0.0046, 7.77e-4, and 5.70e-4 for inf–zc, p1–p2, and inf–p2, respectively). Likewise, the C spikes of sustained optic chiasm axons were broader than those of transient units (Fig. 9D,H). Calculated as p1–p2, the spikes of the sustained units (n = 22) were ∼11% broader than those of the transient units (n = 15), and this difference was statistically significant (p = 0.0412). The spike widths calculated from other time points (inf–zc, inf–p2) were also greater in sustained units than in the transient units, although these differences did not reach statistical significance (p = 0.3391 and 0.1655, respectively).
These results (Fig. 9) show that, although extraretinal axons normalize spike duration within cell types (i.e., as spikes propagate from intraretinal to extraretinal segments), they preserve spike duration differences between cell types (i.e., between sustained and transient cells, and between On and Off cells).
Discussion
Our results identify three properties of spikes in RGC axons. First, short interspike intervals broaden spikes in intraretinal axons but not in extraretinal axons. Second, light onset broadens spikes in intraretinal axons but not in extraretinal axons. Third, axonal spike duration differs between functionally different RGCs both intraretinally and extraretinally. Below, we compare these findings to spikes in other preparations and relate spike durations to signal transfer by RGCs.
Spike broadening and normalization
Spike broadening in the nerve fiber layer and not in the optic chiasm distinguish axons of RGCs and other neurons in two respects. First, unlike axons that can conduct and transmit information with spikes of different shapes (Juusola et al., 2007), individual extraretinal axons eliminate analog differences in spikes and thus reduce the number of spike parameters that visual stimuli can modulate. Second, extraretinal axon segments are as far as several millimeters from the AISs that initiate RGC spikes intraretinally (Wollner and Catterall, 1986; Sekirnjak et al., 2008), and thus are much further from somata than the axonal regions that normalize spike duration in other systems (Williams and Stuart, 1999; Meeks et al., 2005; Khaliq and Raman, 2006; Kole et al., 2007; Foust et al., 2010; Popovic et al., 2011). This separation might be a “best possible” design, in which retinae use unmyelinated RGC axons intraretinally to avoid the opacity of myelin (Straatsma et al., 1979) and then compensate for the slower spike propagation in these fibers by increasing conduction velocity in the optic nerve (Stone and Freeman, 1971) and between the optic nerve and optic tract (Baker and Stryker, 1990).
We did not attempt to identify how short interspike intervals broaden spikes intraretinally, in part because the same intervals do not broaden spikes in optic chiasm or optic tract. We also do not know of any ion channel identifications or current measurements in RGCs that would explain or predict extraretinal spike normalization. However, studies of other neurons suggest possibilities, having found that the kinetics and high densities of Kv1, Kv3, and KCa1.1 channels facilitate Na+ current recovery from inactivation (Kole et al., 2007; Li et al., 2011; Ritzau-Jost et al., 2014) and guard against spike broadening. Increases in Na+ current density, and lowered membrane capacitance in myelinated axons, should also help to curb spike broadening extraretinally.
Generality
The effects we report here of short interspike intervals and light onset on spike width are consistent with the possibility that unmyelinated and myelinated axon segments differ in ways besides conduction velocity. Two questions emerge from these results. First, do spikes broaden in unmyelinated segments of all axons at the interspike intervals examined here (Figs. 3, 4)? This would be of interest to test in several preparations, including axons that are unmyelinated during development and myelinated in adulthood (Foster et al., 1982), adult optic nerves that contain several (and, in a few cases, exclusively) unmyelinated axons (Gaze and Peters, 1961; Holmberg, 1972; Hokoç and Oswaldo-Cruz, 1978), myelinated axons with unmyelinated segments interspersed between their proximal and distal ends (Salami et al., 2003; Tomassy et al., 2014), and axons before and after demyelinating diseases (Ulrich and Groebke-Lorenz, 1983; Utzschneider et al., 1993). Second, do spikes broaden in RGC axon terminals and, if so, does this alter signal transfer? Answering these questions will require methods for testing whether spike duration and amplitude both change (Sabatini and Regehr, 1999) and for assessing the variability of Ca2+ signaling and spike duration at different terminal boutons even in single axon branches (Brenowitz and Regehr, 2007; Rowan et al., 2016), especially given the morphology of RGC axon terminals (Bowling and Michael, 1980; Hammer et al., 2015; Morgan et al., 2016).
A recent study found that the amplitudes of extracellularly recorded spikes are “regularized” as they propagate toward the optic disk in M1-type intrinsically photosensitive RGC axons (Milner and Do, 2017). These results differ from ours in that activity-induced increases in spike width were not noted. Whether spike broadening occurs in M1 RGCs at the high frequencies we examined here remains to be tested.
Spike broadening at light onset
Although light significantly increased the mean spike width in our intraretinal recordings (Fig. 6F), this effect was transient, and some axons showed either no change in width or a decrease in width (Fig. 6E). These properties imply that the light-triggered increases in spike duration were not a nonspecific (e.g., photoelectric) artifact (Moody et al., 1969). Nevertheless, because short interspike intervals increase intraretinal spike duration, light onset increases spike frequency in On cells, and light offset increases spike frequency in Off cells, we were surprised to find that light onset increased spike duration in both On and Off cells (Fig. 6G,H). Because light also broadened spikes regardless of whether the spikes were C or F spikes, and the time course of light-induced broadening was much slower than that of activity-dependent broadening (Figs. 4B, 6F–J), a light-modulated synaptic input to On and Off RGCs seems more likely than short interspike intervals to have caused this broadening. How these inputs operate remains to be identified. However, previous studies have found tyrosine hydroxylase-immunopositive (Kolb et al., 1990) and GABA-immunopositive (Mosinger et al., 1986) neurites in the nerve fiber layer and GABA receptors on RGC axons (Hughes et al., 1989).
Cell-specific spike duration differences between axons
Our results show that axonal spikes of functionally different RGCs differ in duration, and that the cell type-specific differences between nerve fiber layer axons are also found in the optic chiasm (Fig. 9). Our finding that axonal spike widths of sustained RGCs are longer than those of transient RGCs, and that axonal spike widths of On RGCs are longer than those of Off RGCs, agree with differences reported for the somata of these cell types in rats (Wong et al., 2012) and cats (Fukuda and Stone, 1974; O’Brien et al., 2002).
The spike widths of these On and Off cells, and those of the transient and sustained cells, differed by ∼15% (Fig. 9). Spike duration differences of similar magnitude have been found in other neurons, before and after activity-induced broadening, to significantly increase the amplitude of responses to neurotransmitter released by these spikes (Sabatini and Regehr, 1997; Geiger and Jonas, 2000). Thus, broader RGC spikes may transmit signals at higher efficacy than narrower spikes at synapses where single spikes elicit postsynaptic responses (Burke and Sefton, 1966; Cleland and Lee, 1985; Blitz and Regehr, 2003; Gollisch and Meister, 2008).
Do narrower spikes drive other types of synapses better than broader spikes? This would be expected to occur at retinogeniculate synapses that transmit at high efficacy if two or more spikes reach synaptic terminals within short time intervals because RGCs that generate the narrowest spikes can fire repeatedly at the highest rates (O’Brien et al., 2002; Wong et al., 2012; for review, see Bean, 2007) and because synaptic efficacy increases as instantaneous spike frequency increases to >20 Hz (Singer et al., 1972; Mastronarde, 1987; Usrey et al., 1998; Sincich et al., 2007; Weyand, 2007). For example, our finding that Off RGC axons generate narrower spikes than On RGC axons raises the possibility that Off RGCs drive thalamic neurons better than On RGCs. Similarly, our finding that transient RGC axons generate narrower spikes than sustained RGC axons raises the possibility that transient RGCs drive thalamic neurons better than sustained RGCs. The magnitudes of these differences are difficult to predict because Off RGCs have been found to conduct spikes more rapidly than On RGCs (Fukada, 1971; Vaney et al., 1981), transient RGCs have been found to conduct spikes more rapidly than sustained RGCs (Fukada, 1971; Cleland et al., 1971), and spike-induced spike slowing reduces the maximum spike frequency generated by slowly conducting RGC axons more than that of rapidly conducting RGC axons (Partida et al., 2018). However, these trends would increase the difference in signal transfer by Off versus On RGC axons, and by transient versus sustained RGC axons, and in turn contribute to the spread of maximum synaptic efficacies observed to date (Usrey et al., 1998). These differences may be even more subtype specific in the cat, in that Off brisk-transient RGC axons conduct spikes faster than On brisk-transient axons, and Off and On cells of other subtypes (brisk sustained, sluggish) do not differ in conduction velocity (Kirk et al., 1975). Similarly, macaque Off midget RGCs generate narrower spikes than On midget RGCs, and the spike duration of Off and On parasol RGCs do not noticeably differ (Jepson et al., 2013). Notably, these differences parallel the findings that black drives macaque parvocellular V1 cortex better than white, and that there is no preference for black in magnocellular V1 (Xing et al., 2010).
Technical questions
One might ask whether the difference in broadening between intraretinal and extraretinal axons was due to our recording the former ex vivo and the latter in vivo. In one of three optic nerve axons we were able to record from, spikes broadened in response to the same interspike intervals that broadened spikes in the nerve fiber layer (see Results). Also, spikes resisted broadening in a small fraction of the nerve fiber layer axons we recorded from, as in all of the optic chiasm axons we recorded from. Because we found both spike broadening and resistance to spike broadening in the nerve fiber layer and optic nerve, neither seems attributable to our recording conditions. This is consistent with our recent finding that C spikes slow the conduction velocity of F spikes in compound action potential recordings from both ex vivo and in vivo rat optic nerves (Partida et al., 2018). We tried to record from more single optic nerve axons, both in vivo, in isolated portions of optic nerve (Partida et al., 2018), and in an ex vivo retina-attached optic nerve preparation (G. Ogata and A.T. Ishida, unpublished observations) like those developed to study retinal inputs to lateral geniculate nucleus (Mooney et al., 1996) and suprachiasmatic nucleus (Wong et al., 2007). However, none of the electrodes we tested (e.g., glass patch-clamp electrodes, carbon fiber electrodes, various sharp-tipped metal electrodes), including the type we used to record from optic chiasm, allowed us to record from additional optic nerve axons.
One might also ask whether the difference in broadening between intraretinal and extraretinal axons was due to our use of different rat strains for our ex vivo and in vivo recordings. We think this is unlikely because the spikes in these strains are similar in every other property we assessed. Specifically, we found similar C spike durations (Fig. 9) and interspike intervals (Fig. 4), similar cell type-specific differences (Fig. 9), and activity-induced propagation slowing (Partida et al., 2018), in them. We know of no differences in synaptic inputs (Pitcher et al., 2007), neuromodulator responses (Gulledge and Jaffe, 1998), or spike-generating mechanisms that would favor broadening in one strain and lack of broadening in the other. Moreover, similar RGC receptive field sizes (Heine and Passaglia, 2011) and flash electroretinogram kinetics (Polosa et al., 2016) have been found in control Brown Norway and Long–Evans rats.
Lastly, because spikes did not broaden in the optic tract of rhesus macaques (Fig. 5) ventilated with isoflurane and N2O/O2 to maintain anesthesia (Sincich et al., 2007), we would not attribute the lack of spike broadening in optic chiasm to our use of ketamine and xylazine to maintain anesthesia during our in vivo rat recordings. Instead, the results in Figures 2–5 raise the possibility that lack of spike broadening is a general property of the extraretinal portion of mammalian RGC axons.
Acknowledgments
Acknowledgments: We thank D.R. Copenhagen, W.M. Usrey, and K.H. Britten for discussions and comments on the manuscript; D.K. Warland and G.D. Field for discussions; J.S. Werner for radiometric measurements; Y. Hayashida and H. Oi for performing the dissociated soma recordings included here; and L.C. Sincich for access to the optic tract spike recordings described here and permission to present the numbers we extracted from those recordings.
Footnotes
The authors declare no competing financial interests.
This work was supported by National Institutes of Health | National Eye Institute Grants RO1-EY-008120 (to A.T.I.), RO1-EY-027037 (to C.L.P.), and P30-EY-012576 (to J.S. Werner); and by a Plum Foundation Research Award (to A.T.I.).
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.