Upper stimulation threshold for retinal ganglion cell activation

Neural activity of retinal ganglion cells under continuous, dynamically-modulated high frequency electrical stimulation

Objective.Current retinal prosthetics are limited in their ability to precisely control firing patterns of functionally distinct retinal ganglion cell (RGC) types. The aim of this study was to characterise RGC responses to continuous, kilohertz-frequency-varying stimulation to assess its utility in controlling RGC activity.Approach.We usedin vitropatch-clamp experiments to assess electrically-evoked ON and OFF RGC responses to frequency-varying pulse train sequences. In each sequence, the stimulation amplitude was kept constant while the stimulation frequency (0.5-10 kHz) was changed every 40 ms, in either a linearly increasing, linearly decreasing or randomised manner. The stimulation amplitude across sequences was increased from 10 to 300µA.Main results.We found that continuous stimulation without rest periods caused complex and irreproducible stimulus-response relationships, primarily due to strong stimulus-induced response adaptation and influence of the preceding stimulus frequency on the response to a subsequent stimulus. In addition, ON and OFF populations showed different sensitivities to continuous, frequency-varying pulse trains, with OFF cells generally exhibiting more dependency on frequency changes within a sequence. Finally, the ability to maintain spiking behaviour to continuous stimulation in RGCs significantly reduced over longer stimulation durations irrespective of the frequency order.Significance.This study represents an important step in advancing and understanding the utility of continuous frequency modulation in controlling functionally distinct RGCs. Our results indicate that continuous, kHz-frequency-varying stimulation sequences provide very limited control of RGC firing patterns due to inter-dependency between adjacent frequencies and generally, different RGC types do not display different frequency preferences under such stimulation conditions. For future stimulation strategies using kHz frequencies, careful consideration must be given to design appropriate pauses in stimulation, stimulation frequency order and the length of continuous stimulation duration.

Bi-directional electrical recording and stimulation of the intact retina with a screen-printed soft probe: a feasibility study

Introduction

Electrophysiological investigations of intact neural circuits are challenged by the gentle and complex nature of neural tissues. Bi-directional electrophysiological interfacing with the retina, in its intact form, is particularly demanding and currently there is no feasible approach to achieve such investigations. Here we present a feasibility study of a novel soft multi-electrode array suitable for bi-directional electrophysiological study of the intact retina.

Methods

Screen-printed soft electrode arrays were developed and tested. The soft probes were designed to accommodate the curvature of the retina in the eye and offer an opportunity to study the retina in its intact form.

Results

For the first time, we show both electrical recording and stimulation capabilities from the intact retina. In particular, we demonstrate the ability to characterize retina responses to electrical stimulation and reveal stable, direct, and indirect responses compared with ex-vivo conditions.

Discussion

These results demonstrate the unique performances of the new probe while also suggesting that intact retinas retain better stability and robustness than ex-vivo retinas making them more suitable for characterizing retina responses to electrical stimulation.

Optimizing the fabrication of a 3D high-resolution implant for neural stimulation

BACKGROUND

Tissue-integrated micro-electronic devices for neural stimulation hold great potential in restoring the functionality of degenerated organs, specifically, retinal prostheses, which are aimed at vision restoration. The fabrication process of 3D polymer-metal devices with high resolution and a high aspect-ratio (AR) is very complex and faces many challenges that impair its functionality.

APPROACH

Here we describe the optimization of the fabrication process of a bio-functionalized 3D high-resolution 1mm circular subretinal implant composed of SU-8 polymer integrated with dense gold microelectrodes (23μm pitch) passivated with 3D micro-well-like structures (20μm diameter, 3μm resolution). The main challenges were overcome by step-by-step planning and optimization while utilizing a two-step bi-layer lift-off process; bio-functionalization was carried out by N2 plasma treatment and the addition of a bio-adhesion molecule.

MAIN RESULTS

In-vitro and in-vivo investigations, including SEM and FIB cross section examinations, revealed a good structural design, as well as a good long-term integration of the device in the rat sub-retinal space and cell migration into the wells. Moreover, the feasibility of subretinal neural stimulation using the fabricated device was demonstrated in-vitro by electrical activation of rat's retina.

CONCLUSIONS

