Merge branch 'mcc'

This commit is contained in:
ackman678
2014-08-15 07:20:40 -04:00
15 changed files with 134 additions and 70 deletions

BIN
Figure1.png Normal file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.8 MiB

BIN
Figure2.png Normal file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.9 MiB

BIN
Figure3.png Normal file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.3 MiB

BIN
Figure4.png Normal file

Binary file not shown.

After

Width:  |  Height:  |  Size: 350 KiB

BIN
Figure5.png Normal file

Binary file not shown.

After

Width:  |  Height:  |  Size: 450 KiB

View File

@@ -2,7 +2,10 @@
This file contains Supplementary Figures 1- with legends (pdf; x MB)
## Supplementary Movie 1
The movie shows an example of transcranial cortical activity patterns in Rx-Cre:GCaMP3 mice imaged at P3. Movie is a 60 s long recording played back at 10x (50 fps) as dF/F. Field of view is 6 mm x 8 mm (Quicktime; x MB).
The movie shows an example of transcranial cortical activity patterns in Rx-Cre:GCaMP3 mice imaged at P3. Movie is a 60 s long recording played back at 6x (30 fps) as dF/F. Field of view is 6 mm x 8 mm (Quicktime; 4.5 MB).
## Supplementary Movie 2
The movie shows an example of transcranial cortical activity patterns in Rx-Cre:GCaMP3 mice imaged at P8. Movie is a 60 s long recording played back at 10x (50 fps) as dF/F. Field of view is 11 mm x 13 mm (Quicktime; x MB).
The movie shows an example of transcranial cortical activity patterns in Rx-Cre:GCaMP3 mice imaged at P8. Movie is a 36.8 s long recording played back at 6x (30 fps) as dF/F. Field of view is 10 mm x 12 mm (Quicktime; 5 MB).
## Supplementary Movie 3
The movie shows an example of transcranial cortical activity patterns in SNAP25:GCaMP6 mice imaged at P13. Movie is a 50.2 s long recording played back at 6x (30 fps) as dF/F. Field of view is 20 mm x 24 mm (Quicktime; 5.3 MB).

View File

@@ -1,3 +1,6 @@
## Supplementary Figure Legends
**Supplmentary Figure 1.** **Subnetworks in developing isocortex.** **a** Areal trace examples. Synchronized activations between hemispheres indicated by red links.
**supplementaryFig-symmetry.ai** **Cortical domain activity exhibits bilateral symmetry.** **a** Examples of domains exhibiting spatially symmetric activations. Notice most timepoints contain a mixture of symmetric and asymmetric domain activations. **b** Cortical active fraction timecourses for both hemispheres. **c** Hemispheric domain centers of mass for coactive frames in a recording along medial-lateral (ML) and anterior-posterior (AP) extents. Bottom left panels show the periods indicated by black bars at expanded view. Pearson's correlation: ML, p = 1.1591e-28; AP, p = 7.0982e-07. **d** Hemispheric autocorrelation and cross-correlation functions for cortical activity at all and short time lags. Notice the peaks above gaussian distributed noise (blue traces).

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.4 MiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 620 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 712 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 321 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 631 KiB

View File

