Abstract
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.
TABLE 1
| Parameters | Values |
| Temperature | T = 32°C |
| Intracellular axial resistivity | Ra = 110 Ωcm |
| Specific membrane resistance | Rm = 15,000 Ωcm2 |
| Specific membrane capacitance | Cm = 1 μF cm–2 |
Model parameters.
TABLE 2
| Channel type | Conductance in soma (S cm–2) | Conductance in dendrites (S cm–2) |
| gNa | 0.08 | 0.025 |
| gCa | 0 | 0 |
| gK | 0 | 0 |
| gK,A | 0 | 0 |
| gK,Ca | 0 | 0 |
Reference distribution of ion channels in cellular compartments.
TABLE 3
| Branch order | Branch length in sOn-α cells (μ m) | Branch diameter in sOn-α cells (μ m) | Branch length in tOn-small cells (μ m) | Branch diameter in tOn-small cells (μ m) |
| Branch 01 | 14.91 | 0.72 | 11.83 | 0.7 |
| Branch 02 | 22.65 | 0.52 | 7.47 | 0.59 |
| Branch 03 | 36.84 | 0.5 | 5.03 | 0.46 |
| Branch 04 | 30.98 | 0.43 | 5.84 | 0.39 |
| Branch 05 | 40.81 | 0.44 | 5.04 | 0.36 |
| Branch 06 | 32.22 | 0.42 | 4.82 | 0.34 |
| Branch 07 | 53.79 | 0.43 | 6.8 | 0.36 |
| Branch 08 | 38.92 | 0.43 | 5.95 | 0.36 |
| Branch 09 | 18.24 | 0.36 | 5.2 | 0.36 |
| Branch 10 | 5.75 | 0.34 | ||
| Branch 11 | 5.12 | 0.31 | ||
| Branch 12 | 5.62 | 0.34 | ||
| Branch 13 | 5.6 | 0.31 | ||
| Branch 14 | 7.09 | 0.33 | ||
| Branch 15 | 6.62 | 0.33 | ||
| Branch 16 | 7.51 | 0.32 | ||
| Branch 17 | 4.13 | 0.33 | ||
| Branch 18 | 7.51 | 0.28 | ||
| Branch 19 | 7.73 | 0.32 | ||
| Branch 20 | 5.71 | 0.25 |
Dendrite morphology parameters.
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).
FIGURE 1
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.
TABLE 4
| RGC type in lgfbp5 line | IPL depth | RGC types in recent studies | |||||||
| Bae et al., 2018 | Kong et al., 2005 | Farrow et al., 2013 | Sümbül et al., 2014 | Sun et al., 2002 | Goetz et al., 2021 | Baden et al., 2016 | Tran et al., 2019 | ||
| Off (medium soma, bright) | S2 | 4i or 4on | cluster 3 | PV4 | CB2 | RGB1 or outer RGB3 | OFF transient small RF | G9-OFF T alpha mini* | Tbr1-S2, cluster 21 |
| On (large soma, bright) | S3 | 6sn | cluster 2 | PV2 | cluster X | inner RGB3 | ON transient small RF | G18a-ON transient or G23-ON alpha mini* | – |
| Trans. Off -α (largest soma, dim) | S2 | 4ow | – | PV5 | W7a | A2 outer | OFF transient alpha | G8-OFF T alpha | Alpha OFF T, cluster 45 |
Igfbp5-positive RGCs and their presumed morphological and genetic counterparts in earlier mouse studies.
Missing entries indicate that type assignment is unclear.
*As discussed in Tran et al. (2019).
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).
FIGURE 2
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).
FIGURE 3
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).
FIGURE 4
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.
FIGURE 5
FIGURE 6
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.
FIGURE 7
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.
Statements
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
References
1
AbdelsayedM.PetersC.RubenP. (2015). Differential thermosensitivity in mixed syndrome cardiac sodium channel mutants.J. Physiol.5934201–4223. 10.1113/JP270139
2
AndersonJ.LamplI.GillespieD.FersterD. (2000). The contribution of noise to contrast invariance of orientation tuning in cat visual cortex.Science2901968–1972. 10.1126/science.290.5498.1968
3
AsariH.MeisterM. (2012). Divergence of visual channels in the inner retina.Nat. Neurosci.151581–1589. 10.1038/nn.3241
4
AwatramaniG.SlaughterM. (2000). Origin of transient and sustained responses in ganglion cells of the retina.J. Neurosci.207087–7095. 10.1523/jneurosci.20-18-07087.2000
5
BaccusS.MeisterM. (2002). Fast and slow contrast adaptation in retinal circuitry.Neuron36909–919.
6
BadenT.BerensP.BethgeM.EulerT. (2013). Spikes in mammalian bipolar cells support temporal layering of the inner retina.Curr. Biol.2348–52. 10.1016/j.cub.2012.11.006
7
BadenT.BerensP.FrankeK.Román RosónM.BethgeM.EulerT. (2016). The functional diversity of retinal ganglion cells in the mouse.Nature529345–350. 10.1038/nature16468
8
BadenT.EulerT.BerensP. (2020). Understanding the retinal basis of vision across species.Nat. Rev. Neurosci.215–20. 10.1038/s41583-019-0242-1
9
BaeJ.MuS.KimJ.TurnerN.TartavullI.KemnitzN.et al (2018). Digital museum of retinal ganglion cells with dense anatomy and physiology.Cell1731293–1306. 10.1016/j.cell.2018.04.040
10
BarlowH.HillR.LevickW. (1964). Retinal ganglion cells responding selectively to direction and speed of image motion in the rabbit.J. Physiol.173377–407. 10.1113/jphysiol.1964.sp007463
11
BleckertA.SchwartzG.TurnerM.RiekeF.WongR. (2014). Visual space is represented by nonmatching topographies of distinct mouse retinal ganglion cell types.Curr. Biol.24310–315. 10.1016/j.cub.2013.12.020
12
BorstA.EulerT. (2011). Seeing things in motion: models, circuits, and mechanisms.Neuron71974–994. 10.1016/j.neuron.2011.08.031
13
BrancoT.HäusserM. (2010). The single dendritic branch as a fundamental functional unit in the nervous system.Curr. Opin. Neurobiol.20494–502. 10.1016/j.conb.2010.07.009
14
BrombasA.ZhouX.WilliamsS. (2022). Light-evoked dendritic spikes in sustained but not transient rabbit retinal ganglion cells.Neuron1102802–2814. 10.1016/j.neuron.2022.06.009
15
CarterB.BeanB. (2011). Incomplete inactivation and rapid recovery of voltage-dependent sodium channels during high-frequency firing in cerebellar Purkinje neurons.J. Neurophysiol.105860–871. 10.1152/jn.01056.2010
16
ChanderD.ChichilniskyE. (2001). Adaptation to temporal contrast in primate and salamander retina.J. Neurosci.219904–9916.
17
ChangL.BreuningerT.EulerT. (2013). Chromatic coding from cone-type unselective circuits in the mouse retina.Neuron77559–571. 10.1016/j.neuron.2012.12.012
18
ChichilniskyE. (2001). A simple white noise analysis of neuronal light responses.Network12199–213.
19
CoombsJ.van der ListD.WangG.ChalupaL. (2006). Morphological properties of mouse retinal ganglion cells.Neuroscience140123–136. 10.1016/j.neuroscience.2006.02.079
20
CuiY.WangY.ParkS.DembJ.ButtsD. (2016). Divisive suppression explains high-precision firing and contrast adaptation in retinal ganglion cells.bioRxiv[Preprint]. 10.1101/064592
21
CumminsT.WaxmanS. (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.173503–3514.
22
DavisS.LinnC. (2003). Activation of NMDA receptors linked to modulation of voltage-gated ion channels and functional implications.Am. J. Physiol. Cell Physiol.284C757–C768. 10.1152/ajpcell.00252.2002
23
DayanP.AbbottL. (2005). Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems.Cambridge, MA: MIT Press.
24
DeVriesS. (2000). Bipolar cells use kainate and AMPA receptors to filter visual information into separate channels.Neuron28847–856. 10.1016/s0896-6273(00)00158-6
25
DiamondJ. (2017). Inhibitory Interneurons in the Retina: Types, Circuitry, and Function.Annu. Rev. Vis. Sci.31–24. 10.1146/annurev-vision-102016-061345
26
DrinnenbergA.FrankeF.MorikawaR.JüttnerJ.HillierD.HantzP.et al (2018). How diverse retinal functions arise from feedback at the first visual synapse.Neuron99117–134. 10.1016/j.neuron.2018.06.001
27
EulerT.HausseltS.MargolisD.BreuningerT.CastellX.DetwilerP.et al (2009). Eyecup scope–optical recordings of light stimulus-evoked fluorescence signals in the retina.Pflugers Arch.4571393–1414. 10.1007/s00424-008-0603-5
28
FarrowK.TeixeiraM.SzikraT.VineyT.BalintK.YoneharaK.et al (2013). Ambient illumination toggles a neuronal circuit switch in the retina and visual perception at cone threshold.Neuron78325–338. 10.1016/j.neuron.2013.02.014
29
FieldG.GauthierJ.SherA.GreschnerM.MachadoT.JepsonL.et al (2010). Functional connectivity in the retina at the resolution of photoreceptors.Nature467673–677. 10.1038/nature09424
30
FjellJ.Dib-HajjS.FriedK.BlackJ.WaxmanS. (1997). Differential expression of sodium channel genes in retinal ganglion cells.Brain Res. Mol. Brain Res.50197–204.
31
FohlmeisterJ.MillerR. (1997). Impulse encoding mechanisms of ganglion cells in the tiger salamander retina.J. Neurophysiol.781935–1947. 10.1152/jn.1997.78.4.1935
32
FrankeK.BadenT. (2017). General features of inhibition in the inner retina.J. Physiol.5955507–5515. 10.1113/JP273648
33
FrankeK.BerensP.SchubertT.BethgeM.EulerT.BadenT. (2017). Inhibition decorrelates visual feature representations in the inner retina.Nature542439–444. 10.1038/nature21394
34
FriedS.LaskerA.DesaiN.EddingtonD.RizzoJ. (2009). Axonal sodium-channel bands shape the response to electric stimulation in retinal ganglion cells.J. Neurophysiol.1011972–1987. 10.1152/jn.91081.2008
35
GeorgeJ.BadenD.GerwickW.MurrayT. (2012). Bidirectional influence of sodium channel activation on NMDA receptor–dependent cerebrocortical neuron structural plasticity.Proc. Natl. Acad. Sci. U. S. A.10919840–19845. 10.1073/pnas.1212584109
36
GoetzJ.JessenZ.JacobiA.ManiA.CoolerS.GreerD.et al (2021). Unified classification of mouse retinal ganglion cells using function, morphology, and gene expression.bioRxiv[Preprint]. 10.1101/2021.06.10.447922
37
GoldinA. (1999). Diversity of mammalian voltage-gated sodium channels.Ann. N. Y. Acad. Sci.86838–50.
38
GollischT.MeisterM. (2010). Eye smarter than scientists believed: neural computations in circuits of the retina.Neuron65150–164.
39
GongS.ZhengC.DoughtyM.LososK.DidkovskyN.SchambraU.et al (2003). A gene expression atlas of the central nervous system based on bacterial artificial chromosomes.Nature425917–925. 10.1038/nature02033
40
HahinR.CampbellD. (1983). Simple shifts in the voltage dependence of sodium channel gating caused by divalent cations.J. Gen. Physiol.82785–805. 10.1085/jgp.82.6.785
41
HahnJ.MonavarfeshaniA.QiaoM.KaoA.KölschY.KumarA.et al (2023). Evolution of neuronal cell classes and types in the vertebrate retina.bioRxiv[Preprint]. 10.1101/2023.04.07.536039
42
HelmstaedterM.BriggmanK.TuragaS.JainV.SeungH.DenkW. (2013). Connectomic reconstruction of the inner plexiform layer in the mouse retina.Nature500168–174. 10.1038/nature12346
43
HodgkinA.HuxleyA. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve.J. Physiol.117500–544. 10.1113/jphysiol.1952.sp004764
44
HuangS.HongS.De SchutterE. (2015). Non-linear leak currents affect mammalian neuron physiology.Front. Cell Neurosci.9:432. 10.3389/fncel.2015.00432
45
JacobyJ.ZhuY.DeVriesS.SchwartzG. (2015). An Amacrine Cell Circuit for Signaling Steady Illumination in the Retina.Cell Rep.132663–2670. 10.1016/j.celrep.2015.11.062
46
JosephsonI. (1996). Depolarization shifts the voltage dependence of cardiac sodium channel and calcium channel gating charge movements.Pflugers Arch.431895–904. 10.1007/s004240050083
47
KanjhanR.SivyerB. (2010). Two types of ON direction-selective ganglion cells in rabbit retina.Neurosci. Lett.483105–109. 10.1016/j.neulet.2010.07.071
48
KastnerD.BaccusS. (2011). Coordinated dynamic encoding in the retina using opposing forms of plasticity.Nat. Neurosci.141317–1322.
49
KerschensteinerD. (2022). Feature detection by retinal ganglion cells.Annu. Rev. Vis. Sci.8:112009. 10.1146/annurev-vision-100419-112009
50
KimK.RiekeF. (2001). Temporal contrast adaptation in the input and output signals of salamander retinal ganglion cells.J. Neurosci.21287–299.
51
KimK.RiekeF. (2003). Slow Na inactivation and variance adaptation in salamander retinal ganglion cells.J. Neurosci.231506–1516. 10.1523/jneurosci.23-04-01506.2003
52
KongJ.FishD.RockhillR.MaslandR. (2005). Diversity of ganglion cells in the mouse retina: unsupervised morphological classification and its limits.J. Comp. Neurol.489293–310. 10.1002/cne.20631
53
KriegerB.QiaoM.RoussoD.SanesJ.MeisterM. (2017). Four alpha ganglion cell types in mouse retina: Function, structure, and molecular signatures.PLoS One12:e0180091. 10.1371/journal.pone.0180091
54
LancasterB.AdamsP. (1986). Calcium-dependent current generating the afterhyperpolarization of hippocampal neurons.J. Neurophysiol.551268–1282. 10.1152/jn.1986.55.6.1268
55
LevickW. (1967). Receptive fields and trigger features of ganglion cells in the visual streak of the rabbit’s retina.J. Physiol.188285–307.
56
LiM.WestJ.LaiY.ScheuerT.CatterallW. (1992). Functional modulation of brain sodium channels by cAMP-dependent phosphorylation.Neuron81151–1159. 10.1016/0896-6273(92)90135-z
57
LiW.Haji GhaffariD.MisraR.WeilandJ. (2022). Retinal ganglion cell desensitization is mitigated by varying parameter constant excitation pulse trains.Front. Cell Neurosci.16:897146. 10.3389/fncel.2022.897146
58
LinB.MaslandR. (2006). Populations of wide-field amacrine cells in the mouse retina.J. Comp. Neurol.499797–809. 10.1002/cne.21126
59
MajumdarS.WässleH.JusufP.HaverkampS. (2008). Mirror-symmetrical populations of wide-field amacrine cells of the macaque monkey retina.J. Comp. Neurol.50813–27. 10.1002/cne.21666
60
ManiA.SchwartzG. (2017). Circuit mechanisms of a retinal ganglion cell with stimulus-dependent response latency and activation beyond its dendrites.Curr. Biol.27471–482. 10.1016/j.cub.2016.12.033
61
MartiszusB.TsintsadzeT.ChangW.SmithS. (2021). Enhanced excitability of cortical neurons in low-divalent solutions is primarily mediated by altered voltage-dependence of voltage-gated sodium channels.eLife101–26. 10.7554/eLife.67914
62
MaussA.VlasitsA.BorstA.FellerM. (2017). Visual circuits for direction selectivity.Annu. Rev. Neurosci.40211–230. 10.1146/annurev-neuro-072116-031335
63
MobbsP.EverettK.CookA. (1992). Signal shaping by voltage-gated currents in retinal ganglion cells.Brain Res.574217–223. 10.1016/0006-8993(92)90819-u
64
MünchT.da SilveiraR.SiegertS.VineyT.AwatramaniG.RoskaB. (2009). Approach sensitivity in the retina processed by a multifunctional neural circuit.Nat. Neurosci.121308–1316. 10.1038/nn.2389
65
NakaK.RushtonW. (1966). S-potentials from luminosity units in the retina of fish (Cyprinidae).J. Physiol.185587–599.
66
NeumannS.HüserL.OndrekaK.AulerN.HaverkampS. (2016). Cell type-specific bipolar cell input to ganglion cells in the mouse retina.Neuroscience316420–432. 10.1016/j.neuroscience.2015.12.051
67
NikolaevA.LeungK.OdermattB.LagnadoL. (2013). Synaptic mechanisms of adaptation and sensitization in the retina.Nat. Neurosci.16934–941. 10.1038/nn.3408
68
NirenbergS.MeisterM. (1997). The light response of retinal ganglion cells is truncated by a displaced amacrine circuit.Neuron18637–650. 10.1016/s0896-6273(00)80304-9
69
NumannR.CatterallW.ScheuerT. (1991). Functional modulation of brain sodium channels by protein kinase C phosphorylation.Science254115–118. 10.1126/science.1656525
70
O’BrienB.CaldwellJ.EhringG.BumstedO.BrienK.LuoS.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.508940–951. 10.1002/cne.21701
71
O’BrienB.IsayamaT.RichardsonR.BersonD. (2002). Intrinsic physiological properties of cat retinal ganglion cells.J. Physiol.538787–802. 10.1113/jphysiol.2001.013009
72
OeschN.EulerT.TaylorW. (2005). Direction-selective dendritic action potentials in rabbit retina.Neuron47739–750. 10.1016/j.neuron.2005.06.036
73
OeschN.TaylorW. (2010). Tetrodotoxin-resistant sodium channels contribute to directional responses in starburst amacrine cells.PLoS One5:e12447. 10.1371/journal.pone.0012447
74
PangJ.GaoF.WuS. (2003). Light-evoked excitatory and inhibitory synaptic inputs to ON and OFF alpha ganglion cells in the mouse retina.J. Neurosci.236063–6073. 10.1523/jneurosci.23-14-06063.2003
75
ParkS.PottackalJ.Bin KeJ.JunN.RahmaniP.KimI.et al (2018). Convergence and divergence of CRH amacrine cells in mouse retinal circuitry.J. Neurosci.383753–3766. 10.1523/JNEUROSCI.2518-17.2018
76
PeichlL.OttH.BoycottB. (1987). Alpha ganglion cells in mammalian retinae.Proc. R. Soc. Lond. B Biol. Sci.231169–197. 10.1098/rspb.1987.0040
77
Poleg-PolskyA.DiamondJ. (2011). Imperfect space clamp permits electrotonic interactions between inhibitory and excitatory synaptic conductances, distorting voltage clamp recordings.PLoS One6:e19463. 10.1371/journal.pone.0019463
78
RaghuramV.WerginzP.FriedS. (2019). Scaling of the AIS and somatodendritic compartments in α S RGCs.Front. Cell Neurosci.13:436. 10.3389/fncel.2019.00436
79
RanY.HuangZ.BadenT.SchubertT.BaayenH.BerensP.et al (2020). Type-specific dendritic integration in mouse retinal ganglion cells.Nat. Commun.11:2101. 10.1038/s41467-020-15867-9
80
RenganathanM.CumminsT.HormuzdiarW.WaxmanS. (2000). Alpha-SNS produces the slow TTX-resistant sodium current in large cutaneous afferent DRG neurons.J. Neurophysiol.84710–718.
81
RheaumeB.JereenA.BolisettyM.SajidM.YangY.RennaK.et al (2018). Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes.Nat. Commun.9:2759. 10.1038/s41467-018-05134-3
82
RoskaB.WerblinF. (2001). Vertical interactions across ten parallel, stacked representations in the mammalian retina.Nature410583–587. 10.1038/35069068
83
RoussoD.QiaoM.KaganR.YamagataM.PalmiterR.SanesJ. (2016). Two Pairs of ON and OFF retinal ganglion cells are defined by intersectional patterns of transcription factor expression.Cell Rep.151930–1944. 10.1016/j.celrep.2016.04.069
84
SahP.McLachlanE. (1991). Ca(2+)-activated K+ currents underlying the afterhyperpolarization in guinea pig vagal neurons: a role for Ca(2+)-activated Ca2+ release.Neuron7257–264. 10.1016/0896-6273(91)90264-z
85
SawantA.EbbinghausB.BleckertA.GamlinC.YuW.BersonD.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. 10.1016/j.celrep.2021.108858
86
SchachterM.OeschN.SmithR.TaylorW. (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. 10.1371/journal.pcbi.1000899
87
ScheuerT.CatterallW. (2006). Control of neuronal excitability by phosphorylation and dephosphorylation of sodium channels.Biochem. Soc. Trans.341299–1302. 10.1042/BST0341299
88
SchubertT.KerschensteinerD.Eggers, MisgeldT.KerschensteinerM.LichtmanJ.et al (2008). Development of presynaptic inhibition onto retinal bipolar cell axon terminals is subclass-specific.J. Neurophysiol.100304–316. 10.1152/jn.90202.2008
89
SchwartzG.OkawaH.DunnF.MorganJ.KerschensteinerD.WongR.et al (2012). The spatial structure of a nonlinear receptive field.Nat. Neurosci.151572–1580. 10.1038/nn.3225
90
SiegertS.CabuyE.ScherfB.KohlerH.PandaS.LeY.et al (2012). Transcriptional code and disease map for adult retinal cell types.Nat. Neurosci.15S1–S2. 10.1038/nn.3032
91
SiegertS.ScherfB.Del PuntaK.DidkovskyN.HeintzN.RoskaB. (2009). Genetic address book for retinal cell types.Nat. Neurosci.121197–1204. 10.1038/nn.2370
92
SilvaJ. (2014). Slow inactivation of Na(+) channels.Handb. Exp. Pharmacol.22133–49. 10.1007/978-3-642-41588-3_3
93
SmithB.McHughC.HiranoA.BrechaN.BarnesS. (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.432291–2304. 10.1523/JNEUROSCI.1723-22.2023
94
StuartG.SprustonN. (2015). Dendritic integration: 60 years of progress.Nat. Neurosci.181713–1721. 10.1038/nn.4157
95
SümbülU.SongS.McCullochK.BeckerM.LinB.SanesJ.et al (2014). A genetic and computational approach to structurally classify neuronal types.Nat. Commun.5:3512. 10.1038/ncomms4512
96
SunW.LiN.HeS. (2002). Large-scale morphological survey of mouse retinal ganglion cells.J. Comp. Neurol.451115–126. 10.1002/cne.10323
97
TienN.KimT.KerschensteinerD. (2016). Target-Specific Glycinergic Transmission from VGluT3-Expressing Amacrine Cells Shapes Suppressive Contrast Responses in the Retina.Cell Rep.151369–1375. 10.1016/j.celrep.2016.04.025
98
Tikidji-HamburyanA.ReinhardK.SeitterH.HovhannisyanA.ProcykC.AllenA.et al (2014). Retinal output changes qualitatively with every change in ambient illuminance.Nat. Neurosci.10.1038/nn.3891[Epub ahead of print].
99
TranN.ShekharK.WhitneyI.JacobiA.BenharI.HongG.et al (2019). Single-cell profiles of retinal ganglion cells differing in resilience to injury reveal neuroprotective genes.Neuron1041039–1055. 10.1016/j.neuron.2019.11.006
100
Tran-Van-MinhA.CazéR.AbrahamssonT.CathalaL.GutkinB.DiGregorioD. (2015). Contribution of sublinear and supralinear dendritic integration to neuronal computations.Front. Cell Neurosci.9:67. 10.3389/fncel.2015.00067
101
TsaiD.MorleyJ.SuaningG.LovellN. (2011). Sodium channel inactivation reduces retinal ganglion cell responsiveness to repetitive prosthetic stimulation. 2011 5th International IEEE/EMBS Conference on Neural Engineering.NER2011550–553. 10.1109/NER.2011.5910607
102
TurnerM.RiekeF. (2016). Synaptic Rectification Controls Nonlinear Spatial Integration of Natural Visual Inputs.Neuron901257–1271. 10.1016/j.neuron.2016.05.006
103
Van HookM.NawyS.ThoresonW. (2019). Voltage- and calcium-gated ion channels of neurons in the vertebrate retina.Prog. Retinal Eye Res.72:100760. 10.1016/j.preteyeres.2019.05.001
104
Van WartA.TrimmerJ.MatthewsG. (2007). Polarized distribution of ion channels within microdomains of the axon initial segment.J. Compar. Neurol.500339–352. 10.1002/cne.21173
105
van WykM.TaylorW.VaneyD. (2006). Local edge detectors: a substrate for fine spatial vision at low temporal frequencies in rabbit retina.J. Neurosci.2613250–13263. 10.1523/JNEUROSCI.1991-06.2006
106
van WykM.WässleH.TaylorW. (2009). Receptive field properties of ON- and OFF-ganglion cells in the mouse retina.Vis. Neurosci.26297–308. 10.1017/S0952523809990137
107
VelteT.MaslandR. (1999). Action potentials in the dendrites of retinal ganglion cells.J. Neurophysiol.811412–1417. 10.1152/jn.1999.81.3.1412
108
WangY.WeickM.DembJ. (2011). Spectral and temporal sensitivity of cone-mediated responses in mouse retinal ganglion cells.J. Neurosci.317670–7681. 10.1523/JNEUROSCI.0629-11.2011
109
WassleH. (2004). Parallel processing in the mammalian retina.Nat. Rev. Neurosci.5747–757.
110
WeickM.DembJ. (2011). Delayed-Rectifier K Channels Contribute to Contrast Adaptation in Mammalian Retinal Ganglion Cells.Neuron71166–179. 10.1016/j.neuron.2011.04.033
111
WerginzP.RaghuramV.FriedS. (2020). Tailoring of the axon initial segment shapes the conversion of synaptic inputs into spiking output in OFF-α T retinal ganglion cells.Sci. Adv.6eabb6642. 10.1126/sciadv.abb6642
112
WienbarS.SchwartzG. (2022). Differences in spike generation instead of synaptic inputs determine the feature selectivity of two retinal cell types.Neuron1102110–2123. 10.1016/j.neuron.2022.04.012
113
YaoX.CafaroJ.McLaughlinA.PostmaF.PaulD.AwatramaniG.et al (2018). Gap junctions contribute to differential light adaptation across direction-selective retinal ganglion cells.Neuron100216–228. 10.1016/j.neuron.2018.08.021
114
YinL.SmithR.SterlingP.BrainardD. (2009). Physiology and morphology of color-opponent ganglion cells in a retina expressing a dual gradient of S and M opsins.J. Neurosci.292706–2724.
115
YuW.El-DanafR.OkawaH.PacholecJ.MattiU.SchwarzK.et al (2018). Synaptic convergence patterns onto retinal ganglion cells are preserved despite topographic variation in pre- and postsynaptic territories.Cell Rep.252017–2026. 10.1016/j.celrep.2018.10.089
116
YuX. (2006). The role of intracellular sodium in the regulation of NMDA-receptor-mediated channel activity and toxicity.Mol. Neurobiol.3363–80.
117
ZhangC.McCallM. (2012). Receptor targets of amacrine cells.Vis. Neurosci.2911–29. 10.1017/S0952523812000028
118
ZhangY.KimI.SanesJ.MeisterM. (2012). The most numerous ganglion cell type of the mouse retina is a selective feature detector.Proc. Natl. Acad. Sci. U. S. A.109E2391–E2398. 10.1073/pnas.1211547109
Summary
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
Volume
17 - 2023
Edited by
Richard Anthony DeFazio, University of Michigan, United States
Reviewed by
Caroline Adams, University of Pennsylvania, United States
Tomomi Ichinose, Wayne State University, United States
Updates
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.deTimm Schubert, timm.schubert@cin.uni-tuebingen.de
Disclaimer
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.