Chemogenetic stimulation of phrenic motor output and diaphragm activity.

Impaired diaphragm activation contributes to morbidity and mortality in many neurodegenerative diseases and neurologic injuries. We conducted experiments to determine if expression of an excitatory DREADD (designer receptors exclusively activation by designer drugs) in the mid-cervical spinal cord would enable respiratory-related activation of phrenic motoneurons to increase diaphragm activation. Wild type (C57/bl6) and ChAT-Cre mice received bilateral intraspinal (C4) injections of an adeno-associated virus (AAV) encoding the hM3D(Gq) excitatory DREADD. In wild type mice, this produced non-specific DREADD expression throughout the mid-cervical ventral horn. In ChAT-Cre mice, a Cre-dependent viral construct was used to drive DREADD expression in C4 ventral horn motoneurons, targeting the phrenic motoneuron pool. Diaphragm EMG was recorded during spontaneous breathing at 6-8 weeks post-AAV delivery. The selective DREADD ligand JHU37160 (J60) caused a bilateral, sustained (>1 hr) increase in inspiratory EMG bursting in both groups; the relative increase was greater in ChAT-Cre mice. Additional experiments in a ChAT-Cre rat model were conducted to determine if spinal DREADD activation could increase inspiratory tidal volume (VT) during spontaneous breathing without anesthesia. Three to four months after intraspinal (C4) injection of AAV driving Cre-dependent hM3D(Gq) expression, intravenous J60 resulted in a sustained (>30 min) increase in VT assessed using whole-body plethysmography. Subsequently, direct nerve recordings confirmed that J60 evoked a >50% increase in inspiratory phrenic output. The data show that mid-cervical spinal DREADD expression targeting the phrenic motoneuron pool enables ligand-induced, sustained increases in the neural drive to the diaphragm. Further development of this technology may enable application to clinical conditions associated with impaired diaphragm activation and hypoventilation.

Compressive sensing of functional connectivity maps from patterned optogenetic stimulation of neuronal ensembles.

Mapping functional connectivity between neurons is an essential step toward probing the neural computations mediating behavior. Accurately determining synaptic connectivity maps in populations of neurons is challenging in terms of yield, accuracy, and experimental time. Here, we developed a compressive sensing approach to reconstruct synaptic connectivity maps based on random two-photon cell-targeted optogenetic stimulation and membrane voltage readout of many putative postsynaptic neurons. Using a biophysical network model of interconnected populations of excitatory and inhibitory neurons, we characterized mapping recall and precision as a function of network observability, sparsity, number of neurons stimulated, off-target stimulation, synaptic reliability, propagation latency, and network topology. We found that mapping can be achieved with far fewer measurements than the standard pairwise sequential approach, with network sparsity and synaptic reliability serving as primary determinants of the performance. Our results suggest a rapid and efficient method to reconstruct functional connectivity of sparsely connected neuronal networks.

Marked point process representation of oscillatory dynamics underlying working memory.

Objective.Computational models of neural activity at the meso-scale suggest the involvement of discrete oscillatory bursts as constructs of cognitive processing during behavioral tasks. Classical signal processing techniques that attempt to infer neural correlates of behavior from meso-scale activity employ spectral representations of the signal, exploiting power spectral density techniques and time-frequency (T-F) energy distributions to capture band power features. However, such analyses demand more specialized methods that incorporate explicitly the concepts of neurophysiological signal generation and time resolution in the tens of milliseconds. This paper focuses on working memory (WM), a complex cognitive process involved in encoding, storing and retrieving sensory information, which has been shown to be characterized by oscillatory bursts in the beta and gamma band. Employing a generative model for oscillatory dynamics, we present a marked point process (MPP) representation of bursts during memory creation and readout. We show that the markers of the point process quantify specific neural correlates of WM.Approach.We demonstrate our results on field potentials recorded from the prelimbic and secondary motor cortices of three rats while performing a WM task. The generative model for single channel, band-passed traces of field potentials characterizes with high-resolution, the timings and amplitudes of transient neuromodulations in the high gamma (80-150 Hz,γ) and beta (10-30 Hz,ß) bands as an MPP. We use standard hypothesis testing methods on the MPP features to check for significance in encoding of task variables, sensory stimulus and executive control while comparing encoding capabilities of our model with other T-F methods.Main Results.Firstly, the advantages of an MPP approach in deciphering encoding mechanisms at the meso-scale is demonstrated. Secondly, the nature of state encoding by neuromodulatory events is determined. Third, we demonstrate the necessity of a higher time resolution alternative to conventionally employed T-F methods. Finally, our results underscore the novelty in interpreting oscillatory dynamics encompassed by the marked features of the point process.Significance.An MPP representation of meso-scale activity not just enables a rich, high-resolution parameter space for analysis but also presents a novel tool for diverse neural applications.

