Multiparametric Brain Hemodynamics Imaging Using a Combined Ultrafast Ultrasound and Photoacoustic System.

Studying brain-wide hemodynamic responses to different stimuli at high spatiotemporal resolutions can help gain new insights into the mechanisms of neuro- diseases and -disorders. Nonetheless, this task is challenging, primarily due to the complexity of neurovascular coupling, which encompasses interdependent hemodynamic parameters including cerebral blood volume (CBV), cerebral blood flow (CBF), and cerebral oxygen saturation (SO2). The current brain imaging technologies exhibit inherent limitations in resolution, sensitivity, and imaging depth, restricting their capacity to comprehensively capture the intricacies of cerebral functions. To address this, a multimodal functional ultrasound and photoacoustic (fUSPA) imaging platform is reported, which integrates ultrafast ultrasound and multispectral photoacoustic imaging methods in a compact head-mountable device, to quantitatively map individual dynamics of CBV, CBF, and SO2 as well as contrast agent enhanced brain imaging at high spatiotemporal resolutions. Following systematic characterization, the fUSPA system is applied to study brain-wide cerebrovascular reactivity (CVR) at single-vessel resolution via relative changes in CBV, CBF, and SO2 in response to hypercapnia stimulation. These results show that cortical veins and arteries exhibit differences in CVR in the stimulated state and consistent anti-correlation in CBV oscillations during the resting state, demonstrating the multiparametric fUSPA system's unique capabilities in investigating complex mechanisms of brain functions.

Dissecting Multiparametric Cerebral Hemodynamics using Integrated Ultrafast Ultrasound and Multispectral Photoacoustic Imaging.

Understanding brain-wide hemodynamic responses to different stimuli at high spatiotemporal resolutions can help study neuro-disorders and brain functions. However, the existing brain imaging technologies have limited resolution, sensitivity, imaging depth and provide information about only one or two hemodynamic parameters. To address this, we propose a multimodal functional ultrasound and photoacoustic (fUSPA) imaging platform, which integrates ultrafast ultrasound and multispectral photoacoustic imaging methods in a compact head-mountable device, to quantitatively map cerebral blood volume (CBV), cerebral blood flow (CBF), oxygen saturation (SO2) dynamics as well as contrast agent enhanced brain imaging with high spatiotemporal resolutions. After systematic characterization, the fUSPA system was applied to quantitatively study the changes in brain hemodynamics and vascular reactivity at single vessel resolution in response to hypercapnia stimulation. Our results show an overall increase in brain-wide CBV, CBF, and SO2, but regional differences in singular cortical veins and arteries and a reproducible anti-correlation pattern between venous and cortical hemodynamics, demonstrating the capabilities of the fUSPA system for providing multiparametric cerebrovascular information at high-resolution and sensitivity, that can bring insights into the complex mechanisms of neurodiseases.

Inferring Insulin Secretion Rate from Sparse Patient Glucose and Insulin Measures.

The insulin secretion rate (ISR) contains information that can provide a personal, quantitative understanding of endocrine function. If the ISR can be reliably inferred from measurements, it could be used for understanding and clinically diagnosing problems with the glucose regulation system. Objective: This study aims to develop a model-based method for inferring a parametrization of the ISR and related physiological information among people with different glycemic conditions in a robust manner. The developed algorithm is applicable for both dense or sparsely sampled plasma glucose/insulin measurements, where sparseness is defined in terms of sampling time with respect to the fastest time scale of the dynamics. Methods: An algorithm for parametrizing and validating a functional form of the ISR for different compartmental models with unknown but estimable ISR function and absorption/decay rates describing the dynamics of insulin accumulation was developed. The method and modeling applies equally to c-peptide secretion rate (CSR) when c-peptide is measured. Accuracy of fit is reliant on reconstruction error of the measured trajectories, and when c-peptide is measured the relationship between CSR and ISR. The algorithm was applied to data from 17 subjects with normal glucose regulatory systems and 9 subjects with cystic fibrosis related diabetes (CFRD) in which glucose, insulin and c-peptide were measured in course of oral glucose tolerance tests (OGTT). Results: This model-based algorithm inferred parametrization of the ISR and CSR functional with relatively low reconstruction error for 12 of 17 control and 7 of 9 CFRD subjects. We demonstrate that when there are suspect measurements points, the validity of excluding them may be interrogated with this method. Significance: A new estimation method is available to infer the ISR and CSR functional profile along with plasma insulin and c-peptide absorption rates from sparse measurements of insulin, c-peptide, and plasma glucose concentrations. We propose a method to interrogate and exclude potentially erroneous OGTT measurement points based on reconstruction errors.