The reported process and optimization steps described here in detail can aid in designing and fabricating retinal prosthetic devices or similar neural implants.

Effect of Stimulation Current in Transcutaneous Vagus Nerve Stimulation (tVNS): A Study Using Concurrent Magnetoencephalography (MEG)

Transcutaneous vagus nerve stimulation (tVNS) is a non-invasive method of brain stimulation that has been investigated for its use in the clinical treatment of a number of different conditions. There has been little investigation into the stimulation current that is delivered and the effect on individual variability in response to tVNS.Seventeen participants underwent tVNS, and stimulation current was determined based on individual pain threshold. To investigate individual variability, brain dynamics were measured concurrently using magnetoencephalography (MEG) in response to two different stimulation protocols of tVNS. The first protocol consisted of a sequence of equally spaced short (1ms) stimulation pulses applied 24 times per second (24 Hz), and the second consisted of a sequence of 24 pulses per second spaced according to a 6 Hz pulse frequency modulation (PFM). Both stimulation sequences were delivered to the cymba concha in the left ear.The difference in brain responses to the two sequences was initially calculated using a one-sample t-test at the group level, based on z-scoring of the data at the individual level, and no statistically significant differences were observed. Further investigation of individual variability suggested that participants fell into two groups; one that responded more strongly to 24 Hz and one that responded more strongly to the irregular spacing of pulses in the PFM protocol.We tested whether the stimulation current that the participant received could predict how they would respond to the stimulation, but we did not observe any correlation. This supports the literature that suggests that selecting stimulation current based on individual pain threshold is a suitable procedure for tVNS, and higher stimulation intensities does not correspond to stronger brain response. Further investigation into individual variability in response to different frequencies and pulse spacing of tVNS should also be investigated further and may lead to the development of personalised stimulation protocols.Clinical relevance- The stimulation current at which tVNS is delivered does not appear to influence brain response to stimulation, and the value of stimulation current should be selected based on individual participant comfort.

A simple model considering spiking probability during extracellular axon stimulation

The spiking probability of an electrically stimulated axon as a function of stimulus amplitude increases in a sigmoidal dependency from 0 to 1. However, most computer simulation studies for neuroprosthetic applications calculate thresholds for neural targets with a deterministic model and by reducing the sigmoid curve to a step function, they miss an important information about the control signal, namely how the spiking efficiency increases with stimulus intensity. Here, this spiking efficiency is taken into account in a compartment model of the Hodgkin Huxley type where a noise current is added in every compartment with an active membrane. A key parameter of the model is a common factor knoise which defines the ion current fluctuations across the cell membrane for every compartment by its maximum sodium ion conductance. In the standard model Gaussian signals are changed every 2.5 μs as a compromise of accuracy and computational costs. Additionally, a formula for other noise transmission times is presented and numerically tested. Spiking probability as a function of stimulus intensity can be approximated by the cumulative distribution function of the normal distribution with RS = σ/μ. Relative spread RS, introduced by Verveen, is a measure for the spread (normalized by the threshold intensity μ), that decreases inversely with axon diameter. Dynamic range, a related measure used in neuroprosthetic studies, defines the intensity range between 10% and 90% spiking probability. We show that (i) the dynamic range normalized by threshold is 2.56 times RS, (ii) RS increases with electrode—axon distance and (iii) we present knoise values for myelinated and unmyelinated axon models in agreement with recoded RS data. The presented method is applicable for other membrane models and can be extended to whole neurons that are described by multi-compartment models.

A finite element method framework to model extracellular neural stimulation

OBJECTIVE

Increasing complexity in extracellular stimulation experiments and neural implant design also requires realistic computer simulations capable of modeling the neural activity of nerve cells under the influence of an electrical stimulus. Classical model approaches are often based on simplifications, are not able to correctly calculate the electric field generated by complex electrode designs, and do not consider electrical effects of the cell on its surrounding. A more accurate approach is the finite element method (FEM), which provides necessary techniques to solve the Poisson equation for complex geometries under consideration of electrical tissue properties. Especially in situations where neurons experience large and non-symmetric extracellular potential gradients, a FEM solution that implements the cell membrane model can improve the computer simulation results. To investigate the response of neurons in an electric field generated by complex electrode designs, a FEM framework for extracellular stimulation was developed in COMSOL.

