Characterizing Cortex-Wide Dynamics with Wide-Field Calcium Imaging.

The brain functions through coordinated activity among distributed regions. Wide-field calcium imaging, combined with improved genetically encoded calcium indicators, allows sufficient signal-to-noise ratio and spatiotemporal resolution to afford a unique opportunity to capture cortex-wide dynamics on a moment-by-moment basis in behaving animals. Recent applications of this approach have been uncovering cortical dynamics at unprecedented scales during various cognitive processes, ranging from relatively simple sensorimotor integration to more complex decision-making tasks. In this review, we will highlight recent scientific advances enabled by wide-field calcium imaging in behaving mice. We then summarize several technical considerations and future opportunities for wide-field imaging to uncover large-scale circuit dynamics.

Functional network properties derived from wide-field calcium imaging differ with wakefulness and across cell type.

To improve 'bench-to-bedside' translation, it is integral that knowledge flows bidirectionally-from animal models to humans, and vice versa. This requires common analytical frameworks, as well as open software and data sharing practices. We share a new pipeline (and test dataset) for the preprocessing of wide-field optical fluorescence imaging data-an emerging mode applicable in animal models-as well as results from a functional connectivity and graph theory analysis inspired by recent work in the human neuroimaging field. The approach is demonstrated using a dataset comprised of two test-cases: (1) data from animals imaged during awake and anesthetized conditions with excitatory neurons labeled, and (2) data from awake animals with different genetically encoded fluorescent labels that target either excitatory neurons or inhibitory interneuron subtypes. Both seed-based connectivity and graph theory measures (global efficiency, transitivity, modularity, and characteristic path-length) are shown to be useful in quantifying differences between wakefulness states and cell populations. Wakefulness state and cell type show widespread effects on canonical network connectivity with variable frequency band dependence. Differences between excitatory neurons and inhibitory interneurons are observed, with somatostatin expressing inhibitory interneurons emerging as notably dissimilar from parvalbumin and vasoactive polypeptide expressing cells. In sum, we demonstrate that our pipeline can be used to examine brain state and cell-type differences in mesoscale imaging data, aiding translational neuroscience efforts. In line with open science practices, we freely release the pipeline and data to encourage other efforts in the community.

Functional Segregation and Development of Mouse Higher Visual Areas.

Recent studies suggest that higher visual areas (HVAs) in the mouse visual cortex are segregated anatomically into two visual streams, likely analogous to the ventral and dorsal streams in primates. However, HVAs in mice have yet to be characterized functionally. Moreover, it is unknown when the functional segregation of HVAs occurs during development. Here, we investigated spatiotemporal selectivity of HVAs and their development using wide-field calcium imaging. We found that lateral HVAs in the anatomical ventral stream shared similar spatiotemporal selectivity, whereas the spatiotemporal selectivity of anterior and medial HVAs in the anatomical dorsal stream was not uniform and these areas were segregated functionally into multiple groups. This functional segregation of HVAs developed and reached an adult-like pattern ∼10 d after eye opening (EO). These results suggest, not only the functional segregation of ventral and dorsal streams, but also the presence of multiple substreams in the dorsal stream, and indicate that the functional segregation of visual streams occurs gradually after EO.SIGNIFICANCE STATEMENT Investigation of the spatiotemporal selectivity of nine higher visual areas (HVAs) in adult and developing mice revealed that lateral HVAs belonging to the putative ventral stream shared similar spatiotemporal selectivity, whereas the spatiotemporal selectivity of anterior and medial HVAs belonging to the putative dorsal stream was not uniform and these areas were segregated functionally into multiple groups. These results suggest the presence of multiple substreams within the putative dorsal stream for visuospatial processing. Furthermore, we found that initially immature functional segregation among HVAs developed to an adult-like pattern ∼10 d after eye opening. These results provide a foundation for using mouse HVAs as a model to understand parallel processing and its developmental mechanism.

Attention-based CNN-BiLSTM for sleep state classification of spatiotemporal wide-field calcium imaging data.

