- Original Paper
Functional emergence of a column-like architecture in layer 5 of mouse somatosensory cortex in vivo
The Journal of Physiological Sciences volume 69, pages 65–77 (2019)
To investigate how the functional architecture is organized in layer 5 (L5) of the somatosensory cortex of a mouse in vivo, the input–output relationship was investigated using an all-optical approach. The neural activity in L5 was optically recorded using a Ca2+ sensor, R-CaMP2, through a microprism inserted in the cortex under two-photon microscopy, while the L5 was regionally excited using optogenetics. The excitability was spread around the blue-light irradiated region, but the horizontal propagation was limited to within a certain distance (λ < 130 μm from the center of the illumination spot). When two regions were photostimulated with a short interval, the excitability of each cluster was reduced. Therefore, a column-like architecture had functionally emerged with reciprocal inhibition through a minimal number of synaptic relays. This could generate a synchronous output from a region of L5 with simultaneous enhancement of the signal-to-noise ratio by silencing of the neighboring regions.
Generally, the mammalian neocortex consists of six layers, each of which is characterized by its cellular and network organization . Neurons in each layer are also distinguished by their functional connectivity [2, 3]. For example, in the primary somatosensory cortex (S1), the major layer 4 neurons receive thalamic inputs and project locally to other layers such as layer 2/3. On the other hand, the layer 5 (L5) pyramidal neurons produce major outputs to other cortical as well as subcortical regions. In other words, the L5 is a primary cortical layer involved in the integration of various sensory modalities as well as the top-down control of other brain regions.
Some regions of the neocortex are also functionally organized in a columnar fashion. That is, neurons having similar response properties are often perpendicularly arrayed through layers. For example, in the primary visual area (V1), neurons with similar orientational/directional preference are grouped together, forming a column of 0.1–0.5 mm in width, and are segregated in distribution from other neurons of different preferences [4, 5]. During development, it has been computationally predicted and experimentally revealed that the spatiotemporal balance between excitatory and inhibitory connections in the local circuitry is one of the determinants of generating the regularly spaced columnar organization in V1 [6, 7]. Similarly, the functional architecture of the adult neocortex is presumed to be organized by the excitatory-inhibitory balance [8,9,10]. The cortical L5 consists of the excitatory pyramidal neurons and the inhibitory interneurons of various subclasses with a characteristic morphological, molecular and electrophysiological feature [11,12,13]. However, because of the technical difficulties, it remains to be elucidated how the individual neurons are orchestrated to organize the functional architecture in the in vivo cortex .
In the present study, the computational processing in the L5 network of mouse S1 was characterized by the input–output relationship. The neural activity was optically recorded while the L5 was regionally excited using optogenetics. Our all-optical approach revealed the functional emergence of a column-like cluster of active L5 neurons that were surrounded by inactive neurons, and that the column-like clusters inhibited each other. It is suggested that the excitatory projecting neurons and inhibitory interneurons form a circuitry of reciprocal inhibition to sculpture the output signals exiting from the L5.
Materials and methods
All experiments were carried out using wild-type C57BL/6 J mice (female, 8–12 weeks old, 20–30 gBW) in accordance with the animal experiment protocol approved by Tohoku University Committee for Animal Experiments (Approval no. 2017LsA-002) under the guidelines for Animal Experiments and Related Activities of Tohoku University as well as the guiding principles of the Physiological Society of Japan and the National Institutes of Health (NIH), USA. The number of animals in this study was kept to a minimum and, when possible, all animals were anesthetized to minimize their suffering. Animals had access to food and water ad libitum and were kept under a 12-h light–dark cycle.
Adeno-associated virus (AAV) Helper-Free System (Agilent Technologies, Inc., Santa Clara, CA, USA) was used for the generation of the AAV-CaMKII-C1V1-Venus (DJ) viral vector. The purification method was modified from a previously published protocol . Briefly, HEK293 cells were transfected with a pAAV-CaMKII-C1V1-Venus vector plasmid, pHelper and pAAV-RC (serotype DJ; purchased from Cell Biolabs Inc, San Diego, CA, USA), using a standard calcium phosphate method. Three days after transfection, cells were collected and suspended in artificial cerebrospinal fluid (ACSF) containing (in mM) 124 NaCl, 3 KCl, 26 NaHCO3, 2 CaCl2, 1 MgSO4, 1.25 KH2PO4, 10 d-glucose. Following multiple freeze–thaw cycles, the cell lysates were treated with benzonase nuclease (Merck, Darmstadt, Germany) at 37 °C for 30 min, and centrifuged 2 times at 16,000g for 10 min at 4 °C. The supernatant was used as the virus-containing solution. Quantitative PCR was performed to measure the titer of purified virus. Virus aliquots were then stored at − 80 °C until use in the experiment.
The AAV1-CAG-R-CaMP2-WPRE-SV40poly was produced using the baculovirus-Sf9 expression system as previously described . Briefly, the baculoviruses with the AAV construct and those with the RepCap1 helper construct were co-infected into Sf9 cells in a suspension culture. After 2–4 days, the cells were treated with a hypotonic lysis buffer [1% Triton X-100, 10 mM HEPES (pH 8.0), 3 mM NaCl and 0.5 mM MgCl2] for 15 min at room temperature. After inactivating the baculoviruses, the lysate was centrifuged and the supernatant was mixed with AVB Sepharose (GE Healthcare UK Ltd. Little Chalfont, England) for 1 h at room temperature. By applying this mixture to a gravity-flow column (Bio-Rad Laboratories, Inc., Hercules, CA, USA), the viruses were eluted using an elution buffer (50 mM glycine at pH 3.0) and were further concentrated using Amicon Ultra 100 K Centrifugal Concentrators (EMD Millipore-Merck KGaA, Darmstadt, Germany) with PBS. The final virus titer was 1.0 × 1013 gc/ml.
Cortical neurons were isolated from embryonic day-16 ICR mice (Japan SLC Inc., Shizuoka, Japan) using Nerve-Cells Dispersion Solutions (Sumitomo Bakelite, Tokyo, Japan) according to the manufacturer’s instructions, and grown in culture medium (Sumitomo Bakelite) under a 5% CO2 atmosphere at 37 °C, and the viral solution containing AAV-CaMKII-C1V1-Venus viruses was diluted in the culture medium (1/1000, finally 7 × 1010 copies/ml) 7 days in vitro (DIV). The cultured tissues were fixed with 4% paraformaldehyde (room temperature, 30 min) at 18 DIV and served for the immunohistochemistry.
After washing 3 times with phosphate buffer solution (PBS), the specimens were reacted with anti-GFP (rat serum, 1:250, 04404-84, Nacalai) and either anti-pan CaMKII (rabbit serum, 1:250, gift from Dr. Kohji Fukunaga, Tohoku University) or anti-GABA (rabbit serum, 1:250, A2052, Sigma-Aldrich, St. Louis, MO) followed by secondary antibodies conjugated with Alexa Fluor 546 (goat anti-rabbit IgG, 1:250, A11081, Thermo Fisher Scientific, Waltham, MA) and Alexa Fluor 488 (goat anti-rat IgG, 1:250, A11006, Thermo Fisher Scientific), respectively.
To evaluate the co-expression of CaMKII or GABA in the Venus-positive neurons, 10 fields (each 344,000 μm2) were selected randomly and viewed under a conventional epi-fluorescent microscope with 20× objective lens (Axiovert200, Carl Zeiss, Oberkochen, Germany) and the filter system either set at 17 (Alexa Fluor 488: excitation, 485 ± 20 nm; dichroic mirror, 510 nm; emission, 515–565 nm) or set at 20 (Alexa Fluor 546: excitation, 546 ± 12 nm; dichroic mirror, 560 nm; emission, 575–640 nm), and imaged using a cooled CCD camera system (VB-7010, KEYENCE CO, Osaka, Japan).
After administration of the anti-inflammatory agents, dexamethasone sodium phosphate (1.3 mg/kgBW, Banyu Pharmaceutical Co Ltd., Tokyo, Japan), sulfadiazine (24 mg/kgBW, Kyoritsu Seiyaku, Tokyo, Japan), trimethoprim (4.8 mg/kgBW, Kyoritsu Seiyaku) and carprofen (6 mg/kg, Pfizer Inc., New York City, USA), the mice were anesthetized by intraperitoneal injection of a ketamine-xylazine mixture (50 mg/kgBW ketamine, Daiichi Sankyo Co. Ltd., Tokyo, Japan and 10 mg/kgBW xylazine, Sigma-Aldrich, St. Louis, MO, USA). The skull was exposed, attached with a head chamber (CP-2, Narishige Inc., Tokyo, Japan) using glue (Loctite 4011, Henkel Co. Rocky Hill, CT06067, USA), and a hole was drilled (diameter, 3 mm) over the cortex S1 region at 2 mm lateral and 1 mm caudal from the bregma while being cleaned with physiological saline. The viral solution (R-CaMP2: 1.0 × 1013 gc/ml, C1V1: 2.0 × 1013 gc/ml) was injected at 0.5 μl/min using a glass pipette (tip diameter, 20 µm, G-1.5, Narishige) that was siliconized with chlorotri-n-butylsilane (Alfa Aesar, Massachusetts, USA), into the cortex at 500–700 µm below the surface in a region without obviously large blood vessels. The reflector-coated microprism (Micro prism 1 × 1 × 1, Nippon Electric Glass Co. Ltd., Otsu, Japan) was inserted (Fig. 1a, b) in order to see the rostral aspect of the coronary incision (width, 1 ~ 1.5 mm; depth, 1 ~ 1.5 mm) made by a microblade (377615, Beaver-Visitec International, Inc., Waltham, MA 02452 USA) and covered by a round cover glass (diameter, 2.7 mm; Matsunami Glass IND., Ltd., Osaka, Japan) with a head chamber (CP-2, Narishige) that was glued on the skull (Loctite 4011, Henkel). The mice recovered from the anesthesia on a heater (MP-914-NV, Vivaria Co., Osaka, Japan) while being administered intraperitoneally physiological saline containing 5% glucose (0.1 ml) for hydration and nutrition, returned to the cage and kept for 1 month with food and water containing antibacterial agents (0.6%, Shinoral®, Boehringer Ingelheim, Ingelheim, Germany).
To illuminate the targets directly through the objective lens we fabricated an epi-illumination system (ASKA Company, Katou, Japan) consisting of an LED array, a focusing lens, two objective lenses and a dichroic mirror (495 nm, FF495-Di03, Semrock, Inc., Rochester, NY, USA). The photostimulation was made by six blue LEDs (460 ± 10 nm, LXML-PB02, Philips Lumileds Lighting Co. San Jose, CA, USA) arrayed 2 × 3 (height × width) on a circuit board (p-ban.com Co., Tokyo, Japan) (Fig. 1c). Four green LEDs (520 ± 15 nm, LXML-PM01-0100, Philips Lumileds Lighting Co.) were set at the corners for position adjustment. Each LED was analog-regulated using an LED driver (RCD-24-0.70, RECOM Power GmbH., Vienna, Austria; Supplementary Figure 1) and a computer with a homemade driving program (C++). The photostimulation of L5 was made by the 3 blue LEDs in the lower line that were turned on individually for 4 ms with a variable interval between 5 and 9 s. The total number of photostimulations was set at 900. The power density (irradiance) of each LED was directly measured under microscopy by a visible light-sensing thermopile (MIR-101Q, SSC Co., Ltd., Kuwana City, Japan) and was 2.0 mW/mm2 at its focus.
Although the cortical tissue was hard to visualize for individual neurons immediately after surgery, it had cleared up over the following 3 weeks, probably because of the suppression of inflammation and gliosis. One month after surgery, each awake mouse was restrained under the microscope (ECLIPSE FN1, Nikon Co., Tokyo, Japan) with 16 × water-immersion objective lens (CFI75 LWD 16 × W, Nikon: NA, 0.8 and WD, 3.0) using a head-holding device (MAG-1, Narishige). The microscope was focused on the region 500–750 µm under the cortical surface (layer 5, L5) of S1 through the microprism, and the images for R-CaMP2 were acquired using a two-photon microscopy system (A1R MP+, Nikon) with a femtosecond-pulsed laser (Mai Tai® DeepSee™, Spectra-Physics, Newport Co., Santa Clara, CA, USA) at 1000 nm for excitation, and a suppression filter (629 ± 56 nm) for detection (Fig. 1d). The imaging field was 256 × 512 pixels (375 × 750 μm) at 50–150 µm away from the surface of the microprism and sampled at 3.6 fps.
The R-CaMP2 images were analyzed using ImageJ and a homemade program (C++) according to the following procedures. (1) Two-dimensional movement of the image was corrected frame-by-frame. (2) Every single frame capturing the LED-dependent noise at the instance of photostimulation (4 ms), which contained the stray LED light and the R-CaMP2 fluorescent light excited by the LED with/without direct activation of neurons, was removed from the image stack. (3) Circular regions of interest (ROIs) were manually set on the neuron-like cell bodies and their average brightness values were sequentially measured as a function of time. (4) For every 10 frames (duration, 10/3.6 ≈ 2.8 s) the standard deviation (SD) was calculated and its minimum was defined as noise SD (SD0). (5) The baseline fluorescence (Fb) of an ROI was updated with time in a series of records as the average fluorescence of 4 serial frames in which the range was within ± 3SD0. Otherwise, it was not updated. (6) The activated fluorescence (Fa) of the same ROI was defined as the value that exceeded + 6SD0 and its magnitude was expressed as ΔF/F = (Fa − Fb)/Fb. (7) The onset of each Ca2+ spike was determined according to the last baseline point just before the activated points.
All data in the text, figures and tables are expressed as mean ± SEM. The statistical evaluation was conducted using Fisher’s exact test for the frequency, Mann–Whitney U test for unpaired data, t-test for paired data, and one-way ANOVA with a Bonferroni post hoc test for data in order unless otherwise noted. The numbers of samples used for the statistics (n) are shown in the figure legends. It was judged as statistically insignificant when P > 0.05.
In the present study, the adult mice were compelled to express channelrhodopsin and Ca2+ sensor in the somatosensory cortex (S1) using AAV. To photostimulate the local network, a chimeric channelrhodopsin between channelrhodopsin-1 (ChR1) and one of the Volvox-derived channelrhodopsins, C1V1 , was expressed in L5 under the CaMKII promotor. To record the network activity, the L5 neurons expressed R-CaMP2 , one of the RFP-derived Ca2+ sensors, non-selectively under the CAG promotor. The C1V1-expressing neurons were identified immunohistochemically for the primary cultures from mouse cortex using anti-Venus antibodies. Among 111 Venus-positive neurons, 86 were obviously positive with CaMKII (Fig. 2a–c), whereas 4 were positive with GABA among 142 Venus-positive neurons (Fig. 2d–f). Thus, the probability of GABA co-expression was negligible relative to that of CaMKII (Fig. 2g).
The L5 neurons were visualized through a microprism inserted in the somatosensory cortex under two-photon microscopy (Fig. 3a, Supplementary Figure 2 and Supplementary Video 1). That is, neurons deep in the L5 were viewed from the side. On the other hand, the light from an array of blue LED chips was focused on the imaging plane through the objective lens and the microprism. Each LED chip illuminated one of the three circular targets (diameter, ~ 100 μm), A, B or C, which were aligned in the L5 parallel to the cortical surface (Fig. 3a). The LED chips were turned on and off by one of the following modes of spatiotemporal pattern: the “single-mode”, by which each of 3 targets was illuminated one by one in a random order and the “multi-mode”, by which 2 of 3 targets were illuminated simultaneously or sequentially with an interval of either 2, 4 or 12 ms. These two modes were scrambled in order while imaging the R-CaMP2 fluorescence (Fig. 3b).
Response to the single-mode photostimulation
A number of R-CaMP2-positive neurons were identified in an imaging plane set in the L5 of S1 (Fig. 3a). While transient increments of ∆F/F, which were judged as Ca2+ spikes, occurred spontaneously under the awakened condition, they were often evoked immediately after photostimulation (Fig. 3b). For every L5 neuron, the Ca2+ spikes were aligned to the single-mode photostimulation at each of 3 targets, A, B and C, as shown by the sample raster plots in Fig. 3c–e. As shown in the histograms in Fig. 3f–h, the frequency of the Ca2+ spikes was significantly increased in the first scanning frame (~ 280 ms) just after the single-mode photostimulation at target B (P < 0.05, Fisher’s exact test), but not at either target A or C. Therefore, this neuron was judged to be the target B-responsive neuron. A target-responsive neuron was thus defined as a ROI in which the frequency of Ca2+ spikes in response to the stimulation of the target was significantly larger than that in response to either of the other two targets. In summary, of 354 neurons in total (n = 5 animals), 224 were the target-responsive neurons (65 ± 9%, n = 5), 121 (32 ± 7%, n = 5) were non-responsive and 9 (3 ± 2%, n = 5) showed significant decreases in the frequency of Ca2+ spikes after the single-mode photostimulation (P < 0.05, Fisher’s exact test).
The target A-, target B-, and target C-responsive neurons were respectively distributed close to the targets A, B and C, although a small number of neurons were responsive to multiple targets (Fig. 4a–c). Therefore, the excitatory influence of photostimulation, directly or indirectly, diminished with distance along a direction parallel to the cortical surface. To further analyze this, the imaging plane was divided into perpendicular stripes of 25 µm width (Fig. 4d, top) and the response probability of neurons in a stripe was averaged for each target to derive the ensemble response probability (p). As shown in Fig. 4d (bottom), p was highest around at the center of the target and decreased with distance (x). The relationship almost followed a single exponential function (Fig. 4d),
where λ is the length constant. The λ values were similar among targets of photostimulation and among animals in a range of 37–129 µm (63 ± 7 µm, n = 14 targets, 5 animals).
Response to the multi-mode photostimulation
Would the simultaneous photostimulation of multiple targets enhance the responsiveness of a neuron? As shown in the example shown in Fig. 5a, the ensemble response probability of neurons in a perpendicular stripe to simultaneous photostimulation at targets A and B (A + B) tended to be smaller than that to single-mode stimulation at target A or B. Therefore, the photostimulations at target A and B reciprocally inhibit each other. Similar reciprocal inhibition was present between photostimulations at targets B and C (Fig. 5b) or between photostimulations at targets A and C (Fig. 5c). To evaluate the reciprocal inhibition quantitatively, neurons distributed in the region of 6 stripes around each target (range ± 75 µm from the center) were selected and their ensemble response probabilities were compared between single-mode and multi-mode photostimulations (Fig. 5d–f). The response probability to the simultaneous photostimulation was indeed smaller than that to single-mode photostimulation with significant (or nearly significant) difference. In summary, the ensemble response probability to the simultaneous multi-mode photostimulation was 0.16 ± 0.01 (n = 30 cases, 5 animals), which was significantly smaller than that to the single-mode photostimulation, 0.25 ± 0.02 (n = 15 cases, 5 animals) (P < 0.05, Mann–Whitney U test).
The magnitude of reciprocal inhibition was actually dependent on the interval between two photostimulations. To evaluate the effect of target B photostimulation on the target A responsiveness (A/B), for example, the neurons distributed in the region of 6 stripes around the target were again selected, and the relative ensemble probabilities were compared for photostimulations at B before those at A by 12 ms (A/B, − 12), by 4 ms (A/B, − 4), by 2 ms (A/B, − 2), simultaneously at A and B (A/B, 0), at B after A by 2 ms (A/B, 2), by 4 ms (A/B, 4) and by 12 ms (A/B, 12). As shown in Fig. 6a, the relative probability was dependent on the interval in any combination of photostimulations: minimal at 0 ms and almost negligible at ± 12 ms. Similar multi-mode photostimulation was then subdivided into two classes: (1) multi-mode photostimulation of two close targets (A/B, B/A, B/C or C/B), and (2) that of two remote targets (A/C or C/A), and summarized for 5 animals (Fig. 6b). As there was no significant difference between these two classes at any interval (P > 0.2), the effects of the interval were evaluated for all pairs (n = 30, 5 animals) and found to be asymmetrical. For example, the relative probability at − 4 and − 2 was respectively smaller than that at 4 and 2 with significant difference (P < 0.05 and < 0.01). Although the differences between − 12 vs − 4 and 0 vs 2 were significant (P < 0.01), the differences between − 2 vs 0 and 4 vs 12 were insignificant (P > 0.05).
The neuronal network in the brain can be regarded as a kind of multi-dimensional transfer function of time and has been investigated through analyzing the input–output relationship. The optical imaging of neuronal activity using genetically encoded Ca2+ indicators or voltage sensors would be advantageous over other methods in the high spatiotemporal resolutions discriminating individual neurons [19,20,21,22]. As an artificial input, the optical stimulation using optogenetics has advantages over other methods such as electrical stimulation in its high spatiotemporal precision [23,24,25,26,27,28]. However, the power of visible light is massively attenuated with distance because of scattering and absorption by the brain tissue [25, 29]. On the other hand, a strong light would actuate neurons non-selectively in its path. To solve this, two-photon optics has been applied for both photostimulation of target neurons and imaging of neuronal activity in the mouse brain in vivo: L2/3 of somatosensory cortex  and visual cortex . Although the conventional diffraction-limited two-photon laser system has been rather ineffective for exciting neurons due to the small single-channel conductance of channelrhodopsins (40 fS ~ 1.1 pS for ChR2 and ChR1/2 chimera) [32,33,34] and insufficient two-photon absorption volume (~ 2–5 μm3) , attempts have been made to apply the high energy of a laser over the cell of interest for a short duration [28, 36,37,38].
Recently, deep cortical layers have been imaged by combining a microscope with a microprism inserted in the in vivo brain . As a result, the microscopic light is expected to be focused on the vertical visual field directly without mediating tissue while collecting the sample light with minimal diffraction. In the present study, the blue LED light was directly focused in the L5 of mouse cortex through a microprism to excite the neurons that expressed C1V1. Our approach has several advantages and disadvantages. It can excite a homogeneous group of neurons at once because of the localization of a light path within the layer, although is difficult to differentiate individual cells for excitation because of the low spatial resolution of a single photon beam. The on–off timing of irradiation can be precisely regulated for multiple targets when each LED is individually regulated. It can be combined with any microscopic system with minimal modification. However, the possibility should be kept in mind that the two-photon imaging could depolarize the membrane potential of a C1V1-expressing neuron with increased excitability. Although some of the cortical connections were disrupted or disorganized by the chronic insertion of a microprism, the basic properties of the local circuit should remain intact . Since the LED pulse was short in duration (4 ms), it would evoke only a few action potentials in neurons that expressed C1V1 [17, 40]. On the other hand, the Ca2+ spike is assumed to be generated by a burst of action potentials because of the response characteristics of R-CaMP2 . Although we made a ROI on a neuron-like cell body, the measured fluorescence may have been contaminated with that from the neighboring neuropil. This approach could be improved by increasing the resolution of the irradiating apparatus using a digital micromirror device (DMD), liquid crystal display, and liquid crystal on a silicon/spatial light modulator (LCoS/SLM) [41,42,43,44]. For example, point irradiation of blue light on the background of yellow/red light would excite individual cells that express step-function opsins (SFO) in the membrane with high precision and efficiency [45,46,47].
Functional emergence of column-like architecture
The present study is the first all-optical approach to investigate the input–output characteristics of the deep cortical layer circuitry in an awake mouse. The excitability was spread around the blue-light irradiated region, but the horizontal propagation was limited within a certain distance (λ < 130 μm from the center of the spot) and within a single frame of imaging (< 0.3 s). Therefore, a column-like excitatory cluster 100–300 μm in diameter appeared to emerge transiently in response to an input into L5. Some neurons in the cluster may have expressed C1V1 to be excited directly by the blue light irradiation on their soma, dendrites and/or axons, and others should have been excited indirectly through synapses. In addition, a C1V1-expressing neuron could also be synaptically activated to generate Ca2+ spikes even if the light-dependent depolarization of its membrane potential would be subthreshold for generating action potentials. Nevertheless, the column-like cluster, which consists of the directly and indirectly activated neurons, functionally emerged with regional photostimulation of L5 with a spatiotemporal restriction. The size of the column-like cluster is consistent with that of functional columns described previously [3,4,5, 48], and that predicted from the receptive field of individual neurons . It is possible that this column-like excitatory cluster consists of smaller substructures, such as an L5 microcolumn/lattice [49, 50] and cortical microcircuit [51,52,53,54]. The functional organization of these substructures should be revealed in the future by differential optogenetic excitation of neuronal subtypes [55, 56] that have been classified by morphology, projection targets, electrophysiological properties, local excitatory connectivity and long-range inputs : the corticocallosally projecting slender-tufted (L5A) and subcortical targeting thick tufted (L5B) pyramidal neurons [58,59,60,61,62]; the corticopontine (CPn) and crossed corticostriatal (CCS) neurons .
The neocortical circuit is composed of excitatory principal neurons and inhibitory interneurons. As a result, the excitation of L5 pyramidal neurons should inhibit other neurons and vice versa through inhibitory interneurons in the L5 and/or via other layers . It has been postulated that the L5 output signal should be sculptured by the excitatory-inhibitory balance [8, 13]. This notion was confirmed by the present study showing that the column-like excitatory clusters inhibited each other reciprocally with a minimal number of the synaptic relays because the inhibition was maximal when neurons in another target were excited before 0–2 ms. However, there was a negligible possibility that the inhibitory neurons were directly activated by light because of the negligible expression of C1V1 under the CaMKII promotor using AAV. The inhibition propagated horizontally without obvious attenuation: the inhibition between adjacent targets (distance, ~ 150 µm) was similar to that between remote targets (distance, ~ 300 µm). These facts are consistent with that some types of cortical interneurons such as the parvalbumin (PV)-expressing basket cells and chandelier cells have long-range connections with L5 pyramidal neurons [64, 65]. These neurons as well as other types of interneurons, such as the somatostatin (SOM)-expressing Martinotti cells [66, 67], may be involved in the reciprocal inhibition. Once an input is made regionally in the L5, the column-like excitatory cluster should be self-organized through local excitatory connections among L5 pyramidal neurons and sculptured spatially and temporally by the long-range reciprocal inhibitory connections. This could generate a synchronous output from a region of L5 with simultaneous enhancement of the signal-to-noise ratio by the silencing of neighboring regions [6, 68, 69]. It is hypothesized that the long-range inhibitory signals also sculpture other layers such as layer 2/3 to help the functional emergence of a column-like architecture through the layers in a manner dependent on the L5 inputs . In the primary sensory cortex, the thalamic inputs would excite L5 pyramidal neurons directly or indirectly through activation of layer 4 [3, 13, 56]. The regional excitation of L5 would, in turn, conduct the cortical circuitry to generate a column-like functional architecture. This hypothesis could be tested in the future through an all-optical approach in combination with a multilayer imaging system .
Harris KD, Mrsic-Flogel TD (2013) Cortical connectivity and sensory coding. Nature 503:51–58. https://doi.org/10.1038/nature12654
Zeng H, Sanes JR (2017) Neuronal cell-type classification: challenges, opportunities and the path forward. Nat Rev Neurosci 18:530–546. https://doi.org/10.1038/nrn.2017.85
Lübke J, Feldmeyer D (2007) Excitatory signal flow and connectivity in a cortical column: focus on barrel cortex. Brain Struct Funct. 212:3–17. https://doi.org/10.1007/s00429-007-0144-2
Hubel DH, Wiesel TN (1962) Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex. J Physiol 160:106–154. https://doi.org/10.1113/jphysiol.1962.sp006837
Hubel DH, Wiesel TN (1977) Ferrier lecture. Functional architecture of macaque monkey visual cortex. Proc R Soc Lond B Biol Sci 198:1–59. https://doi.org/10.1098/rspb.1977.0085
Miller KD, Keller JB, Stryker MP (1989) Ocular dominance column development: analysis and simulation. Science 245:605–615. https://doi.org/10.1126/science.2762813
Hensch TK, Stryker MP (2004) Columnar architecture sculpted by GABA circuits in developing cat visual cortex. Science 303:1678–1681. https://doi.org/10.1126/science.1091031
Isaacson JS, Scanziani M (2011) How inhibition shapes cortical activity. Neuron 72:231–243. https://doi.org/10.1016/j.neuron.2011.09.027
Lee SH, Kwan AC, Zhang S, Phoumthipphavong V, Flannery JG, Masmanidis SC, Taniguchi H, Huang ZJ, Zhang F, Boyden ES, Deisseroth K, Dan Y (2012) Activation of specific interneurons improves V1 feature selectivity and visual perception. Nature 488:379–383. https://doi.org/10.1038/nature11312
Kim K, Kim JH, Song YH, Lee SH (2017) Functional dissection of inhibitory microcircuits in the visual cortex. Neurosci Res 116:70–76. https://doi.org/10.1016/j.neures.2016.09.003
Markram H, Toledo-Rodriguez M, Wang Y, Gupta A, Silberberg G, Wu C (2004) Interneurons of the neocortical inhibitory system. Nat Rev Neurosci 5:793–807. https://doi.org/10.1038/nrn1519
Kubota Y, Karube F, Nomura M, Kawaguchi Y (2016) The diversity of cortical inhibitory synapses. Front Neural Circuits. 10:27. https://doi.org/10.3389/fncir.2016.00027
Markram H, Muller E, Ramaswamy S, Reimann MW, Abdellah M, Sanchez CA, Ailamaki A, Alonso-Nanclares L, Antille N, Arsever S, Kahou GA, Berger TK, Bilgili A, Buncic N, Chalimourda A, Chindemi G, Courcol JD, Delalondre F, Delattre V, Druckmann S, Dumusc R, Dynes J, Eilemann S, Gal E, Gevaert ME, Ghobril JP, Gidon A, Graham JW, Gupta A, Haenel V, Hay E, Heinis T, Hernando JB, Hines M, Kanari L, Keller D, Kenyon J, Khazen G, Kim Y, King JG, Kisvarday Z, Kumbhar P, Lasserre S, Le Bé JV, Magalhães BR, Merchán-Pérez A, Meystre J, Morrice BR, Muller J, Muñoz-Céspedes A, Muralidhar S, Muthurasa K, Nachbaur D, Newton TH, Nolte M, Ovcharenko A, Palacios J, Pastor L, Perin R, Ranjan R, Riachi I, Rodríguez JR, Riquelme JL, Rössert C, Sfyrakis K, Shi Y, Shillcock JC, Silberberg G, Silva R, Tauheed F, Telefont M, Toledo-Rodriguez M, Tränkler T, Van Geit W, Díaz JV, Walker R, Wang Y, Zaninetta SM, DeFelipe J, Hill SL, Segev I, Schürmann F (2015) Reconstruction and simulation of neocortical microcircuitry. Cell 163:456–492. https://doi.org/10.1016/j.cell.2015.09.029
Reid RC (2012) From functional architecture to functional connectomics. Neuron 75:209–217. https://doi.org/10.1016/j.neuron.2012.06.031
Inutsuka A, Yamashita A, Chowdhury S, Nakai J, Ohkura M, Taguchi T, Yamanaka A (2016) The integrative role of orexin/hypocretin neurons in nociceptive perception and analgesic regulation. Sci Rep 6:29480. https://doi.org/10.1038/srep29480
Kawashima T, Kitamura K, Suzuki K, Nonaka M, Kamijo S, Takemoto-Kimura S, Kano M, Okuno H, Ohki K, Bito H (2013) Functional labeling of neurons and their projections using the synthetic activity-dependent promoter E-SARE. Nat Methods 10:889–895. https://doi.org/10.1038/nmeth.2559
Yizhar O, Fenno LE, Prigge M, Schneider F, Davidson TJ, O’Shea DJ, Sohal VS, Goshen I, Finkelstein J, Paz JT, Stehfest K, Fudim R, Ramakrishnan C, Huguenard JR, Hegemann P, Deisseroth K (2011) Neocortical excitation/inhibition balance in information processing and social dysfunction. Nature 477:171–178. https://doi.org/10.1038/nature10360
Inoue M, Takeuchi A, Horigane S, Ohkura M, Gengyo-Ando K, Fujii H, Kamijo S, Takemoto-Kimura S, Kano M, Nakai J, Kitamura K, Bito H (2015) Rational design of a high-affinity, fast, red calcium indicator R-CaMP2. Nat Methods 12:64–70. https://doi.org/10.1038/nmeth.3185
Nagai T, Horikawa K, Saito K, Matsuda T (2014) Genetically encoded Ca2+ indicators; expanded affinity range, color hue and compatibility with optogenetics. Front Mol Neurosci 7:90. https://doi.org/10.3389/fnmol.2014.00090
Antic SD, Empson RM, Knöpfel T (2016) Voltage imaging to understand connections and functions of neuronal circuits. J Neurophysiol 116:135–152. https://doi.org/10.1152/jn.00226.2016
Storace D, Sepehri Rad M, Kang B, Cohen LB, Hughes T, Baker BJ (2016) Toward better genetically encoded sensors of membrane potential. Trends Neurosci 39:277–289. https://doi.org/10.1016/j.tins.2016.02.005
Sepehri Rad M, Choi Y, Cohen LB, Baker BJ, Zhong S, Storace DA, Braubach OR (2017) Voltage and calcium imaging of brain activity. Biophys J 113:2160–2167. https://doi.org/10.1016/j.bpj.2017.09.040
Boyden ES, Zhang F, Bamberg E, Nagel G, Deisseroth K (2005) Millisecond-timescale, genetically targeted optical control of neural activity. Nat Neurosci 8:1263–1268. https://doi.org/10.1038/nn1525
Ishizuka T, Kakuda M, Araki R, Yawo H (2006) Kinetic evaluation of photosensitivity in genetically engineered neurons expressing green algae light-gated channels. Neurosci Res 54:85–94. https://doi.org/10.1016/j.neures.2005.10.009
Yizhar O, Fenno LE, Davidson TJ, Mogri M, Deisseroth K (2011) Optogenetics in neural systems. Neuron 71:9–34. https://doi.org/10.1016/j.neuron.2011.06.004
Yawo H, Asano T, Sakai S, Ishizuka T (2013) Optogenetic manipulation of neural and non-neural functions. Dev Growth Differ 55:474–490. https://doi.org/10.1111/dgd.12053
Oertner TG, Helmchen F, de Lecea L, Beck H, Konnerth A, Kaupp B, Knöpfel T, Yawo H, Häusser M (2013) Optogenetic analysis of mammalian neural circuits. In: Hegemann P, Sigrist S (eds) Optogenetics, Dahlem Workshop Reports. De Gruyter, Berlin, pp 109–126
Yawo H, Egawa R, Hososhima S, Wen L (2015) General description: future prospects of optogenetics. In: Kandori H, Yawo H, Koizumi A (eds) Optogenetics: light-sensing proteins and their applications. Springer, Tokyo, pp 111–132
Aravanis AM, Wang LP, Zhang F, Meltzer LA, Mogri MZ, Schneider MB, Deisseroth K (2007) An optical neural interface: in vivo control of rodent motor cortex with integrated fiberoptic and optogenetic technology. J Neural Eng 4:S143–S156. https://doi.org/10.1088/1741-2560/4/3/S02
Packer AM, Russell LE, Dalgleish HW, Häusser M (2015) Simultaneous all-optical manipulation and recording of neural circuit activity with cellular resolution in vivo. Nat Methods 12:140–146. https://doi.org/10.1038/nmeth.3217
Carrillo-Reid L, Yang W, Bando Y, Peterka DS, Yuste R (2016) Imprinting and recalling cortical ensembles. Science 353:691–694. https://doi.org/10.1126/science.aaf7560
Nagel G, Szellas T, Huhn W, Kateriya S, Adeishvili N, Berthold P, Ollig D, Hegemann P, Bamberg E (2003) Channelrhodopsin-2, a directly light-gated cation-selective membrane channel. Proc Natl Acad Sci USA 100:13940–13945. https://doi.org/10.1073/pnas.1936192100
Feldbauer K, Zimmermann D, Pintschovius V, Spitz J, Bamann C, Bamberg E (2009) Channelrhodopsin-2 is a leaky proton pump. Proc Natl Acad Sci U S A. 106:12317–12322. https://doi.org/10.1073/pnas.0905852106
Lin JY, Lin MZ, Steinbach P, Tsien RY (2009) Characterization of engineered channelrhodopsin variants with improved properties and kinetics. Biophys J 96:1803–1814. https://doi.org/10.1016/j.bpj.2008.11.034
Papagiakoumou E, Anselmi F, Bègue A, de Sars V, Glückstad J, Isacoff EY, Emiliani V (2010) Scanless two-photon excitation of channelrhodopsin-2. Nat Methods 7:848–854. https://doi.org/10.1038/nmeth.1505
Peron S, Svoboda K (2011) From cudgel to scalpel: toward precise neural control with optogenetics. Nat Methods 8:30–34. https://doi.org/10.1038/nmeth.f.325
Chaigneau E, Ronzitti E, Gajowa MA, Soler-Llavina GJ, Tanese D, Brureau AY, Papagiakoumou E, Zeng H, Emiliani V (2016) Two-photon holographic stimulation of ReaChR. Front Cell Neurosci 10:234. https://doi.org/10.3389/fncel.2016.00234
Ronzitti E, Conti R, Zampini V, Tanese D, Foust AJ, Klapoetke N, Boyden ES, Papagiakoumou E, Emiliani V (2017) Submillisecond optogenetic control of neuronal firing with two-photon holographic photoactivation of Chronos. J Neurosci 37:10679–10689. https://doi.org/10.1523/JNEUROSCI.1246-17.2017
Andermann ML, Gilfoy NB, Goldey GJ, Sachdev RN, Wölfel M, McCormick DA, Reid RC, Levene MJ (2013) Chronic cellular imaging of entire cortical columns in awake mice using microprisms. Neuron 80:900–913. https://doi.org/10.1016/j.neuron.2013.07.052
Hososhima S, Yuasa H, Ishizuka T, Hoque MR, Yamashita T, Yamanaka A, Sugano E, Tomita H, Yawo H (2015) Near-infrared (NIR) up-conversion optogenetics. Sci Rep 5:16533. https://doi.org/10.1038/srep16533
Sakai S, Ueno K, Ishizuka T, Yawo H (2013) Parallel and patterned optogenetic manipulation of neurons in the brain slice using a DMD-based projector. Neurosci Res 75:59–64. https://doi.org/10.1016/j.neures.2012.03.009
Favre-Bulle IA, Preece D, Nieminen TA, Heap LA, Scott EK, Rubinsztein-Dunlop H (2015) Scattering of sculpted light in intact brain tissue, with implications for optogenetics. Sci Rep 5:11501. https://doi.org/10.1038/srep11501
Roy A, Osik JJ, Ritter NJ, Wang S, Shaw JT, Fiser J, Van Hooser SD (2016) Optogenetic spatial and temporal control of cortical circuits on a columnar scale. J Neurophysiol 115:1043–1062. https://doi.org/10.1152/jn.00960.2015
Zhao ZD, Yang WZ, Gao C, Fu X, Zhang W, Zhou Q, Chen W, Ni X, Lin JK, Yang J, Xu XH, Shen WL (2017) A hypothalamic circuit that controls body temperature. Proc Natl Acad Sci USA 114:2042–2047. https://doi.org/10.1073/pnas.1616255114
Berndt A, Yizhar O, Gunaydin LA, Hegemann P, Deisseroth K (2009) Bi-stable neural state switches. Nat Neurosci 12:229–234. https://doi.org/10.1038/nn.2247
Bamann C, Gueta R, Kleinlogel S, Nagel G, Bamberg E (2010) Structural guidance of the photocycle of channelrhodopsin-2 by an interhelical hydrogen bond. Biochemistry 49:267–278. https://doi.org/10.1021/bi901634p
Hososhima S, Sakai S, Ishizuka T, Yawo H (2015) Kinetic evaluation of photosensitivity in bi-stable variants of chimeric channelrhodopsins. PLoS One 10:e0119558. https://doi.org/10.1371/journal.pone.0119558
Mountcastle VB (1997) The columnar organization of the neocortex. Brain 120:701–722
Maruoka H, Kubota K, Kurokawa R, Tsuruno S, Hosoya T (2011) Periodic organization of a major subtype of pyramidal neurons in neocortical layer V. J Neurosci 31:18522–18542. https://doi.org/10.1523/JNEUROSCI.3117-11.2011
Maruoka H, Nakagawa N, Tsuruno S, Sakai S, Yoneda T, Hosoya T (2017) Lattice system of functionally distinct cell types in the neocortex. Science 358:610–615. https://doi.org/10.1126/science.aam6125
Yoshimura Y, Dantzker JL, Callaway EM (2005) Excitatory cortical neurons form fine-scale functional networks. Nature 433:868–873. https://doi.org/10.1038/nature03252
Oberlaender M, Boudewijns ZS, Kleele T, Mansvelder HD, Sakmann B, de Kock CP (2011) Three-dimensional axon morphologies of individual layer 5 neurons indicate cell type-specific intracortical pathways for whisker motion and touch. Proc Natl Acad Sci USA 108:4188–4193. https://doi.org/10.1073/pnas.1100647108
Ko H, Cossell L, Baragli C, Antolik J, Clopath C, Hofer SB, Mrsic-Flogel TD (2013) The emergence of functional microcircuits in visual cortex. Nature 496:96–100. https://doi.org/10.1038/nature12015
Ishikawa AW, Komatsu Y, Yoshimura Y (2014) Experience-dependent emergence of fine-scale networks in visual cortex. J Neurosci 34:12576–12586. https://doi.org/10.1523/JNEUROSCI.1346-14.2014
Narayanan RT, Udvary D, Oberlaender M (2017) Cell type-specific structural organization of the six layers in rat barrel cortex. Front Neuroanat 11:91. https://doi.org/10.3389/fnana.2017.00091
Feldmeyer D (2012) Excitatory neuronal connectivity in the barrel cortex. Front Neuroanat 6:24. https://doi.org/10.3389/fnana.2012.00024
Naka A, Adesnik H (2016) Inhibitory circuits in cortical layer 5. Front Neural Circuits 10:35. https://doi.org/10.3389/fncir.2016.00035
Molnár Z, Cheung AF (2006) Towards the classification of subpopulations of layer V pyramidal projection neurons. Neurosci Res 55:105–115. https://doi.org/10.1016/j.neures.2006.02.008
Molyneaux BJ, Arlotta P, Menezes JR, Macklis JD (2007) Neuronal subtype specification in the cerebral cortex. Nat Rev Neurosci 8:427–437. https://doi.org/10.1038/nrn2151
Le Bé JV, Silberberg G, Wang Y, Markram H (2007) Morphological, electrophysiological, and synaptic properties of corticocallosal pyramidal cells in the neonatal rat neocortex. Cereb Cortex 17:2204–2213. https://doi.org/10.1093/cercor/bhl127
Brown SP, Hestrin S (2009) Intracortical circuits of pyramidal neurons reflect their long-range axonal targets. Nature 457:1133–1136. https://doi.org/10.1038/nature07658
Sato TR, Svoboda K (2010) The functional properties of barrel cortex neurons projecting to the primary motor cortex. J Neurosci 30:4256–4260. https://doi.org/10.1523/JNEUROSCI.3774-09.2010
Morishima M, Kawaguchi Y (2006) Recurrent connection patterns of corticostriatal pyramidal cells in frontal cortex. J Neurosci 26:4394–4405. https://doi.org/10.1523/JNEUROSCI.0252-06.2006
Kawaguchi Y, Kubota Y (1998) Neurochemical features and synaptic connections of large physiologically identified GABAergic cells in the rat frontal cortex. Neuroscience 85:677–701. https://doi.org/10.1016/S0306-4522(97)00685-4
Gupta A, Wang Y, Markram H (2000) Organizing principles for a diversity of GABAergic interneurons and synapses in the neocortex. Science 287:273–278. https://doi.org/10.1126/science.287.5451.273
Wang Y, Toledo-Rodriguez M, Gupta A, Wu C, Silberberg G, Luo J, Markram H (2004) Anatomical, physiological and molecular properties of Martinotti cells in the somatosensory cortex of the juvenile rat. J Physiol 561:65–90. https://doi.org/10.1113/jphysiol.2004.073353
Silberberg G, Markram H (2007) Disynaptic inhibition between neocortical pyramidal cells mediated by Martinotti cells. Neuron 53(5):735–746. https://doi.org/10.1016/j.neuron.2007.02.012
Bush P, Sejnowski T (1996) Inhibition synchronizes sparsely connected cortical neurons within and between columns in realistic network models. J Comput Neurosci 3:91–110 PMID:8840227
Bush PC, Mainen ZF (2015) Columnar architecture improves noise robustness in a model cortical network. PLoS ONE 10:e0119072. https://doi.org/10.1371/journal.pone.0119072
Otsuka T, Kawaguchi Y (2009) Cortical inhibitory cell types differentially form intralaminar and interlaminar subnetworks with excitatory neurons. J Neurosci 29:10533–10540. https://doi.org/10.1523/JNEUROSCI.2219-09.2009
We thank K. Fukunaga for the generous gift of rabbit pan CaMKII antiserum and B. Bell for language assistance. This work was supported by a Grant-in-Aid for Scientific Research (No. 25250001 to HY), a Grant-in-Aid for Scientific Research on Innovative Areas (Adaptive Circuit Shift: No. 15H01413 to HY; Comprehensive Brain Science Network; Platforms for Advanced Technologies and Research Resources: No. 16H06276 to HB) of the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan: http://www.jsps.go.jp/english/e-grants/grants01.html as well as JST, Strategic International Collaborative Research Program, SICORP (to HY): http://www.jst.go.jp/inter/english/sicorp/index.html and CREST (JPMJCR1656 to AY).
Conflict of interest
The authors declare that they have no conflicts of interest.
Statement on the welfare of animals
All applicable international, national, and/or institutional guidelines for the care and use of animals were followed. All procedures performed in studies involving animals were in accordance with the ethical standards of the institution or practice at which the studies were conducted. This article does not contain any studies with human participants performed by any of the authors.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary Figure 2 Cortical organization. The neurons expressing R-CaMP2 in layer 2/3 and 5 were imaged through a microprism under two-photon microscopy. Scale, 100 µm. (PDF 485 kb)
Supplementary Video 1 A sample R-CaMP2 fluorescence image at layer 5 through the prism (the same region shown in Figure 3A). The targets of photostimulation, A, B, and C are overlaid (circles). The illuminated targets are colored in blue in the frame capturing the LED-dependent noise at the instant of photostimulation (AVI 1836 kb)
About this article
Cite this article
Koizumi, K., Inoue, M., Chowdhury, S. et al. Functional emergence of a column-like architecture in layer 5 of mouse somatosensory cortex in vivo. J Physiol Sci 69, 65–77 (2019). https://doi.org/10.1007/s12576-018-0618-4