APPROACH

Methods to implement morphologically- and biophysically-detailed neurons including active Hodgkin-Huxley (HH) cell membrane dynamics as well as the stimulation setup are described in detail. Covered methods are (i) development of cell and electrode geometries including meshing strategies, (ii) assignment of physics for the conducting spaces and the realization of active electrodes, (iii) implementation of the HH model, and (iv) coupling of the physics to get a fully described model.

MAIN RESULTS

Several implementation examples are briefly presented: (i) a full FEM implementation of a HH model cell stimulated with a honeycomb electrode, (ii) the electric field of a cochlear electrode placed inside the cochlea, and (iii) a proof of concept implementation of a detailed double-cable cell membrane model for myelinated nerve fibers.

SIGNIFICANCE

The presented concepts and methods provide basic and advanced techniques to realize a full FEM framework for innovative studies of neural excitation in response to extracellular stimulation.

Block Phenomena During Electric Micro-Stimulation of Pyramidal Cells and Retinal Ganglion Cells

Electric micro-stimulation of the nervous system is a means to restore various body functions. The stimulus amplitude necessary to generate action potentials, the lower threshold (LT), is well characterized for many neuronal populations. However, electric overstimulation above an upper threshold (UT) prevents action potential generation and therefore hinders optimal neuro-rehabilitation. Previous studies demonstrated the impact of the UT in micro-stimulation of retinal ganglion cells (RGCs). The observed phenomenon is mostly explained by (i) reversed sodium ion flow in the soma membrane, and (ii) anodal surround block that hinders spike conduction in strongly hyperpolarized regions of the axon at high stimulus intensities. However, up to now, no detailed study of the nature of these phenomena has been presented, particularly for different cell types. Here, we present computational analyses of LT and UT for layer 5 pyramidal cells (PCs) as well as alpha RGCs. Model neurons were stimulated in close vicinity to the cell body and LTs and UTs as well as the ratio UT/LT were compared. Aside from a simple point source electrode and monophasic stimuli also realistic electrode and pulse configurations were examined. The analysis showed: (i) in RGCs, the soma contributed to action potential initiation and block for small electrode distances, whereas in PCs the soma played no role in LTs or UTs. (ii) In both cell types, action potential always initiated within the axon initial segment at LT. (iii) In contrast to a complete block of spike conductance at UT that occurred in RGCs, an incomplete block of spiking appeared in PC axon collaterals. (iv) PC axon collateral arrangement influenced UTs but had small impact on LTs. (v) Population responses of RGCs change from circular regions of activation to ring-shaped patterns for increasing stimulus amplitude. A better understanding of the stimulation window that can reliably activate target neurons will benefit the future development of neuroprostheses.

M1-Type, but Not M4-Type, Melanopsin Ganglion Cells Are Physiologically Tuned to the Central Circadian Clock

Proper circadian photoentrainment is crucial for the survival of many organisms. In mammals, intrinsically photosensitive retinal ganglion cells (ipRGCs) can use the photopigment melanopsin to sense light independently from rod and cone photoreceptors and send this information to many brain nuclei such as the suprachiasmatic nucleus (SCN), the site of the central circadian pacemaker. Here, we measure ionic currents and develop mathematical models of the electrical activity of two types of ipRGCs: M1, which projects to the SCN, and M4, which does not. We illustrate how their ionic properties differ, mainly how ionic currents generate lower spike rates and depolarization block in M1 ipRGCs. Both M1 and M4 cells have large geometries and project to higher visual centers of the brain via the optic nerve. Using a partial differential equation model, we show how axons of M1 and M4 cells faithfully convey information from the soma to the synapse even when the signal at the soma is attenuated due to depolarization block. Finally, we consider an ionic model of circadian photoentrainment from ipRGCs synapsing on SCN neurons and show how the properties of M1 ipRGCs are tuned to create accurate transmission of visual signals from the retina to the central pacemaker, whereas M4 ipRGCs would not evoke nearly as efficient a postsynaptic response. This work shows how ipRGCs and SCN neurons' electrical activities are tuned to allow for accurate circadian photoentrainment.

Computational challenges and opportunities for a bi-directional artificial retina

