Separability of calcium slow waves and functional connectivity during wake, sleep, and anesthesia

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.

Spatiotemporal relationships between neuronal, metabolic, and hemodynamic signals in the awake and anesthetized mouse brain

Neurovascular coupling (NVC) and neurometabolic coupling (NMC) provide the basis for functional magnetic resonance imaging and positron emission tomography to map brain neurophysiology. While increases in neuronal activity are often accompanied by increases in blood oxygen delivery and oxidative metabolism, these observations are not the rule. This decoupling is important when interpreting brain network organization (e.g., resting-state functional connectivity [RSFC]) because it is unclear whether changes in NMC/NVC affect RSFC measures. We leverage wide-field optical imaging in Thy1-jRGECO1a mice to map cortical calcium activity in pyramidal neurons, flavoprotein autofluorescence (representing oxidative metabolism), and hemodynamic activity during wake and ketamine/xylazine anesthesia. Spontaneous dynamics of all contrasts exhibit patterns consistent with RSFC. NMC/NVC relative to excitatory activity varies over the cortex. Ketamine/xylazine profoundly alters NVC but not NMC. Compared to awake RSFC, ketamine/xylazine affects metabolic-based connectomes moreso than hemodynamic-based measures of RSFC. Anesthesia-related differences in NMC/NVC timing do not appreciably alter RSFC structure.

Mapping brain function in adults and young children during naturalistic viewing with high-density diffuse optical tomography

Human studies of early brain development have been limited by extant neuroimaging methods. MRI scanners present logistical challenges for imaging young children, while alternative modalities like functional near-infrared spectroscopy have traditionally been limited by image quality due to sparse sampling. In addition, conventional tasks for brain mapping elicit low task engagement, high head motion, and considerable participant attrition in pediatric populations. As a result, typical and atypical developmental trajectories of processes such as language acquisition remain understudied during sensitive periods over the first years of life. We evaluate high-density diffuse optical tomography (HD-DOT) imaging combined with movie stimuli for high resolution optical neuroimaging in awake children ranging from 1 to 7 years of age. We built an HD-DOT system with design features geared towards enhancing both image quality and child comfort. Furthermore, we characterized a library of animated movie clips as a stimulus set for brain mapping and we optimized associated data analysis pipelines. Together, these tools could map cortical responses to movies and contained features such as speech in both adults and awake young children. This study lays the groundwork for future research to investigate response variability in larger pediatric samples and atypical trajectories of early brain development in clinical populations.

Multimodal measures of spontaneous brain activity reveal both common and divergent patterns of cortical functional organization

Large-scale functional networks have been characterized in both rodent and human brains, typically by analyzing fMRI-BOLD signals. However, the relationship between fMRI-BOLD and underlying neural activity is complex and incompletely understood, which poses challenges to interpreting network organization obtained using this technique. Additionally, most work has assumed a disjoint functional network organization (i.e., brain regions belong to one and only one network). Here, we employ wide-field Ca2+ imaging simultaneously with fMRI-BOLD in mice expressing GCaMP6f in excitatory neurons. We determine cortical networks discovered by each modality using a mixed-membership algorithm to test the hypothesis that functional networks exhibit overlapping organization. We find that there is considerable network overlap (both modalities) in addition to disjoint organization. Our results show that multiple BOLD networks are detected via Ca2+ signals, and networks determined by low-frequency Ca2+ signals are only modestly more similar to BOLD networks. In addition, the principal gradient of functional connectivity is nearly identical for BOLD and Ca2+ signals. Despite similarities, important differences are also detected across modalities, such as in measures of functional connectivity strength and diversity. In conclusion, Ca2+ imaging uncovers overlapping functional cortical organization in the mouse that reflects several, but not all, properties observed with fMRI-BOLD signals.

Propagating population activity patterns during spontaneous slow waves in the thalamus of rodents