Optimization of an unscented Kalman filter for an embedded platform.

The unscented Kalman filter (UKF) is finding increased application in biological fields. While realizing a complex UKF system in a low-power embedded platform offers many potential benefits including wearability, it also poses significant design challenges. Here we present a method for optimizing a UKF system for realization in an embedded platform. The method seeks to minimize both computation time and error in UKF state reconstruction and forecasting. As a case study, we applied the method to a model for the rat sleep-wake regulatory system in which 432 variants of the UKF over six different variables are considered. The optimization method is divided into three stages that assess computation time, state forecast error, and state reconstruction error. We apply a cost function to variants that pass all three stages to identify a variant that computes 27 times faster than the reference variant and maintains required levels of state estimation and forecasting accuracy. We draw the following insights: 1) process noise provides leeway for simplifying the model and its integration in ways that speed computation time while maintaining state forecasting accuracy, 2) the assimilation of observed data during the UKF correction step provides leeway for simplifying the UKF structure in ways that speed computation time while maintaining state reconstruction accuracy, and 3) the optimization process can be accelerated by decoupling variables that directly impact the underlying model from variables that impact the UKF structure.

Awake mouse brain photoacoustic and optical imaging through a transparent ultrasound cranial window.

Optical resolution photoacoustic microscopy (OR-PAM) can map the cerebral vasculature at capillary-level resolution. However, the OR-PAM setup's bulky imaging head makes awake mouse brain imaging challenging and inhibits its integration with other optical neuroimaging modalities. Moreover, the glass cranial windows used for optical microscopy are unsuitable for OR-PAM due to the acoustic impedance mismatch between the glass plate and the tissue. To overcome these challenges, we propose a lithium niobate based transparent ultrasound transducer (TUT) as a cranial window on a thinned mouse skull. The TUT cranial window simplifies the imaging head considerably due to its dual functionality as an optical window and ultrasound transducer. The window remains stable for six weeks, with no noticeable inflammation and minimal bone regrowth. The TUT window's potential is demonstrated by imaging the awake mouse cerebral vasculature using OR-PAM, intrinsic optical signal imaging, and two-photon microscopy. The TUT cranial window can potentially also be used for ultrasound stimulation and simultaneous multimodal imaging of the awake mouse brain.

Development of Mechanistic Neural Mass (mNM) Models that Link Physiology to Mean-Field Dynamics.

Brain rhythms emerge from the mean-field activity of networks of neurons. There have been many efforts to build mathematical and computational embodiments in the form of discrete cell-group activities-termed neural masses-to understand in particular the origins of evoked potentials, intrinsic patterns of activities such as theta, regulation of sleep, Parkinson's disease related dynamics, and mimic seizure dynamics. As originally utilized, standard neural masses convert input through a sigmoidal function to a firing rate, and firing rate through a synaptic alpha function to other masses. Here we define a process to build mechanistic neural masses (mNMs) as mean-field models of microscopic membrane-type (Hodgkin Huxley type) models of different neuron types that duplicate the stability, firing rate, and associated bifurcations as function of relevant slow variables - such as extracellular potassium - and synaptic current; and whose output is both firing rate and impact on the slow variables - such as transmembrane potassium flux. Small networks composed of just excitatory and inhibitory mNMs demonstrate expected dynamical states including firing, runaway excitation and depolarization block, and these transitions change in biologically observed ways with changes in extracellular potassium and excitatory-inhibitory balance.

Proceedings of the Sleep and Epilepsy Workshop: Section 3 Mortality: Sleep, Night, and SUDEP.

Sudden unexpected death in epilepsy (SUDEP) is the leading cause of death in patients with refractory epilepsy. Likely pathophysiological mechanisms include seizure-induced cardiac and respiratory dysregulation. A frequently identified feature in SUDEP cases is that they occur at night. This raises the question of a role for sleep state in regulating of SUDEP. An association with sleep has been identified in a number of studies with patients and in animal models. The focus of this section of the Sleep and Epilepsy Workshop was on identifying and understanding the role for sleep and time of day in the pathophysiology of SUDEP.

Proceedings of the Sleep and Epilepsy Workshop: Section 1 Decreasing Seizures-Improving Sleep and Seizures, Themes for Future Research.