BACKGROUND: Wide-field calcium imaging (WFCI) with genetically encoded calcium indicators allows for spatiotemporal recordings of neuronal activity in mice. When applied to the study of sleep, WFCI data are manually scored into the sleep states of wakefulness, non-REM (NREM) and REM by use of adjunct EEG and EMG recordings. However, this process is time-consuming, invasive and often suffers from low inter- and intra-rater reliability. Therefore, an automated sleep state classification method that operates on spatiotemporal WFCI data is desired. NEW METHOD: A hybrid network architecture consisting of a convolutional neural network (CNN) to extract spatial features of image frames and a bidirectional long short-term memory network (BiLSTM) with attention mechanism to identify temporal dependencies among different time points was proposed to classify WFCI data into states of wakefulness, NREM and REM sleep. RESULTS: Sleep states were classified with an accuracy of 84 % and Cohen's κ of 0.64. Gradient-weighted class activation maps revealed that the frontal region of the cortex carries more importance when classifying WFCI data into NREM sleep while posterior area contributes most to the identification of wakefulness. The attention scores indicated that the proposed network focuses on short- and long-range temporal dependency in a state-specific manner. COMPARISON WITH EXISTING METHOD: On a held out, repeated 3-hour WFCI recording, the CNN-BiLSTM achieved a κ of 0.67, comparable to a κ of 0.65 corresponding to the human EEG/EMG-based scoring. CONCLUSIONS: The CNN-BiLSTM effectively classifies sleep states from spatiotemporal WFCI data and will enable broader application of WFCI in sleep research.

A modular and adaptable analysis pipeline to compare slow cerebral rhythms across heterogeneous datasets.

Neuroscience is moving toward a more integrative discipline where understanding brain function requires consolidating the accumulated evidence seen across experiments, species, and measurement techniques. A remaining challenge on that path is integrating such heterogeneous data into analysis workflows such that consistent and comparable conclusions can be distilled as an experimental basis for models and theories. Here, we propose a solution in the context of slow-wave activity (<1 Hz), which occurs during unconscious brain states like sleep and general anesthesia and is observed across diverse experimental approaches. We address the issue of integrating and comparing heterogeneous data by conceptualizing a general pipeline design that is adaptable to a variety of inputs and applications. Furthermore, we present the Collaborative Brain Wave Analysis Pipeline (Cobrawap) as a concrete, reusable software implementation to perform broad, detailed, and rigorous comparisons of slow-wave characteristics across multiple, openly available electrocorticography (ECoG) and calcium imaging datasets.

Global and subtype-specific modulation of cortical inhibitory neurons regulated by acetylcholine during motor learning.

Inhibitory neurons (INs) consist of distinct subtypes with unique functions. Previous studies on INs mainly focused on single brain regions, and thus it remains unclear whether the modulation of IN subtypes occurs globally across multiple regions. Here, we monitored the activity of different cortical IN subtypes at both macroscale and microscale in mice learning a lever-press task. Learning evoked a global modulation of IN subtypes throughout the cortex. The initial learning phase involved strong activation of vasoactive intestinal peptide-expressing INs (VIP-INs) and weak activation of somatostatin-expressing INs (SOM-INs). Inactivating VIP-INs increased SOM-IN activity and impaired initial learning. Concurrently, cortical cholinergic inputs from the basal forebrain were initially more active but became less engaged over learning. Manipulation of the cholinergic system impaired motor learning and differentially altered activity of IN subtypes. These results reveal that motor learning involves a global and subtype-specific modulation on cortical INs regulated by the cholinergic system.

Tracking the Effect of Therapy With Single-Trial Based Classification After Stroke.