Slow waves (SWs) represent the most prominent electrophysiological events in the thalamocortical system under anesthesia and during deep sleep. Recent studies have revealed that SWs have complex spatiotemporal dynamics and propagate across neocortical regions. However, it is still unclear whether neuronal activity in the thalamus exhibits similar propagation properties during SWs. Here, we report propagating population activity in the thalamus of ketamine/xylazine-anesthetized rats and mice visualized by high-density silicon probe recordings. In both rodent species, propagation of spontaneous thalamic activity during up-states was most frequently observed in dorsal thalamic nuclei such as the higher order posterior (Po), lateral posterior (LP) or laterodorsal (LD) nuclei. The preferred direction of thalamic activity spreading was along the dorsoventral axis, with over half of the up-states exhibiting a gradual propagation in the ventral-to-dorsal direction. Furthermore, simultaneous neocortical and thalamic recordings collected under anesthesia demonstrated that there is a weak but noticeable interrelation between propagation patterns observed during cortical up-states and those displayed by thalamic population activity. In addition, using chronically implanted silicon probes, we detected propagating activity patterns in the thalamus of naturally sleeping rats during slow-wave sleep. However, in comparison to propagating up-states observed under anesthesia, these propagating patterns were characterized by a reduced rate of occurrence and a faster propagation speed. Our findings suggest that the propagation of spontaneous population activity is an intrinsic property of the thalamocortical network during synchronized brain states such as deep sleep or anesthesia. Additionally, our data implies that the neocortex may have partial control over the formation of propagation patterns within the dorsal thalamus under anesthesia.

Mecp2 deletion results in profound alterations of developmental and adult functional connectivity

As a regressive neurodevelopmental disorder with a well-established genetic cause, Rett syndrome and its Mecp2 loss-of-function mouse model provide an excellent opportunity to define potentially translatable functional signatures of disease progression, as well as offer insight into the role of Mecp2 in functional circuit development. Thus, we applied widefield optical fluorescence imaging to assess mesoscale calcium functional connectivity (FC) in the Mecp2 cortex both at postnatal day (P)35 in development and during the disease-related decline. We found that FC between numerous cortical regions was disrupted in Mecp2 mutant males both in juvenile development and early adulthood. Female Mecp2 mice displayed an increase in homotopic contralateral FC in the motor cortex at P35 but not in adulthood, where instead more posterior parietal regions were implicated. An increase in the amplitude of connection strength, both with more positive correlations and more negative anticorrelations, was observed across the male cortex in numerous functional regions. Widespread rescue of MeCP2 protein in GABAergic neurons rescued none of these functional deficits, nor, surprisingly, the expected male lifespan. Altogether, the female results identify early signs of disease progression, while the results in males indicate MeCP2 protein is required for typical FC in the brain.

Multimodal measures of spontaneous brain activity reveal both common and divergent patterns of cortical functional organization

Large-scale functional networks have been characterized in both rodent and human brains, typically by analyzing fMRI-BOLD signals. However, the relationship between fMRI-BOLD and underlying neural activity is complex and incompletely understood, which poses challenges to interpreting network organization obtained using this technique. Additionally, most work has assumed a disjoint functional network organization (i.e., brain regions belong to one and only one network). Here, we employed wide-field Ca2+ imaging simultaneously with fMRI-BOLD in mice expressing GCaMP6f in excitatory neurons. We determined cortical networks discovered by each modality using a mixed-membership algorithm to test the hypothesis that functional networks are overlapping rather than disjoint. Our results show that multiple BOLD networks are detected via Ca2+ signals; there is considerable network overlap (both modalities); networks determined by low-frequency Ca2+ signals are only modestly more similar to BOLD networks; and, despite similarities, important differences are detected across modalities (e.g., brain region "network diversity"). In conclusion, Ca2+ imaging uncovered overlapping functional cortical organization in the mouse that reflected several, but not all, properties observed with fMRI-BOLD signals.

Transient disruption of functional connectivity and depression of neural fluctuations in a mouse model of acute septic encephalopathy

