GraFT: Graph Filtered Temporal Dictionary Learning for Functional Neural Imaging

Long-wavelength traveling waves of vasomotion modulate the perfusion of cortex

Brain arterioles are active, multicellular complexes whose diameters oscillate at ∼ 0.1 Hz. We assess the physiological impact and spatiotemporal dynamics of vaso-oscillations in the awake mouse. First, vaso-oscillations in penetrating arterioles, which source blood from pial arterioles to the capillary bed, profoundly impact perfusion throughout neocortex. The modulation in flux during resting-state activity exceeds that of stimulus-induced activity. Second, the change in perfusion through arterioles relative to the change in their diameter is weak. This implies that the capillary bed dominates the hydrodynamic resistance of brain vasculature. Lastly, the phase of vaso-oscillations evolves slowly along arterioles, with a wavelength that exceeds the span of the cortical mantle and sufficient variability to establish functional cortical areas as parcels of uniform phase. The phase-gradient supports traveling waves in either direction along both pial and penetrating arterioles. This implies that waves along penetrating arterioles can mix, but not directionally transport, interstitial fluids.

Chronic brain functional ultrasound imaging in freely moving rodents performing cognitive tasks

BACKGROUND

Functional ultrasound imaging (fUS) is an emerging imaging technique that indirectly measures neural activity via changes in blood volume. Chronic fUS imaging during cognitive tasks in freely moving animals faces multiple exceptional challenges: performing large durable craniotomies with chronic implants, designing behavioral experiments matching the hemodynamic timescale, stabilizing the ultrasound probe during freely moving behavior, accurately assessing motion artifacts, and validating that the animal can perform cognitive tasks while tethered.

NEW METHOD

We provide validated solutions for those technical challenges. In addition, we present standardized step-by-step reproducible protocols, procedures, and data processing pipelines. Finally, we present proof-of-concept analysis of brain dynamics during a decision making task.

RESULTS

We obtain stable recordings from which we can robustly decode task variables from fUS data over multiple months. Moreover, we find that brain wide imaging through hemodynamic response is nonlinearly related to cognitive variables, such as task difficulty, as compared to sensory responses previously explored.

COMPARISON WITH EXISTING METHODS

Computational pipelines in fUS are nascent and we present an initial development of a full processing pathway to correct and segment fUS data.

CONCLUSIONS

Our methods provide stable imaging and analysis of behavior with fUS that will enable new experimental paradigms in understanding brain-wide dynamics in naturalistic behaviors.

Identifying Interpretable Latent Factors with Sparse Component Analysis

In many neural populations, the computationally relevant signals are posited to be a set of 'latent factors' - signals shared across many individual neurons. Understanding the relationship between neural activity and behavior requires the identification of factors that reflect distinct computational roles. Methods for identifying such factors typically require supervision, which can be suboptimal if one is unsure how (or whether) factors can be grouped into distinct, meaningful sets. Here, we introduce Sparse Component Analysis (SCA), an unsupervised method that identifies interpretable latent factors. SCA seeks factors that are sparse in time and occupy orthogonal dimensions. With these simple constraints, SCA facilitates surprisingly clear parcellations of neural activity across a range of behaviors. We applied SCA to motor cortex activity from reaching and cycling monkeys, single-trial imaging data from C. elegans, and activity from a multitask artificial network. SCA consistently identified sets of factors that were useful in describing network computations.

maskNMF: A denoise-sparsen-detect approach for extracting neural signals from dense imaging data

A number of calcium imaging methods have been developed to monitor the activity of large populations of neurons. One particularly promising approach, Bessel imaging, captures neural activity from a volume by projecting within the imaged volume onto a single imaging plane, therefore effectively mixing signals and increasing the number of neurons imaged per pixel. These signals must then be computationally demixed to recover the desired neural activity. Unfortunately, currently-available demixing methods can perform poorly in the regime of high imaging density (i.e., many neurons per pixel). In this work we introduce a new pipeline (maskNMF) for demixing dense calcium imaging data. The main idea is to first denoise and temporally sparsen the observed video; this enhances signal strength and reduces spatial overlap significantly. Next we detect neurons in the sparsened video using a neural network trained on a library of neural shapes. These shapes are derived from segmented electron microscopy images input into a Bessel imaging model; therefore no manual selection of "good" neural shapes from the functional data is required here. After cells are detected, we use a constrained non-negative matrix factorization approach to demix the activity, using the detected cells' shapes to initialize the factorization. We test the resulting pipeline on both simulated and real datasets and find that it is able to achieve accurate demixing on denser data than was previously feasible, therefore enabling faithful imaging of larger neural populations. The method also provides good results on more "standard" two-photon imaging data. Finally, because much of the pipeline operates on a significantly compressed version of the raw data and is highly parallelizable, the algorithm is fast, processing large datasets faster than real time.

Voxel-Wise Brain Graphs From Diffusion MRI: Intrinsic Eigenspace Dimensionality and Application to Functional MRI

Goal: Structural brain graphs are conventionally limited to defining nodes as gray matter regions from an atlas, with edges reflecting the density of axonal projections between pairs of nodes. Here we explicitly model the entire set of voxels within a brain mask as nodes of high-resolution, subject-specific graphs. Methods: We define the strength of local voxel-to-voxel connections using diffusion tensors and orientation distribution functions derived from diffusion MRI data. We study the graphs' Laplacian spectral properties on data from the Human Connectome Project. We then assess the extent of inter-subject variability of the Laplacian eigenmodes via a procrustes validation scheme. Finally, we demonstrate the extent to which functional MRI data are shaped by the underlying anatomical structure via graph signal processing. Results: The graph Laplacian eigenmodes manifest highly resolved spatial profiles, reflecting distributed patterns that correspond to major white matter pathways. We show that the intrinsic dimensionality of the eigenspace of such high-resolution graphs is only a mere fraction of the graph dimensions. By projecting task and resting-state data on low-frequency graph Laplacian eigenmodes, we show that brain activity can be well approximated by a small subset of low-frequency components. Conclusions: The proposed graphs open new avenues in studying the brain, be it, by exploring their organisational properties via graph or spectral graph theory, or by treating them as the scaffold on which brain function is observed at the individual level.

Review of data processing of functional optical microscopy for neuroscience

Functional optical imaging in neuroscience is rapidly growing with the development of optical systems and fluorescence indicators. To realize the potential of these massive spatiotemporal datasets for relating neuronal activity to behavior and stimuli and uncovering local circuits in the brain, accurate automated processing is increasingly essential. We cover recent computational developments in the full data processing pipeline of functional optical microscopy for neuroscience data and discuss ongoing and emerging challenges.