Stroke is a debilitating disease that leads, in the 50% of cases, to permanent motor or cognitive impairments. The effectiveness of therapies that promote recovery after stroke depends on indicators of the disease state that can measure the degree of recovery or predict treatment response or both. Here, we propose to use single-trial classification of task dependent neural activity to assess the disease state and track recovery after stroke. We tested this idea on calcium imaging data of the dorsal cortex of healthy, spontaneously recovered and rehabilitated mice while performing a forelimb retraction task. Results show that, at a single-trial level for the three experimental groups, neural activation during the reward pull can be detected with high accuracy with respect to the background activity in all cortical areas of the field of view and this activation is quite stable across trials and subjects of the same group. Moreover, single-trial responses during the reward pull can be used to discriminate between healthy and stroke subjects with areas closer to the injury site displaying higher discrimination capability than areas closer to this site. Finally, a classifier built to discriminate between controls and stroke at the single-trial level can be used to generate an index of the disease state, the therapeutic score, which is validated on the group of rehabilitated mice. In conclusion, task-related neural activity can be used as an indicator of disease state and track recovery without selecting a peculiar feature of the neural responses. This novel method can be used in both the development and assessment of different therapeutic strategies.

Automated sleep state classification of wide-field calcium imaging data via multiplex visibility graphs and deep learning.

BACKGROUND: Wide-field calcium imaging (WFCI) allows for monitoring of cortex-wide neural dynamics in mice. When applied to the study of sleep, WFCI data are manually scored into the sleep states of wakefulness, non-REM (NREM) and REM by use of adjunct EEG and EMG recordings. However, this process is time-consuming and often suffers from low inter- and intra-rater reliability and invasiveness. Therefore, an automated sleep state classification method that operates on WFCI data alone is needed. NEW METHOD: A hybrid, two-step method is proposed. In the first step, spatial-temporal WFCI data is mapped to multiplex visibility graphs (MVGs). Subsequently, a two-dimensional convolutional neural network (2D CNN) is employed on the MVGs to be classified as wakefulness, NREM and REM. RESULTS: Sleep states were classified with an accuracy of 84% and Cohen's κ of 0.67. The method was also effectively applied on a binary classification of wakefulness/sleep (accuracy=0.82, κ = 0.62) and a four-class wakefulness/sleep/anesthesia/movement classification (accuracy=0.74, κ = 0.66). Gradient-weighted class activation maps revealed that the CNN focused on short- and long-term temporal connections of MVGs in a sleep state-specific manner. Sleep state classification performance when using individual brain regions was highest for the posterior area of the cortex and when cortex-wide activity was considered. COMPARISON WITH EXISTING METHOD: On a 3-hour WFCI recording, the MVG-CNN achieved a κ of 0.65, comparable to a κ of 0.60 corresponding to the human EEG/EMG-based scoring. CONCLUSIONS: The hybrid MVG-CNN method accurately classifies sleep states from WFCI data and will enable future sleep-focused studies with WFCI.

Revealing the Cortical Glutamatergic Neural Activity During Burst Suppression by Simultaneous wide Field Calcium Imaging and Electroencephalography in Mice.

Burst suppression (BS) is an electroencephalogram (EEG) pattern in which signals alternates between high-amplitude slow waves (burst waves) and nearly flat low-amplitude waves (suppression waves). In this study, we used wide-field (8.32 mm × 8.32 mm) fluorescent calcium imaging to record the activity of glutamatergic neurons in the parietal and occipital cortex, in conjunction with EEG recordings under BS induced by different anesthetics (sevoflurane, isoflurane, and propofol), to investigate the spatiotemporal pattern of neural activity under BS. The calcium signal of all observed cortices was decreased during the phase of EEG suppression. However, during the phase of EEG burst, the calcium signal in areas of the medial cortex, such as the secondary motor and retrosplenial area, was excited, whereas the signal in areas of the lateral cortex, such as the hindlimb cortex, forelimb cortex, barrel field, and primary visual area, was still suppressed or only weakly excited. Correlation analysis showed a strong correlation between the EEG signal and the calcium signal in the medial cortex under BS (except for propofol induced signals). As the burst-suppression ratio (BSR) increased, the regions with strong correlation coefficients became smaller, but strong correlation coefficients were still noted in the medial cortex. Taken together, our results reveal the landscape of cortical activity underlying BS.

