Auditory nerve fiber excitability for alternative electrode placement in the obstructed human cochlea: electrode insertion in scala vestibuli versus scala tympani.

OBJECTIVE: The cochlear implant (CI) belongs to the most successful neuro-prostheses. Traditionally, the stimulating electrode arrays are inserted into the scala tympani (ST), the lower cochlear cavity, which enables simple surgical access. However, often deep insertion is blocked, e.g., by ossification, and the auditory nerve fibers (ANFs) of lower frequency regions cannot be stimulated causing severe restrictions in speech understanding. As an alternative, the CI can be inserted into the scala vestibuli (SV), the other upper cochlear cavity. APPROACH: In this computational study, the excitability of 25 ANFs are compared for stimulation with ST and SV implants. We employed a 3-dimensional realistic human cochlear model with lateral wall electrodes based on a µ-CT dataset and manually traced fibers. A finite element approach in combination with a compartment model of a spiral ganglion cell was used to simulate monophasic stimulation with anodic (ANO) and cathodic (CAT) pulses of 50 µs. MAIN RESULTS: ANO thresholds are lower in ST (mean/std = µ/σ =189/55 µA) stimulation compared to SV (µ/σ = 323/119 µA) stimulation. Contrary, CAT thresholds are higher for the ST array (µ/σ = 165/42 µA) compared to the SV array (µ/σ = 122/46 µA). The threshold amplitude depends on the specific fiber-electrode spatial relationship, such as lateral distance from the cochlear axis, the angle between electrode and target ANF, and the curvature of the peripheral process. For CAT stimulation the SV electrodes show a higher selectivity leading to less cross-stimulation of additional fibers from different cochlear areas. SIGNIFICANCE: We present a first simulation study with a human cochlear model that investigates an additional CI placement into the SV and its impact on the excitation behavior. Results predict comparable outcomes to ST electrodes which confirms that SV implantation might be an alternative for patients with a highly obstructed ST.

Using Patient-Specific 3D Modeling and Simulations to Optimize Microwave Ablation Therapy for Liver Cancer.

Microwave ablation (MWA) of liver tumors presents challenges like under- and over-ablation, potentially leading to inadequate tumor destruction and damage to healthy tissue. This study aims to develop personalized three-dimensional (3D) models to simulate MWA for liver tumors, incorporating patient-specific characteristics. The primary objective is to validate the predicted ablation zones compared to clinical outcomes, offering insights into MWA before therapy to facilitate accurate treatment planning. Contrast-enhanced CT images from three patients were used to create 3D models. The simulations used coupled electromagnetic wave propagation and bioheat transfer to estimate the temperature distribution, predicting tumor destruction and ablation margins. The findings indicate that prolonged ablation does not significantly improve tumor destruction once an adequate margin is achieved, although it increases tissue damage. There was a substantial overlap between the clinical ablation zones and the predicted ablation zones. For patient 1, the Dice score was 0.73, indicating high accuracy, with a sensitivity of 0.72 and a specificity of 0.76. For patient 2, the Dice score was 0.86, with a sensitivity of 0.79 and a specificity of 0.96. For patient 3, the Dice score was 0.8, with a sensitivity of 0.85 and a specificity of 0.74. Patient-specific 3D models demonstrate potential in accurately predicting ablation zones and optimizing MWA treatment strategies.

Integrating DNA methylation and gene expression data in a single gene network using the iNETgrate package.

