Localized semi-nonnegative matrix factorization (LocaNMF) of widefield calcium imaging data

Widefield calcium imaging enables recording of large-scale neural activity across the mouse dorsal cortex. In order to examine the relationship of these neural signals to the resulting behavior, it is critical to demix the recordings into meaningful spatial and temporal components that can be mapped onto well-defined brain regions. However, no current tools satisfactorily extract the activity of the different brain regions in individual mice in a data-driven manner, while taking into account mouse-specific and preparation-specific differences. Here, we introduce Localized semi-Nonnegative Matrix Factorization (Lo-caNMF), a method that efficiently decomposes widefield video data and allows us to directly compare activity across multiple mice by outputting mouse-specific localized functional regions that are significantly more interpretable than more traditional decomposition techniques. Moreover, it provides a natural subspace to directly compare correlation maps and neural dynamics across different behaviors, mice, and experimental conditions, and enables identification of task- and movement-related brain regions. Author summary While recording from multiple regions of the brain, how does one best incorporate prior information about anatomical regions while accurately representing the data? Here, we introduce Localized semi-NMF (LocaNMF), an algorithm that efficiently decomposes widefield video data into meaningful spatial and temporal components that can be decoded and compared across different behavioral sessions and experimental conditions. Mapping the inferred components onto well-defined brain regions using a widely-used brain atlas provides an interpretable, stable decomposition. LocaNMF allows us to satisfactorily extract the activity of the different brain regions in individual mice in a data-driven manner, while taking into account mouse-specific and preparation-specific differences.

[1]  Xilin Shen,et al.  Simultaneous mesoscopic and two-photon imaging of neuronal activity in cortical circuits , 2018, Nature Methods.

[2]  Amy Hu,et al.  Thy1 transgenic mice expressing the red fluorescent calcium indicator jRGECO1a for neuronal population imaging in vivo , 2018, bioRxiv.

[3]  Hongkui Zeng,et al.  Long-Term Optical Access to an Estimated One Million Neurons in the Live Mouse Cortex. , 2016, Cell reports.

[4]  Allan R. Jones,et al.  A mesoscale connectome of the mouse brain , 2014, Nature.

[5]  Thomas R. Barrick,et al.  Atlas-based segmentation of white matter tracts of the human brain using diffusion tensor tractography and comparison with classical dissection , 2008, NeuroImage.

[6]  Kevin M. Cury,et al.  DeepLabCut: markerless pose estimation of user-defined body parts with deep learning , 2018, Nature Neuroscience.

[7]  Allan R. Jones,et al.  Genome-wide atlas of gene expression in the adult mouse brain , 2007, Nature.

[8]  Karl J. Friston,et al.  Unified segmentation , 2005, NeuroImage.

[9]  Denise M. Piscopo,et al.  Large-scale imaging of cortical dynamics during sensory perception and behavior. , 2016, Journal of neurophysiology.

[10]  Kenneth D. Harris,et al.  Coherent encoding of subjective spatial position in visual cortex and hippocampus , 2018, Nature.

[11]  Alan Urban,et al.  Whole-Brain Functional Ultrasound Imaging Reveals Brain Modules for Visuomotor Integration , 2018, Neuron.

[12]  Dustin Scheinost,et al.  Spanning spatiotemporal scales with simultaneous mesoscopic Ca2+ imaging and functional MRI , 2018, bioRxiv.

[13]  G. Blasdel,et al.  Voltage-sensitive dyes reveal a modular organization in monkey striate cortex , 1986, Nature.

[14]  T. Murphy,et al.  Mesoscale Mapping of Mouse Cortex Reveals Frequency-Dependent Cycling between Distinct Macroscale Functional Modules , 2017, The Journal of Neuroscience.

[15]  Dustin Scheinost,et al.  Simultaneous mesoscopic Ca2+ imaging and fMRI: Neuroimaging spanning spatiotemporal scales , 2018 .

[16]  Daniel Rueckert,et al.  Human brain mapping: A systematic comparison of parcellation methods for the human cerebral cortex , 2017, NeuroImage.

[17]  David N. Thibodeaux,et al.  Wide-field optical mapping of neural activity in awake mice and the importance of hemodynamic correction , 2017 .