Epileptic seizures, sleep, and circadian timing share bilateral interactions, but concerted work to characterize these interactions and to leverage them to the advantage of patients with epilepsy remains in beginning stages. To further the field, a multidisciplinary group of sleep physicians, epileptologists, circadian timing experts, and others met to outline the state of the art, gaps of knowledge, and suggest ways forward in clinical, translational, and basic research. A multidisciplinary panel of experts discussed these interactions, centered on whether improvements in sleep or circadian rhythms improve decrease seizure frequency. In addition, education about sleep was lacking in among patients, their families, and physicians, and that focus on education was an extremely important "low hanging fruit" to harvest. Improvements in monitoring technology, experimental designs sensitive to the rigor required to dissect sleep versus circadian influences, and clinical trials in seizure reduction with sleep improvements were appropriate.

Functional hyperemia drives fluid exchange in the paravascular space.

The brain lacks a conventional lymphatic system to remove metabolic waste. It has been proposed that directional fluid movement through the arteriolar paravascular space (PVS) promotes metabolite clearance. We performed simulations to examine if arteriolar pulsations and dilations can drive directional CSF flow in the PVS and found that arteriolar wall movements do not drive directional CSF flow. We propose an alternative method of metabolite clearance from the PVS, namely fluid exchange between the PVS and the subarachnoid space (SAS). In simulations with compliant brain tissue, arteriolar pulsations did not drive appreciable fluid exchange between the PVS and the SAS. However, when the arteriole dilated, as seen during functional hyperemia, there was a marked exchange of fluid. Simulations suggest that functional hyperemia may serve to increase metabolite clearance from the PVS. We measured blood vessels and brain tissue displacement simultaneously in awake, head-fixed mice using two-photon microscopy. These measurements showed that brain deforms in response to pressure changes in PVS, consistent with our simulations. Our results show that the deformability of the brain tissue needs to be accounted for when studying fluid flow and metabolite transport.

A Brain-Heart Biomarker for Epileptogenesis.

Postinjury epilepsy is an potentially preventable sequela in as many as 20% of patients with brain insults. For these cases biomarkers of epileptogenesis are critical to facilitate identification of patients at high-risk of developing epilepsy and to introduce effective anti-epileptogenic interventions. Here, we demonstrate that delayed brain-heart coincidences serve as a reliable biomarker. In a murine model of post-infection acquired epilepsy, we used long-term simultaneous measurements of the brain activity via electroencephalography and autonomic cardiac activity via electrocardiography, in male mice, to quantitatively track brain-heart interactions during epileptogenesis. We find that abnormal cortical discharges precede abnormal fluctuations in the cardiac rhythm at the resolution of single beat-to-beat intervals. The delayed brain-heart coincidence is detectable as early as the onset of chronic measurements, 2-14 weeks before the first seizure, only in animals that become epileptic, and increases during epileptogenesis. Therefore, delayed brain-heart coincidence serves as a biomarker of epileptogenesis and could be used for phenotyping, diagnostic, and therapeutic purposes.SIGNIFICANCE STATEMENT No biomarker that readily predicts and tracks epileptogenesis currently exists for the wide range of human acquired epilepsies. Here, we used long-term measurements of brain and heart activity in a mouse model of post-infection acquired epilepsy to investigate the potential of brain-heart interaction as a biomarker of epileptogenesis. We found that delayed coincidences from brain to heart can clearly separate the mice that became epileptic from those that did not weeks before development of epilepsy. Our findings allow for phenotyping and tracking of epileptogenesis in this and likely other models of acquired epilepsy. Such capability is critical for efficient adjunctive treatment development and for tracking the efficacy of such treatments.

Control of Spreading Depression with Electrical Fields.

Spreading depression or depolarization is a large-scale pathological brain phenomenon related to migraine, stroke, hemorrhage and traumatic brain injury. Once initiated, spreading depression propagates across gray matter extruding potassium and other active molecules, collapsing the resting membrane electro-chemical gradient of cells leading to spike inactivation and cellular swelling, and propagates independently of synaptic transmission. We demonstrate the modulation, suppression and prevention of spreading depression utilizing applied transcortical DC electric fields in brain slices, measured with intrinsic optical imaging and potassium dye epifluorescence. We experimentally observe a surface-positive electric field induced forcing of spreading depression propagation to locations in cortex deeper than the unmodulated propagation path, whereby further propagation is confined and arrested even after field termination. The opposite surface-negative electric field polarity produces an increase in propagation velocity and a confinement of the wave to more superficial layers of cortex than the unmodulated propagation path. These field polarities are of opposite sign to the polarity that blocks neuronal spiking and seizures, and are consistent with biophysical models of spreading depression. The results demonstrate the potential feasibility of electrical control and prevention of spreading depression.