Analyzing different omics data types independently is often too restrictive to allow for detection of subtle, but consistent, variations that are coherently supported based upon different assays. Integrating multi-omics data in one model can increase statistical power. However, designing such a model is challenging because different omics are measured at different levels. We developed the iNETgrate package ( https://bioconductor.org/packages/iNETgrate/ ) that efficiently integrates transcriptome and DNA methylation data in a single gene network. Applying iNETgrate on five independent datasets improved prognostication compared to common clinical gold standards and a patient similarity network approach.

Effects of Degrees of Degeneration on the Electrical Excitation of Human Spiral Ganglion Neurons Based on a High-Resolution Computer Model.

After hearing loss retrograde degeneration of spiral ganglion neurons (SGNs) has been described. Studies modeling the effects of degeneration mostly omitted peripheral processes (dendrites). Recent experimental observations indicated that degenerating SGNs manifested also a reduced diameter of their dendrites. We simulated populations of 400 SGNs inside a high resolution cochlear model with a cochlear implant, based on µCT scans of a human temporal bone. Cochlear implant stimuli were delivered as biphasic pulses in a monopolar configuration. Three SGN situations were simulated, based on our previous measurements of human SGN dendrites: (A) SGNs with intact dendrites (before degeneration), (B) degenerating SGNs, dendrites with a smaller diameter but original length, (C) degenerating SGNs, dendrites omitted. SGN fibers were mapped to characteristic frequency, and place pitch was estimated from excitation profiles. Results from degenerating SGNs (B, C) were similar. Most action potentials were initiated in the somatic area for all cases (A, B, C), except for areas near stimulating electrodes in the apex with intact SGNs (A), where action potentials were initiated in the distal dendrite. In most cases, degenerating SGNs had lower thresholds than intact SGNs (A) (down to -2 dB). Excitation profiles showed increased ectopic activation, i.e., activation of unintended neuronal regions, as well as similar neuronal regions excited by different apical electrodes, for degenerating SGNs (B, C). The estimated pitch showed cases of pitch reversals in apical electrodes for intact SGNs (A), as well as mostly identical pitches evoked by the four most apical electrodes for degenerating SGNs (B, C). In conclusion, neuronal excitation profiles to electrical stimulation exhibited similar traits in both ways of modeling SGN degeneration. Models showed degeneration of dendrites caused increased ectopic activation, as well as similar excitation profiles and pitch evoked by different apical electrodes. Therefore, insertion of electrodes beyond approximately 450° may not provide any benefit if SGN dendrites are degenerated.

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 (a) development of cell and electrode geometries including meshing strategies, (b) assignment of physics for the conducting spaces and the realization of active electrodes, (c) implementation of the HH model, and (d) coupling of the physics to get a fully described model.Main results.Several implementation examples are briefly presented: (a) a full FEM implementation of a HH model cell stimulated with a honeycomb electrode, (b) the electric field of a cochlear electrode placed inside the cochlea, and (c) 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.

Polarity Sensitivity of Human Auditory Nerve Fibers Based on Pulse Shape, Cochlear Implant Stimulation Strategy and Array.

Neural health is of great interest to determine individual degeneration patterns for improving speech perception in cochlear implant (CI) users. Therefore, in recent years, several studies tried to identify and quantify neural survival in CI users. Among all proposed techniques, polarity sensitivity is a promising way to evaluate the neural status of auditory nerve fibers (ANFs) in CI users. Nevertheless, investigating neural health based on polarity sensitivity is a challenging and complicated task that involves various parameters, and the outcomes of many studies show contradictory results of polarity sensitivity behavior. Our computational study benefits from an accurate three-dimensional finite element model of a human cochlea with realistic human ANFs and determined ANF degeneration pattern of peripheral part with a diminishing of axon diameter and myelination thickness based on degeneration levels. In order to see how different parameters may impact the polarity sensitivity behavior of ANFs, we investigated polarity behavior under the application of symmetric and asymmetric pulse shapes, monopolar and multipolar CI stimulation strategies, and a perimodiolar and lateral CI array system. Our main findings are as follows: (1) action potential (AP) initiation sites occurred mainly in the peripheral site in the lateral system regardless of stimulation strategies, pulse polarities, pulse shapes, cochlear turns, and ANF degeneration levels. However, in the perimodiolar system, AP initiation sites varied between peripheral and central processes, depending on stimulation strategies, pulse shapes, and pulse polarities. (2) In perimodiolar array, clusters formed in threshold values based on cochlear turns and degeneration levels for multipolar strategies only when asymmetric pulses were applied. (3) In the perimodiolar array, a declining trend in polarity (anodic threshold/cathodic threshold) with multipolar strategies was observed between intact or slight degenerated cases and more severe degenerated cases, whereas in the lateral array, cathodic sensitivity was noticed for intact and less degenerated cases and anodic sensitivity for cases with high degrees of degeneration. Our results suggest that a combination of asymmetric pulse shapes, focusing more on multipolar stimulation strategies, as well as considering the distances to the modiolus wall, allows us to distinguish the degeneration patterns of ANFs across the cochlea.

Dendritic Degeneration of Human Auditory Nerve Fibers and Its Impact on the Spiking Pattern Under Regular Conditions and During Cochlear Implant Stimulation.

Due to limitations of human in vivo studies, detailed computational models enable understanding the neural signaling in the degenerated auditory system and cochlear implants (CIs). Four human cochleae were used to quantify hearing levels depending on dendritic changes in diameter and myelination thickness from type I of the auditory nerve fibers (ANFs). Type I neurons transmit the auditory information as spiking pattern from the inner hair cells (IHCs) to the cochlear nucleus. The impact of dendrite diameter and degree of myelination on neural signal transmission was simulated for (1) synaptic excitation via IHCs and (2) stimulation from CI electrodes. An accurate three-dimensional human cochlear geometry, along with 30 auditory pathways, mimicked the CI environment. The excitation properties of electrical potential distribution induced by two CI were analyzed. Main findings: (1) The unimodal distribution of control dendrite diameters becomes multimodal for hearing loss cases; a group of thin dendrites with diameters between 0.3 and 1 µm with a peak at 0.5 µm appeared. (2) Postsynaptic currents from IHCs excite such thin dendrites easier and earlier than under control conditions. However, this advantage is lost as their conduction velocity decreases proportionally with the diameter and causes increased spike latency and jitter in soma and axon. Firing probability reduces through the soma passage due to the low intracellular current flow in thin dendrites during spiking. (3) Compared with dendrite diameter, variations in myelin thickness have a small impact on spiking performance. (4) Contrary to synaptic excitation, CIs cause several spike initiation sites in dendrite, soma region, and axon; moreover, fiber excitability reduces with fiber diameter. In a few cases, where weak stimuli elicit spikes of a target neuron (TN) in the axon, dendrite diameter reduction has no effect. However, in many cases, a spike in a TN is first initiated in the dendrite, and consequently, dendrite degeneration demands an increase in threshold currents. (5) Threshold currents of a TN and co-stimulation of degenerated ANFs in other frequency regions depend on the electrode position, including its distance to the outer wall, the cochlear turn, and the three-dimensional pathway of the TN.

Finite element analysis and three-dimensional reconstruction of tonotopically aligned human auditory fiber pathways: A computational environment for modeling electrical stimulation by a cochlear implant based on micro-CT.

The application of cochlear implants can be studied with computational models. The electrical potential distribution induced by an implanted device is evaluated with a volume conductor model, which is used as input for neuron models to simulate the reaction of cochlear neurons to micro-stimulation. In order to reliably predict the complex excitation profiles it is vital to consider an accurate representation of the human cochlea geometry including detailed three-dimensional pathways of auditory neurons reaching from the organ of Corti through the cochlea-volume. In this study, high-resolution micro-CT imaging (Δx = Δy = Δz = 3 µm) was used to reconstruct the pathways of 30 tonotopically organized nerve fiber bundles, distributed over eight octaves (11500-40 Hz). Results of the computational framework predict: (i) the peripheral process is most sensitive to cathodic stimulation (CAT), (ii) in many cases CAT elicits spikes in the peripheral terminal at threshold but with larger stimuli there is a second spike initiation site within the peripheral process, (iii) anodic stimuli (ANO) can excite the central process even at threshold, (iv) the recruitment of fibers by electrodes located in the narrowing middle- and apical turn is complex and impedes focal excitation of low frequency fibers, (v) degenerated cells which lost the peripheral process are more sensitive to CAT when their somata are totally covered with 2 membranes of a glial cell but they become ANO sensitive when the myelin covering is reduced.