A future artificial retina that can restore high acuity vision in blind people will rely on the capability to both read (observe) and write (control) the spiking activity of neurons using an adaptive, bi-directional and high-resolution device. Although current research is focused on overcoming the technical challenges of building and implanting such a device, exploiting its capabilities to achieve more acute visual perception will also require substantial computational advances. Using high-density large-scale recording and stimulation in the primate retina with an ex vivo multi-electrode array lab prototype, we frame several of the major computational problems, and describe current progress and future opportunities in solving them. First, we identify cell types and locations from spontaneous activity in the blind retina, and then efficiently estimate their visual response properties by using a low-dimensional manifold of inter-retina variability learned from a large experimental dataset. Second, we estimate retinal responses to a large collection of relevant electrical stimuli by passing current patterns through an electrode array, spike sorting the resulting recordings and using the results to develop a model of evoked responses. Third, we reproduce the desired responses for a given visual target by temporally dithering a diverse collection of electrical stimuli within the integration time of the visual system. Together, these novel approaches may substantially enhance artificial vision in a next-generation device.

Stimulation Strategies for Improving the Resolution of Retinal Prostheses

Electrical stimulation using implantable devices with arrays of stimulating electrodes is an emerging therapy for neurological diseases. The performance of these devices depends greatly on their ability to activate populations of neurons with high spatiotemporal resolution. To study electrical stimulation of populations of neurons, retina serves as a useful model because the neural network is arranged in a planar array that is easy to access. Moreover, retinal prostheses are under development to restore vision by replacing the function of damaged light sensitive photoreceptors, which makes retinal research directly relevant for curing blindness. Here we provide a progress review on stimulation strategies developed in recent years to improve the resolution of electrical stimulation in retinal prostheses. We focus on studies performed with explanted retinas, in which electrophysiological techniques are the most advanced. We summarize achievements in improving the spatial and temporal resolution of electrical stimulation of the retina and methods to selectively stimulate neurons with different visual functions. Future directions for retinal prostheses development are also discussed, which could provide insights for other types of neuromodulatory devices in which high-resolution electrical stimulation is required.

Mediating Retinal Ganglion Cell Spike Rates Using High-Frequency Electrical Stimulation

Recent retinal studies have directed more attention to sophisticated stimulation strategies based on high-frequency (>1.0 kHz) electrical stimulation (HFS). In these studies, each retinal ganglion cell (RGC) type demonstrated a characteristic stimulus-strength-dependent response to HFS, offering the intriguing possibility of focally targeting retinal neurons to provide useful visual information by retinal prosthetics. Ionic mechanisms are known to affect the responses of electrogenic cells during electrical stimulation. However, how these mechanisms affect RGC responses is not well understood at present, particularly when applying HFS. Here, we investigate this issue via an in silico model of the RGC. We calibrate and validate the model using an in vitro retinal preparation. An RGC model based on accurate biophysics and realistic representation of cell morphology, was used to investigate how RGCs respond to HFS. The model was able to closely replicate the stimulus-strength-dependent suppression of RGC action potentials observed experimentally. Our results suggest that spike inhibition during HFS is due to local membrane hyperpolarization caused by outward membrane currents near the stimulus electrode. In addition, the extent of HFS-induced inhibition can be largely altered by the intrinsic properties of the inward sodium current. Finally, stimulus-strength-dependent suppression can be modulated by a wide range of stimulation frequencies, under generalized electrode placement conditions. In vitro experiments verified the computational modeling data. This modeling and experimental approach can be extended to further our understanding on the effects of novel stimulus strategies by simulating RGC stimulus-response profiles over a wider range of stimulation frequencies and electrode locations than have previously been explored.

Analysis of upper threshold mechanisms of spherical neurons during extracellular stimulation