Emergent coordination underlying learning to reach to grasp with a brain-machine interface.

The development of coordinated reach-to-grasp movement has been well studied in infants and children. However, the role of motor cortex during this development is unclear because it is difficult to study in humans. We took the approach of using a brain-machine interface (BMI) paradigm in rhesus macaques with prior therapeutic amputations to examine the emergence of novel, coordinated reach to grasp. Previous research has shown that after amputation, the cortical area previously involved in the control of the lost limb undergoes reorganization, but prior BMI work has largely relied on finding neurons that already encode specific movement-related information. In this study, we taught macaques to cortically control a robotic arm and hand through operant conditioning, using neurons that were not explicitly reach or grasp related. Over the course of training, stereotypical patterns emerged and stabilized in the cross-covariance between the reaching and grasping velocity profiles, between pairs of neurons involved in controlling reach and grasp, and to a comparable, but lesser, extent between other stable neurons in the network. In fact, we found evidence of this structured coordination between pairs composed of all combinations of neurons decoding reach or grasp and other stable neurons in the network. The degree of and participation in coordination was highly correlated across all pair types. Our approach provides a unique model for studying the development of novel, coordinated reach-to-grasp movement at the behavioral and cortical levels. NEW & NOTEWORTHY Given that motor cortex undergoes reorganization after amputation, our work focuses on training nonhuman primates with chronic amputations to use neurons that are not reach or grasp related to control a robotic arm to reach to grasp through the use of operant conditioning, mimicking early development. We studied the development of a novel, coordinated behavior at the behavioral and cortical level, and the neural plasticity in M1 associated with learning to use a brain-machine interface.

Changes in cortical network connectivity with long-term brain-machine interface exposure after chronic amputation.

Studies on neural plasticity associated with brain-machine interface (BMI) exposure have primarily documented changes in single neuron activity, and largely in intact subjects. Here, we demonstrate significant changes in ensemble-level functional connectivity among primary motor cortical (MI) neurons of chronically amputated monkeys exposed to control a multiple-degree-of-freedom robot arm. A multi-electrode array was implanted in M1 contralateral or ipsilateral to the amputation in three animals. Two clusters of stably recorded neurons were arbitrarily assigned to control reach and grasp movements, respectively. With exposure, network density increased in a nearly monotonic fashion in the contralateral monkeys, whereas the ipsilateral monkey pruned the existing network before re-forming a denser connectivity. Excitatory connections among neurons within a cluster were denser, whereas inhibitory connections were denser among neurons across the two clusters. These results indicate that cortical network connectivity can be modified with BMI learning, even among neurons that have been chronically de-efferented and de-afferented due to amputation.

Proceedings of the Third Annual Deep Brain Stimulation Think Tank: A Review of Emerging Issues and Technologies.

The proceedings of the 3rd Annual Deep Brain Stimulation Think Tank summarize the most contemporary clinical, electrophysiological, imaging, and computational work on DBS for the treatment of neurological and neuropsychiatric disease. Significant innovations of the past year are emphasized. The Think Tank's contributors represent a unique multidisciplinary ensemble of expert neurologists, neurosurgeons, neuropsychologists, psychiatrists, scientists, engineers, and members of industry. Presentations and discussions covered a broad range of topics, including policy and advocacy considerations for the future of DBS, connectomic approaches to DBS targeting, developments in electrophysiology and related strides toward responsive DBS systems, and recent developments in sensor and device technologies.