The systemDrive: a Multisite, Multiregion Microdrive with Independent Drive Axis Angling for Chronic Multimodal Systems Neuroscience Recordings in Freely Behaving Animals.

A multielectrode system that can address widely separated targets at multiple sites across multiple brain regions with independent implant angling is needed to investigate neural function and signaling in systems and circuits of small animals. Here, we present the systemDrive, a novel multisite, multiregion microdrive that is capable of moving microwire electrode bundles into targets along independent and nonparallel drive trajectories. Our design decouples the stereotaxic surgical placement of individual guide cannulas for each trajectory from the placement of a flexible drive structure. This separation enables placement of many microwire multitrodes along widely spaced and independent drive axes with user-set electrode trajectories and depths from a single microdrive body, and achieves stereotaxic precision with each. The system leverages tight tube-cannula tolerances and geometric constraints on flexible drive axes to ensure concentric alignment of electrode bundles within guide cannulas. Additionally, the headmount and microdrive both have an open-center design to allow for the placement of additional sensing modalities. This design is the first, in the context of small rodent chronic research, to provide the capability to finely position microwires through multiple widely distributed cell groups, each with stereotaxic precision, along arbitrary and nonparallel trajectories that are not restricted to emanate from a single source. We demonstrate the use of the systemDrive in male Long-Evans rats to observe simultaneous single-unit and multiunit activity from multiple widely separated sleep-wake regulatory brainstem cell groups, along with cortical and hippocampal activity, during free behavior over multiple many-day continuous recording periods.

A Murine Model to Study Epilepsy and SUDEP Induced by Malaria Infection.

One of the largest single sources of epilepsy in the world is produced as a neurological sequela in survivors of cerebral malaria. Nevertheless, the pathophysiological mechanisms of such epileptogenesis remain unknown and no adjunctive therapy during cerebral malaria has been shown to reduce the rate of subsequent epilepsy. There is no existing animal model of postmalarial epilepsy. In this technical report we demonstrate the first such animal models. These models were created from multiple mouse and parasite strain combinations, so that the epilepsy observed retained universality with respect to genetic background. We also discovered spontaneous sudden unexpected death in epilepsy (SUDEP) in two of our strain combinations. These models offer a platform to enable new preclinical research into mechanisms and prevention of epilepsy and SUDEP.

The Neural Basis for Sleep Regulation - Data Assimilation from Animal to Model.

Sleep is important for normal brain function, and sleep disruption is comorbid with many neurological diseases. There is a growing mechanistic understanding of the neurological basis for sleep regulation that is beginning to lead to mechanistic mathematically described models. It is our objective to validate the predictive capacity of such models using data assimilation (DA) methods. If such methods are successful, and the models accurately describe enough of the mechanistic functions of the physical system, then they can be used as sophisticated observation systems to reveal both system changes and sources of dysfunction with neurological diseases and identify routes to intervene. Here we report on extensions to our initial efforts [1] at applying unscented Kalman Filter (UKF) to models of sleep regulation on three fronts: tools for multi-parameter fitting; a sophisticated observation model to apply the UKF for observations of behavioral state; and comparison with data recorded from brainstem cell groups thought to regulate sleep.

Statistical evaluation of forecasts.

Reliable forecasts of extreme but rare events, such as earthquakes, financial crashes, and epileptic seizures, would render interventions and precautions possible. Therefore, forecasting methods have been developed which intend to raise an alarm if an extreme event is about to occur. In order to statistically validate the performance of a prediction system, it must be compared to the performance of a random predictor, which raises alarms independent of the events. Such a random predictor can be obtained by bootstrapping or analytically. We propose an analytic statistical framework which, in contrast to conventional methods, allows for validating independently the sensitivity and specificity of a forecasting method. Moreover, our method accounts for the periods during which an event has to remain absent or occur after a respective forecast.

Rapid eye movement sleep and hippocampal theta oscillations precede seizure onset in the tetanus toxin model of temporal lobe epilepsy.

