Automated Detection of Aortic Stenosis Using Machine Learning.

BACKGROUND: Aortic stenosis (AS) is a degenerative valve condition that is underdiagnosed and undertreated. Detection of AS using limited two-dimensional echocardiography could enable screening and improve appropriate referral and treatment of this condition. The aim of this study was to develop methods for automated detection of AS from limited imaging data sets. METHODS: Convolutional neural networks were trained, validated, and tested using limited two-dimensional transthoracic echocardiographic data sets. Networks were developed to accomplish two sequential tasks: (1) view identification and (2) study-level grade of AS. Balanced accuracy and area under the receiver operator curve (AUROC) were the performance metrics used. RESULTS: Annotated images from 577 patients were included. Neural networks were trained on data from 338 patients (average n = 10,253 labeled images), validated on 119 patients (average n = 3,505 labeled images), and performance was assessed on a test set of 120 patients (average n = 3,511 labeled images). Fully automated screening for AS was achieved with an AUROC of 0.96. Networks can distinguish no significant (no, mild, mild to moderate) AS from significant (moderate or severe) AS with an AUROC of 0.86 and between early (mild or mild to moderate AS) and significant (moderate or severe) AS with an AUROC of 0.75. External validation of these networks in a cohort of 8,502 outpatient transthoracic echocardiograms showed that screening for AS can be achieved using parasternal long-axis imaging only with an AUROC of 0.91. CONCLUSION: Fully automated detection of AS using limited two-dimensional data sets is achievable using modern neural networks. These methods lay the groundwork for a novel method for screening for AS.

Designing Patient-Specific Optimal Neurostimulation Patterns for Seizure Suppression.

Neurostimulation is a promising therapy for abating epileptic seizures. However, it is extremely difficult to identify optimal stimulation patterns experimentally. In this study, human recordings are used to develop a functional 24 neuron network statistical model of hippocampal connectivity and dynamics. Spontaneous seizure-like activity is induced in silico in this reconstructed neuronal network. The network is then used as a testbed to design and validate a wide range of neurostimulation patterns. Commonly used periodic trains were not able to permanently abate seizures at any frequency. A simulated annealing global optimization algorithm was then used to identify an optimal stimulation pattern, which successfully abated 92% of seizures. Finally, in a fully responsive, or closed-loop, neurostimulation paradigm, the optimal stimulation successfully prevented the network from entering the seizure state. We propose that the framework presented here for algorithmically identifying patient-specific neurostimulation patterns can greatly increase the efficacy of neurostimulation devices for seizures.

Cannabinoids disrupt memory encoding by functionally isolating hippocampal CA1 from CA3.

Much of the research on cannabinoids (CBs) has focused on their effects at the molecular and synaptic level. However, the effects of CBs on the dynamics of neural circuits remains poorly understood. This study aims to disentangle the effects of CBs on the functional dynamics of the hippocampal Schaffer collateral synapse by using data-driven nonparametric modeling. Multi-unit activity was recorded from rats doing an working memory task in control sessions and under the influence of exogenously administered tetrahydrocannabinol (THC), the primary CB found in marijuana. It was found that THC left firing rate unaltered and only slightly reduced theta oscillations. Multivariate autoregressive models, estimated from spontaneous spiking activity, were then used to describe the dynamical transformation from CA3 to CA1. They revealed that THC served to functionally isolate CA1 from CA3 by reducing feedforward excitation and theta information flow. The functional isolation was compensated by increased feedback excitation within CA1, thus leading to unaltered firing rates. Finally, both of these effects were shown to be correlated with memory impairments in the working memory task. By elucidating the circuit mechanisms of CBs, these results help close the gap in knowledge between the cellular and behavioral effects of CBs.

Nonparametric Model of Smooth Muscle Force Production During Electrical Stimulation.

A nonparametric model of smooth muscle tension response to electrical stimulation was estimated using the Laguerre expansion technique of nonlinear system kernel estimation. The experimental data consisted of force responses of smooth muscle to energy-matched alternating single pulse and burst current stimuli. The burst stimuli led to at least a 10-fold increase in peak force in smooth muscle from Mytilus edulis, despite the constant energy constraint. A linear model did not fit the data. However, a second-order model fit the data accurately, so the higher-order models were not required to fit the data. Results showed that smooth muscle force response is not linearly related to the stimulation power.

Distinguishing cognitive state with multifractal complexity of hippocampal interspike interval sequences.

Fractality, represented as self-similar repeating patterns, is ubiquitous in nature and the brain. Dynamic patterns of hippocampal spike trains are known to exhibit multifractal properties during working memory processing; however, it is unclear whether the multifractal properties inherent to hippocampal spike trains reflect active cognitive processing. To examine this possibility, hippocampal neuronal ensembles were recorded from rats before, during and after a spatial working memory task following administration of tetrahydrocannabinol (THC), a memory-impairing component of cannabis. Multifractal detrended fluctuation analysis was performed on hippocampal interspike interval sequences to determine characteristics of monofractal long-range temporal correlations (LRTCs), quantified by the Hurst exponent, and the degree/magnitude of multifractal complexity, quantified by the width of the singularity spectrum. Our results demonstrate that multifractal firing patterns of hippocampal spike trains are a marker of functional memory processing, as they are more complex during the working memory task and significantly reduced following administration of memory impairing THC doses. Conversely, LRTCs are largest during resting state recordings, therefore reflecting different information compared to multifractality. In order to deepen conceptual understanding of multifractal complexity and LRTCs, these measures were compared to classical methods using hippocampal frequency content and firing variability measures. These results showed that LRTCs, multifractality, and theta rhythm represent independent processes, while delta rhythm correlated with multifractality. Taken together, these results provide a novel perspective on memory function by demonstrating that the multifractal nature of spike trains reflects hippocampal microcircuit activity that can be used to detect and quantify cognitive, physiological, and pathological states.