Neuroplasticity subserving the operation of brain-machine interfaces.

Neuroplasticity is key to the operation of brain machine interfaces (BMIs)-a direct communication pathway between the brain and a man-made computing device. Whereas exogenous BMIs that associate volitional control of brain activity with neurofeedback have been shown to induce long lasting plasticity, endogenous BMIs that use prolonged activity-dependent stimulation--and thus may curtail the time scale that governs natural sensorimotor integration loops--have been shown to induce short lasting plasticity. Here we summarize recent findings from studies using both categories of BMIs, and discuss the fundamental principles that may underlie their operation and the longevity of the plasticity they induce. We draw comparison to plasticity mechanisms known to mediate natural sensorimotor skill learning and discuss principles of homeostatic regulation that may constrain endogenous BMI effects in the adult mammalian brain. We propose that BMIs could be designed to facilitate structural and functional plasticity for the purpose of re-organization of target brain regions and directed augmentation of sensorimotor maps, and suggest possible avenues for future work to maximize their efficacy and viability in clinical applications.

Temporal precision in population-but not individual neuron-dynamics reveals rapid experience-dependent plasticity in the rat barrel cortex.

Cortical reorganization following sensory deprivation is characterized by alterations in the connectivity between neurons encoding spared and deprived cortical inputs. The extent to which this alteration depends on Spike Timing Dependent Plasticity (STDP), however, is largely unknown. We quantified changes in the functional connectivity between layer V neurons in the vibrissal primary somatosensory cortex (vSI) (barrel cortex) of rats following sensory deprivation. One week after chronic implantation of a microelectrode array in vSI, sensory-evoked activity resulting from mechanical deflections of individual whiskers was recorded (control data) after which two whiskers on the contralateral side were paired by sparing them while trimming all other whiskers on the rat's mystacial pad. The rats' environment was then enriched by placing novel objects in the cages to encourage exploratory behavior with the spared whiskers. Sensory-evoked activity in response to individual stimulation of spared whiskers and adjacent re-grown whiskers was then recorded under anesthesia 1-2 days and 6-7 days post-trimming (plasticity data). We analyzed spike trains within 100 ms of stimulus onset and confirmed previously published reports documenting changes in receptive field sizes in the spared whisker barrels. We analyzed the same data using Dynamic Bayesian Networks (DBNs) to infer the functional connectivity between the recorded neurons. We found that DBNs inferred from population responses to stimulation of each of the spared whiskers exhibited graded increase in similarity that was proportional to the pairing duration. A significant early increase in network similarity in the spared-whisker barrels was detected 1-2 days post pairing, but not when single neuron responses were examined during the same period. These results suggest that rapid reorganization of cortical neurons following sensory deprivation may be mediated by an STDP mechanism.

Tracking single units in chronic, large scale, neural recordings for brain machine interface applications.

In the study of population coding in neurobiological systems, tracking unit identity may be critical to assess possible changes in the coding properties of neuronal constituents over prolonged periods of time. Ensuring unit stability is even more critical for reliable neural decoding of motor variables in intra-cortically controlled brain-machine interfaces (BMIs). Variability in intrinsic spike patterns, tuning characteristics, and single-unit identity over chronic use is a major challenge to maintaining this stability, requiring frequent daily calibration of neural decoders in BMI sessions by an experienced human operator. Here, we report on a unit-stability tracking algorithm that efficiently and autonomously identifies putative single-units that are stable across many sessions using a relatively short duration recording interval at the start of each session. The algorithm first builds a database of features extracted from units' average spike waveforms and firing patterns across many days of recording. It then uses these features to decide whether spike occurrences on the same channel on one day belong to the same unit recorded on another day or not. We assessed the overall performance of the algorithm for different choices of features and classifiers trained using human expert judgment, and quantified it as a function of accuracy and execution time. Overall, we found a trade-off between accuracy and execution time with increasing data volumes from chronically implanted rhesus macaques, with an average of 12 s processing time per channel at ~90% classification accuracy. Furthermore, 77% of the resulting putative single-units matched those tracked by human experts. These results demonstrate that over the span of a few months of recordings, automated unit tracking can be performed with high accuracy and used to streamline the calibration phase during BMI sessions. Our findings may be useful to the study of population coding during learning, and to improve the reliability of BMI systems and accelerate their deployment in clinical applications.