@@ -8,11 +8,11 @@
**Calcium signal detection.** Image processing and calcium signal detection was performed using custom software routines written in MATLAB (Mathworks, Natick, MA). The mean pixel intensity at each pixel location, F0 was subtracted and normalized to each frame, Ft of the movie to form a dF/F = (Ft - F0)/F0 array. A background estimate was calculated and subtracted from every frame with a top hat filter using a disk shaped structuring element with radius of 620 µm. Each frame was smoothed with a Gaussian having a standard deviation of 56 µm and a signal intensity threshold was computed using Otsu's method on the histogram of pixel intensities at the 99th percentile from the Sobel gradient transformation of the image array. Calcium signals were automatically segmented as contiguously connected components in space and time using the binary mask for the array from the computed Otsu intensity threshold. Components having an area <50 pixels or a duration of 1 frame were ignored.
**Statistical analysis.** Data sets were analyzed using custom routines written in MATLAB (The Mathworks, Natick, MA) and in R (The R Project for Statistical Computing, http://www.r-project.org). Distribution means were compared using two-sample Student's t-Tests or using ANOVA followed by Tukey's HSD post-hoc test when analyzing the effects of multiple grouping factors (p < 0.01 set as significance). Values are reported as means with the 95% confidence interval or standard error of the mean or medians with the median absolute deviation.
**Statistical analysis.** Data sets were analyzed using custom routines written in MATLAB (The Mathworks, Natick, MA) and in R (The R Project for Statistical Computing, http://www.r-project.org). Distribution means were compared using two-sample Student's t-Tests or using ANOVA followed by Tukey's HSD post-hoc test when analyzing the effects of multiple grouping factors (p < 0.05 set as significance). Values are reported as means with the 95% confidence interval or standard error of the mean or medians with the median absolute deviation.
**Calcium domain analysis.** The mean width in the medial-lateral and height in the rostral-caudal dimensions of the bounding box fitted to each segmented calcium domain signal was taken to be the domain diameter. The number of contiguous frames (bounding box depth) for each segmented calcium domain was taken to be the domain duration. The mean and maximum pixel intensities within each domain were taken as the mean and maximum domain amplitudes. Domains were assigned areal membership by intersection of the domain centroid with a cortical ares's pixel mask. The number of individual domains per recording within a hemisphere or cortical area was taken to be domain frequency.
**Functional correlation analysis.** A binary movie array from all the segmented calcium domain masks for a recording was intersected mask representing different cortical areas. The total number of active pixels per frame expressed as a fraction of possibly active pixels per frame for each cortical area gave active pixel fraction timecourses for each cortical area in each recording. Correlation matrices were calculated for each recording by computing pairwise Pearson's product moment correlation coefficents from the matrix containing the cortical active pixel fraction timecourses. The binarized correlation matrix at *r* > 0.1 was used to form an adjacency matrix with each node representing a cortical area and each edge representing an association between a pair of nodes at weight, *r*. Community structure was detected within each functional association matrix using a fast greedy modularity optimization algorithm [#Clauset:2004] to perform hierarchial clustering using the igraph network analysis software library [#Csardi:2013].
**Functional correlation analysis.** A binary movie array from all the segmented calcium domain masks for a recording was intersected mask representing different cortical areas. The total number of active pixels per frame expressed as a fraction of possibly active pixels per frame for each cortical area gave active pixel fraction timecourses for each cortical area in each recording. Correlation matrices were calculated for each recording by computing pairwise Pearson's product moment correlation coefficents from the matrix containing the cortical active pixel fraction timecourses. The binarized correlation matrix at *r* > 0.15 was used to form an adjacency matrix with each node representing a cortical area and each edge representing an association between a pair of nodes at weight, *r*. Community structure was detected within each functional association matrix using a greedy modularity optimization algorithm [#Newman:2004][#Clauset:2004] to perform hierarchical clustering using the igraph network analysis software library [#Csardi:2013].
[#Ackman:2012]: Ackman, J. B., Burbridge, T. J., and Crair, M. C. (2012). Retinal waves coordinate patterned activity throughout the developing visual system, Nature, 490(7419), 219-25

View File

@@ -116,6 +116,8 @@
[#Espinosa:2012]: Espinosa, J. S. and Stryker, M. P. (2012). Development and plasticity of the primary visual cortex, Neuron, 75(2), 230-49
[#Sanes:1999]: Sanes, J. R. and Lichtman, J. W. (1999). Development of the vertebrate neuromuscular junction, Annu Rev Neurosci, 22(), 389-442
[#Marder:2005]: Marder, E. and Rehm, K. J. (2005). Development of central pattern generating circuits, Curr Opin Neurobiol, 15(1), 86-93
[#Mazzoni:2007]: Mazzoni, A., Broccard, F. D., Garcia-Perez, E., Bonifazi, P., Ruaro, M. E., and Torre, V. (2007). On the dynamics of the spontaneous activity in neuronal networks, PLoS ONE, 2(), e439
@@ -212,8 +214,9 @@
[#Laing:2012]: Laing, R. J., Bock, A. S., Lasiene, J., and Olavarria, J. F. (2012). Role of retinal input on the development of striate-extrastriate patterns of connections in the rat, J Comp Neurol, 520(14), 3256-76
[#Petersson:2003]: Petersson, P., Waldenström, A., Fåhraeus, C., and Schouenborg, J. (2003). Spontaneous muscle twitches during sleep guide spinal self-organization, Nature, 424(6944), 72-5
[#Mohns:2008]: Mohns, E. J. and Blumberg, M. S. (2008). Synchronous bursts of neuronal activity in the developing hippocampus: modulation by active sleep and association with emerging gamma and theta rhythms, J Neurosci, 28(40), 10134-44
@@ -266,6 +269,9 @@ http://vision.ucsd.edu/~pdollar/toolbox/doc/index.html
[#Kobayashi:1963]: Kobayashi, T. (1963). Brain-to-body ratios and time of maturation of the mouse brain, Am J Physiol, 204(), 343-6
[#Newman:2004]: Newman, M. E. J. (2004). Fast algorithm for detecting community structure in networks, Phys. Rev. E, 69(), 066133
[#Bullmore:2009]: Bullmore, E. and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems, Nat Rev Neurosci, 10(3), 186-98
[#Ackman:2014]: Ackman, J. B. and Crair, M. C. (2014). Role of emergent neural activity in visual map development, Curr Opin Neurobiol, 24C(), 166-175

View File

@@ -2,7 +2,7 @@ Author: James B. Ackman
Date: 2014-06-05 00:38:46
Tags: paper, draft, manuscript, literature, research, #results, retinal waves, spontaneous activity, development, calcium domains
# Structured population activity across developing neocortex
# Structured dynamics of neural activity across developing neocortex
**Authors and Affiliated Institutions:**
James B. Ackman¹, Hongkui Zeng², and Michael C. Crair¹
@@ -12,24 +12,26 @@ James B. Ackman¹, Hongkui Zeng², and Michael C. Crair¹
# Summary
The cerebral cortex exhibits spontaneous and sensory evoked patterns of activity during fetal and postnatal development that are crucial for the activity-dependent formation and refinement of neural circuits [#Katz:1996]. Knowing the source and flow of these activity patterns locally and globally is crucial to understanding self-organization in the developing brain. Here we show that neural population activity within newborn mice in vivo is characterized by spatially discrete domains that are coordinated in a state dependent and areal dependent fashion throughout developing neocortex. Whole brain optical recordings from neonatal mice expressing a genetic calcium reporter showed that ongoing activity in the cerebral cortex was characterized by distinct and repetitively active domains measuring hundreds of microns in diameter. Domain activity exhibited mirror-symmetric patterns between the hemispheres, with strong correlations between specific portions of frontal and parietal cortex. Ongoing activity across the cortical hemispheres showed characteristic network architectures with a frontal-motor regions functionally connected to a parietal-sensory areas through secondary motor/cingulate cortex, retrosplenial cortex, and posterior parietal cortex. Furthermore, ongoing activity was regulated by physiological state with cortical regions exhibiting areal dependent coordination of activity with motor behavior differentially during the course of development. This study provides the first comprehensive description of population activity in the developing neocortex at a scope and scale that bridges the microscopic or macroscopic spatiotemporal resolutions provided by traditional neurophysiological or neuroimaging techniques. Mesoscale maps of cortical population dynamics within animal models will be vital to future engineering of repair strategies and brain-machine interfaces for neurodevelopmental disorders.
The cerebral cortex exhibits spontaneous and sensory evoked patterns of activity during development that are vital for the activity-dependent formation and refinement of neural circuits. Identifying the source and flow of these activity patterns locally and globally is vital to understanding self-organization in the developing brain. Here we use whole brain transcranial optical imaging to show that the dynamical patterns of neuronal activity in developing mouse neocortex consists of spatially discrete domains that are coordinated in an age, region, and state- dependent fashion. Ongoing cortical activity displayed mirror-symmetric activation patterns across the cerebral hemispheres and showed characteristic network architectures that were shaped during development, with frontal-parietal areas functionally connected to occipital regions regions through cingulate and motor cortex. This study provides the first broad description of population activity in the developing neocortex at a scope and scale that bridges the microscopic and macroscopic spatiotemporal resolutions provided by traditional neurophysiological or functional neuroimaging techniques. Mesoscale maps of cortical population dynamics within animal models will be crucial for future efforts to understand and treat neurodevelopmental disorders.
<!-- Furthermore, ongoing activity was regulated by physiological state with cortical regions exhibiting areal dependent coordination of activity with motor behavior differentially during the course of development. -->
# Introduction
Brain development requires neural activity and calcium dynamics for establishing proper circuit structure and function. The importance of neural activity in the prenatal and postnatal period can be easily recognized in children exposed to chemical agents affecting neurotransmission during the fetal period that result in severe brain malformations, epilepsy, and mental retardation. Indeed, embryonic limb movements in species ranging from chick to human are thought to be initiated by spontaneous motor neuron activity in the spinal cord and thought to be crucial for activity-dependent development of motor synapses [Schoenberg:2003] [Marder,Lichtmann]. However it is only recently that we have begun to appreciate the underlying patterns of persistent neural activity that exist in the developing brain in vivo. For example, sensori-motor feedback associated with spontaneous movement generated by spinal motor neurons triggers synchronized 'spindle-burst' potentials among cells in somatosensory cortex [Yang:2009][Khazipov:2004a] before the start of locomotion and tactile behavior. Correlated bursts of activity occur in the developing rat hippocampus in vivo [#Leinekugel:2002] [Mohns&Blumberg]. Spontaneous retinal waves drive patterned activation of circuits throughout the immature visual system before the onset of vision [#Ackman:2012] [Hanganu,Colonnese?]. Furthermore, prenatal EEG recordings have demonstrated spindle burst oscillations and slow activity transients in the human infant somatosensory and occipital cortices before birth [#Vanhatalo:2005][#Tolonen:2007]. However, a comprehensive account of the dynamical patterns of persistent activity across the developing neocortex in vivo has not been undertaken, largely because a method to assess neural activity between most cortical areas simultaneously and non-invasively has not been available.
Brain development requires neural activity for establishing proper circuit structure and function [#Katz:1996]. Fetal movements, prenatal electroencephalographic oscillations [#Vanhatalo:2005][#Tolonen:2007], and sensitivity to disruptions in periphreal inputs affecting neurotransmission all underscore the presence and importance of neural activity in the developing brain. Indeed, embryonic limb movements in species ranging from chick to human are thought to be initiated by spontaneous motor neuron activity in the spinal cord and thought to be crucial for activity-dependent development of motor synapses [#Sanes:1999][#Petersson:2003][#Marder:2005]. In the visual system spontaneous waves of activity originating eye, 'retinal waves', have long been studied as a model for activity-dependent circuit development before the start of sensory experience [#Ackman:2014]. However it is only recently that we have begun to appreciate the actual nature of persistent neural activity patterns as they exist in the developing brain in vivo. For example, sensori-motor feedback associated with spontaneous movement generated by spinal motor neurons triggers synchronized 'spindle-burst' potentials among cells in somatosensory cortex [#Khazipov:2004a][#Yang:2009] before the start of locomotion and tactile behavior. Correlated bursts of activity occur among neurons in the developing rat hippocampus in vivo [#Leinekugel:2002][#Mohns:2008]. Spontaneous retinal waves drive patterned activation of circuits throughout the visual system before the onset of vision [#Ackman:2012][#Colonnese:2010] and provide spatiotemporal information suitable for organizing connections within and between different visual areas. However, a comprehensive account of the dynamical patterns of persistent activity across the developing neocortex in vivo has not been undertaken, largely because a method to assess neural activity between cortical areas simultaneously and non-invasively has not been available.
# Results
## Ongoing activity in developing neocortex is characterized by discrete domains
We performed transcranial optical recordings from mice expressing the genetic calcium reporter GCaMP (GCaMP3 or GCaMP6) throughout cortical neurons to assess neural population activity patterns at macroscopic scale (millimeters) and with mesoscopic spatial and temporal resolution (10s of microns and 100s of milliseconds). We performed our recordings in three age groups during the first two postnatal weeks during which the mouse brain develops to >90% of its adult weight [#Kobayashi:1963]: P2-P5, P8-P9, and P12-13.
We performed transcranial optical recordings from mice expressing the genetic calcium reporter GCaMP (GCaMP3 or GCaMP6) throughout cortical neurons to assess neural population activity patterns at macroscopic scale (millimeters) and with mesoscopic spatial and temporal resolution (10s of microns and 100s of milliseconds). We performed our recordings in three age groups throughout the first two postnatal weeks during which the mouse brain attains >90% of its adult weight [#Kobayashi:1963]: P2-P5, P8-P9, and P12-13.
Functional mesoscale optical imaging (fMOI) revealed that supracellular cortical activity patterns were characterized by discrete domains of activation (Fig. 1a-c) [Supplementary Movie 1](../wholeBrain_blob/ackmanWholeBrainGcampP3.mov). These activity domains ranged from 250 - 976 µm in diameter and 0.4 - 2.6 s in duration <!--(10-90th percentiles)-->(Fig. 1e-h) (Table 1.) . The duration of cortical domain activations was not significantly affected by age (F = 0.933, p = 0.428, r^2 = 0.00567) or by hemisphere (F = 0.017, p = 0.900) (P2-5, N = 15653; P8-9, N = 70189; P12-13, N = 120214 domains) (Fig. 1e,f). There was a significant effect of age on the diameter of cortical domain activations (F = 25.788, p = 0.000188, r^2 = 0.1277), but not hemisphere (F = 0.192, p = 0.671808) (Fig. 1g,h). The frequency with which cortical domain activations occurred increased with age (F = 29.562, p = 8.86e-12, r^2 = 0.2535) and did not differ significantly between the hemispheres (F = 0.012, p = 0.911) (P2-5, N = 22; P8-9, N = 30; P12-13, N = 38 movies/hemi) (Fig. i,j) (Table 1).
The neocortex exhibits a characteristic modular organization across the cortical surface such that vertical arrays of cells concerned with specific sensory features are grouped together as columns in a topographic fashion [#Mountcastle:1997]. Most evidence suggests that cortical columns range from 300-600µm diameter, even between species whose brain volumes differ by a factor of 10^3 [#Mountcastle:1997]. Functional mesoscale optical imaging (fMOI) revealed that supracellular cortical activity patterns were characterized by discrete domains of activation (Fig. 1a-c) ([Supplementary Movie 1](../wholeBrain_blob/SupplementaryMovie-P3gcamp3.mov)). These activity domains ranged from 250 - 976 µm in diameter and 0.4 - 2.6 s in duration <!--(10-90th percentiles)-->(Fig. 1e-h) (Table 1). The duration of cortical domain activations was not significantly affected by age (F = 0.933, p = 0.428, r^2 = 0.00567) or by hemisphere (F = 0.017, p = 0.900) (P2-5, N = 15653; P8-9, N = 70189; P12-13, N = 120214 domains) (Fig. 1e,f). There was a significant effect of age on the diameter of cortical domain activations (F = 25.788, p = 0.000188, r^2 = 0.1277), but not hemisphere (F = 0.192, p = 0.671808) (Fig. 1g,h). The frequency with which cortical domain activations occurred increased with age (F = 29.562, p = 8.86e-12, r^2 = 0.2535) and did not differ significantly between the hemispheres (F = 0.012, p = 0.911) (P2-5, N = 22; P8-9, N = 30; P12-13, N = 38 movies/hemi) (Fig. i,j) (Table 1) ([Supplementary Movie 2](../wholeBrain_blob/SupplementaryMovie-P8gcamp3.mov)).
The neocortex exhibits a characteristic modular organization across the cortical surface such that vertical arrays of cells concerned with specific sensory features are grouped together as columns in a topographic fashion [#Mountcastle:1997]. Most evidence suggests that cortical columns range from 300-600µm diameter, even between species whose brain volumes differ by a factor of 10^3 [#Mountcastle:1997]. Its intriguing that we found the size of cortical domains to be centered on this range at early ages, because this is in agreement with previous work showing that population activity in neonatal rat barrel cortex maps onto ontogenetic modules centered on each barrel column [#Yang:2012a] and barrels are an archetypical model for columnar cortical function in rodent. Indeed, we found a cortical area in primary somatosensory cortex at P2-5 where cortical domain activations group into rows and individual modules that match primary barrel cortex structure (Fig. 1c) (Supplementary Fig.). This indicates that early cortical activity in some cortical areas is matched to the size the functional cortical modules that are thought to be the fundamental procdessing unit of the cerebral cortex.
Its intriguing that we found the size of cortical domains to be centered on this range at early ages, because this is in agreement with previous work showing that population activity in neonatal rat barrel cortex maps onto ontogenetic modules centered on each barrel column [#Yang:2012a] and barrels are an archetypical model for columnar cortical function in rodent. Indeed, we found a cortical area in primary somatosensory cortex at P2-5 where cortical domain activations group into rows and individual modules that match primary barrel cortex structure (Fig. 1c) (Supplementary Fig. of zoomed maxproj image). This indicates that early activity in developing cortical areas can be matched to the size the functional columns thought to be the fundamental processing unit of the cerebral cortex.
![ **Figure 1.** Calcium domains throughout neonatal mouse neocortex. **a** Experimental schematic. **b** Single image frame showing calcium domains in both hemispheres at postnatal day 3 (P3) and automatically detected domain masks. **c** Centroid positions for segmented domain masks from a 10 min recording. Points are overlaid on a reference map of primary sensory areas determined by thalamocortical inputs (red outlines). Notice rows of whisker barrels are evident in the structure of domain centroid positions. **d** Functional activity map at P3. Based on pixel activation frequency from all detected domains in a single 10 min recording. Map is overlaid on cortical areal parcellations. Notice localized maxima and minima of functional activity between areas that approximate known anatomical cortical area boundaries and the mirroring of map structure bilaterally. **e** Mean domain duration maps from 3 SNAP25-Ai103 mice. **f** Histograms showing domain durations distributions in the P2-5, P8-9, and P12-13 age groups and by cortical hemisphere (L, R). **g** Mean domain diameter maps from same 3 mice in e. **h** Histograms showing the distributions of domain diameters. **i** Mean domain frequency maps from same 3 mice in e. **j** Boxplot distributions of hemispheric domain frequencies.](figure1.png)
![ **Figure 1.** Calcium domains throughout neonatal mouse neocortex. **a** Experimental schematic. **b** Left panel: Single image frame showing calcium domains in both hemispheres at postnatal day 3 (P3) and the mask of detected domain signals. Middle and right panels: Time projection map from a raw dF/F movie segment and the corresponding map from automatically detected domain masks. Notice the individual domains of activity in the area of barrel cortex (arrow) **c** Centroid positions for segmented domain masks from a 10 min recording. Points are overlaid on a reference map of primary sensory areas determined by thalamocortical inputs (red outlines). Notice rows of whisker barrels are evident in the structure of domain centroid positions (arrow). **d** Functional activity map at P3. Based on pixel activation frequency from all detected domains in a single 10 min recording. Map is overlaid on cortical areal parcellations. Notice localized maxima and minima of functional activity between areas that approximate known anatomical cortical area boundaries and the mirroring of map structure bilaterally. **e** Mean domain duration maps from 3 SNAP25-Ai103 mice. **f** Histograms showing domain durations distributions in the P2-5, P8-9, and P12-13 age groups and by cortical hemisphere (L, R). **g** Mean domain diameter maps from same 3 mice in e. **h** Histograms showing the distributions of domain diameters. **i** Mean domain frequency maps from same 3 mice in e. **j** Boxplot distributions of hemispheric domain frequencies.](Figure1.png)
| | duration (s) | diameter (µm) | frequency (hemisphere-min^-1) |
| ------ | ------------ | ----------------- | ----------------------------- |
@@ -41,72 +43,66 @@ The neocortex exhibits a characteristic modular organization across the cortical
[ **Table 1: Domain statistics**]
We examined how the spatiotemporal properties of cortical domains vary among different cortical regions by parcellating the brain into distinct anatomical boundaries using reference coordinates from a mouse line that expressed a tdtomato reporter in thalamocortical afferents at P7 (Fig. 1c,d) (Supplementary Fig.). Patterns of thalamocortical axon terminals can be used to map out areal boundaries of primary sensory cortical areas [wong riley 1979]. We matched these parcellations to a Allen brain atlas adult mouse reference image and than linearly scaled the cortical area reference boundaries for each animal to maps containing functional boundaries for barrel cortex and visual cortex where spontaneous retinal waves functionally map out developing visual areas [#Ackman:2012] (Fig. 1c-e,g,i).
## Cortical dynamics differs with area and age
Cortical domain frequency among different regions scaled as a function of net cortical area and this association became stronger during the course of development (Fig. 2a). The most frequently cortical regions at each age group when normalized to the amount of total amount of cortical space was the limb/trunk representations in somatosensory cortex (Fig. 1i, Supplementary Fig.). Generally, the frequency of activity was remarkably uniform across cortical areas at each age of development (Supplementary Fig) indicating a homeostatic mechanism regulating global activity levels. The long tails in the domain duration and diameter distributions at P2-5 and P8-9 (Fig. 1f,h) were dominated by retinal wave driven cortical activity in V1 that lasted on the order of seconds to tens of seconds (Fig. 1e, Fig. 2b,c), but also by long lasting wave-like activations occurring in motor cortex (Fig. 1e, Fig. 2b,c). Indeed the cortical regions with the highest wave motion indices were V1 and M1 at P2-5, with V1 continuing to have the highest index at P8-9 and then dropping to mean motion idx level similar to other cortical regions at P12-13. The diameter of domain activation became larger among cortical regions during the second postnatal week including the S1-limb/body regions where at P13 a small subpopulation of events had mean diameters approaching that of the entire hemisphere and a higher wave motion index (Fig. 2d-f) (x% of all events, ~2/10min) [Supplementary Movie 2](../wholeBrain_blob/ackmanWholeBrainImaging-lo.mov) (Fig. 2d). These global population events synchronized activity across cortical areas and had centers of mass that were concentrated near the middle of each hemisphere in the S1-limb/body area.
We examined how the spatiotemporal properties of cortical domains vary among different cortical regions by parcellating the brain into distinct anatomical boundaries using reference coordinates from a mouse line that expressed a tdtomato reporter in thalamocortical afferents at P7 (Fig. 1c,d) (Supplementary Fig.). Patterns of thalamocortical axon terminals outline primary sensory cortical areas [#Lebrand:1996] during mouse postnatal development. We aligned these parcellations to the Allen brain mouse atlas and then scaled the resulting cortical area reference coordinates to match activity maps from each animal containing functional boundaries for barrel cortex and visual cortex where spontaneous retinal waves functionally map out developing visual areas [#Ackman:2012] (Fig. 1c-e,g,i).
Cortical domain frequency among different regions scaled as a function of net cortical area and this association became stronger during the course of development (Fig. 2a). The most frequently active cortical regions at each age group when normalized to the amount of total amount of cortical space was the limb/trunk representations in somatosensory cortex (Fig. 1i, [SupplementaryFig-areal-stats.ai](../wholeBrain_blob/SupplementaryFig-areal-stats.ai)). In contrast to diameter and duration, the frequency and amplitude of cortical domain activity was remarkably uniform across areas at each age of development ([SupplementaryFig-areal-stats.ai](../wholeBrain_blob/SupplementaryFig-areal-stats.ai)) indicating a homeostatic regulation of global activity levels. The long tails in the domain duration and diameter distributions at P2-5 and P8-9 (Fig. 1f,h) were dominated by retinal wave driven cortical activity in V1 that lasted on the order of seconds to tens of seconds (Fig. 1e, Fig. 2b,c), but also by long lasting wave-like activations occurring in motor cortex (Fig. 1e, Fig. 2b,c). Indeed the cortical regions with the highest wave motion indices were V1 and M1 at P2-5, with V1 continuing to have the highest index at P8-9 and then dropping to mean motion idx level similar to other cortical regions at P12-13. The diameter of domain activation became larger among cortical regions during the second postnatal week including the S1-limb/body regions where at P13 a small subpopulation of events had mean diameters approaching that of the entire hemisphere and a higher wave motion index (Fig. 2c-f) (x% of all events, ~2/10min) ([Supplementary Movie 3](../wholeBrain_blob/SupplementaryMovie-P13gcamp6.mov)). These global population events synchronized activity across cortical areas and had centers of mass that were concentrated near the middle of each hemisphere in the S1-limb/body area.
![ **Figure 2.** Spatiotemporal characteristics of cortical domains. **a** Domain frequency as function of cortical area size. **b** Scatterplots of domain diameter and duration. **c** Time projection color maps of waves in visual cortex and motor cortex at P5. **d** Time projection color maps of interareal activations at P13. **e** Scatterplots of wave motion index as function of domain diameter. **f** Mean wave motion index over development.](figure2.png)
![ **Figure 2.** Spatiotemporal properties of cortical domains. **a** Domain frequency as function of cortical area size. **b** Scatterplots of domain diameter and duration. **c** Time projection maps of waves in motor cortex at P3, visual cortex at P5, and occipital-parietal-frontal cortex at P13. **d** Scatterplots of wave motion index as function of domain diameter. **e** Mean wave motion index over development.](Figure2.png)
## Cortical activity is coordinated with motor behavior
## Cortical domain activity is coordinated with motor behavior
Next we assessed mesoscale cortical activity patterns as a function of physiological state and motor behavior. It has previously been demonstrated that general anesthesia abolishes spontaneous retinal wave activity in visual system [#Ackman:2012] and spontaneous activity in entorhinal cortex [#Adelsberger:2005]. We found that during anesthesia induction, there is rapid (<60 s) knock down of cortical activity (Supplementary Movie) (Supplementary Fig) at all ages. While in neonates, no cortical activity was found during general anesthesia, at P12-13 within ~10-20min after general induction we found altered spontaneous patterns, with short duration, large diameter population activities synchronizing multiple cortical regions. (Supplementary Movie). The continued spinal motor activity during early anesthesia and the altered cortical activity patterns that ensue under anesthesia at P12-13 suggest a maturational dependence of isoflurane anesthesia on neural activity that affects brain regions differentially during development.
Next we assessed mesoscale cortical activity patterns as a function of physiological state and motor behavior. It has previously been demonstrated that general anesthesia abolishes spontaneous retinal wave activity in visual system [#Ackman:2012]. We found that during anesthesia induction, there is rapid (<60 s) knock down of cortical activity (Supplementary Movie) (Supplementary Fig) at all ages. While in neonates, no cortical activity was found during general anesthesia in neonates, at P12-13 we found altered spontaneous patterns, with short duration, large diameter population activities synchronizing multiple cortical areas. (Supplementary Movie)
Variation in the strength of correlation between cortical areas and the motor movement signal depended on brain region (p < 2.2e-16, anova) and age (p = 1.627e-05, anova) The first age group in which motor cortex exhibited signficant positive correlation with motor movements was at P12-13 (r=0.06±0.02, p-value = 0.001449, t-test).
![ **Figure 3**. Cortical domains are state dependent. **a** Experimental schematic. Red light illumination measured with a photodiode (PD) was used to monitor motor activity. **b** Cortical activity (active fraction) in each hemisphere after onset of gas anesthetic. **c** Cortical activity and coincident motor activity signals. Gray shading indicates active and quiet motor periods determined by the half-rise and decay times from peak in the low pass filtered motor signal. Active pixel fraction traces for motor (M1,M2), somatosensory (HL,FL,T; barrel), and visual (V1,V2) cortex shown at bottom of panel. Red links show synchronized motor movements and brain activity with different cortical regions. **d** Single frame domain masks for times indicated in **c**. **e** Pixel activation frequency maps during quiet and active motor periods. **f** Cross-correlation functions between cortical regions and motor movement signals. Notice the general negative correlation between motor activity and all cortex activity signals (r = , p = ) and the high positive correlation between motor and S1-limb/body signals (r = 0.3019, p < 2.2e-16).](figure3.png)
We monitored motor movements simultaneously with cortical activity during our fMOI recordings to gain insight into the relationship between motor behavior output to cerebral cortical dynamics during development. The highest levels of synchronized cortical domain activity occurred during periods of relatively sparse motor behavior whereas the lowest levels of synchronized cortical activity occurred during periods of increased motor movement (Fig. 3c-e). Variation in the strength of correlation between cortical areas and the motor movement signal depended on both brain region (p < 2.2e-16, anova) and age (p = 1.627e-05, anova) (Fig 3c-f). Interestingly, the first age group in which motor cortex exhibited signficant positive correlation with motor movements was at P12-13 (r=0.06±0.02, p-value = 0.001449, t-test) (Fig. 3f). We hypothesized that just before eye opening around P11-P13 there will be a shift with significant zero lag or preceding correlation between motor cortex and the motor movement signals perhaps conincciding with teh begining fo ggoal directed behavior. Motor and state dependent behavior surprisingly complex, even in neonates.
## Cortical activity is mirrored between the hemispheres
* Inter hemispheric functional connectivity, importance for autism, schizophrenia. Maybe an activity-dependent mechanism for commisural connectivity.
* olavarria work, evidence for inter hemispheric activity dependence
* [#Hanganu:2006], 30% of spindle bursts correlated across hemispheres
To understand the patterns and how they interact we first looked at correlation between the hemispheres. Cortical activity exhibited high temporal correlation between the hemispheres () . In additon this activity was highly correlated in the spatial dimension. We found that activity was correlated in anterior-posterior and medial-lateral directions. It exhibited mirror symmetric and non-mirror symmetric patterns. For example epochs of time would exhibit high correlation in the medial-lateral dimension or in the rostral-caudal dimension. This strength of correlation temporally and spatially increased between the hemsipheres with a function of age.
<!-- * Each hemisphere 'training' the other one in preparation for behaviorally relevant sensory-motor imitations '[[mirror_neurons]]' hypothesis? -->
![ **Figure 4.** Cortical domain activity exhibits bilateral symmetry. **a** Examples of domains exhibiting spatially symmetric activations. Notice most timepoints contain a mixture of symmetric and asymmetric domain activations. **b** Cortical active fraction timecourses for both hemispheres. **c** Hemispheric domain centers of mass for coactive frames in a recording along medial-lateral (ML) and anterior-posterior (AP) extents. Bottom left panels show the periods indicated by black bars at expanded view. Pearson's correlation: ML, p = 1.1591e-28; AP, p = 7.0982e-07. **d** Hemispheric autocorrelation and cross-correlation functions for cortical activity at all and short time lags. Notice the peaks above gaussian distributed noise (blue traces).](figure4.png)
![ **Figure 3**. Cortical domain activity is state dependent. **a** Experimental schematic. Red light illumination measured with a photodiode (PD) was used to monitor motor movements. **b** Cortical activity (active pixel fraction by hemisphere) and motor movement signal after onset of isoflurane anesthesia at 205 s. **c** Time projection maps (40 s segments) at times indicated in recording from **b**. **d** Cortical activity and coincident motor activity signals. Gray shading indicates active and quiet motor periods determined by the half-rise and decay times from peak in the low pass filtered motor signal. Active pixel fraction traces for motor (M1,M2), somatosensory (HL,FL,T; barrel), and visual (V1,V2) cortex shown at bottom of panel. Red links show synchronized motor movements and brain activity with different cortical regions. **e** Single frame domain masks for times indicated in **d**. **f** Pixel activation frequency maps during quiet and active motor periods. **g** Cross-correlation functions between cortical regions and motor movement signals. Notice the general weak correlation between motor movement and all cortical activity signals (r = , p = ) and the high positive correlation between motor movement and S1-limb/body signals (r = 0.3019, p < 2.2e-16).](Figure3.png)
## Cortical activity comprises distinct subnetworks
## Developing cortical activity consists of distinct subnetworks
The complex activity dynamics among nascent neocortical networks suggests that spatiotemporal correlations between areas may exist that provide information for development of intra- and inter-hemispheric connections. Indeed, recent work has suggested that neural activity is required for the migration of some interneuron subtypes (fishell work; ZJ. Huang work) and the development of callosal connections (kir2.1 eporation paper; olavarria work?). Since the timing of neural activity patterns is thought to be important for various aspects of circuit development it remains crucial to understand the correlational structure of ongoing cortical activity between brain regions. Thus to achieve a better understanding of the early activity patterns that may regulate interactions between cortical regions we first looked at correlation between the hemispheres. Cortical activity exhibited high temporal correlation between the hemispheres ([SupplementaryFig-symmetry.ai](../wholeBrain_blob/SupplementaryFig-symmetry.ai)). In addition this activity was highly correlated in the spatial dimension. We found that activity was correlated in anterior-posterior and medial-lateral directions. It exhibited mirror symmetric and non-mirror symmetric patterns. For example epochs of time would exhibit high correlation in the medial-lateral dimension or in the rostral-caudal dimension. This strength of correlation temporally and spatially increased between the hemsipheres with a function of age.
We then calculated a matrix of pearsons correlation coefficients based on the pixel active fraction timecourses for each pair of parcellations. The resulting assocaition matrix was run through a hierarchal clustering alogtithm to reveal functional modules of of activation. These functional modules typically consisted of 3 distinct subnetworks-- a frontal motor network, a posterior parietal network, a S1-body/limb network, and an auditory A1 network at P12.
We computed a correlation matrix for each recording based on the pixel active fraction timecourses for each pair of cortical parcellations. Community structure among cortical nodes in the pairwise association matrix at P12-13 was detected with hierarchical clustering based on optimization of the graph modularity score, which measures the number of edges that fall within groups minus the number expected by chance [#Newman:2004]. The resulting community dendrogram was used to order the mean correlation matrices for each age group (Fig. 4a) (N=). This revealed a non-random network organization and a grossly similar correlational structure present across development (Fig. 4a). There were 4 primary subnetworks detected at P12-13-- motor-S1-face (green), S1-body (red), medial--visual (blue), and auditory (purple). Indeed the 4 network modules detected at 12-13 display as clusters along the diagonal in the correlation matrix and these are apparent in the earlier age groups as well. However, functional correlation between cortical regions was highly dynamic across development, with both increases and decreases in correlation depending on brain area and age. Remarkably, detection of community structure in each age group independently showed similar vertex memberships in the 3 largest modules, with developmental switches in module membership occurring in PPC and M2, two brain regions known to integrate multimodal sensory and motor planning information (Fig. 4b).
We found many similarities but some striking differences as a function of age.
The strongest correlations at each age typically occurred at symmetrical regions between the hemispheres in RSA, M1, M2, and T (Fig. 4a,b). In primary sensory areas such as A1, V1, and S1-barrel, inter-hemispheric correlation was initially weak and became strong at P12-13. At P8-9, correlations among cortical regions became more spread out, with strong correlations becoming stronger and weak correlations becoming weaker. There was increased positive correlation within modules and stronger negative correlation between modules at P12-13. For example, visual regions exhibited increased negative correlation with the body-parietal subnetwork, A1 showed increased negative correlations with multiple areas, and motor-parietal regions exhibited increased intra-module correlations.
![ **Figure 5.** Subnetworks in developing neocortex. **a** Correlation matrix of domain activity among cortical areas. Colormap indicates Pearson's r correlation coefficient values. **b** Graph of functional connections for r > 0.15. Community structure detected with hierarchical clustering are indicated in 4 colors. Notice that S1-limb/body regions form a distinct subgraph (red). **c** Map of cortical area associations for r > 0.1.](figure5.png)
![ **Figure 4.** Functional architecture of developing neocortex. **a** Group averaged correlation matrices of domain activity among cortical areas. Colormap indicates Pearson's r correlation coefficient values. Dendrogram and node order from community structure detected with hierarchical clustering in the P12-13 group. **b** Map of cortical area associations for r > 0.15. Node colors represent cortical communities detected with clustering within each age group. Edge width indicates the squared connection weight (r^2). Note both similarities in module membership and increased connection strength with age.](Figure4.png)
We analysed the topological properties of developing mouse cortical networks with graph theoretic measures used in fMRI studies of large-scale functional connectivity [#Bullmore:2009]. Graphs of mean functional connectivity illustrated decreased randomness and tighter clustering among cortical modules during network development (Fig. 5a). The global clustering coefficient-- calculated as the mean of local clustering coefficients that reflect how highly interconnected each node's neighbors are-- significantly increased during development (Fig. 5c). In contrast, the average network diameter and shortest path changed little over the course of development. Both of these measures were significantly different than that of equivalent random networks at each age group. The mean degree (number of edges per node) was higher at P12-13 and the degree distribution shifted during development to match the node order seen at P12-13 (Fig. 5b) (Supplementary Fig graph metrics). The mean node strength--the sum of all edge weights per node-- was higher at P12-13 and the rank order of node strengths also changed to that of the P12-13 distribution during development. The 4 cortical areas having both the highest degree and node strength were M1, M2, PPC, and V2M (Fig. 5b).
Measures of node centrality were computed to identify potential network hubs. Betweeness centrality measures the fraction of all shortest paths in a network that pass through a node, therefore identifing important throughputs in the network. M1, M2, PPC, and RSA had the highest betweeness centrality scores, which generally decreased during development. Eigenvector centrality is proportional to the sum of centralities for a node's connections, giving high scores to nodes that are linked to many other highly connected nodes. High scores (near the top of the linear fit) indicated hub nodes which included M1 and M2. Deviations from betweeness-eigenvector linearity can indicate whether a node is potentially more important as a network throughput (higher betweeness) or a driver node (higher eigenvector centrality). These results indicate that greater local connectivity together with higher connection strengths are key features in development of global network architecture of the cerebral cortex. Furthermore, the developing mouse cerebral cortex shares features of network topology that are similar to those seen in resting states networks in human infant. Networks having both high clustering and short average path lengths are key features that distinguish small world networks from random or lattice based networks-- , the clustering coefficient and average path length in the mean graphs from each age group. While the average path length decreased only slightly over the course of development, the clustering coefficient increased significantly with age.
![ **Figure 5.** Dynamics of functional connectivity in developing neocortex. **a** Graph of functional connections for r > 0.15. Node colors represent cortical communities detected with clustering within each age group. **b** Boxplots of degree (number of links) and node strength (sum of connection weights) by cortical area. The distributions become increasingly ordered like the P12-13 group with age. **c** Boxplots of clustering coefficient and average path length by recording. **d** Scatterplots of mean network centrality scores by cortical area.](Figure5.png)
# Conclusions
* Neural population activity constitutes discrete spatial and temporal activations among developing cortical areas
* Cortical activity exhibits symmetrical spatial and temporal activations across the hemispheres
* Cortical activity is coordinated with motor behavioral state in an areal dependent fashion
* Developing cerebral networks comprise distinct functional modules among cortical areas
We have provided the large-scale account of spatially discrete neural population activity among the developing cortical hemispheres in vivo. In contrast to classical in vitro imaging studies (Yuste, Katz)
* **Ongoing activity in developing cortex is not random** -- it is coordinated in space and time within and between the hemispheres among cortical areas
* Functional mesoscale imaging technique as a template for assessing altered functional dynamics in models for neurological disorders
We found that cortical activity was coordinated with motor behavioral state in an areal dependent fashion. This is consistent with reports of twitch activated LFPs called spindle bursts in S1 and M1 during development (Khazipov 2004; Yang & luhmann, Luhman JNS 2014). The high spatial resolution of our fMOI recordings provide the first evidence that
Cortical activity exhibited symmetrical spatial and temporal activations across the hemispheres. This is consistent with the stronger interhemsipheric resting state connectivity between homotopic cortical regions in the human infant. Increased intra-hemispheric connectivity and modularity is thought to be key features of functional connectivity development in large scale brain networks.
Developing cerebral networks comprised distinct functional modules among cortical areas. This indicated that functional network identities may be shaped early in development. The inter-module dynamics we found may have important relevance for cortical plasticity seen after sensory deprivation or traumatic injury early in development-- such as that seen in enucleated monkeys (H. kennedy experiments), children born blind or deaf, or epilepsy patients undergoing hemispherectomies. Given that altered inter- and intra- hemispheric functional connectivity is thought to be relevant in autism and schizophrenia extension of the work performed in this study to mouse models for these neurological diseases will be interesting.
This study demonstrates that ongoing activity in developing cortex is not random -- it is coordinated in space and time between hemispheric networks in the neocortex. Furthermore, functional mesoscale optical imaging will be useful in assessing potentially altered functional connectivity dynamics in animal models for neurological disorders.
# Methods Summary
Anesthetized Rx-Cre:GCaMP3 or SNAP25-GCaMP6 mice between postnatal day 2 to 13 (P2-P13) were were prepared for transcranial optical imaging. Calcium imaging was performed in vivo using wide-field epifluoresence microsopy using a DC-Hg2+ lamp, 1x macro objective, and pco.edge sCMOS camera after a 1 hour recovery period from general anesthesia. Calcium transients and waves were detected using custom MATLAB routines.
Anesthetized Rx-Cre:GCaMP3 or SNAP25-GCaMP6 mice between postnatal day 2 to 13 (P2-P13) were were prepared for transcranial optical imaging. Calcium imaging was performed in vivo using wide-field epifluoresence microsopy using a DC-Hg2+ lamp, 1x macro objective, and pco.edge sCMOS camera after a 1 hour recovery period from general anesthesia. Automated image segmentation and calcium event detection was performed using custom MATLAB routines.
**Full methods** and any associated references are available in the online version of the paper at www.nature.com/nature
**Supplementary Information** is linked to the online version of the paper at www.nature.com/nature.
**Acknowledgements** We thank Y. Zhang for technical support. We would like to thank members of the Crair lab for valuable comments on the manuscript. This work was supported by NIH grants P30 EY000785, R01 EY015788 to M.C.C. M.C.C. also thanks the family of William Ziegler III for their support.
**Acknowledgements** We thank Y. Zhang for technical support. We would like to thank members of the Crair lab for valuable comments on the manuscript. This work was supported by NIH Grants RR19895, RR029676-01 for the Yale University Biomedical High Performance Computing Center and NIH grants P30 EY000785, R01 EY015788 to M.C.C. M.C.C. also thanks the family of William Ziegler III for their support.
**Author Contributions** J.B.A. and M.C.C. designed the experiments. J.B.A. performed in vivo imaging experiments, wrote the image processing and data analysis code, and analyzed the recordings. H.Z. created the GCaMP3 and GCaMP6 mouse lines. J.B.A. and M.C.C. wrote the manuscript.
@@ -123,6 +119,8 @@ Anesthetized Rx-Cre:GCaMP3 or SNAP25-GCaMP6 mice between postnatal day 2 to 13 (
<!--Figure 1 metadata
* neonate_ms_fig.png
* binary masks: Screen_Shot_2013-03-29_at_12.06.25_PM_crop.png, ..._crop1.png, ..._crop2.png
* 133.6s timeColorMap projection: 120518_07-fr698-1365-20140710-222230.tif
* Segmentation mask of 133.6s timeColorMap projection: 120518_07-fr698-1365-20140715-230320-cadj.png
* parcellation map
* update 2013-10-03 11:07:29:
* 120518_07_2013-09-11-225029_d2rImageCoords20130930-144657.ai
@@ -132,17 +130,28 @@ Anesthetized Rx-Cre:GCaMP3 or SNAP25-GCaMP6 mice between postnatal day 2 to 13 (
* 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022_d2r_20130930-124950_eps.png
* [x] add a domain centroid size/duration map similar to: ![](../figures/Screen_Shot_2013-04-03_at_8.42.49_AM.png)
* ![](../figures/Screen_Shot_2013-04-03_at_10.04.36_AM.png)
* domainFreq map:
* domainFreq map: 20140613-084032_ActivityMapFigRawProj-domainFreqRdBu.eps, P5mapDomainFreq-crop.tif, P8mapDomainFreq-crop.tif, P13mapDomainFreq-crop.tif
* domainFreq boxplot:
* domainDur map: 20140613-082954_ActivityMapFigRawProj-domainDurRdBu.eps
* domainDur map: 20140613-082954_ActivityMapFigRawProj-domainDurRdBu.eps, P5mapDur-crop.tif, P8mapDur-crop.tif, P13mapDur-crop.tif
* domainDur cdf:
* domainDiam map: 20140613-083156_ActivityMapFigRawProj-domainDiamRdBu.eps, P5mapDiam-crop.tif, P8mapDiam-crop.tif, P13mapDiam-crop.tif
-->
<!-- Figure 2 metadata
* Scatterplot mean freq: 140704-143614-freq_min-scatter.pdf
* Scatterplot domain duration-diameter timeline: 140702-141423-durDiam-scatter-img.pdf
* P3 M1 wave time projection, lomag: 140218_13-fr2564-2601-20140716-144024.tif
* P5 V1 wave time projection, himag: 140328_10-fr2016-2153-20140710-092344.tif
* P13 global wave time projection, lomag: ~~140509_07-fr1725-2004-20140717-172658.tif~~
* P13 global wave time projection, lomag: 140509_22-fr484-609-20140725-155542.tif
* wave motion index scatter: 140702-152109-rhoDiam-scatter.png
* Wave motion index timeline: 140702-162029-waveMotionIdx-jitterTimeline.pdf
-->
<!--Figure 3 metadata
* Isoflurane traces:
* P8 isoflurane active px fraction and motor movement signal traces: 140331_11_20140425-171130_d2ractvFraction20140425-171354-20140425-171355.ai
* P8 isoflurane time color projection sequence: 140331_11-fr775-975-20140710-163609.tif, 140331_11-fr1800-2000-20140710-162448.tif, 140331_11-fr2800-3000-20140710-162206.tif
* Cortical and photodiode signals and moving averages color coded at diff lags: 120518_07_2013-09-11-225029_d2r_motorSignalFiltFilt_fig.eps
@@ -157,27 +166,70 @@ Anesthetized Rx-Cre:GCaMP3 or SNAP25-GCaMP6 mice between postnatal day 2 to 13 (
-->
<!--Figure 4 metadata
* binary mask snapshots, cropped from screen shots in [[2013-04-19_analysis]]
* Screen_Shot_2013-04-19_at_8.26.00_AM_fr1786.png
* Screen_Shot_2013-04-19_at_8.27.49_AM_fr2134.png
* Screen_Shot_2013-04-19_at_8.30.27_AM_fr759.png
* Screen_Shot_2013-04-19_at_8.30.51_AM_fr373.png
* Screen_Shot_2013-04-19_at_8.38.54_AM_fr177.png
* Temporal correlation of activity between the hemispheres and preceding motor activation:
* ![](../figures/Screen_Shot_2013-04-30_at_3.02.20_PM.png)
* hemisphere active fraction traces: Screen_Shot_2013-04-08_at_8.47.19_AM.png
* activefraction hemis AP and ML all: ![](../figures/Screen_Shot_2013-04-23_at_8.45.18_AM.png) | 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022_d2ractiveFractionPixelLocaCorr20130423-094506.eps
* activefraction hemis AP and ML segment: ![](../figures/Screen_Shot_2013-04-23_at_8.46.27_AM.png)
* activefraction hemis AP and ML segment: ![](../figures/Screen_Shot_2013-04-23_at_8.51.55_AM.png)
### Cortical activity correlated between the hemispheres and is periodic
* hemi auto and xcorr:
* 2500fr lags: 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022activeFraction20130408-143100.eps
* 250fr lags: 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022activeFraction20130408-151655.eps
* 1500fr lags: 120518_07_2013-10-18_AllgoodactiveFraction20131023-145023.eps
* corr matrix: MeanCorrMatrix-age-2014-06-02.ai, (140602-113013-age.g-groupCorrMatrix.pdf, 140602-112947-dendr-ageP12-13.pdf)
* 140808-093608-graph-squaredEdgeScale-auto.pdf
~~* corr graph spatial layout: 140602-100049-P12-13_0.15.pdf ~~
-->
<!--Figure 5 metadata
* corr matrix: MeanCorrMatrix-age-2014-06-02.ai, (140602-113013-age.g-groupCorrMatrix.pdf, 140602-112947-dendr-ageP12-13.pdf)
* corr graph force layout: 140602-101610-P12-13_0.15.pdf
* corr graph spatial layout: 140602-100049-P12-13_0.15.pdf
* 140807-121720-graph-squaredEdgeScale.pdf
~~* corr graph force layout: 140602-101610-P12-13_0.15.pdf~~
~~140806-174912-modularity-boxplot.pdf~~
~~140806-174917-diam-boxplot.pdf~~
140806-174920-clusterCoeffGlobal-boxplot.pdf
140806-174915-pathLength-boxplot.pdf
* 140804-111208-eigCentr-scatterplot.pdf, r>0.15
~~* 140804-075217-degree-hist.pdf, binwidth=2, r>0.15~~
-->
<!-- Supplementary Movies
* supplementaryMovie-P3gcamp3.mov: wholeBrain-shortAlpha-lomed.mov, (120518_07.tif) 10 s long playback, 30fps, = 300fr = 60 s real time
* supplementaryMovie-P8gcamp3.mov: 131208_06_std_lomed-all.mov, 6 s long playback, 30fps, 6s*30fps+3fr = 184fr = 36.8 s real time
* supplementaryMovie-P13gcamp6.mov: 140509_22_fr484-735-lo-all-trans.mov, 8 s long playback, 30fps, 8*30+11fr=. 50.2 s real time.
-->
<!-- Supplementary Figs
* supplementaryFig-domain-stats.ai:
* supplementaryFig-areal-stats.ai:
* 140702-204319-freq_min-boxplot.pdf
* 140703-083800-freq_min-boxplot.pdf
* 140606-091326-diamnodeT-boxplot.png
* 140606-091332-MeanIntensitynodeT-boxplot.png
* 140606-091338-MaxIntensitynodeT-boxplot.png
* 140606-091345-Duration_snodeT-boxplot.png
* 140606-143236-pairwisetPvalMatrix.pdf
* supplementaryFig-montage.ai:
* 140509_22-fr484-735-20140729-102651-montage.png, 140509_22-fr484-735-20140729-102651-montage.tif
* supplementaryFig-corticaltracesCorr.ai:
* P5: 140328_10_20140420-111218_d2rActvRaster20140420-111804.eps, 3000fr
* P8: 140331_05_20140425-160657_d2rActvRaster20140425-161244.eps, 3000fr
* P13: 140509_22_20140522-074327_d2rActvRaster20140522-075627.eps, 3000fr
* supplementaryFig-symmetry.ai:
* binary mask snapshots, cropped from screen shots in [[2013-04-19_analysis]]
* Screen_Shot_2013-04-19_at_8.26.00_AM_fr1786.png
* Screen_Shot_2013-04-19_at_8.27.49_AM_fr2134.png
* Screen_Shot_2013-04-19_at_8.30.27_AM_fr759.png
* Screen_Shot_2013-04-19_at_8.30.51_AM_fr373.png
* Screen_Shot_2013-04-19_at_8.38.54_AM_fr177.png
* Temporal correlation of activity between the hemispheres and preceding motor activation:
* ![](../figures/Screen_Shot_2013-04-30_at_3.02.20_PM.png)
* hemisphere active fraction traces: Screen_Shot_2013-04-08_at_8.47.19_AM.png
* activefraction hemis AP and ML all: ![](../figures/Screen_Shot_2013-04-23_at_8.45.18_AM.png) | 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022_d2ractiveFractionPixelLocaCorr20130423-094506.eps
* activefraction hemis AP and ML segment: ![](../figures/Screen_Shot_2013-04-23_at_8.46.27_AM.png)
* activefraction hemis AP and ML segment: ![](../figures/Screen_Shot_2013-04-23_at_8.51.55_AM.png)
### Cortical activity correlated between the hemispheres and is periodic
* hemi auto and xcorr:
* 2500fr lags: 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022activeFraction20130408-143100.eps
* 250fr lags: 120518_07_connComponents_BkgndSubtr-60px_noWatershed-20130327-151022activeFraction20130408-151655.eps
* 1500fr lags: 120518_07_2013-10-18_AllgoodactiveFraction20131023-145023.eps
-->