- 1Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
- 2Werner Reichardt Centre for Integrative Neuroscience (CIN), University of Tübingen, Tübingen, Germany
- 3Key Laboratory of Primate Neurobiology, Institute of Neuroscience, CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai, China
- 4Key Laboratory of Preclinical Study for New Drugs of Gansu Province, and Institute of Physiology, School of Basic Medical Sciences, Lanzhou University, Lanzhou, China
- 5Max-Planck-Institute for Brain Research, Frankfurt am Main, Germany
- 6Department of Computational Neuroethology, Max Planck Institute for Neurobiology of Behavior – Caesar, Bonn, Germany
In the vertebrate retina, several dozens of parallel channels relay information about the visual world to the brain. These channels are represented by the different types of retinal ganglion cells (RGCs), whose responses are rendered selective for distinct sets of visual features by various mechanisms. These mechanisms can be roughly grouped into synaptic interactions and cell-intrinsic mechanisms, with the latter including dendritic morphology as well as ion channel complement and distribution. Here, we investigate how strongly ion channel complement can shape RGC output by comparing two mouse RGC types, the well-described ON alpha cell and a little-studied ON cell that is EGFP-labelled in the Igfbp5 mouse line and displays an unusual selectivity for stimuli with high contrast. Using patch-clamp recordings and computational modelling, we show that a higher activation threshold and a pronounced slow inactivation of the voltage-gated Na+ channels contribute to the distinct contrast tuning and transient responses in ON Igfbp5 RGCs, respectively. In contrast, such a mechanism could not be observed in ON alpha cells. This study provides an example for the powerful role that the last stage of retinal processing can play in shaping RGC responses.
1 Introduction
Retinal ganglion cells (RGCs) and their upstream circuits detect and encode specific visual features and relay this information along parallel pathways to higher visual centres in the brain [reviewed in (Kerschensteiner, 2022)]. Functional, anatomical, and genetic evidence (Baden et al., 2016; Bae et al., 2018; Rheaume et al., 2018; Goetz et al., 2021) support the presence of at least 40 RGC types in the mouse retina.
Retinal ganglion cells receive their excitatory drive mostly from bipolar cells (BCs), which relay the photoreceptor signal to the inner retina, and their inhibitory input from amacrine cells (ACs) [reviewed in (Diamond, 2017)]. As the dendrites of different RGC types arborize at distinct inner plexiform layer (IPL) depths (Helmstaedter et al., 2013; Sümbül et al., 2014; Bae et al., 2018), they pick up inputs from distinct sets of BC and AC types (Field et al., 2010; Helmstaedter et al., 2013). The selective connectivity with presynaptic neurons in the IPL is considered the foundation of the feature-selectivity of RGC pathways [e.g., (Roska and Werblin, 2001; Helmstaedter et al., 2013; Neumann et al., 2016; Yu et al., 2018)].
The response properties of an RGC type are typically determined by a hierarchy of mechanisms. For instance, the temporal response of an RGC to a light-step is shaped by presynaptic circuit components, such as glutamate receptor kinetics along the BC-RGC pathway (Awatramani and Slaughter, 2000; DeVries, 2000; Turner and Rieke, 2016; Yu et al., 2018), as well as by AC input (Nirenberg and Meister, 1997; Asari and Meister, 2012; Nikolaev et al., 2013; Jacoby et al., 2015; Cui et al., 2016; Franke et al., 2017). Moreover, the spatial receptive field (RF) is jointly formed by horizontal cells (Drinnenberg et al., 2018) and ACs (Diamond, 2017) in the outer and inner retina, respectively. In particular, AC circuits are very versatile “function modifiers”: For example, in the highly contrast-sensitive On alpha RGCs (Krieger et al., 2017), different AC circuits converge to provide On and Off inhibition to balance tonic excitatory drive from BCs (Park et al., 2018; Sawant et al., 2021); in On delayed RGCs, they provide a fast, excitatory surround through disinhibition (Mani and Schwartz, 2017).
In addition, RGC types differ in their expression of ion channels (Siegert et al., 2012; Rheaume et al., 2018; Tran et al., 2019) and how they are distributed across the cell [reviewed in (Van Hook et al., 2019)]. This ion channel complement, in combination with the RGC’s specific dendritic geometry, determines how synaptic input is integrated (Schachter et al., 2010; Poleg-Polsky and Diamond, 2011; Ran et al., 2020) and how the resulting signal is translated by the cell’s spike generator into action potentials on the optic nerve (Mobbs et al., 1992; Kim and Rieke, 2001, 2003; Raghuram et al., 2019; Werginz et al., 2020; Wienbar and Schwartz, 2022). Hence, RGCs “themselves” can significantly transform the input they receive from their circuits and thereby shape the retina’s output to the brain [reviewed in (Branco and Häusser, 2010; Stuart and Spruston, 2015; Tran-Van-Minh et al., 2015)].
The contribution of the aforementioned mechanisms to the RGC output varies among RGC types, yielding the diversity of feature representations, such as local edges (Levick, 1967; van Wyk et al., 2006; Zhang et al., 2012), approaching objects (Münch et al., 2009), “uniformity” (Levick, 1967; Jacoby et al., 2015; Tien et al., 2016), and direction of motion (Barlow et al., 1964). For this, intrinsic properties may play a pivotal role: For instance, On-Off direction-selective (DS) RGCs employ a range of mechanisms, including directionally tuned inhibitory and likely excitatory input from asymmetrically wired ACs, but their dendritic morphology and active channel distribution substantially contribute to these cells’ DS tuning [reviewed in (Borst and Euler, 2011; Mauss et al., 2017)].
Here, we study the role of a mechanism at the very end of the signal-shaping hierarchy, the encoding of the RGC membrane voltage into spike trains. We performed electrical single-cell recordings from a transient On RGC type that is EGFP-labelled in the transgenic Igfbp5 mouse line (Siegert et al., 2009), which likely corresponds to G18a [“ON trans”; (Baden et al., 2016)], “ON transient small RF” (Goetz et al., 2021), and “6 sn” (Bae et al., 2018). These cells display an unusual non-linear response behaviour in that they are sharply tuned to light stimuli with a high contrast while ignoring stimuli with low contrast (high contrast selectivity). Based on data from single-cell current analysis and computational modelling, we suggest that desensitisation of the spike generator of Igfbp5-positive transient On small (tOn-small) RGCs significantly contributes to both the transience of their light response and their selectivity for high contrasts.
2 Materials and methods
2.1 Animals and tissue preparation
The Tg(Igfbp5-EGFP)JE168Gsat transgenic mice of either sex were obtained from the Mutant Mouse Regional Resource Centre (MMRRC; University of California, Davis, CA, USA). In this transgenic mouse line, the EGFP reporter gene, followed by a polyadenylation sequence, was inserted into the BAC clone RP24-159O10 at the initiating ATG codon of the first coding exon of the Igfbp5 gene so that EGFP expression is driven by the regulatory sequences of the BAC gene (Gong et al., 2003). The resulting modified BAC (BX1812) was used to generate this transgenic mouse line (The Gene Expression Nervous System Atlas [GENSAT] Project, The Rockefeller University, New York, NY, USA), to which we refer in the following as “Igfbp5” line.
All procedures were performed in accordance with the law on animal protection issued by the German Federal Government (Tierschutzgesetz) and approved by the institutional animal welfare committee of the University of Tübingen or the MPI for Brain Research, Frankfurt/M. Mice of both genders (4–8 weeks of age) were housed under a standard 12 h day/night rhythm. For electrical recordings, mice were dark-adapted for ≥1 h before the experiment. The animals were then anesthetized with isoflurane (Baxter) and killed by cervical dislocation.
For immunostainings, the mouse eyes were dissected in cold 0.1 M phosphate buffer (PB), pH 7.4 and the posterior eyecups were immersion fixed in 4% paraformaldehyde (PFA) in PB for 15–30 min at room temperature. Following fixation, retinas were dissected from the eyecup, cryo-protected in graded sucrose solutions (10, 20, 30% w/v), and stored at –20°C in 30% sucrose until use. Retinal pieces were sectioned vertically at 16–20 μm using a cryostat.
For electrical recordings, the mouse eyes were enucleated and hemisected in carboxygenated (95% O2, 5% CO2) artificial cerebral spinal fluid (ACSF) solution, which contained (in mM): 125 NaCl, 2.5 KCl, 2 CaCl2, 1 MgCl2, 1.25 NaH2PO4, 26 NaHCO3, 0.5 L-glutamine, and 20 glucose; maintained at pH 7.4. After removal of the vitreous body, each retina was flat-mounted onto an Anodisc (#13, 0.2 um pore size, GE Healthcare, Maidstone, UK) with the ganglion cell layer facing up. The Anodisc with the tissue was then transferred to the recording chamber of a two-photon (2P) microscope or a Zeiss Axioscope (see below), where it was continuously perfused with carboxygenated ACSF. When using the two-photon microscope, ACSF contained 0.5–1 μM Sulforhodamine 101 (SR101, Invitrogen Steinheim, Germany) to reveal blood vessels and any damaged cells in the red fluorescence channel (see below). All procedures were carried out under very dim red (>650 nm) light.
2.2 Antibodies and immunohistochemistry
Cholinergic amacrine cells were labelled with goat anti-choline acetyltransferase (ChAT, 1:200; Chemicon), GABAergic amacrine cells with rabbit anti-GABA (1:2000; Sigma). A rat anti-glycine antibody labelled all glycinergic amacrine cells and ON cone bipolar cells (1:1000, kindly provided by David Pow, Brisbane, QLD, Australia). Rabbit anti-GFP was used to increase the EGFP signal (1:2000, Molecular Probes). The neurofilament marker SMI32 (mouse monoclonal, 1:1000; Sternberger Monoclonals) and a goat anti-osteopontin antibody (OPN, 1:1000; R&D Systems) were used to label alpha RGCs. Immunocytochemical labelling was performed using the indirect fluorescence method. Sections were incubated overnight with primary antibodies in 3% normal donkey serum (NDS), 0.5% Triton X-100, and 0.02% sodium azide in PB. After washing in PB, secondary antibodies were applied for 1 h. These were conjugated either to Cy3 (Dianova), or Alexa Fluor 488 (Invitrogen).
Confocal images were taken by using a Zeiss LSM 5 Pascal confocal microscope equipped with an argon and a HeNe laser. Images were taken with a Plan-Neofluar 40x/1.3 objective. Figures represent projections calculated from stacks of images with the LSM software or ImageJ (W.S. Rasband).1 Brightness and contrast of the final images were adjusted using Adobe Photoshop.
2.3 Single cell injection with DiI
For dye injections of BCs, enucleated eyes were transferred to oxygenated Ames medium (Sigma-Aldrich, Taufkirchen, Germany) and opened by an encircling cut. The retinas were dissected and embedded in 2% low melting agar (2-hydroxymethyl agarose, Sigma Aldrich), mounted on a vibratome (DSK Microslicer, DTK-1000, Ted Pella, Inc), and cut into 150 μm sections. After another 10 min in Ames Medium, sections were fixed in 4% PFA in PB at 4°C for 15 min. For injections with the fluorescent lipophilic tracer DiI (Molecular Probes) sharp microelectrodes were pulled from borosilicate glass tubing (Hilgenberg, Malsfeld, Germany) and filled with 0.5% DiI solution in 100% ethanol. The dye was injected into EGFP-labelled bipolar cells with 1 nA positive current for 3 min. For DiI injections of ACs and RGCs, retinal whole mounts were fixed in 4% PFA in PB for 15–30 min. After fixation, the tissue was kept in PB at 4°C overnight for the efficient diffusion of dye into the fine dendrites and axons. Filled cells were imaged with a confocal microscope.
2.4 Two-photon microscopy
We used a MOM-type two-photon microscope (designed by W. Denk, MPImF, Heidelberg; purchased from Science Products/Sutter Instruments, Novato, CA, USA). Both design and procedures were described previously (Euler et al., 2009). The system was equipped with a mode-locked Ti:Sapphire laser (MaiTai-HP DeepSee, Newport Spectra-Physics, Germany) tuned to 927 nm, two detection channels for fluorescence imaging (red, HQ 622 BP36; green, D 535 BP 50, or 520 BP 39; AHF, Tübingen, Germany) and a 20x objective (XLUMPlanFL, 0.95 NA, Olympus, Hamburg, Germany). The red fluorescence channel was used to visualize the retinal structure using SR101 staining (see above), the green channel to target EGFP-labelled Igfbp5 RGCs.
2.5 Electrophysiology
In the Igfbp5 line, two types of RGCs were targeted for single-cell electrical recordings: Non-fluorescent sOn-α cells, identified by their very large polygonal somata, and EGFP-labelled tOn-small RGCs, which possess medium-sized somata (∼15 μm in diameter). The identity of the recorded RGC type was confirmed after the recording based on the dendritic morphology. To avoid confounds related to morphological variability across the retina (Bleckert et al., 2014), the cells were collected (and recorded) from the ventral retina (∼0.8 mm from the optic disc).
For recording of light responses (current-clamp), electrodes (with resistances of 5–10 MΩ) contained (in mM): 120 K-gluconate, 5 NaCl, 10 KCl, 1 MgCl2, 1 EGTA, 10 HEPES, 2 Mg-ATP, and 0.5 Tris-GTP, adjusted to pH 7.2 using 1 M KOH. In addition, 4% Neurobiotin (Molecular Probes, Eugene, ORE, USA) and 0.2 mM SR101 were added to reveal the cell’s dendritic morphology. Membrane voltage was corrected for a liquid junction potential of ∼14 mV (Chang et al., 2013). Data were acquired using an Axoclamp-900A amplifier (Molecular Devices GmbH, San Jose, CA, USA) and digitized at 10 kHz. Experiments were carried out at 37°C. After the electrical recordings, the tissue was fixed, and cells were visualized by overnight incubation in 1:1000 Streptavidin-Alexa Fluor 594 (Invitrogen, Darmstadt, Germany). The RGC’s morphology (as image stack) was documented, and Z-projection images were made using ImageJ.
For voltage-clamp recordings, CdCl2 (100 μM) was added to the ACSF to block currents through voltage-gated Ca2+ channels. Electrodes (with resistances of 5–8 MΩ) contained (in mM): 120 Cs-gluconate, 1 CaCl2, 1 MgCl2, 10 Na-HEPES, 11 EGTA, 10 TEA-Cl, Sulforhodamine B (0.005%), adjusted to pH 7.2 with CsOH. Liquid junction potentials of 15 mV were corrected before the measurement with the pipette offset function of the amplifier (Schubert et al., 2008). Typically, the series resistance was around 7–18 MΩ for sON-a RGCs (13.4 ± 4.2 MΩ, mean ± SD, n = 18) and tOn-small RGCs (13.1 ± 3.4 MΩ, mean ± SD, n = 20). Cell capacitance was not compensated. Seal resistances >2 GΩ were routinely obtained. Cells were voltage-clamped at −90 mV and the different step protocols applied. All experiments were carried out at room temperature (20–22°C). Data were acquired using an Axopatch-200B amplifier (Molecular Devices, San Jose, CA, USA), digitized at 10 kHz using the pClamp software (Molecular Devices), and Bessel-filtered at 2 kHz.
2.6 Light stimulation
For light stimulation, we used a small reflective liquid-crystal-on-Silicon (LCoS) display (i-glasses; EST), coupled in the microscope’s optical path. The display was alternately illuminated by two bandpass-filtered (blue, 400 BP 10; green, 578 BP 10; AHF) LEDs, projecting spatio-temporally structured stimuli through the objective lens onto the retina (Euler et al., 2009). Stimulus intensity was measured using a calibrated photometer (Model 842-PE, 200–1100 nm, Newport) set to the respective centre wavelength of the LED filters. Cone photo-isomerization rates were calculated as described previously (Chang et al., 2013); all stimuli featured equal photo-isomerization rates for M- and S-opsin. Five stimulus protocols were used:
(a) a 200 μm-diameter bright spot flashed for 1 s,
(b) a 300 × 1000 μm bright bar moving at 500 μm/s in 8 directions,
(c) a 1 Hz sinusoidally modulated bright spot of varying diameter at 57% mean contrast,
(d) a 1 Hz sinusoidally modulated spot of varying contrasts, with the optimal spot diameter determined by stimulus (c),
(e) a flickering (white noise) spot, with the optimal spot diameter determined by stimulus (c); spot intensity was randomly chosen from a binary distribution for each frame (80 Hz refresh rate) and the mean intensity equal to background intensity.
For sinusoidally modulated stimuli (c, d), contrast was defined as the ratio of s.d. and mean intensity. The retina was always illuminated with a constant background in the (low) photopic range. In case of the protocols a and b, background illumination generated photo-isomerization rates (in 104⋅P*s–1 per photoreceptor) of 1.2, 0.9, and 2.2, and in case of protocols c to e, 2.3, 2.2, and 4.8 for M-opsin, S-opsin and rhodopsin, respectively.
2.7 Data analysis
All electrophysiological data were analysed off-line using custom MATLAB scripts (Mathworks, Ismaning, Germany). To analyse tOn-small RGC responses to moving bars (protocol b), we defined the stimulus direction that generated most spikes (calculating the vector sum of the total spiking responses in 8 directions) as “preferred direction,” and calculated a direction selective index (DSi) as follows:
with preferred direction-response RP and null ( = opposite) direction response RN.
For the white noise stimulus (protocol e), we used the linear-non-linear (LN) cascade model described earlier (Chichilnisky, 2001; Kim and Rieke, 2001; Baccus and Meister, 2002; Field et al., 2010; Wang et al., 2011) to interpret RGC responses. This model consists of a linear filter that determines the cell’s temporal, chromatic and spatial sensitivities, as well as a “static” non-linearity that converts the filtered stimulus into a firing rate. In the time domain, the linear filter is proportional to the spike-triggered average stimulus (STA, the average stimulus preceding each spike) (Chichilnisky, 2001). Therefore, for LN models with identical linear filter but different non-linearities, spike-triggered average stimuli are identical up to a scale factor. We estimated the “static” non-linearity from the non-linear relationship between linear filtered responses (“generator signal”) and real responses.
Due to the limited refresh rate of white noise stimulus (80 Hz), we defined number of spikes in a time bin of △t = 12.5 ms (1/80 Hz) as the cell’s response, and STA was computed as the cross-correlation between response and stimulus. For graded signal, the cell’s response was defined as relative average voltage within each time bin (by subtracting the minimum voltage of the recorded trace).
We used the Difference-of-Gaussians (DOG) receptive field model in two dimensions to describe the spatial structure of the recorded RGCs. We assumed that the RF of each cell could be approximated by the weighted difference of two concentric Gaussian, one for the RF centre and one for the RF surround. A non-linear function was included into the model to convert this weighted sum into the firing rate (Yin et al., 2009). An RGC’s response R to a spot stimulus centred in the RF and with diameter D and contrast c was calculated as follows:
with the non-linear function f, amplitude (AC) and spatial extent (σC) of the RF centre, as well as amplitude (AC) and spatial extent (σS) of the RF surround.
The non-linear function f was determined using the same RGC’s responses to spots with fixed diameter but varying contrasts, fitted by the Naka-Rushton Equation (Naka and Rushton, 1966):
The other four parameters (AC, AS, σC, σS) were determined by fitting the model to area summation data using a numerical search.
2.8 Computational modelling
To explain our experimental data, we constructed a one-compartment Hodgkin-Huxley model (Dayan and Abbott, 2005):
with the specific membrane capacitance cm = 10 nF/mm2, the voltage across the cell membrane (V) in mV, the maximal conductance for transient Na+ current (), delayed-rectifier K+ current and a leakage current ( = 1.2 mS/mm2; = 0.05 mS/mm2; = 0.003 mS/mm2), the reversal potentials for the three currents (ENa = 50 mV; EK = −76 mV; EL = −70 mV), the area of the cell membrane (A = 0.0013 mm2), the injected current (Ie), and a noise component (In) with a 1/f power spectrum (pink noise). Ie is the estimated response to a sequence of binary white noise, using the linear filter and static non-linearity determined from graded voltage responses. In brief, the linear temporal filter was estimated by computing the spike-triggered average (STA). The inner product between the STA and the actual stimuli yielded a generator signal, and the relationship between the actual response and the generator signal was used to estimate the static non-linearity (Chichilnisky, 2001). Ie was computed by first convolving the light stimuli with the linear filter and then passing the results through the non-linearity. Ie were current steps with varying amplitudes. The gating variables were updated as following:
As rate constants for m, h, and n we used (Hodgkin and Huxley, 1952):
After each spike, the voltage dependences of the rate constants (α and β) for m and h were increased by 1.55 mV for the sensitising model (i.e., tOn-small RGCs) and by 0.01 mV for the non-desensitising case (i.e., sON-α RGCs). For example, after an action potential, changed with dV = 1.55 mV for the desensitising model and 0.01 mV for the non-desensitising model. The recovery from this shift in voltage dependence followed an exponential function with a time constant of 5 s.
2.9 Space-clamp modelling
Based on the differences in neuronal morphology, we built a ball-and-stick model for each cell type (Ran et al., 2020). To get precise measurements of dendritic lengths and radii for sOn-α and tOn-small cells, we extracted this information from published morphologies of 8 w (n = 4) and 6 sn (n = 6) cells, reconstructed from a mouse retina electron microscopy data set.2 Soma diameters were estimated from our data [∼20 μm for sOn-α (8 w) cells and ∼15 μm for tOn-small (6 sn) cells]. We used the median of these values in our model, which was implemented in NEURON.3 Each dendritic portion was further divided into multiple segments with a maximal length of 7 μm. The 1D model can be characterized by the cable equation:
where V is the voltage across the cell membrane, x is the distance along the cable, d is the dendritic diameter, ra is the intracellular resistivity, Cm is the specific membrane capacitance, Icomp is the compensated current to clamp the voltage at the soma according to the command voltage, Iion represents the sum of the sodium current and the leak current, the current dynamics are described following Fohlmeister and Miller (1997) as:
The elicited potential of each dendritic segment was recorded. Model parameters, channel conductances, and dendritic parameters are shown in Tables 1–3, respectively.
2.10 Statistics
Data presentation and statistical analysis were performed using Matlab or in the R programming language. The Wilcoxon signed-rank test was used to determine statistical significance between different conditions. Significance was defined as[: p > 0.05 = non-significant (ns), p < 0.05 = *, p < 0.01 = **]. Mean values in text and figures are given as mean ± standard error of the mean (SEM).
3 Results
3.1 EGFP-expressing neurons in the Igfbp5 mouse retina
Transgenic mice, in which specific RGC types are fluorescently labelled, greatly facilitate investigating RGC function (e.g., Münch et al., 2009; Bleckert et al., 2014; Rousso et al., 2016; Yao et al., 2018). While screening the Gene Expression Nervous System Atlas (GENSAT) database of transgenic mice for selective lines (Siegert et al., 2009), we were struck by the very distinctive pattern of retinal EGFP expression in the Igfbp5 (insulin-like growth factor-binding protein 5) mouse line (Figures 1A, B). In the Igfbp5 retina, the EGFP is expressed in two prominent, equally thick bands of processes along IPL sublaminae 2 and 3 (Figure 1A and Supplementary Figures 1A–C), extending just in between the two choline-acetyltransferase (ChAT) -positive bands (Supplementary Figures 1D, E). This pattern is reminiscent of the glypho (glycogen phosphorylase) staining in the macaque monkey retina (Majumdar et al., 2008). We detected EGFP in a subset of BCs and ACs in the inner nuclear layer (INL; Figure 1A; for details, see Supplementary Figures 1A–D), as well as in some RGCs and displaced ACs in the ganglion cell layer (GCL; Figure 1B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g001.jpg)
Figure 1. EGFP-labelled cells in the Igfbp5 transgenic mouse line. (A) EGFP expression in a vertical section of the Igfbp5 retina, with EGFP expressed in bipolar cells (BCs, asterisks) and amacrine cells (ACs) in the inner nuclear layer (INL), as well as in displaced amacrine cells and retinal ganglion cells (RGCs) in the ganglion cell layer (GCL). EGFP-labelled processes extend along sublaminae 2 and 3 of the inner plexiform layer (IPL). (B) Whole-mount showing mosaic of EGFP-expressing somata in the GCL. Examples of putative displaced ACs (red arrows), Off RGCs with medium-sized (white arrows) and large (yellow arrow) somata, and On RGCs (circles) indicated. Examples of EGFP-positive RGCs: tOff-α RGC (C), and tOn-small RGC (D) filled with Neurobiotin during electrical recording. Below each panel, z-projections of the cell’s dendritic arborisation. Scale bars: (A), 20 μm; (B), 25 μm; (C) and (D), 50 μm.
In the GCL, four cell types could be distinguished based on their soma size and labelling intensity (Figure 1B): very bright cells with small somata (soma diameter: 9.1 ± 0.5 μm, mean ± s.d., n = 29), dimmer cells with medium-sized (12.0 ± 0.5 μm, n = 24) or larger somata (14.8 ± 0.9 μm, n = 34), and very dim cells with large somata (17.5 ± 0.7 μm, n = 14).
Dye-injections (n = 40) revealed that the small, brightly labelled cells were likely displaced wide-field ACs (Lin and Masland, 2006): They featured 6–10 primary dendrites that sometimes bifurcated close to the soma (∼14 dendrites/cell in total), rarely crossed each other and narrowly stratified in IPL sublamina 3 (Supplementary Figures 1F, 2A). The dendritic fields of neighbouring cells overlapped frequently, resulting in dense retinal coverage (Supplementary Figures 2B, C). These putative ACs closely resembled glypho-positive On ACs in the macaque (Majumdar et al., 2008).
The three other lgfbp5-positive cells in the GCL were monostratified RGCs (Figures 1C, D and Supplementary Figures 3A–C) with their dendrites in either one of the two EGFP bands (cf. Table 4). The cells with the largest, dimmest somata were transient Off-alpha (tOff-α) RGCs (Figures 1B, C and Supplementary Figure 3B; van Wyk et al., 2009; Krieger et al., 2017; Ran et al., 2020), as confirmed by immunolabeling (Supplementary Figures 3E–H) for markers such as SMI32 (Supplementary Figures 3E, F), which strongly labels sustained On alpha (sOn-α) and tOff-α cells in mouse (Coombs et al., 2006; Bleckert et al., 2014). The remaining two labelled cell populations had dendritic arbour diameters of around 200 μm.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-t004.jpg)
Table 4. Igfbp5-positive RGCs and their presumed morphological and genetic counterparts in earlier mouse studies.
The lgfbp5-positive cells with medium-sized somata were presumably Off cells, as their dendrites stratified in sublamina 2; their morphology (Supplementary Figure 3A) resembled that of “4i” or “4on” cells described in Bae et al. (2018), cluster 3 (Kong et al., 2005), PV4 (Farrow et al., 2013), CB2 (Sümbül et al., 2014), and RGB1 or “outer” RGB3 cells (Sun et al., 2002).
The mouse lgfbp5-positive RGCs with the larger somata (Figures 1D, 2A, C and Supplementary Figure 3C) were presumably On cells as they stratified in sublamina 3; their morphology resembled that of “6sn” (Bae et al., 2018), “ON transient small RF” (Goetz et al., 2021), cluster 2 (Kong et al., 2005), PV2 (Farrow et al., 2013), “cluster X” (Sümbül et al., 2014), and “inner” RGB3 cells (Sun et al., 2002). In the transcriptomics study by Tran et al. (2019), this RGC type remained elusive. For simplicity, we refer to the lgfbp5-positive On RGCs in the following as transient On small (tOn-small) cells. Mouse transient On-alpha (tOn-α) RGCs stratify at the same level as the tOn-small cells but have a slightly larger dendritic tree (Goetz et al., 2021) and are osteopontin-positive (Krieger et al., 2017). We think that the tOn-small RGCs are good candidates for being homologous to primate parasol cells (for discussion see also Hahn et al., 2023; Supplementary Figure 3I).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g002.jpg)
Figure 2. Igfbp5-positive tOn-small ganglion cells differ morphologically from sustained sOn-α cells. Examples of dye-injected RGCs: three tOn-small (A) and three sOn-α cells (B). To avoid confounds related to morphological variability across the retina (Bleckert et al., 2014), the cells were collected (and recorded) from the ventral retina (∼0.8 mm from the optic disc). (C) Soma diameter vs. dendritic field diameter (left) and vs. number of dendritic tips (right) for the two RGC types, with distribution histograms at the sides. Scale bar: (A,B), 100 μm. For a more detailed overview of Igfbp5-positive RGCs and comparison with primate retina see Supplementary Figure 3.
3.2 tOn-small RGC responses are highly transient
When we characterized light-evoked signals in tOn-small cells using two-photon-guided electrical recordings (see section “2 Materials and methods”), we were intrigued by their very transient responses—with the spike rate increasing almost instantaneously at light-onset but then quickly dropping to zero within ∼200 ms and no response at light-offset (Supplementary Figure 4A). As such transient responses were described in a type of direction-selective On RGC in rabbit (Kanjhan and Sivyer, 2010), we recorded tOn-small cell responses to a moving bar stimulus, but did not find any substantial directional tuning in these cells (Supplementary Figure 4B; DSi = 0.043 ± 0.019, n = 5 cells).
In the following, we studied the mechanisms underlying the characteristic transient responses of tOn-small cells in the mouse retina (Figure 2A). For comparison, we recorded the well-described sOn-α RGCs (Peichl et al., 1987), which are known for their sustained responses and their high sensitivity for small contrast changes (high contrast sensitivity) (Pang et al., 2003; van Wyk et al., 2009; Schwartz et al., 2012; Bleckert et al., 2014; Krieger et al., 2017; Figure 2B). Morphologically, the two RGC types could be distinguished easily, with tOn-small cells having smaller somata and dendritic fields, but more dendritic tips than sOn-α cells (Figure 2C).
3.3 tOn-small cells encode selectively high-contrast signals
First, we studied the spatial RF organisation of tOn-small cells using spot stimuli with varying diameters centred on the soma and sinusoidally modulated at 1 Hz (Figures 3A, B), in comparison to sOn-α cells. Using spectral analysis, we then calculated the amplitude of the cells’ fundamental response component (see section “2 Materials and methods”). With increasing spot diameter, amplitudes first increased and then declined (Figure 3B), indicative of centre-surround antagonism. To quantify the spatial RFs, we used a Difference-of-Gaussians (DOG) model to interpret the recorded area summation data (see section “2 Materials and methods”). Consistent with their larger dendritic field diameters (Figure 2C), sOn-α cells had larger RF centres than tOn-small cells (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g003.jpg)
Figure 3. tOn-small ganglion cells have smaller receptive fields (RFs) and higher contrast thresholds compared to sOn-α cells. (A) Electrically recorded voltage responses of a tOn-small (blue) and a sOn-α (black) to a 1-s flashed spot (200 μm in diameter). (B) Fundamental (F1) spiking response (F1 Freq in Hz) as a function of spot diameter (1 Hz, 57% mean contrast) for the same tOn-small cell as in (A) (dotted blue line). Data fit with Difference-of-Gaussians (DOG) receptive field model (solid blue line). (C) Mean and SEM of estimated RF profiles, using difference-of-Gaussians (DOG) RF model for n = 7 tOn-small and n = 6 sOn-α cells (see section “2 Materials and methods”). (D) Responses of the same cells as in (A) to a 1-Hz sinusoidal stimulus (270 μm in diameter). (E) Representative responses of a tOn-small cell for different contrasts. (F) F1 Freq as a function of stimulus contrast (D1; 270 μm diameter spot; data fit with Naka-Rushton equation; see section “2 Materials and methods”) for the same tOn-small cell as in (E). (G) Response-contrast functions for tOn-small cells (blue, spikes, n = 7 cells; cyan, voltage, n = 3 cells). (H) Normalized response-contrast curves for an additional set of sOn-α cells (black, spikes, n = 6 cells; grey, voltage, n = 6 cells; Naka-Rushton fits). (I) Normalized response-contrast curves [dataset from (G,H); Naka-Rushton fits]. Symbols and error bars represent mean and SEM, respectively. For responses of tOn-small cells to moving bar stimuli see Supplementary Figure 4.
In the DOG model, we included a static-non-linearity to convert the weighted stimulus into a firing rate. Because the area summation data by itself was not sufficient to characterize this non-linearity, we also presented spots with fixed diameter but with different contrasts (Figure 3D). Here, we used the spot diameter (270 μm) that elicited the maximal spiking response. Other than sOn-α cells, tOn-small cells were insensitive to weak stimuli: Sinusoidally modulated spots with contrasts lower than 40% hardly evoked any spiking response (Figures 3E, F).
We also analysed the graded voltage responses recorded in the whole-cell current-clamp mode (with spikes digitally removed, see section “2 Materials and methods”) and found that the strong threshold-like behaviour was not reflected in the graded response: tOn-small cells showed a detectable depolarization even to the smallest contrast tested (5%, Figure 3G, cyan circles vs. blue circles), resulting in a contrast sensitivity curve similar to that for the sOn-α cells’ spiking response (Figure 3H; black squares). Note that sOn-α cells displayed a very similar contrast sensitivity curve for spikes and graded voltage (Figure 3H; filled black vs. open grey squares). Together, this suggests that the conversion of voltage into a spiking response differs between tOn-small and sOn-α cells (Figure 3I). A possible explanation may be a larger difference between resting potential and spiking threshold in tOn-small vs. sOn-α cells. While the resting potential (Vrest) of tOn-small cells (−66.8 ± 3.4 mV; n = 6) indeed was slightly more hyperpolarized than that of sOn-α cells (−64.8 ± 3.3 mV; n = 6), the difference seems too small to explain the observed difference in contrast sensitivity alone. Note that the estimate of the resting potentials may not be accurate due to the perturbation of the intracellular milieu in whole-cell experiments. We tried to mitigate this problem by reading the resting potential as soon as the whole-cell recording is established. We also estimated the difference between the spiking threshold and the resting potential and found that the results of the tOn-small cells (12.2 ± 3.5 mV; n = 3) were only slightly higher than those of sOn-α cells (9.3 ± 2.3 mV; n = 4).
The difference in contrast sensitivity between spiking and graded voltage responses in tOn-small cells is striking. Related to this was our finding that when stimulated with a bright flash, the cells displayed a sustained graded voltage response—in contrast to the transient nature of their spiking response (Figure 3A). We therefore studied temporal processing in tOn-small (and, for comparison, sOn-α) cells with a homogeneous white-noise stimulus (see section “2 Materials and methods”), using again the optimal spot diameter. First, we calculated the spike triggered average (STA) as an estimate of the cell’s temporal linear filter (Chichilnisky, 2001). tOn-small cells had highly biphasic filters, with an upward-pointing lobe (On-lobe) close to the time of spike and an almost similarly strong lobe pointing downward (Off-lobe) further away (Figure 4A1; blue curve). We quantified the filter’s biphasic nature by calculating the ratio between Off- and On-lobe peak (Chander and Chichilnisky, 2001), yielding a biphasicity index (Bi) of 0.60 ± 0.19 (n = 6) for the spike response of tOn-small cells. In contrast, the linear filter estimated from their graded voltage response (Bi = 0.28 ± 0.17, n = 2) and the STA-derived linear filters of sOn-α cells (Bi = 0.15 ± 0.09, n = 5) were more monophasic (Figure 4A1; dashed cyan and black curve, respectively).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g004.jpg)
Figure 4. Spike desensitisation in tOn-small ganglion cells accounts for their biphasic temporal filter and a high contrast-selective threshold. (A) Temporal linear filter (A1) and static non-linearity (A2) estimated using a white noise stimulus (see section “2 Materials and methods”) for tOn-small (blue, spikes, n = 6 cells; cyan, graded voltage responses, n = 2 cells) and sOn-α cells (black, spikes, n = 5 cells). (B) Electrically recorded voltage responses of a tOn-small cell [(B1), bottom, resting potential Vrest = –66 mV] in response to current injections of varying amplitude [(B1), top]. (B2), Firing rate of the same cell right after current injection onset (0–0.1 s, red) and at steady-state (1–3 s, brown), as functions of injected current [for time windows, see (B1)]. (B3) Relationship between spike desensitisation and injected current for tOn-small (blue, spikes; cyan, graded voltage responses, n = 3) and sOn-α cells (black, spikes, n = 4). Spike desensitisation was quantified as transience index (Tispike1−(r1−3s/r0−0.1s) or TiV1−(V1−3s/V0−0.1s), with spike rate r and baseline-corrected graded potential V, respectively). Symbols and error bars represent mean and SEM, respectively. (C) Temporal linear filter (C1) and static non-linearity (C2) estimated using simulated responses of two modelled RGCs (left: with strong spike desensitisation; right: with weak spike desensitisation) to a white noise stimulus (for details, see text). (D) Voltage responses of the two modelled cells to current injection of 0.1 nA (D1) and their firing rates as functions of injected current [(D2), analogous to recording in (B2)]. (E) Voltage recording of a sOn-α cell (Vrest = –64 mV) in response to current injections of 0.1 nA (E1) and its firing rate as a function of injected current [(E2), analogous to (B2)].
The difference between linear filters estimated from spike vs. voltage responses indicates that the non-linearity underlying spike generation in tOn-small cells is not a “trivial” static non-linearity, i.e., a simple non-linear function whose output only depends on the present value of input. Because the linear filter estimation is expected to be resistant to a static non-linearity (Chichilnisky, 2001), the involvement of a static non-linearity should result in similar filters—independent of the response modality used (spikes vs. graded voltage). To characterize the non-linearity in tOn-small cells, we calculated the relationship between the linear response (“generator signal” = the convolution of estimated linear filter and stimulus) and the recorded response to the white noise stimuli. The resulting non-linear function for the spike responses of tOn-small cells had a much higher threshold than that for their voltage responses and the spike responses of sOn-α cells (Figure 4A2).
Many response features of RGCs, including temporal properties and contrasts sensitivity, are shaped by synaptic interactions (see sections also “1 Introduction and 4 Discussion”). However, the clear differences in transiency and contrast sensitivity between spiking and graded voltage responses observed in tOn-small (but not sOn-α) RGCs point at a contribution of the tOn-small cells’ intrinsic spike generator (Li et al., 2022). Therefore, in the following we focused on this final step in retinal signal transformation.
3.4 Spike generation in tOn-small RGCs is readily desensitised by constant current injection
The comparison between spike and graded voltage responses revealed two distinct features of spike generation in tOn-small cells: (i) the high-threshold non-linearity and (ii) the pronounced biphasicity of the linear filter. The first feature could arise simply from a spike threshold that is substantially higher than the resting potential. The second feature implies more complicated non-linearities, such as spike desensitisation (Goldin, 1999). We tested the spike generator by injecting constant current of different amplitudes into the cells while recording their voltage responses (Figure 4B1). We found that for tOn-small cells, injection of a large current elicited a burst of spikes at injection onset but then spiking stopped rapidly (Figure 4B1, dark blue trace). The initial spike frequency (time window: 0–0.1 s) scaled almost linearly with the amplitude of injected current (Figure 4B2, red symbols), whereas the steady-state spike frequency remained as low as ∼10 spikes/s or less (Figure 4B2, orange symbols). As a result, transiency of the spike response increased with stimulus strength, as indicated by the increasing ratio between initial spiking response and steady state (Figure 4B3). This increase was absent or less pronounced in graded voltage responses of tOn-small cells and spiking responses of sOn-α cells, respectively. Together, this points at strong desensitisation of the spike generator in tOn-small cells.
3.5 Spike desensitisation model mimics response properties of tOn-small RGCs
We constructed a simple model to explore how varying levels of spike desensitisation affect the signal encoding properties of desensitising and non-desensitizing cells, such as tOn-small and sOn-α cells, respectively (see section “2 Materials and methods”). We modelled responses to the white-noise stimulus by calculating the amplitude of the input current to the cell using a simple linear-non-linear (LN) model. First, for the desensitising case (Figures 4C1, C2; left), we estimated the linear filter and non-linearity from graded voltage responses measured in tOn-small RGCs (Figures 4A1, A2). Using the voltage output of the modelled cell, we then estimated linear filter and non-linearity for the spike responses. We found that the spike response-derived linear filter was more biphasic than that based on the graded voltage response, very similar to what we measured in tOn-small cells (cf. Figures 4A1, A2). For the non-desensitising case (Figures 4C1, C2, right), the two linear filters were nearly identical. In both model cells, the non-linearity for spike responses exhibited a higher threshold than that for graded voltage responses, however, only the threshold in the desensitising cell was higher than the linear response to 0% contrast (generator signal = 0; Figure 4C2; cf. panel Figure 4A2). Therefore, low-contrast stimuli (generator signal ∼0) evoked only sub-threshold responses in the desensitising model cell—very similar to the responses recorded in tOn-small cells (cf. Figure 3E).
We also used current steps with different amplitudes as model input (Figure 4D). As expected, the desensitising model cell displayed little or no spike response at the steady state (1–3 s) of the current step (Figure 4D, left), very similar to the responses recorded in tOn-small RGCs (Figure 4B1). In contrast, the steady-state response of the non-desensitising model cell was nearly as strong as its response at the current injection onset (0–0.1 s; r2−3sr0−0.1s) for each current amplitude (Figure 4D, right), very similar to the recorded response of sOn-α cells (Figure 4E).
In conclusion, our model supports that spike desensitisation, and a high spiking threshold may contribute to two distinct response features of tOn-small RGCs: their strong transience and their selectivity for high contrasts. Note that our simple desensitising model generates responses very similar to those of tOn-small cells, whereas the match between the non-desensitising model and sOn-α cells was not as good: The non-linearity is “steeper” in the non-desensitising model cell than that estimated from sOn-α cell recordings (Figure 4C2 right vs. Figure 4A2). This discrepancy suggests that a sOn-α cell model requires additional properties, such as, for example, a higher noise level in the membrane potential, which can result in a smoother non-linearity (Anderson et al., 2000). Although our model is rather simple and does not take the detailed desensitization states of the voltage-gated channels into account, it supports that desensitization of voltage-gated sodium channels (VGSCs) can explain the differences we observed in two types of RGCs.
3.6 Voltage-gated Na+ channels in tOn-small RGCs are high-threshold activated
Our modelling predicts that the Na+ channel desensitization observed in tOn-small cells may result from differences in VGSC activation, inactivation and/or recovery from inactivation (Carter and Bean, 2011). To test this experimentally, we recorded Na+ currents in these cells and, for comparison, in sOn-α cells (Figures 5, 6). To isolate VGSCs, we blocked voltage-gated K+ channels with Cs+ and TEA in the electrode solution and voltage-gated Ca2+ channels with Cd2+ in the extracellular solution (see section “2 Materials and methods”). With voltage-step protocols to characterize VGSC activation (Figure 5A1) and inactivation (Figure 5A2) we found VGSCs in tOn-small cells to activate at higher potentials than those in sOn-α cells (Figure 5A3; tOn-small, −51.1 ± 3.9 mV, n = 6; sOn-α, −67.5 ± 3.2 mV, n = 7; at 5% activation estimated from sigmoidal fit), providing a likely explanation why a substantial depolarisation (approx. >15 mV from VRest) is required to trigger spikes in tOn-small cells. This higher activation threshold is also consistent with the tOn-small cells’ selectivity for high contrasts. However, the inactivation profiles were almost identical (Figure 5A3) for the used step protocol, suggesting that the VGSCs in the two cell types share similar steady-state inactivation properties.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g005.jpg)
Figure 5. Voltage-gated Na+ channels in tOn-small cells desensitise strongly. (A) Voltage-clamp recording of a tON-small cell to step protocols probing activation (A1) and inactivation (A2) of voltage-gated Na+ channels [VGSCs; (A1), from –90 mV for 1 s to different voltages (–90 to –10 mV); (A2), from different potentials (–90 to –10 mV) for 1 s to 0 mV; ΔV = 10 mV for both protocols]. (A3) Normalized conductance as functions of step voltage (Vstep) for channel activation (squares) and inactivation (circles) in tOn-small (blue, n = 6) and sOn-α cells (black, n = 7). (B) Probing inactivation of VGSCs with 3-ms test pulses (to 0 mV) delivered before (P1) and after (P2) the 1-s inactivating steps [(B1), top; –90 to 0 V, ΔV = 10 mV, followed by a hyperpolarising step for 20 ms to –90 mV]. Na+ currents elicited by the 2nd test pulse (P2) in a tOn-small cell [(B1), bottom]. (B2) Ratio of Na+ current amplitudes elicited by the two test pulses [cf. (B1)] plotted against inactivating voltage. (C) Probing “inactivation sensitivity” using sequences of 10-ms depolarizing voltage pulses (–70 to 0 mV). (C1) Na2+ currents elicited in a tOn-small cell by a 20-Hz sequence. (C2) Peak Na+ currents as a function of pulse number for 20 (Δt = 50 ms) and 8 Hz (Δt = 125 ms blue circles and black squares for tOn-small and sOn-α cells, respectively). Peak Na+ currents were normalized to the amplitude evoked by the first pulse. Data from tOn-small (circles, n = 7) and sOn-α cells (squares, n = 5) for 20 Hz pulses, and tOn-small (circles, n = 6) and sOn-α cells (squares, n = 3) for 8 Hz pulses. Wilcoxon rank sum test was used for determining statistical significance (*p < 0.05, **p < 0.01) between different conditions in (A3,B2,C2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g006.jpg)
Figure 6. Steady-state slow inactivation of VGSCs shapes tOn-small RGC responses. (A) Same protocol as in Figure 5B, but with a 100-ms hyperpolarizing pre-pulse (VPP = –90 mV) inserted between inactivation step and test pulse P2; P1 not shown (A1). Normalized conductance as a function of inactivation step voltage for tOn-small (filled circles) and sOn-α RGCs [open squares; (A2)]. (B) Like in (A) but with pre-pulse duration varied (Δtpp = 0, 10, 20, 100, and 500 ms) for a single inactivation step voltage (Vinact = –50 mV). (C) Like in (A) but with fixed Δtpp = 10 ms and Vinact = –40 mV, while varying pre-pulse voltage (VPP = –90 to –50 mV). (D) Same protocol as in Figure 5B but with the duration of the inactivation step (Δtinact_step) varied (D1). Test pulse current ratio (I(P2)/I(P1)) as a function of Vinact for Δtinact_step = 1 (D2), 2 (D3), and 4 s (D4). Wilcoxon rank sum test was used for determining statistical significance (* = p < 0.05) between different conditions in (A–D).
3.7 Voltage-gated Na+ channels in tOn-small RGCs need stronger hyperpolarization and longer times to recover from inactivation
It is known that VGSCs can undergo inactivation at different timescales and that their recovery from different states of inactivation may differ in its time- and voltage-dependence (Tsai et al., 2011; Martiszus et al., 2021). Hence, we next wanted to look at these aspects in the two cell types. To this end, we probed the cells with two test pulses (P1, P2), with P2 following a 1 s “inactivating” voltage-step and a short 20 ms hyperpolarizing pre-pulse (Figure 5B1). While the previous protocol (Figure 5A2) tested the combined effects of fast and slow inactivation, the insertion of a pre-pulse enables VGSCs to recover from fast inactivation and, hence, allowed us to probe the effect of slow inactivation (Kim and Rieke, 2003; Silva, 2014). In the presence of this pre-pulse, the inactivating step reduced the test pulse response in both RGC types, but much more so in tOn-small cells (Figure 5B2). This suggests that the VGSCs’ slow inactivation differed between the two cell types with respect to voltage-dependence, with the VGSCs in tOn-small displaying consistently stronger inactivation for step voltages ≥−70 mV.
To probe the time-dependence of recovery from inactivation, we first applied a pulse-train protocol with two different frequencies (8 and 20 Hz, Figures 5C1, C2), mimicking the activation by spikes (Kim and Rieke, 2003). We quantified VGSC inactivation by calculating the ratio between the Na+ current activated by the ith pulse (INa(i)) and that activated by the 1st pulse of a train (INa(1)). We found that INa ratio (INa(i)/INa(1)) in tOn-small cells dramatically dropped at 20 Hz already for the 2nd pulse to ∼75% and for the 10th pulse to less than 50%, whereas in sOn-α cells, INa ratio was only slightly reduced by consecutive pulses (Figure 5C2). This indicates that VGSCs in tOn-small cells may need substantially more time to recover from inactivation compared to those in sOn-α cells.
3.8 Steady-state slow inactivation of VGSCs shapes tOn-small RGC responses
To study time- and voltage-dependence of RGC responses in our voltage-clamp experiments more closely, we went back to the two-pulse protocol (Figure 6A1). First, we varied the duration of the hyperpolarizing pre-pulse (ΔtPP) between inactivation step and second test pulse (P2) to determine the time-dependence of recovery from slow inactivation (Figure 6B). With ΔtPP0, inactivation was virtually the same for both cell types (cf. Figures 5A2, A3), but became significantly different for ΔtPP20ms before starting to approach similar values for ΔtPP500ms (Figure 6B2), which supports a difference in time-dependence of recovery from slow inactivation. Next, we kept the duration of the pre-pulse constant but varied its voltage (VPP) to test the voltage-dependence (Figure 6C1). We found tOn-small cells recovered only to approx. half the levels compared to sOn-α cells for all tested VPP values (Figure 6C2), arguing for an additional difference in voltage-dependence of recovery. Finally, we changed the duration of the inactivating step (Δtinact_step of 1, 2, or 4 s; Figure 6D1) to make sure the difference between the cell types were not due to insufficient inactivation of the sOn-α cells; this seemed to have not been the case, as for all tested inactivation durations, tOn-small cells consistently showed stronger inactivation than sOn-α cells (Figures 6D2–D4).
Taken together, our data suggests that the differences observed in the spiking behaviour of tOn-small vs. sOn-α cells at least partially arise from VGSCs in tOn-small cells: These VGSCs (i) activate at a higher threshold and require (ii) stronger hyperpolarization and (iii) more time to recover from steady-state slow inactivation. These results support that idea that properties of the spike generator importantly contribute to the temporal coding properties observed in tOn-small RGCs.
3.9 Differences between tON-small and sOn-α cell responses are unlikely due to space-clamp errors
While we isolated the recorded cells from Ca2+-dependent synaptic input and blocked voltage-gated K+ channels in our voltage-clamp experiments, we cannot exclude that the cells’ morphologies resulted in different space-clamp situations that affected our results. We therefore assessed the effect of morphological differences with computational methods (Poleg-Polsky and Diamond, 2011). We built a ball-and-stick model for both sOn-α and tOn-small cells (see section “2 Materials and methods”; Figure 7) that captured dendritic features, such as diameter and segment length (Fohlmeister and Miller, 1997). We tested the model for different reversal potentials for Na+ (RevNa) and unspecific leak currents (Revleak) and applied a step protocol comprising inactivating voltage steps and the P2 test pulse from Figure 5B1. The goal was to test if voltage-clamp errors may have contributed to the differences we observed between the two RGC types in Figures 5B, C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1337768/fncel-17-1337768-HTML/image_m/fncel-17-1337768-g007.jpg)
Figure 7. Modelling of space-clamp properties in tON-small and sOn-α retinal ganglion cells. (A) Ball-and-stick model showing membrane depolarisations at soma (left) and distal dendritic tips (right) elicited with 3-ms test pulses (P2, –90 to 0 mV) delivered after 1-s inactivating steps (–90 to 0 mV, ΔV = 10 mV, followed by a hyperpolarising step for 20 ms to –90 mV). Membrane potentials elicited by inactivating steps and 3-ms test pulse in a tOn-small cell (top row) and a sOn-α cell (bottom row) are shown. Reversal potentials (Revs) of the “ball-and-stick” model were set to RevNa = +35 mV and Revleak = –40 mV. (B) As in (A) but model parameters set to RevNa = +90 mV and Revleak = –40 mV. recov, recovery period. (C) As in (A) but with model parameters set to RevNa = +35 mV and Revleak = –18 mV.
The specific dendritic morphology may affect the local ion reversal potential because intracellular ion concentrations in high-resistance distal tips are likely to be strongly influenced by ion flux through local channels (vs. diffusion of intracellular solution from the electrode tip at the soma along the dendrite). Therefore, the actual local reversal potential for a particular ion is difficult to estimate but likely different from the theoretical value. We therefore tested for RevNa = +90 mV (the theoretical values defined by our solutions) and RevNa = +35 mV (assuming a higher dendritic Na+ concentration). We found that for RevNa = +35 mV, the voltage drop after the depolarization step had a longer time constant in sOn-α vs. tON-small cells, but the voltage recovered to baseline before the P2 test pulse (within less than 20 ms) (Figure 7A). Importantly, the modelled voltage amplitudes were very similar for the two RGC types. Therefore, it is unlikely that the smaller current amplitudes observed in tOn-small cells (Figure 5B2) are due to space-clamp issues. For RevNa = +90 mV, the dendritic voltage amplitudes for the P2 test pulse were even higher in tOn-small cells than in sOn-α cells (Figure 7B), which are expected to result in higher currents in tOn-small cells, but which was opposite to what we measured experimentally (Figure 5B2).
Similar to RevNa, also non-voltage-sensitive leak currents of a neuron can affect the space-clamp properties (Huang et al., 2015). Thus, we modelled the impact of the leak current (Revleak) for the condition with the higher dendritic Na+ concentration (RevNa = +35 mV) (Figures 7A, C). We estimated Revleak for two extreme ratios of Cl– and Na+ conductances: Revleak = −40 mV (Cl–:Na+ conductance = 12:1) and Revleak = −18 mV (Cl–:Na+ conductance = 3:1). For both Revleak values, the voltage amplitudes in the dendritic tips were virtually identical (Figures 7A, C), which suggests that the leak currents were unlikely to substantially affect our measurements.
In summary, our modelling results suggest that while space-clamp problems may indeed have affected our Na+ current measurements, they also indicate that sOn-α RGCs were expected to be more strongly affected than tOn-small RGCs (i.e., time constant)—which does not match our experimental observations. Therefore, we consider it unlikely that space-clamp issues caused the observed difference in Na+ channel inactivation between sOn-α and tOn-small RGCs.
4 Discussion
Here we studied the underlying cell-intrinsic mechanisms that encode membrane voltage into spike patterns in two distinct types of RGC in the mouse retina. Recordings of light- and current-elicited voltage responses showed differences in the relationship of membrane voltage and elicited spikes between tOn-small and On-α cells: Unlike in sOn-α RGCs, spiking but not the voltage signal in tOn-small RGCs was very transient and high-contrast selective. Modelling and voltage-clamp recordings revealed that this response pattern of tOn-small cells is largely shaped by stronger desensitisation of the cell’s spike generator, i.e., its VGSCs. Our results support that notion that—in complement with upstream mechanisms and as the last step of the signal retinal processing chain—the intrinsic properties of an RGC can importantly shape the retina’s output to the brain.
4.1 Multiple mechanisms can shape RGC response kinetics
Retinal ganglion cell response properties can be shaped at different stages of retinal signal processing. First, RGCs can directly inherit diverse temporal response properties from their excitatory presynaptic partners, the BCs. It has been shown that transient BCs, such as types 5t, 5o, 5i, and XBC, mostly stratify in the central bulk of the IPL (Baden et al., 2013), therefore RGCs may become transient simply by having their dendrites tap into the respective IPL layers to collect transient excitatory input. Interestingly, the tOn-small cells likely correspond to the “6 sn” cells classified in an EM dataset (Bae et al., 2018), which stratify slightly more toward the centre of the IPL than the sOn-α cells, and thus can form synapses with more transient BC types. This may contribute to transient responses in tOn-small RGCs and sustained responses in sOn-α cells. Second, inhibitory AC circuits can sharpen the time course of RGC responses either by shaping the signals at BC axon terminals or RGC dendrites [reviewed in (Zhang and McCall, 2012; Diamond, 2017)]. Third, intrinsic RGC properties, such as kinetics of postsynaptic glutamate receptors, dendritic morphology and active dendritic channels, define the computations executed of distinct RGC types (Brombas et al., 2022). Forth, RGC responses can be shaped by properties of the intrinsic spike generator. For example, in suppressed-by-contrast RGCs, a depolarization block, resulting from a low VGSC conductance, short axonal initial segment (AIS), and selective expression of Nav isoforms (see below), defines the cells’ response (Wienbar and Schwartz, 2022). Another example would be direct (neuro)modulation of Nav function by reactive oxygen species (ROS), as it was recently shown for sustained RGCs (Smith et al., 2023). In complement, our study suggests that the functional properties of VGSCs expressed in tOn-small cells are one of the key mechanisms that shape these cells’ responses. Further experiments are needed to reveal how synaptic and intrinsic features interact with each other to shape RGC output (Li et al., 2022). In particular, it would be interesting to unravel how excitatory and inhibitory synaptic input sets the operational range of VGSCs in RGCs (Davis and Linn, 2003) and how, in turn, VGSC activity—particularly in the dendrites—modulates integration of neurotransmitter input (Yu, 2006; George et al., 2012; Ran et al., 2020; Brombas et al., 2022).
4.2 Spike generator desensitisation and functional consequences
Our results suggest that desensitisation of the spike generator, likely based on slow inactivation of the VGSCs themselves, shapes spike coding in tOn-small cells in two ways: their responses are more transient and their sensitivity for low contrast stimuli is decreased compared with sOn-α cells. In previous studies, it was shown that slow inactivation of VGSCs participates in the adaptation of RGCs to temporal contrasts (Kim and Rieke, 2003; Weick and Demb, 2011). Another study investigating contrast adaptation in RGCs revealed two types of “plasticity” mechanisms: adaptation and sensitisation (Kastner and Baccus, 2011). Interestingly, the plasticity of the adapting type supports encoding of strong stimuli (i.e., high-contrast), whereas plasticity of the sensitising type promotes encoding of weak stimuli (i.e., low-contrast)—very similar to what we found in tOn-small (high contrast-selective) and sOn-α cells (also sensitive to low-contrast stimuli). Therefore, it is likely that desensitisation of the spike generator plays a functionally crucial role in at least three different aspects of RGC signalling: temporal kinetics (transient vs. sustained), contrast adaptation, and limiting the dynamic range for contrast encoding.
In our model, we simulated spike generator desensitization by changing the voltage dependence of the gating functions in a history-dependent manner (see section “2 Materials and methods”). This can be thought as equivalent to slowing the recovery of the VGSCs from inactivation. Also, our voltage-clamp data points at the VGSCs to be involved in defining the response properties of tOn-small cells: compared to sOn-α cells, the Na+ currents in tOn-small cells exhibited a higher threshold and recovered more slowly from inactivation. Therefore, that tOn-small and sOn-α cells likely differ in their VGSC complement.
Using RT-PCR, four VGSC α-subunit isoforms have been identified in rodent RGCs: Nav1.1, 1.2, 1.3, and 1.6 (Fjell et al., 1997). Among other things, these isoforms differ in their persistent current [reviewed in (Goldin, 1999)]: For example, for more positive membrane potentials, the persistent current of Nav1.6 increases, whereas that of Nav1.1 decreases. A persistent current that increases with depolarization “pulls” the membrane potential toward the spiking threshold and, thus, may favour sustained spiking. Indeed, a recent study provided evidence that TTX-sensitive Nav1.6 channels are dominantly expressed in mouse sustained OFF alpha RGCs, whereas suppressed-by-contrast RGCs likely express a different isoform (Wienbar and Schwartz, 2022). Thus, the expression ratio of VGSC isoforms affects the time course of an RGC’s spiking response. In addition, VGSCs are strongly modulated by accessory ß-subunits (reviewed in Goldin, 1999), which can, for example, slow or accelerate inactivation, and shift voltage-dependence. Therefore, not only VGSC properties but also differential ß-subunit expression may also contribute to the differences we observed between tOn-small and sOn-α cells. Like many other proteins in neurons, the intracellular domains of VGSCs are a substrate for phosphorylation, and thus, VGSCs can be functionally modulated. Depending on the phosphorylation/dephosphorylation state, amplitude and inactivation of currents through VGSCs can significantly vary (Numann et al., 1991; Li et al., 1992; Scheuer and Catterall, 2006). In complement, “extracellular” factors such as temperature and concentration of divalent cations in the extracellular space can modulate the voltage dependence of VGSCs (Hahin and Campbell, 1983; Josephson, 1996; Abdelsayed et al., 2015).
In addition to the Nav1.6 channels, which can be selectively blocked with TTX, it was shown that a small subset of RGCs additionally expresses an unusual TTX-insensitive VGSC α-subunit isoform in the soma and the proximal dendrites: Nav1.8 (O’Brien et al., 2008). These cells had large somata and were neurofilament-positive, suggesting that alpha (but not tOn-small) cells express Nav1.8. In contrast to the above VGSC isoforms, Nav1.8 mediates currents that exhibit very little subthreshold inactivation and recover much faster from inactivation (Cummins and Waxman, 1997; Renganathan et al., 2000). It is possible that channels with such properties contribute to the high contrast sensitivity—for instance, by enabling dendritic spiking (Velte and Masland, 1999)—and the sustained spiking response of sOn-α cells. Nevertheless, we did not observe differences in inactivation between tOn-small and sOn-α cells which argues against a prominent role of Nav1.8 in sOn-α cells—at least under our experimental conditions. Future pharmacological experiments will unravel how TTX-sensitive and TTX-resistant VGCSs together generate response patterns in individual RGCs.
Finally, VGSCs are precisely arranged in type-specific subdomains along the axon initial segment (Van Wart et al., 2007). These VGSC “bands” appear to vary in length and location (relative to the soma) between RGC types (Fried et al., 2009; Wienbar and Schwartz, 2022). It was suggested that this organization shapes the response properties of the RGC spike generator, including the activation threshold (Fried et al., 2009). Whether the VGSC bands along the axon initial segment of tOn-small and sON-α cells differ remains to be investigated.
4.3 Diverse retinal computations are performed by intrinsic properties of RGCs
A great deal of neural computations taking place in the retina ultimately segregates features of the visual stimulus into diverse parallel channels to the brain (reviewed in Wassle, 2004; Gollisch and Meister, 2010; Baden et al., 2020; Kerschensteiner, 2022). Complex interactions involving different sets of interneurons in the two plexiform layers of the retina contribute to this parallel feature extraction (reviewed in Diamond, 2017; Franke and Baden, 2017). Our study highlights that a mechanism at the last step of retinal processing—intrinsic properties of the RGCs—can strongly participate in forming their responses. This is in line with several previous studies showing that spike generation in RGCs is rather diverse in its location and dynamics (see also above): For example, it was reported that Na+-based spikes can originate from both the dendrites and the soma region of RGCs (Velte and Masland, 1999; Oesch et al., 2005) and dendritic spikes play an important role in sharpening the directional tuning of DS RGCs (Velte and Masland, 1999; Oesch and Taylor, 2010; Schachter et al., 2010). Moreover, RGCs display very different levels of spike frequency adaptation in response to current injection (O’Brien et al., 2002). Depending on the somatic or dendritic location of the Nav, the integration properties of RGCs vary, and the cell’s sensitivity can be optimized and tuned for the preferred visual stimulus (Brombas et al., 2022). These studies, together with the present one, support the notion that RGCs are not the passive integrators of upstream signals—as they are still depicted in textbooks—but play an active role in retinal processing.
Our voltage-clamp experiments suggest that the slow VGSC desensitisation in tOn-small RGCs is one of the mechanisms shaping the cells’ distinct light responses. We cannot answer the question, though, where the relevant VGSCs are located: at the soma or along the dendrites. Yet, the subcellular VGSCs distribution pattern shapes how a cell processes synaptic input (Ran et al., 2020; Brombas et al., 2022). The relationship between VGSC channel distribution and kinetics for the overall cell’s light responses needs to be addressed in future experiments. In addition to VGSCs, calcium-gated potassium channels are thought to cause hyperpolarized potentials after depolarization (Lancaster and Adams, 1986; Sah and McLachlan, 1991), and may therefore also contribute to the differential desensitization observed in the two RGC types. Future work is needed to systematically investigate this possibility.
Still, the question arises what the advantage of feature extraction at the very last level of retinal processing may be. Considering the “feature” that is extracted by the tOn-small RGC—selective high contrast signals—one may speculate that this RGC type exclusively relays very strong and reliable events to higher visual areas. Extracting this feature late in the retinal network may bear the advantage that signal-to-noise is at its maximum at the RGC level, providing a suitable high-quality signal to cut off low-contrast responses. Additionally, a feature extraction mechanism intrinsic to RGCs would be largely independent from the highly dynamic inner retinal circuits and could provide “stability” to the output signal. In light of evidence suggesting that the retinal code changes over the full range of illumination conditions (Tikidji-Hamburyan et al., 2014), a spike generator with properties like the one in tOn-small RGCs may provide a “stable” element for encoding high contrast signals across different visual environments.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by the Institutional Animal Welfare Committee of the University of Tübingen and the MPI for Brain Research, Frankfurt/M. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
LC: Conceptualization, Data curation, Formal analysis, Methodology, Visualization, Writing – original draft. YR: Data curation, Formal analysis, Writing – review and editing. MY: Data curation, Formal analysis, Visualization, Writing – review and editing. OA: Conceptualization, Data curation, Visualization, Writing – review and editing, Formal analysis, Investigation. EB: Data curation, Formal analysis, Writing – review and editing. LH: Data curation, Formal analysis, Writing – review and editing. SH: Conceptualization, Data curation, Formal analysis, Funding acquisition, Project administration, Supervision, Visualization, Writing – original draft. TE: Conceptualization, Formal analysis, Funding acquisition, Project administration, Supervision, Validation, Visualization, Writing – original draft. TS: Conceptualization, Data curation, Visualization, Writing – original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by Deutsche Forschungsgemeinschaft (DFG) (HA 5277/3-1, EU 42/4-1 and FOR 701).
Acknowledgments
We thank Tom Baden, Philipp Berens, and Katrin Franke for helpful discussions, Gordon Eske for technical assistance, and Ziwei Huang for extracting cell morphologies. We acknowledge support from the Open Access Publication Fund of the University of Tübingen.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2023.1337768/full#supplementary-material
Footnotes
References
Abdelsayed, M., Peters, C., and Ruben, P. (2015). Differential thermosensitivity in mixed syndrome cardiac sodium channel mutants. J. Physiol. 593, 4201–4223. doi: 10.1113/JP270139
Anderson, J., Lampl, I., Gillespie, D., and Ferster, D. (2000). The contribution of noise to contrast invariance of orientation tuning in cat visual cortex. Science 290, 1968–1972. doi: 10.1126/science.290.5498.1968
Asari, H., and Meister, M. (2012). Divergence of visual channels in the inner retina. Nat. Neurosci. 15, 1581–1589. doi: 10.1038/nn.3241
Awatramani, G., and Slaughter, M. (2000). Origin of transient and sustained responses in ganglion cells of the retina. J. Neurosci. 20, 7087–7095. doi: 10.1523/jneurosci.20-18-07087.2000
Baccus, S., and Meister, M. (2002). Fast and slow contrast adaptation in retinal circuitry. Neuron 36, 909–919.
Baden, T., Berens, P., Bethge, M., and Euler, T. (2013). Spikes in mammalian bipolar cells support temporal layering of the inner retina. Curr. Biol. 23, 48–52. doi: 10.1016/j.cub.2012.11.006
Baden, T., Berens, P., Franke, K., Román Rosón, M., Bethge, M., and Euler, T. (2016). The functional diversity of retinal ganglion cells in the mouse. Nature 529, 345–350. doi: 10.1038/nature16468
Baden, T., Euler, T., and Berens, P. (2020). Understanding the retinal basis of vision across species. Nat. Rev. Neurosci. 21, 5–20. doi: 10.1038/s41583-019-0242-1
Bae, J., Mu, S., Kim, J., Turner, N., Tartavull, I., Kemnitz, N., et al. (2018). Digital museum of retinal ganglion cells with dense anatomy and physiology. Cell 173, 1293–1306. doi: 10.1016/j.cell.2018.04.040
Barlow, H., Hill, R., and Levick, W. (1964). Retinal ganglion cells responding selectively to direction and speed of image motion in the rabbit. J. Physiol. 173, 377–407. doi: 10.1113/jphysiol.1964.sp007463
Bleckert, A., Schwartz, G., Turner, M., Rieke, F., and Wong, R. (2014). Visual space is represented by nonmatching topographies of distinct mouse retinal ganglion cell types. Curr. Biol. 24, 310–315. doi: 10.1016/j.cub.2013.12.020
Borst, A., and Euler, T. (2011). Seeing things in motion: models, circuits, and mechanisms. Neuron 71, 974–994. doi: 10.1016/j.neuron.2011.08.031
Branco, T., and Häusser, M. (2010). The single dendritic branch as a fundamental functional unit in the nervous system. Curr. Opin. Neurobiol. 20, 494–502. doi: 10.1016/j.conb.2010.07.009
Brombas, A., Zhou, X., and Williams, S. (2022). Light-evoked dendritic spikes in sustained but not transient rabbit retinal ganglion cells. Neuron 110, 2802–2814. doi: 10.1016/j.neuron.2022.06.009
Carter, B., and Bean, B. (2011). Incomplete inactivation and rapid recovery of voltage-dependent sodium channels during high-frequency firing in cerebellar Purkinje neurons. J. Neurophysiol. 105, 860–871. doi: 10.1152/jn.01056.2010
Chander, D., and Chichilnisky, E. (2001). Adaptation to temporal contrast in primate and salamander retina. J. Neurosci. 21, 9904–9916.
Chang, L., Breuninger, T., and Euler, T. (2013). Chromatic coding from cone-type unselective circuits in the mouse retina. Neuron 77, 559–571. doi: 10.1016/j.neuron.2012.12.012
Chichilnisky, E. (2001). A simple white noise analysis of neuronal light responses. Network 12, 199–213.
Coombs, J., van der List, D., Wang, G., and Chalupa, L. (2006). Morphological properties of mouse retinal ganglion cells. Neuroscience 140, 123–136. doi: 10.1016/j.neuroscience.2006.02.079
Cui, Y., Wang, Y., Park, S., Demb, J., and Butts, D. (2016). Divisive suppression explains high-precision firing and contrast adaptation in retinal ganglion cells. bioRxiv [Preprint]. doi: 10.1101/064592
Cummins, T., and Waxman, S. (1997). Downregulation of tetrodotoxin-resistant sodium currents and upregulation of a rapidly repriming tetrodotoxin-sensitive sodium current in small spinal sensory neurons after nerve injury. J. Neurosci. 17, 3503–3514.
Davis, S., and Linn, C. (2003). Activation of NMDA receptors linked to modulation of voltage-gated ion channels and functional implications. Am. J. Physiol. Cell Physiol. 284, C757–C768. doi: 10.1152/ajpcell.00252.2002
Dayan, P., and Abbott, L. (2005). Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems. Cambridge, MA: MIT Press.
DeVries, S. (2000). Bipolar cells use kainate and AMPA receptors to filter visual information into separate channels. Neuron 28, 847–856. doi: 10.1016/s0896-6273(00)00158-6
Diamond, J. (2017). Inhibitory Interneurons in the Retina: Types, Circuitry, and Function. Annu. Rev. Vis. Sci. 3, 1–24. doi: 10.1146/annurev-vision-102016-061345
Drinnenberg, A., Franke, F., Morikawa, R., Jüttner, J., Hillier, D., Hantz, P., et al. (2018). How diverse retinal functions arise from feedback at the first visual synapse. Neuron 99, 117–134. doi: 10.1016/j.neuron.2018.06.001
Euler, T., Hausselt, S., Margolis, D., Breuninger, T., Castell, X., Detwiler, P., et al. (2009). Eyecup scope–optical recordings of light stimulus-evoked fluorescence signals in the retina. Pflugers Arch. 457, 1393–1414. doi: 10.1007/s00424-008-0603-5
Farrow, K., Teixeira, M., Szikra, T., Viney, T., Balint, K., Yonehara, K., et al. (2013). Ambient illumination toggles a neuronal circuit switch in the retina and visual perception at cone threshold. Neuron 78, 325–338. doi: 10.1016/j.neuron.2013.02.014
Field, G., Gauthier, J., Sher, A., Greschner, M., Machado, T., Jepson, L., et al. (2010). Functional connectivity in the retina at the resolution of photoreceptors. Nature 467, 673–677. doi: 10.1038/nature09424
Fjell, J., Dib-Hajj, S., Fried, K., Black, J., and Waxman, S. (1997). Differential expression of sodium channel genes in retinal ganglion cells. Brain Res. Mol. Brain Res. 50, 197–204.
Fohlmeister, J., and Miller, R. (1997). Impulse encoding mechanisms of ganglion cells in the tiger salamander retina. J. Neurophysiol. 78, 1935–1947. doi: 10.1152/jn.1997.78.4.1935
Franke, K., and Baden, T. (2017). General features of inhibition in the inner retina. J. Physiol. 595, 5507–5515. doi: 10.1113/JP273648
Franke, K., Berens, P., Schubert, T., Bethge, M., Euler, T., and Baden, T. (2017). Inhibition decorrelates visual feature representations in the inner retina. Nature 542, 439–444. doi: 10.1038/nature21394
Fried, S., Lasker, A., Desai, N., Eddington, D., and Rizzo, J. (2009). Axonal sodium-channel bands shape the response to electric stimulation in retinal ganglion cells. J. Neurophysiol. 101, 1972–1987. doi: 10.1152/jn.91081.2008
George, J., Baden, D., Gerwick, W., and Murray, T. (2012). Bidirectional influence of sodium channel activation on NMDA receptor–dependent cerebrocortical neuron structural plasticity. Proc. Natl. Acad. Sci. U. S. A. 109, 19840–19845. doi: 10.1073/pnas.1212584109
Goetz, J., Jessen, Z., Jacobi, A., Mani, A., Cooler, S., Greer, D., et al. (2021). Unified classification of mouse retinal ganglion cells using function, morphology, and gene expression. bioRxiv [Preprint]. doi: 10.1101/2021.06.10.447922
Goldin, A. (1999). Diversity of mammalian voltage-gated sodium channels. Ann. N. Y. Acad. Sci. 868, 38–50.
Gollisch, T., and Meister, M. (2010). Eye smarter than scientists believed: neural computations in circuits of the retina. Neuron 65, 150–164.
Gong, S., Zheng, C., Doughty, M., Losos, K., Didkovsky, N., Schambra, U., et al. (2003). A gene expression atlas of the central nervous system based on bacterial artificial chromosomes. Nature 425, 917–925. doi: 10.1038/nature02033
Hahin, R., and Campbell, D. (1983). Simple shifts in the voltage dependence of sodium channel gating caused by divalent cations. J. Gen. Physiol. 82, 785–805. doi: 10.1085/jgp.82.6.785
Hahn, J., Monavarfeshani, A., Qiao, M., Kao, A., Kölsch, Y., Kumar, A., et al. (2023). Evolution of neuronal cell classes and types in the vertebrate retina. bioRxiv [Preprint]. doi: 10.1101/2023.04.07.536039
Helmstaedter, M., Briggman, K., Turaga, S., Jain, V., Seung, H., and Denk, W. (2013). Connectomic reconstruction of the inner plexiform layer in the mouse retina. Nature 500, 168–174. doi: 10.1038/nature12346
Hodgkin, A., and Huxley, A. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117, 500–544. doi: 10.1113/jphysiol.1952.sp004764
Huang, S., Hong, S., and De Schutter, E. (2015). Non-linear leak currents affect mammalian neuron physiology. Front. Cell Neurosci. 9:432. doi: 10.3389/fncel.2015.00432
Jacoby, J., Zhu, Y., DeVries, S., and Schwartz, G. (2015). An Amacrine Cell Circuit for Signaling Steady Illumination in the Retina. Cell Rep. 13, 2663–2670. doi: 10.1016/j.celrep.2015.11.062
Josephson, I. (1996). Depolarization shifts the voltage dependence of cardiac sodium channel and calcium channel gating charge movements. Pflugers Arch. 431, 895–904. doi: 10.1007/s004240050083
Kanjhan, R., and Sivyer, B. (2010). Two types of ON direction-selective ganglion cells in rabbit retina. Neurosci. Lett. 483, 105–109. doi: 10.1016/j.neulet.2010.07.071
Kastner, D., and Baccus, S. (2011). Coordinated dynamic encoding in the retina using opposing forms of plasticity. Nat. Neurosci. 14, 1317–1322.
Kerschensteiner, D. (2022). Feature detection by retinal ganglion cells. Annu. Rev. Vis. Sci. 8:112009. doi: 10.1146/annurev-vision-100419-112009
Kim, K., and Rieke, F. (2001). Temporal contrast adaptation in the input and output signals of salamander retinal ganglion cells. J. Neurosci. 21, 287–299.
Kim, K., and Rieke, F. (2003). Slow Na inactivation and variance adaptation in salamander retinal ganglion cells. J. Neurosci. 23, 1506–1516. doi: 10.1523/jneurosci.23-04-01506.2003
Kong, J., Fish, D., Rockhill, R., and Masland, R. (2005). Diversity of ganglion cells in the mouse retina: unsupervised morphological classification and its limits. J. Comp. Neurol. 489, 293–310. doi: 10.1002/cne.20631
Krieger, B., Qiao, M., Rousso, D., Sanes, J., and Meister, M. (2017). Four alpha ganglion cell types in mouse retina: Function, structure, and molecular signatures. PLoS One 12:e0180091. doi: 10.1371/journal.pone.0180091
Lancaster, B., and Adams, P. (1986). Calcium-dependent current generating the afterhyperpolarization of hippocampal neurons. J. Neurophysiol. 55, 1268–1282. doi: 10.1152/jn.1986.55.6.1268
Levick, W. (1967). Receptive fields and trigger features of ganglion cells in the visual streak of the rabbit’s retina. J. Physiol. 188, 285–307.
Li, M., West, J., Lai, Y., Scheuer, T., and Catterall, W. (1992). Functional modulation of brain sodium channels by cAMP-dependent phosphorylation. Neuron 8, 1151–1159. doi: 10.1016/0896-6273(92)90135-z
Li, W., Haji Ghaffari, D., Misra, R., and Weiland, J. (2022). Retinal ganglion cell desensitization is mitigated by varying parameter constant excitation pulse trains. Front. Cell Neurosci. 16:897146. doi: 10.3389/fncel.2022.897146
Lin, B., and Masland, R. (2006). Populations of wide-field amacrine cells in the mouse retina. J. Comp. Neurol. 499, 797–809. doi: 10.1002/cne.21126
Majumdar, S., Wässle, H., Jusuf, P., and Haverkamp, S. (2008). Mirror-symmetrical populations of wide-field amacrine cells of the macaque monkey retina. J. Comp. Neurol. 508, 13–27. doi: 10.1002/cne.21666
Mani, A., and Schwartz, G. (2017). Circuit mechanisms of a retinal ganglion cell with stimulus-dependent response latency and activation beyond its dendrites. Curr. Biol. 27, 471–482. doi: 10.1016/j.cub.2016.12.033
Martiszus, B., Tsintsadze, T., Chang, W., and Smith, S. (2021). Enhanced excitability of cortical neurons in low-divalent solutions is primarily mediated by altered voltage-dependence of voltage-gated sodium channels. eLife 10, 1–26. doi: 10.7554/eLife.67914
Mauss, A., Vlasits, A., Borst, A., and Feller, M. (2017). Visual circuits for direction selectivity. Annu. Rev. Neurosci. 40, 211–230. doi: 10.1146/annurev-neuro-072116-031335
Mobbs, P., Everett, K., and Cook, A. (1992). Signal shaping by voltage-gated currents in retinal ganglion cells. Brain Res. 574, 217–223. doi: 10.1016/0006-8993(92)90819-u
Münch, T., da Silveira, R., Siegert, S., Viney, T., Awatramani, G., and Roska, B. (2009). Approach sensitivity in the retina processed by a multifunctional neural circuit. Nat. Neurosci. 12, 1308–1316. doi: 10.1038/nn.2389
Naka, K., and Rushton, W. (1966). S-potentials from luminosity units in the retina of fish (Cyprinidae). J. Physiol. 185, 587–599.
Neumann, S., Hüser, L., Ondreka, K., Auler, N., and Haverkamp, S. (2016). Cell type-specific bipolar cell input to ganglion cells in the mouse retina. Neuroscience 316, 420–432. doi: 10.1016/j.neuroscience.2015.12.051
Nikolaev, A., Leung, K., Odermatt, B., and Lagnado, L. (2013). Synaptic mechanisms of adaptation and sensitization in the retina. Nat. Neurosci. 16, 934–941. doi: 10.1038/nn.3408
Nirenberg, S., and Meister, M. (1997). The light response of retinal ganglion cells is truncated by a displaced amacrine circuit. Neuron 18, 637–650. doi: 10.1016/s0896-6273(00)80304-9
Numann, R., Catterall, W., and Scheuer, T. (1991). Functional modulation of brain sodium channels by protein kinase C phosphorylation. Science 254, 115–118. doi: 10.1126/science.1656525
O’Brien, B., Caldwell, J., Ehring, G., Bumsted, O., Brien, K., Luo, S., et al. (2008). Tetrodotoxin-resistant voltage-gated sodium channels Na(v)1.8 and Na(v)1.9 are expressed in the retina. J. Comp. Neurol. 508, 940–951. doi: 10.1002/cne.21701
O’Brien, B., Isayama, T., Richardson, R., and Berson, D. (2002). Intrinsic physiological properties of cat retinal ganglion cells. J. Physiol. 538, 787–802. doi: 10.1113/jphysiol.2001.013009
Oesch, N., Euler, T., and Taylor, W. (2005). Direction-selective dendritic action potentials in rabbit retina. Neuron 47, 739–750. doi: 10.1016/j.neuron.2005.06.036
Oesch, N., and Taylor, W. (2010). Tetrodotoxin-resistant sodium channels contribute to directional responses in starburst amacrine cells. PLoS One 5:e12447. doi: 10.1371/journal.pone.0012447
Pang, J., Gao, F., and Wu, S. (2003). Light-evoked excitatory and inhibitory synaptic inputs to ON and OFF alpha ganglion cells in the mouse retina. J. Neurosci. 23, 6063–6073. doi: 10.1523/jneurosci.23-14-06063.2003
Park, S., Pottackal, J., Bin Ke, J., Jun, N., Rahmani, P., Kim, I., et al. (2018). Convergence and divergence of CRH amacrine cells in mouse retinal circuitry. J. Neurosci. 38, 3753–3766. doi: 10.1523/JNEUROSCI.2518-17.2018
Peichl, L., Ott, H., and Boycott, B. (1987). Alpha ganglion cells in mammalian retinae. Proc. R. Soc. Lond. B Biol. Sci. 231, 169–197. doi: 10.1098/rspb.1987.0040
Poleg-Polsky, A., and Diamond, J. (2011). Imperfect space clamp permits electrotonic interactions between inhibitory and excitatory synaptic conductances, distorting voltage clamp recordings. PLoS One 6:e19463. doi: 10.1371/journal.pone.0019463
Raghuram, V., Werginz, P., and Fried, S. (2019). Scaling of the AIS and somatodendritic compartments in α S RGCs. Front. Cell Neurosci. 13:436. doi: 10.3389/fncel.2019.00436
Ran, Y., Huang, Z., Baden, T., Schubert, T., Baayen, H., Berens, P., et al. (2020). Type-specific dendritic integration in mouse retinal ganglion cells. Nat. Commun. 11:2101. doi: 10.1038/s41467-020-15867-9
Renganathan, M., Cummins, T., Hormuzdiar, W., and Waxman, S. (2000). Alpha-SNS produces the slow TTX-resistant sodium current in large cutaneous afferent DRG neurons. J. Neurophysiol. 84, 710–718.
Rheaume, B., Jereen, A., Bolisetty, M., Sajid, M., Yang, Y., Renna, K., et al. (2018). Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes. Nat. Commun. 9:2759. doi: 10.1038/s41467-018-05134-3
Roska, B., and Werblin, F. (2001). Vertical interactions across ten parallel, stacked representations in the mammalian retina. Nature 410, 583–587. doi: 10.1038/35069068
Rousso, D., Qiao, M., Kagan, R., Yamagata, M., Palmiter, R., and Sanes, J. (2016). Two Pairs of ON and OFF retinal ganglion cells are defined by intersectional patterns of transcription factor expression. Cell Rep. 15, 1930–1944. doi: 10.1016/j.celrep.2016.04.069
Sah, P., and McLachlan, E. (1991). Ca(2+)-activated K+ currents underlying the afterhyperpolarization in guinea pig vagal neurons: a role for Ca(2+)-activated Ca2+ release. Neuron 7, 257–264. doi: 10.1016/0896-6273(91)90264-z
Sawant, A., Ebbinghaus, B., Bleckert, A., Gamlin, C., Yu, W., Berson, D., et al. (2021). Organization and emergence of a mixed GABA-glycine retinal circuit that provides inhibition to mouse ON-sustained alpha retinal ganglion cells. Cell Rep. 34:108858. doi: 10.1016/j.celrep.2021.108858
Schachter, M., Oesch, N., Smith, R., and Taylor, W. (2010). Dendritic spikes amplify the synaptic signal to enhance detection of motion in a simulation of the direction-selective ganglion cell. PLoS Comput. Biol. 6:1000899. doi: 10.1371/journal.pcbi.1000899
Scheuer, T., and Catterall, W. (2006). Control of neuronal excitability by phosphorylation and dephosphorylation of sodium channels. Biochem. Soc. Trans. 34, 1299–1302. doi: 10.1042/BST0341299
Schubert, T., Kerschensteiner, D., Eggers, Misgeld, T., Kerschensteiner, M., Lichtman, J., et al. (2008). Development of presynaptic inhibition onto retinal bipolar cell axon terminals is subclass-specific. J. Neurophysiol. 100, 304–316. doi: 10.1152/jn.90202.2008
Schwartz, G., Okawa, H., Dunn, F., Morgan, J., Kerschensteiner, D., Wong, R., et al. (2012). The spatial structure of a nonlinear receptive field. Nat. Neurosci. 15, 1572–1580. doi: 10.1038/nn.3225
Siegert, S., Cabuy, E., Scherf, B., Kohler, H., Panda, S., Le, Y., et al. (2012). Transcriptional code and disease map for adult retinal cell types. Nat. Neurosci. 15, S1–S2. doi: 10.1038/nn.3032
Siegert, S., Scherf, B., Del Punta, K., Didkovsky, N., Heintz, N., and Roska, B. (2009). Genetic address book for retinal cell types. Nat. Neurosci. 12, 1197–1204. doi: 10.1038/nn.2370
Silva, J. (2014). Slow inactivation of Na(+) channels. Handb. Exp. Pharmacol. 221, 33–49. doi: 10.1007/978-3-642-41588-3_3
Smith, B., McHugh, C., Hirano, A., Brecha, N., and Barnes, S. (2023). Transient and sustained ganglion cell light responses are differentially modulated by intrinsically produced reactive oxygen species acting upon specific voltage-gated Na+ channel isoforms. J. Neurosci. 43, 2291–2304. doi: 10.1523/JNEUROSCI.1723-22.2023
Stuart, G., and Spruston, N. (2015). Dendritic integration: 60 years of progress. Nat. Neurosci. 18, 1713–1721. doi: 10.1038/nn.4157
Sümbül, U., Song, S., McCulloch, K., Becker, M., Lin, B., Sanes, J., et al. (2014). A genetic and computational approach to structurally classify neuronal types. Nat. Commun. 5:3512. doi: 10.1038/ncomms4512
Sun, W., Li, N., and He, S. (2002). Large-scale morphological survey of mouse retinal ganglion cells. J. Comp. Neurol.451, 115–126. doi: 10.1002/cne.10323
Tien, N., Kim, T., and Kerschensteiner, D. (2016). Target-Specific Glycinergic Transmission from VGluT3-Expressing Amacrine Cells Shapes Suppressive Contrast Responses in the Retina. Cell Rep. 15, 1369–1375. doi: 10.1016/j.celrep.2016.04.025
Tikidji-Hamburyan, A., Reinhard, K., Seitter, H., Hovhannisyan, A., Procyk, C., Allen, A., et al. (2014). Retinal output changes qualitatively with every change in ambient illuminance. Nat. Neurosci. doi: 10.1038/nn.3891 [Epub ahead of print].
Tran, N., Shekhar, K., Whitney, I., Jacobi, A., Benhar, I., Hong, G., et al. (2019). Single-cell profiles of retinal ganglion cells differing in resilience to injury reveal neuroprotective genes. Neuron 104, 1039–1055. doi: 10.1016/j.neuron.2019.11.006
Tran-Van-Minh, A., Cazé, R., Abrahamsson, T., Cathala, L., Gutkin, B., and DiGregorio, D. (2015). Contribution of sublinear and supralinear dendritic integration to neuronal computations. Front. Cell Neurosci. 9:67. doi: 10.3389/fncel.2015.00067
Tsai, D., Morley, J., Suaning, G., and Lovell, N. (2011). Sodium channel inactivation reduces retinal ganglion cell responsiveness to repetitive prosthetic stimulation. 2011 5th International IEEE/EMBS Conference on Neural Engineering. NER 2011, 550–553. doi: 10.1109/NER.2011.5910607
Turner, M., and Rieke, F. (2016). Synaptic Rectification Controls Nonlinear Spatial Integration of Natural Visual Inputs. Neuron 90, 1257–1271. doi: 10.1016/j.neuron.2016.05.006
Van Hook, M., Nawy, S., and Thoreson, W. (2019). Voltage- and calcium-gated ion channels of neurons in the vertebrate retina. Prog. Retinal Eye Res. 72:100760. doi: 10.1016/j.preteyeres.2019.05.001
Van Wart, A., Trimmer, J., and Matthews, G. (2007). Polarized distribution of ion channels within microdomains of the axon initial segment. J. Compar. Neurol. 500, 339–352. doi: 10.1002/cne.21173
van Wyk, M., Taylor, W., and Vaney, D. (2006). Local edge detectors: a substrate for fine spatial vision at low temporal frequencies in rabbit retina. J. Neurosci. 26, 13250–13263. doi: 10.1523/JNEUROSCI.1991-06.2006
van Wyk, M., Wässle, H., and Taylor, W. (2009). Receptive field properties of ON- and OFF-ganglion cells in the mouse retina. Vis. Neurosci. 26, 297–308. doi: 10.1017/S0952523809990137
Velte, T., and Masland, R. (1999). Action potentials in the dendrites of retinal ganglion cells. J. Neurophysiol. 81, 1412–1417. doi: 10.1152/jn.1999.81.3.1412
Wang, Y., Weick, M., and Demb, J. (2011). Spectral and temporal sensitivity of cone-mediated responses in mouse retinal ganglion cells. J. Neurosci. 31, 7670–7681. doi: 10.1523/JNEUROSCI.0629-11.2011
Weick, M., and Demb, J. (2011). Delayed-Rectifier K Channels Contribute to Contrast Adaptation in Mammalian Retinal Ganglion Cells. Neuron 71, 166–179. doi: 10.1016/j.neuron.2011.04.033
Werginz, P., Raghuram, V., and Fried, S. (2020). Tailoring of the axon initial segment shapes the conversion of synaptic inputs into spiking output in OFF-α T retinal ganglion cells. Sci. Adv. 6, eabb6642. doi: 10.1126/sciadv.abb6642
Wienbar, S., and Schwartz, G. (2022). Differences in spike generation instead of synaptic inputs determine the feature selectivity of two retinal cell types. Neuron 110, 2110–2123. doi: 10.1016/j.neuron.2022.04.012
Yao, X., Cafaro, J., McLaughlin, A., Postma, F., Paul, D., Awatramani, G., et al. (2018). Gap junctions contribute to differential light adaptation across direction-selective retinal ganglion cells. Neuron 100, 216–228. doi: 10.1016/j.neuron.2018.08.021
Yin, L., Smith, R., Sterling, P., and Brainard, D. (2009). Physiology and morphology of color-opponent ganglion cells in a retina expressing a dual gradient of S and M opsins. J. Neurosci. 29, 2706–2724.
Yu, W., El-Danaf, R., Okawa, H., Pacholec, J., Matti, U., Schwarz, K., et al. (2018). Synaptic convergence patterns onto retinal ganglion cells are preserved despite topographic variation in pre- and postsynaptic territories. Cell Rep. 25, 2017–2026. doi: 10.1016/j.celrep.2018.10.089
Yu, X. (2006). The role of intracellular sodium in the regulation of NMDA-receptor-mediated channel activity and toxicity. Mol. Neurobiol. 33, 63–80.
Zhang, C., and McCall, M. (2012). Receptor targets of amacrine cells. Vis. Neurosci. 29, 11–29. doi: 10.1017/S0952523812000028
Keywords: retina, action potential, sodium channel, light response, retinal ganglion cell, spike generator, mouse
Citation: Chang L, Ran Y, Yang M, Auferkorte O, Butz E, Hüser L, Haverkamp S, Euler T and Schubert T (2024) Spike desensitisation as a mechanism for high-contrast selectivity in retinal ganglion cells. Front. Cell. Neurosci. 17:1337768. doi: 10.3389/fncel.2023.1337768
Received: 13 November 2023; Accepted: 19 December 2023;
Published: 10 January 2024.
Edited by:
Richard Anthony DeFazio, University of Michigan, United StatesReviewed by:
Caroline Adams, University of Pennsylvania, United StatesTomomi Ichinose, Wayne State University, United States
Copyright © 2024 Chang, Ran, Yang, Auferkorte, Butz, Hüser, Haverkamp, Euler and Schubert. 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: Thomas Euler, thomas.euler@cin.uni-tuebingen.de; Timm Schubert, timm.schubert@cin.uni-tuebingen.de