Hippocampal closed-loop modeling and implications for seizure stimulation design.

OBJECTIVE: Traditional hippocampal modeling has focused on the series of feedforward synapses known as the trisynaptic pathway. However, feedback connections from CA1 back to the hippocampus through the entorhinal cortex (EC) actually make the hippocampus a closed-loop system. By constructing a functional closed-loop model of the hippocampus, one may learn how both physiological and epileptic oscillations emerge and design efficient neurostimulation patterns to abate such oscillations. APPROACH: Point process input-output models where estimated from recorded rodent hippocampal data to describe the nonlinear dynamical transformation from CA3 → CA1, via the schaffer-collateral synapse, and CA1 → CA3 via the EC. Each Volterra-like subsystem was composed of linear dynamics (principal dynamic modes) followed by static nonlinearities. The two subsystems were then wired together to produce the full closed-loop model of the hippocampus. MAIN RESULTS: Closed-loop connectivity was found to be necessary for the emergence of theta resonances as seen in recorded data, thus validating the model. The model was then used to identify frequency parameters for the design of neurostimulation patterns to abate seizures. SIGNIFICANCE: Deep-brain stimulation (DBS) is a new and promising therapy for intractable seizures. Currently, there is no efficient way to determine optimal frequency parameters for DBS, or even whether periodic or broadband stimuli are optimal. Data-based computational models have the potential to be used as a testbed for designing optimal DBS patterns for individual patients. However, in order for these models to be successful they must incorporate the complex closed-loop structure of the seizure focus. This study serves as a proof-of-concept of using such models to design efficient personalized DBS patterns for epilepsy.

Understanding spike-triggered covariance using Wiener theory for receptive field identification.

Receptive field identification is a vital problem in sensory neurophysiology and vision. Much research has been done in identifying the receptive fields of nonlinear neurons whose firing rate is determined by the nonlinear interactions of a small number of linear filters. Despite more advanced methods that have been proposed, spike-triggered covariance (STC) continues to be the most widely used method in such situations due to its simplicity and intuitiveness. Although the connection between STC and Wiener/Volterra kernels has often been mentioned in the literature, this relationship has never been explicitly derived. Here we derive this relationship and show that the STC matrix is actually a modified version of the second-order Wiener kernel, which incorporates the input autocorrelation and mixes first- and second-order dynamics. It is then shown how, with little modification of the STC method, the Wiener kernels may be obtained and, from them, the principal dynamic modes, a set of compact and efficient linear filters that essentially combine the spike-triggered average and STC matrix and generalize to systems with both continuous and point-process outputs. Finally, using Wiener theory, we show how these obtained filters may be corrected when they were estimated using correlated inputs. Our correction technique is shown to be superior to those commonly used in the literature for both correlated Gaussian images and natural images.

System identification of point-process neural systems using probability based Volterra kernels.

BACKGROUND: Neural information processing involves a series of nonlinear dynamical input/output transformations between the spike trains of neurons/neuronal ensembles. Understanding and quantifying these transformations is critical both for understanding neural physiology such as short-term potentiation and for developing cognitive neural prosthetics. NEW METHOD: A novel method for estimating Volterra kernels for systems with point-process inputs and outputs is developed based on elementary probability theory. These Probability Based Volterra (PBV) kernels essentially describe the probability of an output spike given q input spikes at various lags t1, t2, …, tq. RESULTS: The PBV kernels are used to estimate both synthetic systems where ground truth is available and data from the CA3 and CA1 regions rodent hippocampus. The PBV kernels give excellent predictive results in both cases. Furthermore, they are shown to be quite robust to noise and to have good convergence and overfitting properties. Through a slight modification, the PBV kernels are shown to also deal well with correlated point-process inputs. COMPARISON WITH EXISTING METHODS: The PBV kernels were compared with kernels estimated through least squares estimation (LSE) and through the Laguerre expansion technique (LET). The LSE kernels were shown to fair very poorly with real data due to the large amount of input noise. Although the LET kernels gave the best predictive results in all cases, they require prior parameter estimation. It was shown how the PBV and LET methods can be combined synergistically to maximize performance. CONCLUSIONS: The PBV kernels provide a novel and intuitive method of characterizing point-process input-output nonlinear systems.

Model-based asessment of an in-vivo predictive relationship from CA1 to CA3 in the rodent hippocampus.

Although an anatomical connection from CA1 to CA3 via the Entorhinal Cortex (EC) and through backprojecting interneurons has long been known it exist, it has never been examined quantitatively on the single neuron level, in the in-vivo nonpatholgical, nonperturbed brain. Here, single spike activity was recorded using a multi-electrode array from the CA3 and CA1 areas of the rodent hippocampus (N = 7) during a behavioral task. The predictive power from CA3→CA1 and CA1→CA3 was examined by constructing Multivariate Autoregressive (MVAR) models from recorded neurons in both directions. All nonsignificant inputs and models were identified and removed by means of Monte Carlo simulation methods. It was found that 121/166 (73 %) CA3→CA1 models and 96/145 (66 %) CA1→CA3 models had significant predictive power, thus confirming a predictive 'Granger' causal relationship from CA1 to CA3. This relationship is thought to be caused by a combination of truly causal connections such as the CA1→EC→CA3 pathway and common inputs such as those from the Septum. All MVAR models were then examined in the frequency domain and it was found that CA3 kernels had significantly more power in the theta and beta range than those of CA1, confirming CA3's role as an endogenous hippocampal pacemaker.