A fully automated rodent conditioning protocol for sensorimotor integration and cognitive control experiments.

Rodents have been traditionally used as a standard animal model in laboratory experiments involving a myriad of sensory, cognitive, and motor tasks. Higher cognitive functions that require precise control over sensorimotor responses such as decision-making and attentional modulation, however, are typically assessed in nonhuman primates. Despite the richness of primate behavior that allows multiple variants of these functions to be studied, the rodent model remains an attractive, cost-effective alternative to primate models. Furthermore, the ability to fully automate operant conditioning in rodents adds unique advantages over the labor intensive training of nonhuman primates while studying a broad range of these complex functions. Here, we introduce a protocol for operantly conditioning rats on performing working memory tasks. During critical epochs of the task, the protocol ensures that the animal's overt movement is minimized by requiring the animal to 'fixate' until a Go cue is delivered, akin to nonhuman primate experimental design. A simple two alternative forced choice task is implemented to demonstrate the performance. We discuss the application of this paradigm to other tasks.

Sorting and tracking neuronal spikes via simple thresholding.

A fundamental goal in systems neuroscience is to assess the individual as well as the synergistic roles of single neurons in a recorded ensemble as they relate to an observed behavior. A mandatory step to achieve this goal is to sort spikes in an extracellularly recorded mixture that belong to individual neurons through feature extraction and clustering techniques. Here, we propose an approach for approximating the often nonlinear and time varying decision boundaries between spike-derived feature classes based on a simple, yet optimal thresholding mechanism. Because thresholding is a binary classifier, we show that the complex nonlinear decision boundaries required for spike class discrimination can be achieved by adequately fusing a set of weak binary classifiers. The thresholds for these binary classifiers are adaptively estimated through a learning algorithm that maximizes the separability between the sparsely represented classes. Based on our previous work, the approach substantially reduces the computational complexity of extracting, aligning and sorting multiple single unit activity early in the data stream. Here, we also show its ability to track changes in spike features over extended periods of time, making it highly suitable for basic neuroscience studies as well as for implementation in miniaturized, fully implantable electronics in brain-machine interface applications.

Operant conditioning of a multiple degree-of-freedom brain-machine interface in a primate model of amputation.

Operant conditioning with biofeedback has been shown to be an effective method to modify neural activity to generate goal-directed actions in a brain-machine interface. It is particularly useful when neural activity cannot be mathematically mapped to motor actions of the actual body such as in the case of amputation. Here, we implement an operant conditioning approach with visual feedback in which an amputated monkey is trained to control a multiple degree-of-freedom robot to perform a reach-to-grasp behavior. A key innovation is that each controlled dimension represents a behaviorally relevant synergy among a set of joint degrees-of-freedom. We present a number of behavioral metrics by which to assess improvements in BMI control with exposure to the system. The use of non-human primates with chronic amputation is arguably the most clinically-relevant model of human amputation that could have direct implications for developing a neural prosthesis to treat humans with missing upper limbs.

Optimal space-time precoding of artificial sensory feedback through mutichannel microstimulation in bi-directional brain-machine interfaces.

Brain-machine interfaces (BMIs) aim to restore lost sensorimotor and cognitive function in subjects with severe neurological deficits. In particular, lost somatosensory function may be restored by artificially evoking patterns of neural activity through microstimulation to induce perception of tactile and proprioceptive feedback to the brain about the state of the limb. Despite an early proof of concept that subjects could learn to discriminate a limited vocabulary of intracortical microstimulation (ICMS) patterns that instruct the subject about the state of the limb, the dynamics of a moving limb are unlikely to be perceived by an arbitrarily-selected, discrete set of static microstimulation patterns, raising questions about the generalization and the scalability of this approach. In this work, we propose a microstimulation protocol intended to activate optimally the ascending somatosensory pathway. The optimization is achieved through a space-time precoder that maximizes the mutual information between the sensory feedback indicating the limb state and the cortical neural response evoked by thalamic microstimulation. Using a simplified multi-input multi-output model of the thalamocortical pathway, we show that this optimal precoder can deliver information more efficiently in the presence of noise compared to suboptimal precoders that do not account for the afferent pathway structure and/or cortical states. These results are expected to enhance the way microstimulation is used to induce somatosensory perception during sensorimotor control of artificial devices or paralyzed limbs.

