Functional connectivity structure of cortical calcium dynamics in anesthetized and awake mice

The interplay between hemodynamic-based markers of cortical activity (e.g. fMRI and optical intrinsic signal imaging), which are an indirect and relatively slow report of neural activity, and underlying synaptic electrical and metabolic activity through neurovascular coupling is a topic of ongoing research and debate. As application of resting state functional connectivity measures is extended further into topics such as brain development, aging and disease, the importance of understanding the fundamental physiological basis for functional connectivity will grow. Here we extend functional connectivity analysis from hemodynamic- to calcium-based imaging. Transgenic mice (n = 7) expressing a fluorescent calcium indicator (GCaMP6) driven by the Thy1 promoter in glutamatergic neurons were imaged transcranially in both anesthetized (using ketamine/xylazine) and awake states. Sequential LED illumination (λ = 454, 523, 595, 640nm) enabled concurrent imaging of both GCaMP6 fluorescence emission (corrected for hemoglobin absorption) and hemodynamics. Functional connectivity network maps were constructed for infraslow (0.009-0.08Hz), intermediate (0.08-0.4Hz), and high (0.4-4.0Hz) frequency bands. At infraslow and intermediate frequencies, commonly used in BOLD fMRI and fcOIS studies of functional connectivity and implicated in neurovascular coupling mechanisms, GCaMP6 and HbO2 functional connectivity structures were in high agreement, both qualitatively and also quantitatively through a measure of spatial similarity. The spontaneous dynamics of both contrasts had the highest correlation when the GCaMP6 signal was delayed with a ~0.6-1.5s temporal offset. Within the higher-frequency delta band, sensitive to slow wave sleep oscillations in non-REM sleep and anesthesia, we evaluate the speed with which the connectivity analysis stabilized and found that the functional connectivity maps captured putative network structure within time window lengths as short as 30 seconds. Homotopic GCaMP6 functional connectivity maps at 0.4-4.0Hz in the anesthetized states show a striking correlated and anti-correlated structure along the anterior to posterior axis. This structure is potentially explained in part by observed propagation of delta-band activity from frontal somatomotor regions to visuoparietal areas. During awake imaging, this spatio-temporal quality is altered, and a more complex and detailed functional connectivity structure is observed. The combined calcium/hemoglobin imaging technique described here will enable the dissociation of changes in ionic and hemodynamic functional structure and neurovascular coupling and provide a framework for subsequent studies of neurological disease such as stroke.