Abstract
The laminar structure and conserved cellular organization of mouse visual cortex provide a useful model to determine the mechanisms underlying the development of visual system function. However, the normal development of many receptive field properties has not yet been thoroughly quantified, particularly with respect to layer identity and in the absence of anesthesia. Here, we use multisite electrophysiological recording in the awake mouse across an extended period of development, starting at eye opening, to measure receptive field properties and behavioral-state modulation of responsiveness. We find selective responses for orientation, direction, and spatial frequency at eye opening, which are similar across cortical layers. After this initial similarity, we observe layer-specific maturation of orientation selectivity, direction selectivity, and linearity over the following week. Developmental increases in selectivity are most robust and similar between layers 2–4, whereas layers 5 and 6 undergo distinct refinement patterns. Finally, we studied layer-specific behavioral-state modulation of cortical activity and observed a striking reorganization in the effects of running on response gain. During week 1 after eye opening, running increases responsiveness in layers 4 and 5, whereas in adulthood, the effects of running are most pronounced in layer 2/3. Together, these data demonstrate that response selectivity is present in all layers of the primary visual cortex (V1) at eye opening in the awake mouse and identify the features of basic V1 function that are further shaped over this early developmental period in a layer-specific manner.
Introduction
Studies of mouse primary visual cortex can serve as a powerful tool to causally link developmental mechanisms to normal cellular and circuit function in the adult visual system. The majority of this work has focused primarily on understanding the mechanisms underlying retinotopy and ocular dominance (OD) plasticity (for review, see Huberman et al., 2008; Morishita and Hensch, 2008; Cang and Feldheim, 2013; Triplett, 2014). However, the visual system must also appropriately wire up circuits to support the emergence of many basic receptive field (RF) properties, including orientation selectivity (OS), direction selectivity (DS), spatial frequency (SF) preference, and linearity. Thus far, only a handful of studies have quantified the early maturation of a subset of these properties in mouse V1 (Fagiolini et al., 2003; Rochefort et al., 2011; Wang et al., 2010, 2013; Ko et al., 2013, 2014). No single study has systematically quantified the simultaneous development of these RF properties in the awake mouse as a function of cortical layer in monocular V1. Such a study would serve as a foundation for understanding the molecular and cellular developmental mechanisms that give rise to these fundamental features of visual system function.
The visual cortex of mice, as well as other mammals, comprises multiple layers, each of which is defined by a distinct set of functional properties, subcortical inputs, and intracortical connectivity (Ringach et al., 2002; Niell and Stryker, 2008; Clemens et al., 2012). Furthermore, selective responses, such as DS and OS, emerge in multiple layers, and there is evidence in mammals that this selectivity arises within each layer via distinct developmental mechanisms (Trachtenberg et al., 2000; Trachtenberg and Stryker, 2001; Clemens et al., 2012). It remains unknown whether distinct mechanisms contribute to the development of basic RF properties and V1 function across cortical layers in the mouse. Careful examination of the development of the spatial–temporal structure of the RF and selectivity of neurons throughout V1 as a function of layer may provide the information necessary to begin to address this outstanding question.
Here, we systematically quantify RF structure and a core set of V1 response properties as a function of layer in awake, quiescent mice over development between eye opening and adulthood. This work is further distinguished from existing studies in that we additionally measure the emergence of mature patterns of behavioral-state modulation of response gain and local field potential (LFP) power (Niell and Stryker, 2010). Overall, we find that all measured RF properties are present in each cortical layer of V1 at eye opening. We also observe that the spatial structure of the RF of linear responses in V1 at eye opening closely resembles those measured in the adult. However, significant refinement of selectivity is seen specifically in layers 2–4 in the week after eye opening, whereas the majority of response selectivity observed in layers 5 and 6 does not change significantly after eye opening. Interestingly, although behavioral state affects response gain already at eye opening, the laminar profile of this effect switches over development.
In general, we find evidence for layer-specific developmental differences in both RF properties and behavioral-state modulation of V1 responsiveness. DS, RF shape, response gain, and LFP modulation significantly develop in a layer-specific manner over early development from eye opening. These findings argue that, despite a substantial degree of cortical function already established at eye opening, there are important layer-specific differences in the development of monocular V1 function from eye opening through adulthood.
Materials and Methods
All studies were conducted with approved protocols from the University of Oregon Institutional Animal Care and Use Committees, in compliance with National Institutes of Health guidelines for the care and use of experimental animals.
In vivo physiology.
Recordings were performed on awake female and male C57BL/6J mice, ages P12 up to 2.5 months of age. Animals were maintained on a 14/10 h light/dark cycle. All recordings were performed 3–4 h before the start of their normal dark (awake) cycle. This facilitated comparable epochs of moving and stationary behaviors at all ages.
We observed that eye opening could occur between P12 and P14, and therefore pups in this range were monitored every 12 h to determine the precise time of eye opening. Pups within 0–24 h of eye opening were pooled in the EO1 (day 1 of eye opening) developmental group and similarly for EO3 (3–4 d after eye opening) and EO7 (7–8 d after eye opening).
Before recording, all animals underwent one 10 min handling session before head-plate implant and two 10 min sessions after implant. This handling was followed up with one 10 min “practice” session in which the animal was attached to the spherical treadmill recording apparatus described previously (Niell and Stryker, 2010) and initially developed by Dombeck et al. (2007) (Fig. 1A).
To facilitate head fixation during recording, custom stainless steel head-plate implants were cemented to the skull between 2 and 7 d before the recording session as described previously (Niell and Stryker, 2010). Briefly, adults and pups were anesthetized with isoflurane in oxygen (3% induction, 1.0–1.5% maintenance), warmed with a heating pad at 37°C (adult) to 39°C (pups), and given subcutaneous injections of 0.1 ml (pups) or 0.25 ml (adults) of lactated Ringer's solution and carprofen (5 mg/kg). In all animals, the scalp and fascia from bregma to just behind lambda were removed, and the skull was covered with a thin layer of cyanoacrylate (VetBond; WPI) before attaching the head plate with dental acrylic (Ortho-JET; Lang Dental Manufacturing). The well of the head plate was filled with silicone elastomer (Kwik-Sil; WPI) to protect the skull before recordings. Head plates were stable on pups for up to 3 d, after which skull growth and their mother's grooming behavior decreased stability of the implant.
On the day of recording, the animals were anesthetized to perform a craniotomy over the visual cortex. The craniotomies were ∼1 mm in diameter and centered at 2.5 mm lateral of midline and 1 mm anterior of the posterior suture. The brain surface was covered in 1.5% agarose (Sigma-Aldrich) in sterile saline and then capped with silicone elastomer during a 3–4 h recovery. To begin the recording session, the animal was placed in the head-plate holder on the free-floating ball and allowed to habituate for 15 min in low light. After removing the protective agarose and silicone plug, the ground wire was set, and a fresh layer of 1.5% agarose in saline was applied to the well. A multisite electrode (A2X32-5mm-25-200-177-A64; Neuronexus Technologies), coated with DiO (Invitrogen) to allow post hoc track recovery, was then inserted through the craniotomy and overlying agarose. We used a custom 64-channel electrode arranged in a two-shank, linear configuration, with each shank having 32 sites separated by 25 μm. The electrode was placed at an angle of 45° relative to the cortical surface and inserted to an appropriate depth. After insertion, agarose was added to stabilize the electrode, and it was allowed to settle for 40 min. Targeting the penetrations in this way ensured sampling of all cortical layers across groups. Only one penetration per animal was made, but up to two recording sessions per animal were performed in approximately two-thirds of the animals by advancing the electrode until all sites had moved past previous recording depths. All units stably isolated over the recording were included in subsequent analysis, with specific response properties calculated offline.
Recordings sessions typically lasted 1.5 h. Stimulus presentation timing was optimized previously to ensure that each stimulus was presented five to seven times in each behavioral state (Niell and Stryker, 2010). After the recording session, mice were deeply anesthetized in 3% isoflurane and killed via cervical dislocation. The brains were then fixed whole in 4% PFA (Electron Microscopy Sciences) overnight at room temperature. Brains were subsequently sectioned and mounted in Flouromount G with DAPI (Southern Biotechnology). We imaged the DiO electrode tracks on a Zeiss Axio Imager 2 to confirm that penetrations were specific to monocular V1 (Antonini et al., 1999; Altamura et al., 2007; Niell and Stryker, 2008) and to determine the location of the electrode tips relative to the cortical layer.
Visual stimuli.
Visual stimuli were presented as described previously (Niell and Stryker, 2008). Briefly, stimuli were generated in MATLAB (MathWorks) using the Psychophysics Toolbox extension (Brainard, 1997; Pelli, 1997) and displayed with gamma correction on a monitor (Planar SA2311W, 30 × 50 cm, 60 Hz refresh rate). The screen was placed 25 cm from the mouse's eye, subtending ∼60–75° visual space. Monitor mean luminance was measured to be 70 cd/m2. The monitor was raised or lowered up to 10 cm, depending on the hand-mapped RF location of multiunit activity.
For LFP measurements, we presented full-length drifting bars (5° width, 30°/s velocity, eight directions, two contrasts). To characterize selectivity of neural responses with single-unit recordings, we presented drifting sinusoidal gratings of 1.5 s duration at 100% contrast, with temporal frequency of 2 Hz, SF of 0.01, 0.02, 0.04, 0.08, 0.16, 0.32, and 0 cycles/° (cpd) in 12 evenly spaced directions. Gratings with a temporal frequency of 2 Hz were chosen to compare responses with those from the study by Niell and Stryker (2008), which showed a preferred temporal frequency of 2 Hz in adults. However, we also confirmed that there were no differences in the OS index (OSI), DS index (DSI), or SF preference when temporal frequency was varied between 1, 2, and 4 Hz in layers 2–5 of juveniles (OSI: 0.65 ± 0.08, 0.58 ± 0.07, 0.66 ± 0.07, p = 0.830, 1 vs 2 Hz, p = 0.674, 2 vs 4 Hz; DSI: 0.13 ± 0.05, 0.18 ± 0.05, 0.18 ± 0.06, p = 0.684, 1 vs 2 Hz, p = 0.865, 2 vs 4 Hz; and SF preference: 0.05 ± 0.02, 0.04 ± 0.01, 0.03 ± 0.01, p = 0.923, 1 vs 2 Hz, p = 0.657, 2 vs 4 Hz, n = 30, Wilcoxon's rank-sum test, Bonferroni's correction). The stimulus conditions were randomly interleaved, and a gray blank condition (mean luminance) was included to estimate the spontaneous firing rate. To identify cortical layers via current source density (CSD) analysis, as in the study by Niell and Stryker (2008), we presented a contrast-reversing square checkerboard (0.04 cpd, square-wave reversing at 0.5 Hz). The spatial RF was estimated by the spike-triggered average (STA) in response to Gaussian 1/f noise movies as described previously.
Data acquisition and analysis.
Signals were acquired using a System 3 workstation (Tucker Davis Technologies) and analyzed with custom routines written in MATLAB.
As in the study by Niell and Stryker (2010), movement signals from the optical mice were acquired at up to 300 Hz and integrated at 100 ms intervals (Mx310; Logitech). We used these measurements to calculate the net physical displacement of the top surface of the ball and calculate the average speed of the animal during a stimulus presentation. Stationary trials were classified as ball speeds <0.7 cm/s and moving >0.7 cm/s for all age groups. The 0.7 cm/s separates the central peak at ∼0 cm/s (stationary) from the wide distribution of moving speeds for all developmental groups tested, as well as adults. This threshold is slightly lower than that used by Niell and Stryker (2010) (1 cm/s) to account for overall lower speeds in young mice.
For LFP analysis, the extracellular signal was filtered from 1 to 300 Hz and sampled at 1.5 kHz. The power spectrum was computed using multi-taper estimation in MATLAB with the Chronux package (Mitra and Pesaran, 1999; Mitra and Bokil, 2007), using a 3 s sliding window and three to five tapers. Spectra were normalized for presentation by applying a 1/f correction (Sirota et al., 2008).
For single-unit activity, the extracellular signal was filtered from 0.7 to 7 kHz and sampled at 25 kHz. Spiking events were detected online by voltage threshold crossing, and a 1 ms waveform sample on four adjacent recording surfaces constituting a “virtual” tetrode was acquired around the time of threshold crossing. Single-unit clustering and spike waveform analysis was performed as described previously (Niell and Stryker, 2008) with a combination of custom software in MATLAB and Klusta-Kwik (Harris et al., 2000). Quality of separation was determined based on the following criteria: (1) contained <0.1% of spikes within a 1.0 ms interspike interval; (2) at least one feature of a given spike had to form a visually identifiable well segregated cluster in a projection onto a two-dimension feature subspace, confirmed by the computed L ratio (Schmitzer-Torbert et al., 2005); (3) the peak amplitude of a unit remained stable over the entire recording session; and (4) the peak amplitude of a unit had a signal-to-noise ratio (SNR) of >2 (mean amplitude at the peak is two times greater than the SD of the amplitude). In fact, the median SNR for adult units and EO juveniles was 5.95 ± 0.09 and 5.29 ± 0.23, respectively. Representative isolated units along with their SNR are shown in Figure 3A. Typical recordings yielded ∼10–12 single units across each shank of the electrode. Most units appeared predominantly on a single recording site, although units often contributed a signal below the voltage threshold on neighboring sites, which allowed improved unit discrimination.
Units were classified as narrow or broad spiking based on properties of their average waveforms at the electrode site with largest amplitude. As detailed in the study by Niell and Stryker (2008), two parameters—(1) height of the positive peak relative to the initial negative trough and (2) time from the minimum of the initial trough to maximum of the following peak—were sufficient to generate two linearly separable clusters corresponding to narrow-spiking (putative inhibitory) and broad-spiking (putative excitatory) neurons. These clusters were separated using K means.
Layer identity was derived independently and corroborated from both CSD analysis (Figure 2B) and recovery of electrode tracks labeled by DiO in histological sections (Fig. 2A). In both cases, all sites were given a specific layer assignment based on their position along the probe relative to the probe tip and the geometry of the penetration angle.
The average spontaneous rate for each unit was calculated by averaging the rate over all blank condition presentations. Responses at each orientation and SF were calculated by averaging the spike rate during the 1.5 s presentation of stimulus type and subtracting the spontaneous rate during gray screen presentations. We first determined the neurons that were visually responsive by calculating the t statistic (mean evoked firing divided by SE) across trials for the optimal stimulus. Units with t > 2 were considered reliably visually responsive. To ensure robust assessment of RF and tuning properties, we imposed an additional more stringent requirement for those units entered into all subsequent analyses: that units have an average peak evoked firing rate, across trials of the optimal stimulus, of >2 Hz.
Tuning curves were analyzed following the methods of Niell and Stryker (2008). Briefly, the preferred orientation was determined by averaging the response across all SFs and calculating half the complex phase of the value. The orientation tuning curve was constructed for the SF that gave peak response at this orientation. Given this fixed preferred orientation, θpref, the tuning curve was fitted as the sum of two Gaussians centered on θpref and θpref + π, of different amplitudes A1 and A2 but equal width σ, with a constant baseline B. From this fit, we calculated two metrics: (1) an OSI representing the ratio of the tuned versus untuned component of the response; and (2) the width of the tuned component. The OSI was calculated as the depth of modulation from the preferred orientation to its orthogonal orientation θortho = θpref + π/2, as (Rpref − Rortho)/(Rpref + Rortho). Tuning width was the half-width at half-maximum of the fit above the baseline, Rortho. We also included the global measures of orientation and direction selectivity as represented by 1 − Circular variance (1 − CirVar) and 1 − direction space circular variance (DirCirVar). These values were computed as in the study by Mazurek et al. (2014), after Ringach et al. (2002) and Batschelet (1981):
Statistics.
Comparisons of medians and variation around the median are performed via Wilcoxon's rank-sum test unless noted otherwise. SEM is obtained via bootstrapping (data sampled 1000×). Differences in frequency of categorical observations between groups are determined via χ2, and differences in distributions of continuous, quantitative data between groups are determined via Kolmogrov–Smirnov (KS) test. To test for the modality of a distribution, we used the Hartigan's dip test that probes for deviations from unimodality (Hartigan and Hartigan, 1985). The Lillifor test is used to determine distribution skew. All data represent data from all animals in each group: EO1, N = 13, n = 227; EO3, N = 11, n = 295; EO7, N = 11, n = 225; and adult, N = 10, n = 283. When necessary, Holm–Bonferroni test is applied to correct for multiple comparisons of data. For the comparisons of ratios, SEMs are computed via a bootstrapping procedure, and significance was determined via Wilcoxon's signed-rank tests.
Results
Our goal was to quantify the early development of neuronal response properties and behavioral-state modulation as a function of layer in monocular V1 of the awake mouse. To accomplish this goal, we measured V1 function over several stages of development from the onset of eye opening to adulthood in awake mice (Fig. 1A). We used a custom multisite silicon electrode to record neuronal responses while delivering stimuli designed to measure canonical cortical responses (Fig. 1B–E; Niell and Stryker, 2008). Recordings were obtained between the first day of eye opening (EO1) and ∼2 months of age (adult). Initial experiments specifically compared the EO1 group (n = 227 units in 13 animals) with the adult group (n = 283 units in 11 animals) to first understand which visual responses significantly changed over development and whether we could detect layer-specific developmental differences. In follow-up experiments, we compared these measures with two intermediate developmental time points, EO3 (N = 11) and EO7 (N = 10) to determine the time course over which significant changes occurred.
Laminar location, cell-type identity, and responsiveness
To identify the laminar location of recording sites within V1, we recovered electrode tracks through histology and reconstructed site localization based on the geometry of the probe relative to the penetration angle at the surface of cortex (Fig. 2A; see Materials and Methods). Localization to monocular V1 was confirmed simultaneously. CSD analysis (Fig. 2B) and qualitative observations of previously identified layer-specific activity patterns in the adult confirmed the site localization identified via histology (Niell and Stryker, 2008).
Well isolated single units (Fig. 3A; see Materials and Methods) were then classified based on their spike waveform parameters. As shown previously (Niell and Stryker, 2008), clustering reliably separated two distinct classes of waveforms: narrow-spiking and broad-spiking waveforms (Fig. 3B, red and black, respectively). However, we observed no narrow-spike waveforms among units recorded at EO1 (Fig. 3C). In our adult group, ∼11% (30 of 283 units) of visually responsive units had classical, narrow waveforms. The narrow waveforms did not make up a significant percentage of the responsive population compared with the adult group until 7 d after eye opening (EO3 fraction of population with narrow waveform: 1%, n = 130, χ = 8.72, p = 0.0031; EO7 fraction with narrow waveform: 12%, n = 121, χ = 0.0821, p = 0.774).
To characterize unit responsiveness over development, we calculated an SNR for the firing rate response to repeated presentations of the preferred stimulus (Fig. 4E; see Materials and Methods). The fraction of cells that were deemed responsive by this statistical measure increases over development in layers 4–6 and slightly decreases in layers 2/3 between EO1 and adulthood. Although this measure determines which units have a statistically significant visual response, responses meeting this criterion could still be weak and not suitable for calculating tuning curves and RFs. Therefore, we used a second criterion—a firing rate of at least 2 Hz to the optimal stimulus—to select units for subsequent analysis. The distribution of cells meeting this criterion is shown in Figure 4F.
Visually evoked firing rates in layers 2/3, 5, and 6 were higher in adults than at eye opening (Fig. 4A,B). Additionally, the spontaneous firing rate specifically in layers 2/3 decreased over development (Fig. 4C,D). Thus, overall visual responsiveness of cortex, both in terms of the number of responsive neurons and the peak firing rate, increases after eye opening in a layer-specific manner.
RF structure
Previous studies have shown that simple cells in V1 have a characteristic spatial RF structure that is well described by a subset of two-dimensional Gabor functions. Notably, this RF organization is shared across mice, cats, and monkeys (Jones and Palmer, 1987; Ringach, 2002; Niell and Stryker, 2008). To determine when this structure is established over development, we obtained the STA of the response of each unit to stochastic noise (Fig. 1E; see Materials and Methods) and fit these STAs to a two-dimensional Gabor function (Ringach, 2002; Niell and Stryker, 2008; Ko et al., 2013). STAs in which >60% of the signal variance was accounted for by the Gabor fit were included for analysis. We compared STA structure using the scale-free parameters nx and ny, which express the width and length of the fitted Gabor function in terms of the number of cycles of the underlying sinusoidal grating (Ringach, 2002).
RF structure, as measured by STA, was established at eye opening in all cortical layers. We found no significant differences in the fraction of units with STAs well described by Gabor functions over development (EO1: 49 ± 8%, n = 164; adult: 50 ± 7%, n = 191; p = 0.0939, χ2 test). In contrast, the median nx and ny values for the EO1 group were significantly smaller than the adults (Fig. 5A), indicating decreased complexity of spatial RF structure at eye opening.
When measures of nx and ny were segregated by layer, significant differences were only observed in layer 5 after correction for multiple comparisons (Fig. 5E,F). Overall median nx and ny values reached adult-like levels by EO3 (Fig. 5G). Although we found these differences in RF organization over development, the median size of the RFs overall, as measured by the Gaussian envelope of the Gabor function, was not significantly different between the EO1 and adult groups (data not shown). Thus, although total RF size did not change over development, the spatial structure within this area did refine, with a slight shift toward more subregions and a significant elongation of the RF.
Layer-specific maturation of orientation and direction selectivity
We next characterized the selectivity of V1 neurons for specific visual properties throughout the time course of development. At eye opening, drifting sinusoidal gratings elicited orientation- and direction-selective responses (Fig. 6A), consistent with previous findings that these fundamental properties are established before visual experience (Fagiolini et al., 2003; Rochefort et al., 2011; Ko et al., 2013). However, there are important differences in OS between EO1 animals and adults (Fig. 6, compare A, B). Notably, OS significantly increases over development as indicated by developmental changes in the median OSI and median orientation tuning width of responsive units (Fig. 6C–F; Table 1). We find the same developmental increase (approximately threefold) in OS by comparing a global measure of OS, 1 − CirVar (Mazurek et al., 2014; Fig. 7A; Table 1).
Consistent with previous findings in the anesthetized adult mouse (Niell and Stryker, 2008), we find that OS is higher in layers 2–4 relative to layer 5 and inhibitory units (Fig. 6C,D; Table 1). Interestingly, this difference in OS between excitatory cells from layers 2–6 is not yet manifested at EO1 (Fig. 6C,D; Table 1). Another notable aspect of the development of OS was that units at EO1 exhibited a bias for orientations along the cardinal axes (Fig. 6G; Table 1). This bias, which was observed across layers, was significantly reduced by adulthood (Fig. 6H; Table 1).
In examining intermediate time points, we found that OSI and OS tuning width change most rapidly over week 1 of eye opening, peaking at adult-like levels by EO7 for excitatory responses (Fig. 7A,B; Table 1). The OSI for inhibitory responses also did not change significantly between EO7 and adulthood (EO7, 0.38 ± 0.09; adult, 0.41 ± 0.06; n.s., Wilcoxon's rank-sum test; Table 1). In contrast, the bias toward cardinal axes persisted throughout week 1 after eye opening, suggesting relatively late maturation of this feature of selectivity (Fig. 7C). Overall, this suggests that OS is primarily refined in week 1 after eye opening in layers 2–5 but that changes must continue to balance the representation of all orientations.
We next characterized DS, the preference for one particular direction of motion for a given oriented grating. Although recent work has shown a lack of significant developmental changes in DS in superficial layers of V1 in anesthetized mice (Rochefort et al., 2011), this has not yet been examined in awake mice or in deeper layers. Our experiments reveal a significant change in DS in layers 2/3 that occur specifically between EO1 and EO7 (Fig. 8C; Table 1). Unexpectedly, we find that DS significantly decreases between EO7 and adulthood in layers 2/3. In comparison, we find no significant differences in DS over any stage of development in layer 4 (Fig. 8D). These specific observations of DS also hold true when measuring 1 − DirCirVar (Fig. 8C,D; Table 1). Together, our data argue that there are important differences in the development of DS in the mouse visual cortex. However, future work will require both layer-specific examination of this development and the inclusion of data from 1 week after eye opening.
SF preference and response linearity significantly mature after eye opening
One hallmark of selectivity in mammalian V1 is preference for particular SFs (SF preference). Furthermore, the range of SFs over which neurons respond sets constraints on cortical visual acuity. Similar to OS and DS, SF preference across a range of SFs is observed in each layer at eye opening (Fig. 9). However, the distribution of SF preferences in EO1 juveniles was skewed toward lower SFs and was significantly different from the distribution of SF preference in the adult (EO1 vs adult; n = 164 of 220; p = 0.0017, KS test; Fig. 9A). Despite this difference in the distribution, the median SF preference was not significantly different between EO1 juveniles and adults in any layer after correction for multiple comparisons (Fig. 9B; Table 1).
To probe the differences more closely, we examined the populations responding to either full-field flicker or high SFs. We observe a decrease in the population responding to full-field stimuli in layers 2–4 between EO1 juveniles and adults (Fig. 9C). This change is accompanied by the emergence of a significant population of responses to high SF gratings (>0.25 cpd) in all layers (Fig. 9D). These layer-specific shifts in the population response underlie the overall shift in SF preferences seen in V1 development. Interestingly, we also find that this developmental progression results in layer-specific differences in responses to low and high SFs in the adult. In layers 2–4, there are significantly fewer responses to full-field flashes relative to layers 5 and 6, whereas responses to high-frequency gratings are also significantly higher in layers 2/3 relative to all other layers in the adult. These developmental changes in SF preference occurred steadily between eye opening and adulthood (Fig. 9E,F; Table 1). Concurrent with changes in the preferred SF over development, SF tuning width significantly narrowed within layers 2/3 (Fig. 10A–C).
V1 neurons were characterized originally as simple or complex based on their ability to be defined by segregated ON and OFF subregions (Hubel and Wiesel, 1962, 1968). This aspect of RF organization, reflecting the linearity of spatial summation, has been shown to vary across the layers of cortex and is often measured by the periodicity of response to drifting gratings (F1/F0 ratio; Skottun et al., 1991; Mechler and Ringach, 2002). Therefore, we measured developmental changes in linearity of responses to drifting gratings over development. As predicted by our ability to measure linear STAs, we find linear responses to gratings (F1/F0 > 1) already established at the time of eye opening (Fig. 11). Furthermore, response linearity is already significantly enhanced in layer 2/3 relative to layers 4 and 5 at EO1 (39.58 ± 3.17 vs 22.22 ± 3.48%, layers 2/3 and 4, respectively; p < 0.05, χ2 test; Fig. 11A). In adults, we observe a larger proportion of linear responses in layers 2–4 relative to layer 5 (47.22 ± 4.08 vs 24.07 ± 1.51%, layers 4 and 5, respectively; p < 0.01, χ2 test; Fig. 11A), consistent with previous studies of anesthetized mouse V1 (Niell and Stryker, 2008). Given this layer difference, we specifically compared developmental changes in layers 2–4 between age groups. We found that the overall proportion of linear responses increases between EO1 and the adult in layers 2–4 (p < 0.05, KS test; Fig. 11B; Table 1). Furthermore, linearity increases over each stage of development examined here (Fig. 11C; Table 1). Thus, together with the STA data, we find that RFs do have distinct ON and OFF regions at eye opening but that there is an increase in the spatial segregation of these regions over an extended period of development after eye opening.
Modulation of visual responses by behavioral state refines over week 1 after eye opening
Previous studies of V1 function in the awake mouse have shown that cortical activity is modulated as a function of behavioral state (Niell and Stryker, 2010). In particular, locomotion induces an increase in visually evoked neuronal firing rate and enhances the amplitude of a high-frequency oscillation in the gamma range observed in the LFP. Because these features of cortical activity may reflect the output of important computations that affect visual processing (Heeger, 1992; Vidyasagar et al., 1996; Sompolinsky and Shapley, 1997; Cardin et al., 2009; Katzner et al., 2011), we monitored their development from eye opening. Importantly, recent studies have found that mouse pups exhibit spontaneous bouts of quadruped locomotion by P10 (Clarke and Still, 2001; Dehorter et al., 2011). However, several aspects of the gait and velocity during locomotion do not fully mature until P24 (Clarke and Still, 2001). Consistent with these findings, we find that P12–P14 mouse pups exhibit clear shifts in locomotive state on the trackball during recordings (Figs. 12A, 13A, bottom gray trace). Although we observed a similar threshold for locomotion in EO1 juveniles compared with adults (Fig. 12A), their speeds during locomotive bouts were significantly slower than adults (EO1, 4.74 ± 0.01 cm/s vs adults, 7.57 ± 0.02 cm/s; p < 0.001, KS test).
A previous study primarily investigated locomotor induced modulation of firing rates in layer 2/3 (Niell and Stryker, 2010). Here we determined the time line for the emergence of these activities across all cortical layers in V1. In the adults, we find that firing rate is significantly modulated in layers 2–5. However, the strength of modulation is strongest for those units in layer 2/3 (ratios, 1.98 ± 0.21 vs 1.31 ± 0.19, layers 2/3 and 4 in the adult, respectively; Wilcoxon's rank-sum test, p < 0.05; Fig. 12B). Surprisingly, we find dramatic layer-specific differences in response gain between eye opening and adulthood (Fig. 12B). The increase in visually evoked firing rate during locomotion is absent in layer 2/3 at eye opening (1.12 ± 0.15 vs 1; n.s., Wilcoxon's signed-rank test; Fig. 12B; Table 1). Instead, locomotion induces increases in response gain specifically in layers 4 and 5 (1.91 ± 0.24 and 1.98 ± 0.23 vs 1, layers 4 and 5, respectively; Wilcoxon's signed-rank test, p < 0.05; Fig. 12B; Table 1). The juvenile laminar pattern of behavioral-state modulation of gain reverses between EO7 and adulthood (Fig. 12C; Table 1), resulting in the mature pattern of large gain changes seen in layer 2/3.
We also observe significant maturation of locomotor-induced changes in LFP power between eye opening and the adult (Fig. 13). In adults, locomotion induces a decrease in power at frequencies between 7 and 10 Hz and an increase in a narrow band of LFP power between 50 and 60 Hz (gamma), consistent with previous reports (Fig. 13C,D; Niell and Stryker, 2010). Several features of this characteristic activity are not observed at EO1 during locomotive bouts (Fig. 13, compare A, C). Instead, the most prominent feature induced by locomotion in juveniles is a broad increase in power over a lower-frequency range of 20–40 Hz. This suggests that the characteristic modulations of LFP activity are immature at eye opening.
Given our findings for layer-specific differences in firing rate modulation over development, we further analyzed LFP activity as a function of layer (Fig. 14). In the adult, the locomotion-induced increase in power (at 50–60 Hz) was observed in layers 2–5 but sharply attenuated in layer 6 (Fig. 14B). In contrast, at eye opening, the locomotion-induced increase (at 20–40 Hz) is evident only in layers 4 and 5 (Fig. 14A), consistent with the site of firing rate changes. Therefore, there are distinct patterns of locomotor-induced changes in the LFP power over development that vary as a function of layer. When we analyzed the peak frequency within the range increased by locomotion (30–60 Hz), we did not observe the fully mature pattern of frequency modulation until after EO7 (Fig. 14C). Similarly, the modulation of the amplitude of the peak signal between 20 and 60 Hz matured significantly between EO7 and adulthood (Fig. 14D). Thus, although cortical state, as assayed by firing rate changes and oscillatory activity, is already regulated by behavior at eye opening, the pattern of activity changes over development in a layer-specific manner.
Discussion
A timeline for the development of RF properties across layers of mouse V1
We quantified the normal development of OS, DS, SF preference, and linearity as a function of cortical layer in the quiescent, awake mouse. We find that a significant degree of response selectivity is established at eye opening across all layers of monocular V1. However, several measures of selectivity significantly increase over week 1 after eye opening, and this development is particularly prominent in layers 2–4. We also provide novel evidence that the development of layers 5 and 6 diverges from layers 2–4. This suggests that the maturation of basic response selectivity in layers 5 and 6 may be mostly complete at eye opening. Interestingly, features such as a developmental reduction in bias for selectivity along the cardinal axes and enhanced responses to high SF gratings (>0.25 cpd) are observed in nearly every layer and continue to mature after week 1 of eye opening. We observe that these early developmental changes are tightly coupled between cortical layers that share common features of selectivity in adulthood.
Our finding for an early increase in OS (between P14 and P22) is consistent with previous studies examining superficial, monocular V1 (Rochefort et al., 2011; Ko et al., 2013). Interestingly, this increase in OS in monocular V1 precedes that seen in binocular cortex because one study found that the OSI increases for contralateral responses significantly between P21 (∼1 week after eye opening) and adulthood (Wang et al., 2013). The same study also reports significant development of binocular matching between simple cells that take places between P21 and P23. Here, we find a significant increase in linearity in monocular V1 over week 1 after eye opening (P12–P21) that directly precedes orientation matching of binocular simple cells. Thus, our study provides an important link in relating the developmental trajectories of V1 function in the monocular and binocular zones. In addition, our finding for a developmental change in DS conflicts with previous observations (Rochefort et al., 2011). Differences in developmental time points examined, alertness, or variability by layer may account for this particular discrepancy, highlighting the importance in generating a more complete characterization of mouse visual cortex development.
Overall, our specific findings compliment other recent work on OS and the development of V1 function in the mouse and suggest that week 1 after eye opening (P14–P21) is particularly important in the initial refinement of many basic RF properties across several layers of V1 (Fagiolini et al., 2003; Wang et al., 2010, 2013; Rochefort et al., 2011; Ko et al., 2014). In addition, we also find that some developmental changes occur in monocular V1 over a time period that coincides with the onset of the classical critical period for OD plasticity and orientation preference matching (Hensch et al., 1998; Antonini et al., 1999; Wang et al., 2010). Importantly, this overall developmental trajectory mirrors classical studies in other models of visual system development. In cats, monkeys, and ferrets, selectivity for orientation is established at eye opening but undergoes significant refinement over a defined period thereafter (Wiesel and Hubel, 1963; Blakemore and Van Sluyters, 1975; Frégnac and Imbert, 1978; Imbert et al., 1978; Chapman and Stryker, 1993; Krug et al., 2001; White et al., 2001). However, the mechanisms that underlie this initial establishment and refinement of RF selectivity remain unclear and may vary by species.
The possible role of selective subcortical input in shaping cortical response selectivity
Recent evidence has suggested a dedicated circuit that links direction- and orientation-selective retinal ganglion cell responses to superficial layers of V1 (Cruz-Martín et al., 2014). It has been proposed that this subcortical selectivity may contribute substantially to the selectivity that is observed in mature V1 (Tan et al., 2011; Clemens et al., 2012; Cruz-Martín et al., 2014). However, there are distinct hallmarks of this subcortical selectivity that differentiate it from the type of selectivity that dominates in V1. In particular, OS/DS cells in the LGN, which are a minority of the population, are generally nonlinear (ON/OFF) with a concentric RF (Wang et al., 2010; Piscopo et al., 2013; Zhao et al., 2013a,b; Sarnaik et al., 2014). In contrast, orientation-selective cells in superficial mouse V1 are generally linear with elongated ON and OFF subregions aligned along their preferred orientation (Niell and Stryker, 2008; Bonin et al., 2011). Notably, the selective responses that we observe at eye opening have these canonical oriented RFs (Fig. 5). This suggests that the majority of neurons in V1 are not simply recapitulating subcortical selectivity in either the adult or the visually naive juvenile.
However, we do find an interesting parallel between the features of selectivity represented in the retina and the cortex at eye opening. The majority of OS and DS neurons present at eye opening in the cortex are tuned along the cardinal axes. This feature of DS tuning in the early developing mouse cortex was reported previously by Rochefort et al. (2011). Interestingly, this bias in tuning is a property of DS responses in the retina and LGN (Marshel et al., 2012; Yonehara et al., 2013) and is established before eye opening in the retina (Elstrott et al., 2008; Chen et al., 2009). Thus, although cortical neurons may not directly inherit this tuning from subcortical projections, it raises the possibility that these DS/OS inputs may provide a bias in the initial establishment of RFs.
Maturation of behavioral-state modulation of the cortical circuit
Locomotion has been shown to exert a powerful effect on cortical processing in the mouse (Niell and Stryker, 2010; Andermann et al., 2011; Keller et al., 2012; Ayaz et al., 2013). Recent studies have suggested that cholinergic input to a disinhibitory circuit mediated by vasointestinal peptide-expressing interneurons may at least partially contribute to this activity (Fu et al., 2014; Lee et al., 2014). However, other modulatory pathways are likely required as well (Polack et al., 2013). Surprisingly, we find that the laminar profile of this effect changes between week 1 after eye opening and adulthood. Early in development, locomotion primarily modulates the firing rate of layers 4 and 5, whereas in adulthood, gain effects are predominantly in layer 2/3. This shift could be attributable to refinement of the modulatory inputs to cortex, changes in the intracortical synaptic connectivity mediating the effects, or both. Interestingly, this behavioral-state effect has been shown to regulate plasticity in the adult (Kaneko and Stryker, 2014). Therefore, it is possible that this could represent a shift in the sites of plasticity within the cortical circuit over development.
At the same time, immature patterns of behavioral modulation of coherent LFP are apparent across layers at EO1 and also do not reach maturity until after EO7 (Figs. 13, 14). Overall, we observe gamma-like modulation responses specifically in layers 4 and 5 at eye opening, whereas strong modulatory effects on LFP are instead observed in layers 2–5 in the adult. Our finding that behavioral modulation of gamma occurs specifically in layers 4 and 5 at eye opening is consistent with recent work that has reported the existence of a “developmental” gamma oscillation in somatosensory cortex that couples thalamocortical circuits and is generated within the thalamus (Minlebaev et al., 2011).
The development of inhibition and V1 function
Several important computational roles for broadly tuned feedforward inhibition in V1 have been proposed, including in the development of cortical circuits (Vidyasagar et al., 1996; Sompolinsky and Shapley, 1997; Ferster and Miller, 2000; Bartos et al., 2007; Runyan et al., 2010). In particular, the tuning response of inhibitory interneurons in mouse V1 broadens over the same developmental period that RF property selectivity increases in excitatory neurons (Kuhlman et al., 2011; Li et al., 2012). To assess the relationship between specific classes of inhibition and V1 development across cortical layers, we monitored the emergence and development of narrow waveforms from eye opening. However, we were unable to detect significant numbers of this response type before EO7. By EO7, we are able to detect a significant number of narrow waveforms and find that their OS is already broadly tuned, consistent with the time course observed in previous studies of the broadening of parvalbumin-expressing interneuron tuning over early development (Kuhlman et al., 2011; Li et al., 2012).
In summary, we provide a comprehensive picture of the progressive establishment of multiple key features of V1 function in the mouse as a function of cortical layer. This work will spur future studies to determine the molecular mechanisms underlying the development of cortical visual processing with greater spatial and temporal precision.
Footnotes
This work was supported by National Institutes of Health Grants 1RO1EY023337 and 1DP2EY023190 (C.M.N.) and 1F32EY024179 (J.L.H.). We thank Drs. Sunil Gandhi, Denise Piscopo, and Wayne Tschetter for comments on this manuscript and members of the Niell laboratory for helpful discussions.
The authors declare no competing financial interests.
- Correspondence should be addressed to Cristopher Niell, Institute of Neuroscience, 1254 University of Oregon, Eugene, OR 97403. cniell{at}uoregon.edu