A Fully Implantable, Programmable and Multimodal Neuroprocessor for Wireless, Cortically Controlled Brain-Machine Interface Applications.

Reliability, scalability and clinical viability are of utmost importance in the design of wireless Brain Machine Interface systems (BMIs). This paper reports on the design and implementation of a neuroprocessor for conditioning raw extracellular neural signals recorded through microelectrode arrays chronically implanted in the brain of awake behaving rats. The neuroprocessor design exploits a sparse representation of the neural signals to combat the limited wireless telemetry bandwidth. We demonstrate a multimodal processing capability (monitoring, compression, and spike sorting) inherent in the neuroprocessor to support a wide range of scenarios in real experimental conditions. A wireless transmission link with rate-dependent compression strategy is shown to preserve information fidelity in the neural data. At 32 channels, the neuroprocessor has been fully implemented on a 5mm×5mm nano-FPGA, and the prototyping resulted in 5.19 mW power consumption, bringing its performance within the power-size constraints for clinical use. The optimal design for compression and sorting performance was evaluated for multiple sampling frequencies, wavelet basis choice and power consumption.

NeuroQuest: a comprehensive analysis tool for extracellular neural ensemble recordings.

Analyzing the massive amounts of neural data collected using microelectrodes to extract biologically relevant information is a major challenge. Many scientific findings rest on the ability to overcome these challenges and to standardize experimental analysis across labs. This can be facilitated in part through comprehensive, efficient and practical software tools disseminated to the community at large. We have developed a comprehensive, MATLAB-based software package - entitled NeuroQuest - that bundles together a number of advanced neural signal processing algorithms in a user-friendly environment. Results demonstrate the efficiency and reliability of the software compared to other software packages, and versatility over a wide range of experimental conditions.

Neural feedback for instantaneous spatiotemporal modulation of afferent pathways in bi-directional brain-machine interfaces.

In bi-directional brain-machine interfaces (BMIs), precisely controlling the delivery of microstimulation, both in space and in time, is critical to continuously modulate the neural activity patterns that carry information about the state of the brain-actuated device to sensory areas in the brain. In this paper, we investigate the use of neural feedback to control the spatiotemporal firing patterns of neural ensembles in a model of the thalamocortical pathway. Control of pyramidal (PY) cells in the primary somatosensory cortex (S1) is achieved based on microstimulation of thalamic relay cells through multiple-input multiple-output (MIMO) feedback controllers. This closed loop feedback control mechanism is achieved by simultaneously varying the stimulation parameters across multiple stimulation electrodes in the thalamic circuit based on continuous monitoring of the difference between reference patterns and the evoked responses of the cortical PY cells. We demonstrate that it is feasible to achieve a desired level of performance by controlling the firing activity pattern of a few "key" neural elements in the network. Our results suggest that neural feedback could be an effective method to facilitate the delivery of information to the cortex to substitute lost sensory inputs in cortically controlled BMIs.

Millisecond-timescale local network coding in the rat primary somatosensory cortex.