Septic encephalopathy leads to major and costly burdens for a large percentage of admitted hospital patients. Elderly patients are at an increased risk, especially those with dementia. Current treatments are aimed at sedation to combat mental status changes and are not aimed at the underlying cause of encephalopathy. Indeed, the underlying pathology linking together peripheral infection and altered neural function has not been established, largely because good, acutely accessible readouts of encephalopathy in animal models do not exist. Behavioral testing in animals lasts multiple days, outlasting the time frame of acute encephalopathy. Here, we propose optical fluorescent imaging of neural functional connectivity (FC) as a readout of encephalopathy in a mouse model of acute sepsis. Imaging and basic behavioral assessment were performed at baseline, Hr8, Hr24, and Hr72 following injection of either lipopolysaccharide or phosphate buffered saline. Neural FC strength decreased at Hr8 and returned to baseline by Hr72 in motor, somatosensory, parietal, and visual cortical regions. Additionally, neural fluctuations transiently declined at Hr8 and returned to baseline by Hr72. Both FC strength and fluctuation tone correlated with neuroscore indicating this imaging methodology is a sensitive and acute readout of encephalopathy.

Simulations approaching data: cortical slow waves in inferred models of the whole hemisphere of mouse

The development of novel techniques to record wide-field brain activity enables estimation of data-driven models from thousands of recording channels and hence across large regions of cortex. These in turn improve our understanding of the modulation of brain states and the richness of traveling waves dynamics. Here, we infer data-driven models from high-resolution in-vivo recordings of mouse brain obtained from wide-field calcium imaging. We then assimilate experimental and simulated data through the characterization of the spatio-temporal features of cortical waves in experimental recordings. Inference is built in two steps: an inner loop that optimizes a mean-field model by likelihood maximization, and an outer loop that optimizes a periodic neuro-modulation via direct comparison of observables that characterize cortical slow waves. The model reproduces most of the features of the non-stationary and non-linear dynamics present in the high-resolution in-vivo recordings of the mouse brain. The proposed approach offers new methods of characterizing and understanding cortical waves for experimental and computational neuroscientists.

Open-source statistical and data processing tools for wide-field optical imaging data in mice

Significance

Wide-field optical imaging (WOI) can produce concurrent hemodynamic and cell-specific calcium recordings across the entire cerebral cortex in animal models. There have been multiple studies using WOI to image mouse models with various environmental or genetic manipulations to understand various diseases. Despite the utility of pursuing mouse WOI alongside human functional magnetic resonance imaging (fMRI), and the multitude of analysis toolboxes in the fMRI literature, there is not an available open-source, user-friendly data processing and statistical analysis toolbox for WOI data.

Aim

To assemble a MATLAB toolbox for processing WOI data, as described and adapted to combine techniques from multiple WOI groups and fMRI.

Approach

We outline our MATLAB toolbox on GitHub with multiple data analysis packages and translate a commonly used statistical approach from the fMRI literature to the WOI data. To illustrate the utility of our MATLAB toolbox, we demonstrate the ability of the processing and analysis framework to detect a well-established deficit in a mouse model of stroke and plot activation areas during an electrical paw stimulus experiment.

Results

Our processing toolbox and statistical methods isolate a somatosensory-based deficit 3 days following photothrombotic stroke and cleanly localize sensory stimulus activations.

Conclusions

The toolbox presented here details an open-source, user-friendly compilation of WOI processing tools with statistical methods to apply to any biological question investigated with WOI techniques.

Wide-Field Optical Imaging in Mouse Models of Ischemic Stroke

Functional neuroimaging is a powerful tool for evaluating how local and global brain circuits evolve after focal ischemia and how these changes relate to functional recovery. For example, acutely after stroke, changes in functional brain organization relate to initial deficit and are predictive of recovery potential. During recovery, the reemergence and restoration of connections lost due to stroke correlate with recovery of function. Thus, information gleaned from functional neuroimaging can be used as a proxy for behavior and inform on the efficacy of interventional strategies designed to affect plasticity mechanisms after injury. And because these findings are consistently observed across species, bridge measurements can be made in animal models to enrich findings in human stroke populations. In mice, genetic engineering techniques have provided several new opportunities for extending optical neuroimaging methods to more direct measures of neuronal activity. These developments are especially useful in the context of stroke where neurovascular coupling can be altered, potentially limiting imaging measures based on hemodynamic activity alone. This chapter is designed to give an overview of functional wide-field optical imaging (WFOI) for applications in rodent models of stroke, primarily in the mouse. The goal is to provide a protocol for laboratories that want to incorporate an affordable functional neuroimaging assay into their current research thrusts, but perhaps lack the background knowledge or equipment for developing a new arm of research in their lab. Within, we offer a comprehensive guide developing and applying WFOI technology with the hope of facilitating accessibility of neuroimaging technology to other researchers in the stroke field.