Improved understanding of the interaction between state of vigilance (SOV) and seizure onset has therapeutic potential. Six rats received injections of tetanus toxin (TeTX) in the ventral hippocampus that resulted in chronic spontaneous seizures. The distribution of SOV before 486 seizures was analyzed for a total of 19 d of recording. Rapid eye movement sleep (REM) and exploratory wake, both of which express prominent hippocampal theta rhythm, preceded 47 and 34%, for a total of 81%, of all seizures. Nonrapid eye movement sleep (NREM) and nonexploratory wake, neither of which expresses prominent theta, preceded 6.8 and 13% of seizures. We demonstrate that identification of SOV yields significant differentiation of seizure susceptibilities, with the instantaneous seizure rate during REM nearly 10 times higher than baseline and the rate for NREM less than half of baseline. Survival analysis indicated a shorter duration of preseizure REM bouts, with a maximum transition to seizure at ∼90 s after the onset of REM. This study provides the first analysis of a correlation between SOV and seizure onset in the TeTX model of temporal lobe epilepsy, as well as the first demonstration that hippocampal theta rhythms associated with natural behavioral states can serve a seizure-promoting role. Our findings are in contrast with previous studies suggesting that the correlations between SOV and seizures are primarily governed by circadian oscillations and the notion that hippocampal theta rhythms inhibit seizures. The documentation of significant SOV-dependent seizure susceptibilities indicates the potential utility of SOV and its time course in seizure prediction and control.

Frequency dependence of behavioral modulation by hippocampal electrical stimulation.

Electrical stimulation offers the potential to develop novel strategies for the treatment of refractory medial temporal lobe epilepsy. In particular, direct electrical stimulation of the hippocampus presents the opportunity to modulate pathological dynamics at the ictal focus, although the neuroanatomical substrate of this region renders it susceptible to altering cognition and affective processing as a side effect. We investigated the effects of three electrical stimulation paradigms on separate groups of freely moving rats (sham, 8-Hz and 40-Hz sine-wave stimulation of the ventral/intermediate hippocampus, where 8- and 40-Hz stimulation were chosen to mimic naturally occurring hippocampal oscillations). Animals exhibited attenuated locomotor and exploratory activity upon stimulation at 40 Hz, but not at sham or 8-Hz stimulation. Such behavioral modifications were characterized by a significant reduction in rearing frequency, together with increased freezing behavior. Logistic regression analysis linked the observed changes in animal locomotion to 40-Hz electrical stimulation independently of time-related variables occurring during testing. Spectral analysis, conducted to monitor the electrophysiological profile in the CA1 area of the dorsal hippocampus, showed a significant reduction in peak theta frequency, together with reduced theta power in the 40-Hz vs. the sham stimulation animal group, independent of locomotion speed (theta range: 4-12 Hz). These findings contribute to the development of novel and safe medical protocols by indicating a strategy to constrain or optimize parameters in direct hippocampal electrical stimulation.

Challenges in wearable personal health monitoring systems.

Wearable sensors give the users convenience in daily health monitoring, though several challenges in such sensor systems should be overcome. This paper discusses the challenges in wearable health monitoring sensors and solutions for multi-modal and multi-functional wrist-worn devices based on novel circuit design techniques to reject DC offset. This paper also presents a novel sophisticated algorithm to reject motion artifacts. The system has the capability to simultaneously acquire several biomedical signals (i.e. electrocardiogram, PPG, and body-electrode impedance). The system can also help patients who want to monitor their psychological signals to mitigate health risks.

Modulation of hippocampal rhythms by subthreshold electric fields and network topology.

Theta (4-12 Hz) and gamma (30-80 Hz) rhythms are considered important for cortical and hippocampal function. Although several neuron types are implicated in rhythmogenesis, the exact cellular mechanisms remain unknown. Subthreshold electric fields provide a flexible, area-specific tool to modulate neural activity and directly test functional hypotheses. Here we present experimental and computational evidence of the interplay among hippocampal synaptic circuitry, neuronal morphology, external electric fields, and network activity. Electrophysiological data are used to constrain and validate an anatomically and biophysically realistic model of area CA1 containing pyramidal cells and two interneuron types: dendritic- and perisomatic-targeting. We report two lines of results: addressing the network structure capable of generating theta-modulated gamma rhythms, and demonstrating electric field effects on those rhythms. First, theta-modulated gamma rhythms require specific inhibitory connectivity. In one configuration, GABAergic axo-dendritic feedback on pyramidal cells is only effective in proximal but not distal layers. An alternative configuration requires two distinct perisomatic interneuron classes, one exclusively receiving excitatory contacts, the other additionally targeted by inhibition. These observations suggest novel roles for particular classes of oriens and basket cells. The second major finding is that subthreshold electric fields robustly alter the balance between different rhythms. Independent of network configuration, positive electric fields decrease, while negative fields increase the theta/gamma ratio. Moreover, electric fields differentially affect average theta frequency depending on specific synaptic connectivity. These results support the testable prediction that subthreshold electric fields can alter hippocampal rhythms, suggesting new approaches to explore their cognitive functions and underlying circuitry.