Unsupervised Bayesian Ising Approximation for decoding neural activity and other biological dictionaries

  1. Damián G Hernández  Is a corresponding author
  2. Samuel J Sober  Is a corresponding author
  3. Ilya Nemenman  Is a corresponding author
  1. Centro Atómico Bariloche and Instituto Balseiro, Argentina
  2. Emory University, United States

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.

The following previously published data sets were used

Article and author information

Author details

  1. Damián G Hernández

    Department of Medical Physics, Centro Atómico Bariloche and Instituto Balseiro, Bariloche, Argentina
    For correspondence
    damian.g.h.l@gmail.com
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8995-7495
  2. Samuel J Sober

    Department of Biology, Emory University, Atlanta, United States
    For correspondence
    samuel.j.sober@emory.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1140-7469
  3. Ilya Nemenman

    Department of Physics, Emory University, Atlanta, United States
    For correspondence
    ilya.nemenman@emory.edu
    Competing interests
    The authors declare that no competing interests exist.

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

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Damián G Hernández
  2. Samuel J Sober
  3. Ilya Nemenman
(2022)
Unsupervised Bayesian Ising Approximation for decoding neural activity and other biological dictionaries
eLife 11:e68192.
https://doi.org/10.7554/eLife.68192

Share this article

https://doi.org/10.7554/eLife.68192

Further reading

    1. Neuroscience
    Kaspar E Vogt, Ashwinikumar Kulkarni ... Robert W Greene
    Research Article

    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.

    1. Neuroscience
    Hans Auer, Donna Gift Cabalo ... Jessica Royer
    Research Article

    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.