Transcranial cortex-wide Ca2+ imaging for the functional mapping of cortical dynamics

Visualization and tracking of the information flow in the broader brain area are essential because nerve cells make a vast network in the brain. Fluorescence Ca²⁺ imaging is a simultaneous visualization of brain cell activities in a wide area. Instead of classical chemical indicators, developing various types of transgenic animals that express Ca²⁺-sensitive fluorescent proteins enables us to observe brain activities in living animals at a larger scale for a long time. Multiple kinds of literature have reported that transcranial imaging of such transgenic animals is practical for monitoring the wide-field information flow across the broad brain regions, although it has a lower spatial resolution. Notably, this technique is helpful for the initial evaluation of cortical function in disease models. This review will introduce fully intact transcranial macroscopic imaging and cortex-wide Ca²⁺ imaging as practical applications.

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.

REM sleep is associated with distinct global cortical dynamics and controlled by occipital cortex

The cerebral cortex is spontaneously active during sleep, yet it is unclear how this global cortical activity is spatiotemporally organized, and whether such activity not only reflects sleep states but also contributes to sleep state switching. Here we report that cortex-wide calcium imaging in mice revealed distinct sleep stage-dependent spatiotemporal patterns of global cortical activity, and modulation of such patterns could regulate sleep state switching. In particular, elevated activation in the occipital cortical regions (including the retrosplenial cortex and visual areas) became dominant during rapid-eye-movement (REM) sleep. Furthermore, such pontogeniculooccipital (PGO) wave-like activity was associated with transitions to REM sleep, and optogenetic inhibition of occipital activity strongly promoted deep sleep by suppressing the NREM-to-REM transition. Thus, whereas subcortical networks are critical for initiating and maintaining sleep and wakefulness states, distinct global cortical activity also plays an active role in controlling sleep states.

Validation of Deep Learning-based Sleep State Classification

Deep learning methods have been developed to classify sleep states of mouse electroencephalogram (EEG) and electromyogram (EMG) recordings with accuracy reported as high as 97%. However, when applied to independent datasets, with a variety of experimental and recording conditions, sleep state classification accuracy often drops due to distributional shift. Mixture z-scoring, a pre-processing standardization of EEG/EMG signals, has been suggested to account for these variations. This study sought to validate mixture z-scoring in combination with a deep learning method on an independent dataset. The open-source software Accusleep, which implements mixture z-scoring in combination with deep learning via a convolutional neural network, was used to classify sleep states in 12, three-hour EEG/EMG recordings from mice sleeping in a head-fixed position. Mixture z-scoring with deep learning classified sleep states on two independent recordings with 85-92% accuracy and a Cohen's κ of 0.66-0.71. These results validate mixture z-scoring in combination with deep learning to classify sleep states with the potential for widespread use.

Normal aging in mice is associated with a global reduction in cortical spectral power and network-specific declines in functional connectivity

Normal aging is associated with a variety of neurologic changes including declines in cognition, memory, and motor activity. These declines correlate with neuronal changes in synaptic structure and function. Degradation of brain network activity and connectivity represents a likely mediator of age-related functional deterioration resulting from these neuronal changes. Human studies have demonstrated both general decreases in spontaneous cortical activity and disruption of cortical networks with aging. Current techniques used to study cerebral network activity are hampered either by limited spatial resolution (e.g. electroencephalography, EEG) or limited temporal resolution (e.g., functional magnetic resonance imaging, fMRI). Here we utilize mesoscale imaging of neuronal activity in Thy1-GCaMP6f mice to characterize neuronal network changes in aging with high spatial resolution across a wide frequency range. We show that while evoked activity is unchanged with aging, spontaneous neuronal activity decreases across a wide frequency range (0.01-4 Hz) involving all regions of the cortex. In contrast to this global reduction in cortical power, we found that aging is associated with functional connectivity (FC) deterioration of select networks including somatomotor, cingulate, and retrosplenial nodes. These changes are corroborated by reductions in homotopic FC and node degree within somatomotor and visual cortices. Finally, we found that whole-cortex delta power and delta band node degree correlate with exploratory activity in young but not aged animals. Together these data suggest that aging is associated with global declines in spontaneous cortical activity and focal deterioration of network connectivity, and that these reductions may be associated with age-related behavioral declines.