Correlation among neocortical neurons is thought to play an indispensable role in mediating sensory processing of external stimuli. The role of temporal precision in this correlation has been hypothesized to enhance information flow along sensory pathways. Its role in mediating the integration of information at the output of these pathways, however, remains poorly understood. Here, we examined spike timing correlation between simultaneously recorded layer V neurons within and across columns of the primary somatosensory cortex of anesthetized rats during unilateral whisker stimulation. We used bayesian statistics and information theory to quantify the causal influence between the recorded cells with millisecond precision. For each stimulated whisker, we inferred stable, whisker-specific, dynamic bayesian networks over many repeated trials, with network similarity of 83.3±6% within whisker, compared to only 50.3±18% across whiskers. These networks further provided information about whisker identity that was approximately 6 times higher than what was provided by the latency to first spike and 13 times higher than what was provided by the spike count of individual neurons examined separately. Furthermore, prediction of individual neurons' precise firing conditioned on knowledge of putative pre-synaptic cell firing was 3 times higher than predictions conditioned on stimulus onset alone. Taken together, these results suggest the presence of a temporally precise network coding mechanism that integrates information across neighboring columns within layer V about vibrissa position and whisking kinetics to mediate whisker movement by motor areas innervated by layer V.

Model-based analysis and control of a network of basal ganglia spiking neurons in the normal and parkinsonian states.

Controlling the spatiotemporal firing pattern of an intricately connected network of neurons through microstimulation is highly desirable in many applications. We investigated in this paper the feasibility of using a model-based approach to the analysis and control of a basal ganglia (BG) network model of Hodgkin-Huxley (HH) spiking neurons through microstimulation. Detailed analysis of this network model suggests that it can reproduce the experimentally observed characteristics of BG neurons under a normal and a pathological Parkinsonian state. A simplified neuronal firing rate model, identified from the detailed HH network model, is shown to capture the essential network dynamics. Mathematical analysis of the simplified model reveals the presence of a systematic relationship between the network's structure and its dynamic response to spatiotemporally patterned microstimulation. We show that both the network synaptic organization and the local mechanism of microstimulation can impose tight constraints on the possible spatiotemporal firing patterns that can be generated by the microstimulated network, which may hinder the effectiveness of microstimulation to achieve a desired objective under certain conditions. Finally, we demonstrate that the feedback control design aided by the mathematical analysis of the simplified model is indeed effective in driving the BG network in the normal and Parskinsonian states to follow a prescribed spatiotemporal firing pattern. We further show that the rhythmic/oscillatory patterns that characterize a dopamine-depleted BG network can be suppressed as a direct consequence of controlling the spatiotemporal pattern of a subpopulation of the output Globus Pallidus internalis (GPi) neurons in the network. This work may provide plausible explanations for the mechanisms underlying the therapeutic effects of deep brain stimulation (DBS) in Parkinson's disease and pave the way towards a model-based, network level analysis and closed-loop control and optimization of DBS parameters, among many other applications.

Using brain-computer interfaces to induce neural plasticity and restore function.

Analyzing neural signals and providing feedback in realtime is one of the core characteristics of a brain-computer interface (BCI). As this feature may be employed to induce neural plasticity, utilizing BCI technology for therapeutic purposes is increasingly gaining popularity in the BCI community. In this paper, we discuss the state-of-the-art of research on this topic, address the principles of and challenges in inducing neural plasticity by means of a BCI, and delineate the problems of study design and outcome evaluation arising in this context. We conclude with a list of open questions and recommendations for future research in this field.

Evaluation of a closed loop inductive power transmission system on an awake behaving animal subject.

This paper presents in vivo experimental results for a closed loop wireless power transmission system to implantable devices on an awake behaving animal subject. In this system, wireless power transmission takes place across an inductive link, controlled by a commercial off-the-shelf (COTS) radio frequency identification (RFID) transceiver (TRF7960) operating at 13.56 MHz. Induced voltage on the implantable secondary coil is rectified, digitized by a 10-bit analog to digital converter, and transmitted back to the primary via back telemetry. Transmitter (Tx) and receiver (Rx) circuitry were mounted on the back of an adult rat with a nominal distance of ~7 mm between their coils. Our experiments showed that the closed loop system was able to maintain the Rx supply voltage at the designated 3.8 V despite changes in the coils' relative distance and alignment due to animal movements. The Tx power consumption changed between 410 ~ 560 mW in order to deliver 27 mW to the receiver. The open loop system, on the other hand, showed undesired changes in the Rx supply voltage while the Tx power consumption was constant at 660 mW.