[18]  B. N. Chatterji,et al.  An FFT-based technique for translation, rotation, and scale-invariant image registration , 1996, IEEE Trans. Image Process..

[19]  Shigeyoshi Itohara,et al.  Excitatory Neuronal Hubs Configure Multisensory Integration of Slow Waves in Association Cortex. , 2018, Cell reports.

[20]  D. McVea,et al.  Mirrored Bilateral Slow-Wave Cortical Activity within Local Circuits Revealed by Fast Bihemispheric Voltage-Sensitive Dye Imaging in Anesthetized and Awake Mice , 2010, The Journal of Neuroscience.

[21]  Mariel G Kozberg,et al.  Resting-state hemodynamics are spatiotemporally coupled to synchronized and symmetric neural activity in excitatory neurons , 2016, Proceedings of the National Academy of Sciences.

[22]  Nicholas A. Steinmetz,et al.  Distinct contributions of mouse cortical areas to visual discrimination , 2018 .

[23]  Kenneth D. Harris,et al.  A perceptual decision requires sensory but not action coding in mouse cortex , 2020 .

[24]  Matthias Kaschube,et al.  Distributed network interactions and their emergence in developing neocortex , 2018, Nature Neuroscience.

[25]  E. Hillman,et al.  Hyperspectral in vivo two-photon microscopy of intrinsic contrast. , 2008, Optics letters.

[26]  Ying Ma,et al.  Wide-field optical mapping of neural activity and brain haemodynamics: considerations and novel approaches , 2016, Philosophical Transactions of the Royal Society B: Biological Sciences.

[27]  Rafael Yuste,et al.  Multi-scale approaches for high-speed imaging and analysis of large neural populations , 2016 .

[28]  Matthew T. Kaufman,et al.  Single-trial neural dynamics are dominated by richly varied movements , 2019, Nature Neuroscience.

[29]  T. Wiesel,et al.  Functional architecture of cortex revealed by optical imaging of intrinsic signals , 1986, Nature.

[30]  Andrzej Cichocki,et al.  Nonnegative Matrix and Tensor Factorization T , 2007 .

[31]  Matthew R Whiteway,et al.  Parallel Processing of Sound Dynamics across Mouse Auditory Cortex via Spatially Patterned Thalamic Inputs and Distinct Areal Intracortical Circuits. , 2019, Cell reports.

[32]  Jesper Andersson,et al.  A multi-modal parcellation of human cerebral cortex , 2016, Nature.

[33]  Cornelia I Bargmann,et al.  The Brain Research Through Advancing Innovative Neurotechnologies (BRAIN) initiative and neurology. , 2014, JAMA neurology.

[34]  Mark W. Woolrich,et al.  Network modelling methods for FMRI , 2011, NeuroImage.

[35]  Naoshige Uchida,et al.  Demixed principal component analysis of neural population data , 2014, eLife.

[36]  Matteo Carandini,et al.  Cortical State Fluctuations during Sensory Decision Making , 2018, Current Biology.

[37]  M. Carandini,et al.  Focal cortical seizures start as standing waves and propagate respecting homotopic connectivity , 2016, Nature Communications.

[38]  Yong Fan,et al.  Large-scale sparse functional networks from resting state fMRI , 2017, NeuroImage.

[39]  Liren Zhu,et al.  High‐resolution deep functional imaging of the whole mouse brain by photoacoustic computed tomography in vivo , 2018, Journal of biophotonics.

[40]  William E. Allen,et al.  Global Representations of Goal-Directed Behavior in Distinct Cell Types of Mouse Neocortex , 2017, Neuron.

[41]  D. Tank,et al.  Imaging Cortical Dynamics in GCaMP Transgenic Rats with a Head-Mounted Widefield Macroscope , 2018, Neuron.

[42]  Ying Ma,et al.  Penalized matrix decomposition for denoising, compression, and improved demixing of functional imaging data , 2018, bioRxiv.

[43]  M. V. D. Heuvel,et al.  Exploring the brain network: A review on resting-state fMRI functional connectivity , 2010, European Neuropsychopharmacology.

[44]  Timothy O. Laumann,et al.  Functional Network Organization of the Human Brain , 2011, Neuron.

[45]  David Pfau,et al.  Simultaneous Denoising, Deconvolution, and Demixing of Calcium Imaging Data , 2016, Neuron.