Selective blockade of rat brain T-type calcium channels provides insights on neurophysiological basis of arousal dependent resting state functional magnetic resonance imaging signals

A number of studies point to slow (0.1–2 Hz) brain rhythms as the basis for the resting-state functional magnetic resonance imaging (rsfMRI) signal. Slow waves exist in the absence of stimulation, propagate across the cortex, and are strongly modulated by vigilance similar to large portions of the rsfMRI signal. However, it is not clear if slow rhythms serve as the basis of all neural activity reflected in rsfMRI signals, or just the vigilance-dependent components. The rsfMRI data exhibit quasi-periodic patterns (QPPs) that appear to increase in strength with decreasing vigilance and propagate across the brain similar to slow rhythms. These QPPs can complicate the estimation of functional connectivity (FC) via rsfMRI, either by existing as unmodeled signal or by inducing additional wide-spread correlation between voxel-time courses of functionally connected brain regions. In this study, we examined the relationship between cortical slow rhythms and the rsfMRI signal, using a well-established pharmacological model of slow wave suppression. Suppression of cortical slow rhythms led to significant reduction in the amplitude of QPPs but increased rsfMRI measures of intrinsic FC in rats. The results suggest that cortical slow rhythms serve as the basis of only the vigilance-dependent components (e.g., QPPs) of rsfMRI signals. Further attenuation of these non-specific signals enhances delineation of brain functional networks.

SIRT1 mediates hypoxic postconditioning- and resveratrol-induced protection against functional connectivity deficits after subarachnoid hemorrhage

Functional connectivity (FC) is a sensitive metric that provides a readout of whole cortex coordinate neural activity in a mouse model. We examine the impact of experimental SAH modeled through endovascular perforation, and the effectiveness of subsequent treatment on FC, through three key questions: 1) Does the endovascular perforation model of SAH induce deficits in FC; 2) Does exposure to hypoxic conditioning provide protection against these FC deficits and, if so, is this neurovascular protection SIRT1-mediated; and 3) does treatment with the SIRT1 activator resveratrol alone provide protection against these FC deficits? Cranial windows were adhered on skull-intact mice that were then subjected to either sham or SAH surgery and either left untreated or treated with hypoxic post-conditioning (with or without EX527) or resveratrol for 3 days. Mice were imaged 3 days post-SAH/sham surgery, temporally aligned with the onset of major SAH sequela in mice. Here we show that the endovascular perforation model of SAH induces global and network-specific deficits in FC by day 3, corresponding with the time frame of DCI in mice. Hypoxic conditioning provides SIRT1-mediated protection against these network-specific FC deficits post-SAH, as does treatment with resveratrol. Conditioning-based strategies provide multifaceted neurovascular protection in experimental SAH.

Homotopic contralesional excitation suppresses spontaneous circuit repair and global network reconnections following ischemic stroke

Understanding circuit-level manipulations that affect the brain's capacity for plasticity will inform the design of targeted interventions that enhance recovery after stroke. Following stroke, increased contralesional activity (e.g. use of the unaffected limb) can negatively influence recovery, but it is unknown which specific neural connections exert this influence, and to what extent increased contralesional activity affects systems- and molecular-level biomarkers of recovery. Here, we combine optogenetic photostimulation with optical intrinsic signal imaging to examine how contralesional excitatory activity affects cortical remodeling after stroke in mice. Following photothrombosis of left primary somatosensory forepaw (S1FP) cortex, mice either recovered spontaneously or received chronic optogenetic excitation of right S1FP over the course of 4 weeks. Contralesional excitation suppressed perilesional S1FP remapping and was associated with abnormal patterns of stimulus-evoked activity in the unaffected limb. This maneuver also prevented the restoration of resting-state functional connectivity (RSFC) within the S1FP network, RSFC in several networks functionally distinct from somatomotor regions, and resulted in persistent limb-use asymmetry. In stimulated mice, perilesional tissue exhibited transcriptional changes in several genes relevant for recovery. Our results suggest that contralesional excitation impedes local and global circuit reconnection through suppression of cortical activity and several neuroplasticity-related genes after stroke, and highlight the importance of site selection for targeted therapeutic interventions after focal ischemia.

