Unsupervised Bayesian Ising Approximation for decoding neural activity and other biological dictionaries
Abstract
The problem of deciphering how low-level patterns (action potentials in the brain, amino acids in a protein, etc.) drive high-level biological features (sensorimotor behavior, enzymatic function) represents the central challenge of quantitative biology. The lack of general methods for doing so from the size of datasets that can be collected experimentally severely limits our understanding of the biological world. For example, in neuroscience, some sensory and motor codes have been shown to consist of precisely timed multi-spike patterns. However, the combinatorial complexity of such pattern codes have precluded development of methods for their comprehensive analysis. Thus, just as it is hard to predict a protein’s function based on its sequence, we still do not understand how to accurately predict an organism's behavior based on neural activity. Here we introduce the unsupervised Bayesian Ising Approximation (uBIA) for solving this class of problems. We demonstrate its utility in an application to neural data, detecting precisely timed spike patterns that code for specific motor behaviors in a songbird vocal system. In data recorded during singing from neurons in a vocal control region, our method detects such codewords with an arbitrary number of spikes, does so from small data sets, and accounts for dependencies in occurrences of codewords. Detecting such comprehensive motor control dictionaries can improve our understanding of skilled motor control and the neural bases of sensorimotor learning in animals. To further illustrate the utility of uBIA, used it to identify the distinct sets of activity patterns that encode vocal motor exploration versus typical song production. Crucially, our method can be used not only for analysis of neural systems, but also for understanding the structure of correlations in other biological and nonbiological datasets.
Data availability
The software implementation of uBIA is available from https://github.com/dghernandez/decomotor. The data used in this work is availablefrom https://figshare.com/articles/Songbird_premotor_dictionaries/10315844.
-
Bengalese finch: Spike timings and acoustic measurements for all casesPLoS Biol, doi.org/10.1371/journal.pbio.1002018.s004.
Article and author information
Author details
Funding
National Institutes of Health (R01-EB022872)
- Damián G Hernández
- Samuel J Sober
- Ilya Nemenman
Simons Foundation (Simons Investigator in MPS)
- Ilya Nemenman
National Institutes of Health (R01-NS084844)
- Samuel J Sober
- Ilya Nemenman
National Institutes of Health (R01-NS099375)
- Damián G Hernández
- Samuel J Sober
- Ilya Nemenman
National Science Foundation (BCS-1822677 (CRCNS Program))
- Samuel J Sober
- Ilya Nemenman
KITP ((NSF) PHY-1748958)
- Ilya Nemenman
KITP ((NIH) R25GM067110)
- Ilya Nemenman
KITP ((Gordon and Betty Moore) 2919.01)
- Ilya Nemenman
Aspen Center for Physics ((NSF) PHY-1607611)
- Samuel J Sober
- Ilya Nemenman
Simons Foundation (The Simons-Emory International Consortium on Motor Control)
- Samuel J Sober
- Ilya Nemenman
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Hernández et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,145
- views
-
- 226
- downloads
-
- 8
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Sleep loss increases AMPA-synaptic strength and number in the neocortex. However, this is only part of the synaptic sleep loss response. We report an increased AMPA/NMDA EPSC ratio in frontal-cortical pyramidal neurons of layers 2–3. Silent synapses are absent, decreasing the plastic potential to convert silent NMDA to active AMPA synapses. These sleep loss changes are recovered by sleep. Sleep genes are enriched for synaptic shaping cellular components controlling glutamate synapse phenotype, overlap with autism risk genes, and are primarily observed in excitatory pyramidal neurons projecting intra-telencephalically. These genes are enriched with genes controlled by the transcription factor, MEF2c, and its repressor, HDAC4. Sleep genes can thus provide a framework within which motor learning and training occur mediated by the sleep-dependent oscillation of glutamate-synaptic phenotypes.
-
- Neuroscience
The amygdala is a subcortical region in the mesiotemporal lobe that plays a key role in emotional and sensory functions. Conventional neuroimaging experiments treat this structure as a single, uniform entity, but there is ample histological evidence for subregional heterogeneity in microstructure and function. The current study characterized subregional structure-function coupling in the human amygdala, integrating post-mortem histology and in vivo MRI at ultra-high fields. Core to our work was a novel neuroinformatics approach that leveraged multiscale texture analysis as well as non-linear dimensionality reduction techniques to identify salient dimensions of microstructural variation in a 3D post-mortem histological reconstruction of the human amygdala. We observed two axes of subregional variation in this region, describing inferior-superior as well as mediolateral trends in microstructural differentiation that in part recapitulated established atlases of amygdala subnuclei. Translating our approach to in vivo MRI data acquired at 7 Tesla, we could demonstrate the generalizability of these spatial trends across 10 healthy adults. We then cross-referenced microstructural axes with functional blood-oxygen-level dependent (BOLD) signal analysis obtained during task-free conditions, and revealed a close association of structural axes with macroscale functional network embedding, notably the temporo-limbic, default mode, and sensory-motor networks. Our novel multiscale approach consolidates descriptions of amygdala anatomy and function obtained from histological and in vivo imaging techniques.