Exceeding a certain stimulation strength can prevent the generation of somatic action potentials, as has been demonstrated in vitro with extracellularly stimulated dorsal root ganglion cells as well as retinal ganglion cells. This phenomenon, termed upper threshold, is currently thought to be a consequence of sodium current reversal in strongly depolarized regions. Here we analyze the contribution of membrane kinetics, using spherical model neurons that are stimulated externally with a microelectrode, in more detail. During extracellular pulse application, the electric field depolarizes one part and hyperpolarizes the other part of the cell. Strong transmembrane currents are generated only in the active depolarized region, changing the overall polarization level. The asymmetric membrane voltage distribution caused by the stimulus strongly influences the cell's behavior during and even after the stimulus. Effects on membrane voltage and transmembrane currents during and after the stimulus are shown and discussed in detail. Aside from the sodium current reversal, two more key mechanisms were identified in causing the upper threshold: strong potassium currents and inactivation of sodium channels. The contributions of the mechanisms involved strongly depend on cell properties, stimulus parameters, and other factors such as temperature. The conclusions presented here are based on several retinal ganglion cell models of the Fohlmeister group, a model with original Hodgkin-Huxley membrane, and a pyramidal cell model. NEW & NOTEWORTHY The upper threshold phenomenon in extracellular stimulation is analyzed in detail for spherical cells. Three main mechanisms were identified that prevent the generation of action potentials at high stimulation strengths: 1) strong potassium currents, 2) inactivating sodium ion channels, and 3) sodium current reversal. Ion channel kinetics in retinal ganglion cells, pyramidal cells, and the original Hodgkin-Huxley model were investigated under the influence of an extracellular stimulus.

Retinotopic Responses in the Visual Cortex Elicited by Epiretinal Electrical Stimulation in Normal and Retinal Degenerate Rats

Purpose

Electronic retinal prostheses restore vision in people with outer retinal degeneration by electrically stimulating the inner retina. We characterized visual cortex electrophysiologic response elicited by electrical stimulation of retina in normally sighted and retinal degenerate rats.

Methods

Nine normally sighted Long Evans and 11 S334ter line 3 retinal degenerate (rd) rats were used to map cortical responses elicited by epiretinal electrical stimulation in four quadrants of the retina. Six normal and six rd rats were used to compare the dendritic spine density of neurons in the visual cortex.

Results

The rd rats required higher stimulus amplitudes to elicit responses in the visual cortex. The cortical electrically evoked responses (EERs) for both healthy and rd rats show a dose-response characteristic with respect to the stimulus amplitude. The EER maps in healthy rats show retinotopic organization. For rd rats, cortical retinotopy is not well preserved. The neurons in the visual cortex of rd rats show a 10% higher dendritic spine density than in the healthy rats.

Conclusions

Cortical activity maps, produced when epiretinal stimulation is applied to quadrants of the retina, exhibit retinotopy in normal but not rd rats. This is likely due to a combination of degeneration of the retina and increased stimulus thresholds in rd, which broadens the activated area of the retina.

Translational Relevance

Loss of retinotopy is evident in rd rats. If a similar loss of retinotopy is present in humans, retinal prostheses design must include flexibility to account for patient specific variability.

Irregularly timed electrical pulses reduce adaptation of retinal ganglion cells

OBJECTIVE

Retinal prostheses aim to provide visual percepts to blind people affected by diseases caused by photoreceptor degeneration. One of the main challenges presented by current devices is neural adaptation in the retina, which is believed to be the cause of fading-an effect where artificially produced percepts disappear over a short period of time, despite continuous stimulation of the retina. We aim to understand the neural adaptation generated in retinal ganglion cells (RGCs) during electrical stimulation.

APPROACH

Current visual prostheses use electrical pulses with fixed frequencies and amplitudes modulated over hundreds of milliseconds to stimulate the retina. However, in nature, neuronal spiking occurs with stochastic timing, hence the information received naturally from other neurons by RGCs is irregularly timed. We used a single epiretinal electrode to stimulate and compare rat RGC responses to stimulus trains of biphasic pulses delivered at regular and random inter-pulse intervals (IPI), the latter taken from an exponential distribution.

MAIN RESULTS

Our observations suggest that stimulation with random IPIs result in lower adaptation rates than stimulation with constant IPIs at frequencies of 50 Hz and 200 Hz. We also found a high proportion of lower amplitude action potentials, or spikelets. The spikelets were more prominent at high stimulation frequencies (50 Hz and 200 Hz) and were less susceptible to adaptation, but it was not clear if they propagated along the axon.

SIGNIFICANCE

Using random IPI stimulation in retinal prostheses reduces the decay of RGCs and this could potentially reduce fading of electrically induced visual perception.