A Multivariate Functional Connectivity Approach to Mapping Brain Networks and Imputing Neural Activity in Mice

Temporal correlation analysis of spontaneous brain activity (e.g., Pearson "functional connectivity," FC) has provided insights into the functional organization of the human brain. However, bivariate analysis techniques such as this are often susceptible to confounding physiological processes (e.g., sleep, Mayer-waves, breathing, motion), which makes it difficult to accurately map connectivity in health and disease as these physiological processes affect FC. In contrast, a multivariate approach to imputing individual neural networks from spontaneous neuroimaging data could be influential to our conceptual understanding of FC and provide performance advantages. Therefore, we analyzed neural calcium imaging data from Thy1-GCaMP6f mice while either awake, asleep, anesthetized, during low and high bouts of motion, or before and after photothrombotic stroke. A linear support vector regression approach was used to determine the optimal weights for integrating the signals from the remaining pixels to accurately predict neural activity in a region of interest (ROI). The resultant weight maps for each ROI were interpreted as multivariate functional connectivity (MFC), resembled anatomical connectivity, and demonstrated a sparser set of strong focused positive connections than traditional FC. While global variations in data have large effects on standard correlation FC analysis, the MFC mapping methods were mostly impervious. Lastly, MFC analysis provided a more powerful connectivity deficit detection following stroke compared to traditional FC.

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.

Neurophotonic tools for microscopic measurements and manipulation: status report

Neurophotonics was launched in 2014 coinciding with the launch of the BRAIN Initiative focused on development of technologies for advancement of neuroscience. For the last seven years, Neurophotonics' agenda has been well aligned with this focus on neurotechnologies featuring new optical methods and tools applicable to brain studies. While the BRAIN Initiative 2.0 is pivoting towards applications of these novel tools in the quest to understand the brain, this status report reviews an extensive and diverse toolkit of novel methods to explore brain function that have emerged from the BRAIN Initiative and related large-scale efforts for measurement and manipulation of brain structure and function. Here, we focus on neurophotonic tools mostly applicable to animal studies. A companion report, scheduled to appear later this year, will cover diffuse optical imaging methods applicable to noninvasive human studies. For each domain, we outline the current state-of-the-art of the respective technologies, identify the areas where innovation is needed, and provide an outlook for the future directions.

Altered dorsal functional connectivity after post-weaning social isolation and resocialization in mice

BACKGROUND

Social isolation (SI) leads to various mental health disorders. Despite abundant studies on behavioral and neurobiological changes induced by post-weaning SI, the characterization of its imaging correlates, such as resting-state functional connectivity (RSFC), is critically lacking. In addition, the effects of resocialization after isolation remain unclear. Therefore, this study aimed to explore the effects of 1) SI on cortical functional connectivity and 2) subsequent resocialization on behavior and functional connectivity.

METHODS

Behavioral tests were conducted to validate the post-weaning SI mouse model, which is isolated during the juvenile period. Wide-field optical mapping was performed to observe both neuronal and hemodynamic signals in the cortex under anesthesia. Using seed-based and graph theoretical analyses, RSFC was analyzed. SI mice were then resocialized and the array of behavior and imaging tests was conducted.

RESULTS

Behaviorally, SI mice showed elevated anxiety, social preference, and aggression. RSFC analyses using the seed-based approach revealed decreased cortical functional connectivity in SI mice, especially in the frontal region. Graph network analyses demonstrated significant reduction in network segregation measures. After resocialization, mice exhibited recovered anxiogenic and aggressive behavior, but RSFC data did not show significant changes.

CONCLUSIONS

We observed an overall decrease in functional connectivity in SI mice. Moreover, resocialization restored the disruptions in behavioral patterns but functional connectivity was not recovered. To our knowledge, this is the first study to report that, despite the recovering tendencies of behavior in resocialized mice, similar changes in RSFC were not observed. This suggests that disruptions in functional connectivity caused by social isolation remain as long-term sequelae.

Combined resting state-fMRI and calcium recordings show stable brain states for task-induced fMRI in mice under combined ISO/MED anesthesia