Mesoscale cortical dynamics reflect the interaction of sensory evidence and temporal expectation during perceptual decision-making.

How sensory evidence is transformed across multiple brain regions to influence behavior remains poorly understood. We trained mice in a visual change detection task designed to separate the covert antecedents of choices from activity associated with their execution. Wide-field calcium imaging across the dorsal cortex revealed fundamentally different dynamics of activity underlying these processes. Although signals related to execution of choice were widespread, fluctuations in sensory evidence in the absence of overt motor responses triggered a confined activity cascade, beginning with transient modulation of visual cortex and followed by sustained recruitment of the secondary and primary motor cortex. Activation of the motor cortex by sensory evidence was modulated by animals' expectation of when the stimulus was likely to change. These results reveal distinct activation timescales of specific cortical areas by sensory evidence during decision-making and show that recruitment of the motor cortex depends on the interaction of sensory evidence and temporal expectation.

Sensory and Behavioral Components of Neocortical Signal Flow in Discrimination Tasks with Short-Term Memory.

In the neocortex, each sensory modality engages distinct sensory areas that route information to association areas. Where signal flow converges for maintaining information in short-term memory and how behavior may influence signal routing remain open questions. Using wide-field calcium imaging, we compared cortex-wide neuronal activity in layer 2/3 for mice trained in auditory and tactile tasks with delayed response. In both tasks, mice were either active or passive during stimulus presentation, moving their body or sitting quietly. Irrespective of behavioral strategy, auditory and tactile stimulation activated distinct subdivisions of the posterior parietal cortex, anterior area A and rostrolateral area RL, which held stimulus-related information necessary for the respective tasks. In the delay period, in contrast, behavioral strategy rather than sensory modality determined short-term memory location, with activity converging frontomedially in active trials and posterolaterally in passive trials. Our results suggest behavior-dependent routing of sensory-driven cortical signals flow from modality-specific posterior parietal cortex (PPC) subdivisions to higher association areas.

Excitatory Neuronal Hubs Configure Multisensory Integration of Slow Waves in Association Cortex.

Multisensory integration (MSI) is a fundamental emergent property of the mammalian brain. During MSI, perceptual information encoded in patterned activity is processed in multimodal association cortex. The systems-level neuronal dynamics that coordinate MSI, however, are unknown. Here, we demonstrate intrinsic hub-like network activity in the association cortex that regulates MSI. We engineered calcium reporter mouse lines based on the fluorescence resonance energy transfer sensor yellow cameleon (YC2.60) expressed in excitatory or inhibitory neurons. In medial and parietal association cortex, we observed spontaneous slow waves that self-organized into hubs defined by long-range excitatory and local inhibitory circuits. Unlike directional source/sink-like flows in sensory areas, medial/parietal excitatory and inhibitory hubs had net-zero balanced inputs. Remarkably, multisensory stimulation triggered rapid phase-locking mainly of excitatory hub activity persisting for seconds after the stimulus offset. Therefore, association cortex tends to form balanced excitatory networks that configure slow-wave phase-locking for MSI. VIDEO ABSTRACT.

Transformation of Cortex-wide Emergent Properties during Motor Learning.

Learning involves a transformation of brain-wide operation dynamics. However, our understanding of learning-related changes in macroscopic dynamics is limited. Here, we monitored cortex-wide activity of the mouse brain using wide-field calcium imaging while the mouse learned a motor task over weeks. Over learning, the sequential activity across cortical modules became temporally more compressed, and its trial-by-trial variability decreased. Moreover, a new flow of activity emerged during learning, originating from premotor cortex (M2), and M2 became predictive of the activity of many other modules. Inactivation experiments showed that M2 is critical for the post-learning dynamics in the cortex-wide activity. Furthermore, two-photon calcium imaging revealed that M2 ensemble activity also showed earlier activity onset and reduced variability with learning, which was accompanied by changes in the activity-movement relationship. These results reveal newly emergent properties of macroscopic cortical dynamics during motor learning and highlight the importance of M2 in controlling learned movements.