For fMRI in animal models, the combination of low-dose anesthetic, isoflurane (ISO), and the sedative medetomidine (MED) has recently become an advocated regimen to achieve stable neuronal states and brain networks in rats that are required for reliable task-induced BOLD fMRI. However, in mice the temporal stability of neuronal states and networks in resting-state (rs)-fMRI experiments during the combined ISO/MED regimen has not been systematically investigated. Using a multimodal approach with optical calcium (Ca²⁺) recordings and rs-fMRI, we investigated cortical neuronal/astrocytic Ca²⁺activity states and brain networks at multiple time points while switching from anesthesia with 1% ISO to a combined ISO/MED regimen. We found that cortical activity states reached a steady-state 45 min following start of MED infusion as indicated by stable Ca²⁺ transients. Similarly, rs-networks were not statistically different between anesthesia with ISO and the combined ISO/MED regimen 45 and 100 min after start of MED. Importantly, during the transition time we identified changed rs-network signatures that likely reflect the different mode of action of the respective anesthetic; these included a dose-dependent increase in cortico-cortical functional connectivity (FC) presumably caused by reduction of ISO concentration and decreased FC in subcortical arousal nuclei due to MED infusion. Furthermore, we report detection of visual stimulation-induced BOLD fMRI during the stable ISO/MED neuronal state 45 min after induction. Based on our findings, we recommend a 45-minute waiting period after switching from ISO anesthesia to the combined ISO/MED regimen before performing rs- or task-induced fMRI experiments.

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.

Cortical propagation tracks functional recovery after stroke

Stroke is a debilitating condition affecting millions of people worldwide. The development of improved rehabilitation therapies rests on finding biomarkers suitable for tracking functional damage and recovery. To achieve this goal, we perform a spatiotemporal analysis of cortical activity obtained by wide-field calcium images in mice before and after stroke. We compare spontaneous recovery with three different post-stroke rehabilitation paradigms, motor training alone, pharmacological contralesional inactivation and both combined. We identify three novel indicators that are able to track how movement-evoked global activation patterns are impaired by stroke and evolve during rehabilitation: the duration, the smoothness, and the angle of individual propagation events. Results show that, compared to pre-stroke conditions, propagation of cortical activity in the subacute phase right after stroke is slowed down and more irregular. When comparing rehabilitation paradigms, we find that mice treated with both motor training and pharmacological intervention, the only group associated with generalized recovery, manifest new propagation patterns, that are even faster and smoother than before the stroke. In conclusion, our new spatiotemporal propagation indicators could represent promising biomarkers that are able to uncover neural correlates not only of motor deficits caused by stroke but also of functional recovery during rehabilitation. In turn, these insights could pave the way towards more targeted post-stroke therapies.

Millions of people worldwide suffer from long-lasting motor deficits caused by stroke. Very recently, the two basic therapeutic approaches, motor training and pharmacological intervention, have been combined in order to achieve a more efficient functional recovery. In this study, we analyze the neurophysiological activity in the brain of mice observed with in vivo calcium imaging before and after the induction of a stroke. We use a newly developed universal approach based on the temporal sequence of local activation in different brain regions to quantify three properties of global propagation patterns: duration, smoothness and angle. These innovative spatiotemporal propagation indicators allow us to track damage and functional recovery following stroke and to quantify the relative success of motor training, pharmacological inactivation, and a combination of both, compared to spontaneous recovery. We show that all three treatments reverse the alterations observed during the subacute phase right after stroke. We also find that combining motor training and pharmacological intervention does not restore pre-stroke features but rather leads to the emergence of new propagation patterns that, surprisingly, are even faster and smoother than the pre-stroke patterns.

Characterization of Developmental Changes in Spontaneous Electrical Activity of Medial Superior Olivary Neurons Before Hearing Onset With a Combination of Injectable and Volatile Anesthesia

In this work the impact of two widely used anesthetics on the electrical activity of auditory brainstem neurons was studied during postnatal development. Spontaneous electrical activity in neonate rats of either sex was analyzed through a ventral craniotomy in mechanically ventilated pups to carry out patch clamp and multi-electrode electrophysiology recordings in the medial region of the superior olivary complex (SOC) between birth (postnatal day 0, P0) and P12. Recordings were obtained in pups anesthetized with the injectable mix of ketamine/xylazine (K/X mix), with the volatile anesthetic isoflurane (ISO), or in pups anesthetized with K/X mix that were also exposed to ISO. The results of patch clamp recordings demonstrate for the first time that olivary and periolivary neurons in the medial region of the SOC fire bursts of action potentials. The results of multielectrode recordings suggest that the firing pattern of single units recorded in K/X mix is similar to that recorded in ISO anesthetized rat pups. Taken together, the results of this study provide a framework to use injectable and volatile anesthetics for future studies to obtain functional information on the activity of medial superior olivary neurons in vivo.

Peripheral sensory stimulation elicits global slow waves by recruiting somatosensory cortex bilaterally

Slow waves (SWs) are globally propagating, low-frequency (0.5- to 4-Hz) oscillations that are prominent during sleep and anesthesia. SWs are essential to neural plasticity and memory. However, much remains unknown about the mechanisms coordinating SW propagation at the macroscale. To assess SWs in the context of macroscale networks, we recorded cortical activity in awake and ketamine/xylazine-anesthetized mice using widefield optical imaging with fluorescent calcium indicator GCaMP6f. We demonstrate that unilateral somatosensory stimulation evokes bilateral waves that travel across the cortex with state-dependent trajectories. Under anesthesia, we observe that rhythmic stimuli elicit globally resonant, front-to-back propagating SWs. Finally, photothrombotic lesions of S1 show that somatosensory-evoked global SWs depend on bilateral recruitment of homotopic primary somatosensory cortices. Specifically, unilateral lesions of S1 disrupt somatosensory-evoked global SW initiation from either hemisphere, while spontaneous SWs are largely unchanged. These results show that evoked SWs may be triggered by bilateral activation of specific, homotopically connected cortical networks.

The (Un)Conscious Mouse as a Model for Human Brain Functions: Key Principles of Anesthesia and Their Impact on Translational Neuroimaging

In recent years, technical and procedural advances have brought functional magnetic resonance imaging (fMRI) to the field of murine neuroscience. Due to its unique capacity to measure functional activity non-invasively, across the entire brain, fMRI allows for the direct comparison of large-scale murine and human brain functions. This opens an avenue for bidirectional translational strategies to address fundamental questions ranging from neurological disorders to the nature of consciousness. The key challenges of murine fMRI are: (1) to generate and maintain functional brain states that approximate those of calm and relaxed human volunteers, while (2) preserving neurovascular coupling and physiological baseline conditions. Low-dose anesthetic protocols are commonly applied in murine functional brain studies to prevent stress and facilitate a calm and relaxed condition among animals. Yet, current mono-anesthesia has been shown to impair neural transmission and hemodynamic integrity. By linking the current state of murine electrophysiology, Ca2+ imaging and fMRI of anesthetic effects to findings from human studies, this systematic review proposes general principles to design, apply and monitor anesthetic protocols in a more sophisticated way. The further development of balanced multimodal anesthesia, combining two or more drugs with complementary modes of action helps to shape and maintain specific brain states and relevant aspects of murine physiology. Functional connectivity and its dynamic repertoire as assessed by fMRI can be used to make inferences about cortical states and provide additional information about whole-brain functional dynamics. Based on this, a simple and comprehensive functional neurosignature pattern can be determined for use in defining brain states and anesthetic depth in rest and in response to stimuli. Such a signature can be evaluated and shared between labs to indicate the brain state of a mouse during experiments, an important step toward translating findings across species.

Brain segmentation, spatial censoring, and averaging techniques for optical functional connectivity imaging in mice

Resting-state functional connectivity analysis using optical neuroimaging holds the potential to be a powerful bridge between mouse models of disease and clinical neurologic monitoring. However, analysis techniques specific to optical methods are rudimentary, and algorithms from magnetic resonance imaging are not always applicable to optics. We have developed visual processing tools to increase data quality, improve brain segmentation, and average across sessions with better field-of-view. We demonstrate improved performance using resting-state optical intrinsic signal from normal mice. The proposed methods increase the amount of usable data from neuroimaging studies, improve image fidelity, and should be translatable to human optical neuroimaging systems.