Layer 1 of somatosensory cortex: an important site for input to a tiny cortical compartment

Abstract Neocortical layer 1 has been proposed to be at the center for top-down and bottom-up integration. It is a locus for interactions between long-range inputs, layer 1 interneurons, and apical tuft dendrites of pyramidal neurons. While input to layer 1 has been studied intensively, the level and effect of input to this layer has still not been completely characterized. Here we examined the input to layer 1 of mouse somatosensory cortex with retrograde tracing and optogenetics. Our assays reveal that local input to layer 1 is predominantly from layers 2/3 and 5 pyramidal neurons and interneurons, and that subtypes of local layers 5 and 6b neurons project to layer 1 with different probabilities. Long-range input from sensory-motor cortices to layer 1 of somatosensory cortex arose predominantly from layers 2/3 neurons. Our optogenetic experiments showed that intra-telencephalic layer 5 pyramidal neurons drive layer 1 interneurons but have no effect locally on layer 5 apical tuft dendrites. Dual retrograde tracing revealed that a fraction of local and long-range neurons was both presynaptic to layer 5 neurons and projected to layer 1. Our work highlights the prominent role of local inputs to layer 1 and shows the potential for complex interactions between long-range and local inputs, which are both in position to modify the output of somatosensory cortex.


Introduction
Neocortical layer (L) 1, long an enigma (Hubel 1982;Burkhalter 1989;Marín-Padilla 1998;G ȃm ȃnuţ et al. 2018;Ibrahim et al. 2020), has recently become central to ideas about consciousness and perception (Aru et al. 2020;Gidon et al. 2020;Guest et al. 2021;Ibrahim et al. 2021;Larkum 2013;Schroeder et al. 2022;Suzuki and Larkum 2020;Takahashi et al. 2016Takahashi et al. , 2020)).These ideas seem astonishing for such a cell sparse, tiny layer, located at the surface of cortex, spanning only 100 microns in width in mice, and containing only GABAergic interneurons (Rudy et al. 2011;Schuman et al. 2021).While L1 is sparse in cell bodies, it is filled with processes, both dendrites of a variety of local pyramidal neurons and interneurons, and axons of a variety of inputs (Cauller 1995;DeFelipe and Fariñas 1992;DeFelipe 1997;Gabbott and Somogyi 1986;Ibrahim et al. 2020;Ito et al. 1998;Karimi et al. 2020;Larkum 2013;Lee et al. 2010;Marín-Padilla 1998;Mason and Larkman 1990;Schuman et al. 2019).A consistent theme for L1 in all cortices is that L2/3 and L5 pyramidal neurons have dendrites in L1, which makes this layer a potential locus for contextual interactions between the dendrites of the pyramidal neurons and long-range feedback inputs (Rudy et al. 2011;Cruikshank et al. 2012;Abs et al. 2018;Doron et al. 2020;Anastasiades et al. 2021;Guest et al. 2021;Schuman et al. 2021).This organization of inputs and apical tuft dendrites has led to the proposal that the apical tuft dendrites in L1 are a key element for sensory perception and learning.
One long-standing-related idea is that primary sensory cortical areas generate feedforward input to higher-order cortical areas, to the middle layers of cortex, which, in turn, provide feedback to lower order cortical areas, to the outer layers-including L1of cortex (Rockland and Pandya 1979;Felleman and Van Essen 1991;Markov et al. 2014).This work suggested specific laminar profiles of cortico-cortical connectivity with a class of longrange inputs arising exclusively from infragranular layers, and a second class from both infragranular and supragranular layers (Felleman and Van Essen 1991).These ideas have been synthesized into a hypothesis with the pyramidal neuron and its apical tuft dendrites as a key locus for the contextual interaction between feedforward and feedback input (Larkum 2013).According to this hypothesis, when a pyramidal neuron receives strong feedforward input targeting its somatic region, action potentials are triggered that backpropagate into the apical dendrite; if feedback inputs activate the apical tufts at the same time (Takahashi et al. 2020), this leads to long lasting calcium potentials in the apical dendrite and bursts of somatic action potentials (Doron et al. 2020;Guest et al. 2021).
Earlier work with the retrograde tracers fast blue (fb) and diamidino yellow has shown that long-range cortico-cortical projections are an important source of input to L1 of somatosensory cortex (Cauller 1995;Cauller et al. 1998).Single-cell fills of motor cortical and thalamocortical neurons also reveal that long-range axons target L1 (Veinante and Deschênes 2003;

Ethics statement
All experiments were conducted under the license G0278/16 in accordance with the guidelines of animal welfare of the Charité Universitätsmedizin Berlin and the local authorities, the "Landesamt für Gesundheit und Soziales."

Rabies virus production
All viruses for rabies tracing were produced in the Viral Core Facility of the Charité, Universitätsmedizin Berlin (vcf.charite.de).The production was previously described by the Callaway lab (Osakada and Callaway 2013).In brief, we transfected B7GG cells by Lipofectamine 3000 (Thermo Fisher) with the rabies virus genomic vectors that contained either mCherry (pSADdeltaG-F3mcherry, Addgene plasmid #32634), or GFP and Synaptophysin-RFP cDNA (pRVdG-N-P-M-EGFP-SynPhRFP-L, Addgene plasmid #52483), and with these additional plasmids: pcDNA-SADB19N Addgene plasmid #32630, pcDNA-SADB19P Addgene plasmid #32631, pcDNA-SADB19L Addgene plasmid #32632, pcDNA-SADB19G Addgene plasmid #32633.We collected the supernatant over several days and re-transfected the recovered virus in B7GG cells for a final collection step.For additional pseudotyping, the rabies with the envelope protein EnvA of the Avian Sarcoma and Leukosis virus and the obtained supernatant containing non pseudotyped viruses were applied onto BHK-EnvA cells.After 3 to 5 days, EnvA-pseudotyped rabies viruses containing supernatant were collected, filtered, and concentrated by ultracentrifugation.Rabies virus titer was determined by infecting a serial dilution of the virus in HEK293T-TVA cells.All components for rabies virus production were gifts from Edward Callaway (Salk Institute, San Diego, CA, United States).
To allow expression in Cre-positive starter cells, we applied a previously described AAV vector that contained a Cre-dependent expression cassette of a nuclear GFP, EnvA interacting cognate avian viral TVA receptor, and an optimized rabies G protein (pAAV-Syn-Flex-nGToG-WPRE3; Choi and Callaway 2011;Sun et al. 2014;Kim et al. 2016;Zolnik et al. 2020).To distinguish between rabies GFP expression and GFP expression from AAV infected starter cells, we exchanged the nuclear GFP with mCerulean3 (addgene plasmid #54730; Markwardt et al. 2011) to form pAAV-Syn-Flex-Ce3ToG-WPRE3.After sequence verification, AAV was produced at the Viral Core Facility of the Charité-Universitätsmedizin Berlin.

Preparation of mice for surgeries
We deeply anesthetized the mice (n = 42) according to their weight with a mixture of ketamine (100 mg/kg) and xylazine (10 mg/kg).After full anesthesia, confirmed by toe pinch for the absence of ref lexes, mice were placed into a stereotaxic frame with nonpuncture ear bars and a nose clamp (Kopf Stereotaxic device, California, United States, Inc.).For analgesics, we applied 100 μL of lidocaine under the scalp before surgery, and administered Carprofen (5 mg/kg) and Buprenorphine (5 mg/kg) intraperitoneal for analgesic posttreatment after surgery.Body temperature was maintained at 37 • C, using a heating pad during the surgery and while mice recovered from the operation.

Fast blue application
Fast blue (fb) was applied on L1 of S1 cortex, where it is taken up for retrograde transport along the axons (Kuypers et al. 1980;Keizer et al. 1983).We made a craniotomy of ∼ 0.5 mm 2 over S1 cortex (AP: −1.2 to −1.7 mm, lateral: 2.7 to 3.3 mm) and without removing the dura mater applied 2 μL (1% in H 2 O) of the retrograde, blue f luorescent neuronal tracer fb (Polyscience Inc., 17740-1, 1 mg) onto the surface of L1 for 5 min.The f luid that remained on the brain was washed off with PBS.After the procedure, the skull was cleaned, the wound edge sutured, and the mouse was put back into its home cage for recovery.Seven days later, brains were removed and processed for analysis as described below.

Injections of retrobeads
Injections of retrobeads (Lumaf lour) were made in S1 cortex in L1 at a depth of ∼ 150 to 200 microns from pia (200 nl, 20 nl/min).The depth was chosen to be just under L1 in order to take advantage of backf low along the injection track.After the procedure, the skull was cleaned, the wound edge sutured, and the mouse was put back into its home cage for recovery.Seven days later, brains were removed and processed for analysis as described below.

Stereotaxic viral injections
For stereotaxic injections, an incision was made in the scalp, the skull was cleaned with PBS, and a small craniotomy was drilled above the injection location using a dental drill.The brain tissue was kept moist by applying sterile PBS at regular intervals.For injections, we tip-filled pipettes using negative pressure, and injected with pipettes of a tip size of ∼4 μm.Viral injections-200 nl under constant positive pressure with a f low rate of 20 nl/min-were made using the following stereotaxic coordinates for S1 cortex: AP −1.5 mm, lateral: 2.7 mm, ventral: 0.8 mm; for POm: AP: −1.5 mm, lateral: 1.5 mm, ventral: 3.0 mm.
To induce Cre-specific expression of channelrhodopsin, an AAV1-Flex-ChR2-eYFP (1 × 10 13 vg/mL) was injected into S1 of Tlx3-Cre and Sim1-Cre mice.To induce rabies virus expression, we first injected 200 nl AAV-Syn-Flex-nGToG-WPRE3 (8.1 × 10 11 GC/mL) into L5a and L5b and allowed 2 to 3 weeks for expression, followed by injection of 200 nl of RABVSADB19dG-mCherry (1.4 × 10 8 IU/mL).When combining rabies virus injection with fb application, we applied fb as on the day of rabies virus injection and waited 7 days for expression.We analyzed every second section.In rabies virus injections, GFP-positive neurons were distributed from 600 to 1,100 μm across the injection site.

Definition of cell types for recording
We recorded from eYFP-negative neurons in L5 at the site of the viral injection (where expression was maximal).We classified neurons as L5b pyramidal cells if they had a single apical dendrite that extended to L1, their somata were about 800 μm from pia, and they had electrophysiological characteristics of PT neurons, such as large sag current upon hyperpolarization and low input resistance.L5a pyramidal neurons had somata less than about 700 μm from pia, and a thin apical dendrite and small tuft in L1.Neurons with somata in L1 with non-spiny dendrites were classified as L1 interneurons.

Channelrhodopsin-based circuit mapping (sCRACM)
sCRAM was used to test the effect of Tlx3-Cre or Sim1-Cre input to L1 interneurons and to examine connections between L5 pyramidal neurons.We injected AAV1-Flex-ChR2-eYFP (1 × 10 13 vg/mL) into S1 of Tlx3-Cre and Sim1-Cre mice.Following 3 to 5 weeks of expression, we sliced the brains and recorded from neurons in L1 and L5 in the region of the densest ChR2-eYFP label.ChR2 photocurrents and EPSPs were distinguishable based on the response time during optical stimulation.Recordings of synaptic input were done in the presence of tetrodotoxin (TTX) and 4-aminopyradine (4-AP) to ensure that the connections being activated were monosynaptic.

Histology and Nissl stain
Mice were deeply anesthetized with isof lurane (Baxter Inc.), and transcardially perfused with 0.1% PBS for 3 min, and subsequently with 4% paraformaldehyde (PFA) in 0.1 M phosphate buffer (PB, pH 7.2) for 7 min.Brains were carefully removed from the skull and post-fixed in 4% PFA overnight at 4 • C. For sectioning with a vibratome, brains were embedded in 3% agarose; for sectioning with a microtome, brains were soaked in 30% sucrose overnight.All brains were sectioned at 70 μm.Every second section was mounted on slides and cover slipped.The remaining sections were stored in cryo-protect medium (30% Ethylene Glycol, 30% Glycerol in PBS) at −20 • C. We performed nissl stain with Neu-roTrace 500/525 Fluorescent Nissl Stain (Thermo Fisher, N21480) according to the protocol: Brain sections were washed in PBS 3 times, and then incubated in nissl solution (1:500) for 20 min.After the labeling step, brain sections were washed 3 times with PBS and mounted in glycerol (80%, 2.5% DAPCO in PBS).

Definition of application and injection sites
For histological analysis, the application site for fb and the injection sites for retrobeads and rabies were defined by visual examination as the brain section with the densest labeling.Note that fb could spread from the application site at barrel cortex into lateral hindlimb/trunk S1 cortex).For analysis, we used every second brain section and stored alternate sections for additional examination.

Cell counting and definition of cortical areas and cortical layers
Cell bodies were manually counted using the cell counter and analyze particles plugins from Fiji/ImageJ.Labeling directly at the center of the application site can be dense and individual neurons difficult to count; therefore, we counted local-labeled neurons at the center of the application site and at 140 microns on each side of the center of the application site where the fb label was less dense.For retrobeads, counts were made only at the injection site because label of tracer was less dense.
Cortical layers in the brain areas were identified according to the Allen brain reference atlas and in bins of 100/200 microns for each layer.In S1 cortex, layers were identified according to these criteria, in microns from pia: L1 1 to 100 μm, L2/3 100 to 300 μm, L4 300 to 400 μm, L5a 400 to 500 μm, L5b 500 to 700 μm, L6a 700 to 900 μm.Layers for M1 and M2 cortex were identified according to these criteria, in microns from pia: L1 1 to 100 μm, L2/3 100 to 300 μm, L5a 300 to 400 μm, L5b 400 to 600 μm, L6a 600 to 800 μm.For S2 and visual cortices (including V1, V2L), we set the criteria for layers, in microns from pia: L1 1 to 100 μm, L2/3 100 to 250 μm, L4 250 to 350 μm, L5a 350 to 450 μm, L5b 450 to 600 μm, L6a 600 to 800 μm.In all these areas, L6b was also defined as the layer within 100 microns above the white matter.In our scheme, L5a was 100 microns and L5b was 200 microns.Because laminar definitions are to some extent arbitrary (Oberlaender et al. 2011), to examine how our definition of layers affects our results, in a subset of brains, we also counted labeled neurons in a 30-micron strip in the middle of each 100-micron bin.

Statistics
Statistical analyses were carried out with graphpad/Prism.Data are shown as mean ± standard error of the mean (SEM).P-values are indicated in the text: * P = 0.05 to 0.01, * * P = 0.01 to 0.001, * * * P = 0.001 to 0.0001; * * * * P = < 0.0001.For Fig. 3, a Fisher's exact test and Student's t-test were used.For comparison between fb application and retrobead injection, the Kruskal-Wallis test was used.For all other figures, a 1-way ANOVA and Bonferroni post hoc test were used.

Local input to L1 of S1 cortex
To assess input to L1 of primary somatosensory cortex, the retrograde tracer fb was applied on the cortical surface (Cauller et al. 1998;Clancy and Cauller 1999;Keizer et al. 1983).Seven days after application of the tracer, labeled neurons were evident at the application site and in sites distributed throughout the brain (Fig. 1A).The laminar pattern of the local uptake was specific and in agreement with the earlier qualitative observations in rat somatosensory cortex (Cauller et al. 1998;Clancy and Cauller 1999).Because label under the application site can be dense and individually labeled neurons difficult to distinguish from background, we assessed counts both at the center of the application site and 140 microns from the center of the application site.At 140 microns from the center of the application site, we found the highest uptake of fb in L2/3 (40% of all fb-labeled neurons in n = 6 sections from n = 6 mice) and in L5 (L5a: 18%, L5b: 11%).L4 and L6a had a smaller percentage of labeled neurons (L4: 5%, L6a: 1%), and L6b had 13% of fb-labeled neurons.The proportion of total fb-labeled neurons in L1 was 13%.When we compared the counts 140 microns from the center of the application site to counts directly under the center, we found that the laminar pattern of local-labeled neurons was similar (n = 6 mice, 3,096 ± 459 fb-labeled neurons per S1 cortex, 1-way ANOVA * * * * P < 0.0001; Fig. 1B, Supplementary Table 1a, and Supplementary Figs.1a to c and 2a).
The percentage of neurons labeled with fb in each layer was assessed by combining nissl staining with fb application.The number of double-labeled (nissl and fb) neurons in a 30-micron bin in the middle of each defined layer was counted and the laminar pattern assessed.These analyses showed that 45% of all neurons in L2/3, 33% of all neurons in L5a, and 20% of all neurons in L5b were double labeled.Few neurons in L4 (4%) and L6a (0.4%), and 32% of L6b neurons were double labeled.Finally, 43% of neurons in L1 were double labeled (n = 6 mice, n = 3 nissl + fb-stained sections, 1-way ANOVA * * P < 0.01; Fig. 1C and Supplementary Table 1a).The laminar pattern of fb with relation to nissl stain was similar when we used 100 micron bins to define each layer (Supplementary Fig. 1D).Next, we assessed input to L1 with retrobeads injected into cortex.When retrobeads were injected in L1, we obtained a similar laminar pattern as with application of fb on the surface of cortex (n = 3 mice, 2,490 ± 374 retrobead-labeled neurons per S1 cortex; Fig. 1D and E, Supplementary Fig. 2B, and Supplementary Table 1a).Locally, 45% of all retrobead-labeled neurons were in L2/3, followed by 12% of L5 neurons (L5a: 12%, L5b: 11%).We found little uptake of retrobeads in L4 (0.1%), no label in L6a (0.3%), 16% of neurons labeled in L6b, and 15% in L1.We found no significant differences in the laminar distribution of fb-labeled versus retrobead-labeled neurons at the application/injection site (Supplementary Table 1a).From these experiments, we conclude that input to L1 has a specific local laminar pattern that can be assessed by fb or retrobeads, with the majority of local projections arising from L2/3, L5, and L6b.

Long-range input to L1
To estimate the laminar distribution of long-range input to L1 of S1 cortex, we counted and categorized fb-labeled neurons in a variety of cortical areas (Fig. 2, Supplementary Fig. 3, and Supplementary Table 1b).
Fast blue-labeled neurons were also found in subcortical areas, POm, zona incerta, midline thalamic nuclei, basal forebrain, and   brainstem.Subcortically we only counted the number of fblabeled neurons in the thalamus and zona incerta (Fig. 2G).We made no attempt to determine the percentage of labeled neurons in each nucleus (i.e. the ratio of fb-labeled neurons to nissl-stained cells).
When fb was injected into cortex-label was distributed from L1 to L6 at the injection site-the laminar pattern of long-range input to S1 was statistically, significantly different from that observed with fb application on the cortical surface (Supplementary Table 1a and c).The cortical areas containing retrogradely labeled neurons were the same when fb was injected into cortex or applied on the cortical surface.But in contralateral S1, injections of fb revealed more labeled neurons in L6a that were not evident when fb was applied on the cortical surface.Note that whether fb was applied on cortex or injected in cortex, most of the retrogradely labeled neurons in contralateral S1 were in L2/3.In S2, the largest number of fb-labeled neurons after injection of fb was found in L2/3 and L6a (1-way ANOVA, * P < 0.05).In motor cortices, fb neurons were distributed across all layers (Supplementary Fig. 4A to F).These results suggest that longrange projections to L1 consist of a specific set of neurons located in L2/3 and L5.In contrast, long-range projections to deeper parts of S1 engage more neurons in L6a in both contralateral S1 and ipsilateral S2.
The injections of fb also served as a control.Because VPm thalamic axons target L4 and can even target pyramidal cell dendrites in L2/3 of cortex, label in VPm or VPl thalamus is only expected if fb applied on cortical surface seeps into L2/3 or L4, or after injection of fb into cortex (Killackey and Ebner 1973;Lu and Lin 1993;Sermet et al. 2019;Guest et al. 2021).In experiments with fb application on cortical surface, we therefore excluded brains that showed fb label in VPm and VPl because label at these sites would indicate that fb had seeped into L2/3 or L4 at the application site (Supplementary Fig. 1E and F).Application of fb in Scnn1a-Cre mice (Madisen et al. 2010) confirmed that fb was not labeling cells in L4 (Supplementary Fig. 1G).These results confirm that local laminar distribution of fb label was significantly different when fb was applied on L1 compared with when fb was injected into cortex (Supplementary Table 1a and c).

Functional effects of L5 monosynaptic connections in L1
The local pattern of fb distribution in S1 cortex suggests that a substantial input to L1 arises from local L2/3 and L5 neurons.The synaptic connectivity between L2/3 neurons and neurons in L1 has been examined previously (Wozny and Williams 2011), so here we focused on the effect of L5 inputs to L1.We examined the synaptic connectivity between L5 pyramidal neurons and L1 interneurons as well as between L5 pyramidal neurons themselves (via connections within L1).We patched L1 interneurons or L5 pyramidal neurons in Tlx3-Cre or Sim1-Cre mice expressing ChR2.To examine (i) input to L1 interneurons, we photo-stimulated over L1 while recording from L1 interneurons using normal ASCF solution without TTX and 4-AP.To examine (ii) connectivity between L5 pyramidal neurons, in the same experiment, we then photo-stimulated at different points along the somato-dendritic axis of L5 pyramidal neurons in the presence of TTX and 4-AP.This approach blocks di-synaptic transmission (Petreanu et al. 2007).We expected the IT input to generate a larger PSP than the PT input, in part because more IT neurons project to L1.The PSP parameters-amplitude, latency, and short-term plasticity (at 10 Hz)-evoked in interneurons and in L5 pyramidal neurons-were measured but because of the use of TTX and 4-AP, which can alter normal release characteristics, these were not characterized in depth.

Effect of L5 input on L1 interneurons
We recorded from interneurons in L1 in Tlx3-Cre or Sim1-Cre slices while photo-stimulating in L1.This experiment was performed without TTX and 4-AP (Fig. 3A).Photo-stimulation of the axons in L1 showed that IT neurons generate a powerful PSP in L1 interneurons and that this input was significantly stronger than the PT input to the interneurons (Fig. 3B and C; n = 5 neurons, Tlx3 to L1 interneuron amplitude: 25.5 ± 4.9 mV; Sim1 to L1 interneuron amplitude: 0.9 ± 0.2 mV, * * * * P < 0.0001, Student's t-test).Also, note that 5 out of 5 interneurons received input from IT, whereas only 2 out of 5 interneurons received input from PT.The input from IT to L1 interneurons had depressing synaptic dynamics at 10 Hz stimulation (Fig. 3D).

IT-PT interactions with other L5 pyramidal neurons
Next, we examined whether IT neurons could drive non-Tlx3 neurons in L5b, and whether PT neurons could drive non-Sim1 neurons in L5a (Fig. 3E and F, schematic top).In order to restrict the transmitter release to the photo-stimulated area, and to prevent di-synaptic activation, here we applied TTX and 4-AP to the bath solution.The entire somato-dendritic axis of the filled neurons was systematically photo-stimulated by moving the optogenetic light in 50 μm steps, while simultaneously measuring the amplitude of the evoked EPSP.Activation of the Tlx3-Cre axons did not evoke any EPSPs in L5b tufts in L1.Similarly, activation of Sim1-Cre axons did not evoke EPSPs in non-Sim-L5a tufts in L1 (Fig. 3E, bottom graph).The Tlx3 neurons did, however, evoke a significant depolarization (while photostimulating) at the oblique dendrites and to a lesser extent at the basal dendrites and soma in L5b (14.97 ± 1.56 mV).The Sim1-Cre axons, on the other hand, had a negligible effect on non-Sim L5a neurons across the entire somato-dendritic axis (2.30 ± 0.46 mV; Fig. 3F, bottom graph; Ramaswamy et al. 2012).
In summary, we show that the L5 Tlx3 IT neurons strongly drive interneurons in L1 but have a negligible effect on the apical tuft dendrites in L1 of non-Tlx3 L5b neurons.Layer 5 Sim1 PT neurons do not drive interneurons in L1 or the tuft dendrites of L5 non-Sim1 neurons in L1.This suggests that 2 major classes of L5 pyramidal neurons do not interact with each other directly within L1.

Glutamate uncaging
To confirm that the input to the apical tuft dendrite could evoke a synaptic response in our experimental conditions, we performed glutamate uncaging along the dendritic axis while recording from neurons in L5a and L5b (Supplementary Fig. 5A and B).Our results indicate that the effect of localized input at the tuft dendrite can be measured at the soma.Glutamate uncaging evokes measurable EPSPs in both cell types at the soma, when stimulating a local spot in the apical tuft dendrite for 30 ms in the presence of TTX/4-AP.Our work and previous work (Hooks et al. 2011;Passlick and Ellis-Davies 2019) indicate that our methods using sCRACM are potentially sensitive enough to detect inputs to the apical tuft dendrites of L5 cells.

Laminar profile of L5, L6b, and inhibitory neurons
We first assessed the percentages of layer-specific cell types for each Ai9-reporter line (Supplementary Fig. 6A and B and Supplementary Table 2a).Tlx3-Cre neurons were exclusively found in L5, primarily in L5a.They constitute 42% of nissl-stained neurons in L5, with 76% in L5a, and the remaining 24% in L5b.By comparison, Sim1-Cre neurons were distributed from L4 to L6, with 54% in L5b, and 28% in L5a and the rest (18%) in other layers.Overall, Sim1-Cre neurons constituted 18% of all neurons in L5.Note that although tdTom-positive neurons are found across L5 in both Cre-lines Tlx3-Cre and Sim1-Cre, the majority of tdTompositive neurons for Tlx3 were in L5a, and the majority tdTompositive neurons in Sim1 were in L5b.Drd1-positive neurons were found in L6a and L6b, with the majority in L6b.
All Ctgf neurons were located in L6b, but Ctgf neurons constituted only 15% of all neurons in L6b.By comparison, Drd1apositive neurons were distributed in L6a and L6b, with 76% of all Drd1a neurons in L6b.Drd1a-positive neurons were 43% of all neurons L6b (Supplementary Fig. 6C and D).
High magnification images of Sim1-Cre-Ai9 show fb/tdTom processes and neurons in L1 and L5b.These images show that adjacent tdTom-positive neurons that have dendrites extending toward L1 do not necessarily take up fb, i.e. their somas are not uniformly fb positive (Supplementary Fig. 6H).These images argue for a specific retrograde uptake of fb.

Local input from classes of L5, L6b, and inhibitory neurons to L1
Local projections from layers 5, 6b, and from inhibitory neurons to L1.To examine the IT and PT connectivity to L1, we applied fb on L1 in Tlx3-Cre-Ai9 and in Sim1-Cre-Ai9 mice and counted doublelabeled neurons in those lines (Fig. 4A and B, Supplementary Fig. 6I, and Supplementary Table 2c and e).In ipsilateral S1, 29% of all IT neurons in L5 were fb labeled, with the majority located in L5a.The PT neurons projecting to L1 showed a different profile, with 20% of all PT neurons double labeled for fb ipsilaterally, with the majority of L1 projecting neurons located in L5b.In L6b, a larger proportion, 52%, of Ctgf neurons projected to L1, whereas only 13% of Drd1a neurons projected to L1 (Fig. 4C and D and Supplementary Table 2c).Fast blue-labeled Sst neurons were most abundant in L2/3 and L5a, fb-labeled Vip neurons were most abundant in L2/3, and fb-labeled Pv neurons were most abundant in L5b (Fig. 4E to G and Supplementary Table 2d).

Cell-type-specific long-range projections to S1 L1
IT and PT neurons in contralateral S1, ipsilateral motor cortices, and S2 cortex were labeled with fb (Fig. 4H and I and Supplementary Table 2c).For contralateral S1, 8% of IT neurons and only 1% of PT neurons were labeled.The percentage of double-labeled neurons in motor cortices was similar for IT and PT neurons (∼10%).In S2, 11% of IT neurons, and 8% of PT neurons were double labeled.Thus, IT neurons project contralaterally to L1, and slightly more IT than PT neurons in motor and S2 cortices have long-range projections to L1.Overall, 1.5% of the IT neurons in perirhinal cortex project to contralateral S1 L1.No Ctgf or Drd1apositive fb double-labeled neurons were found in the contralateral hemisphere (Fig. 4K).In Sst projecting inhibitory neurons, there were also a small number of long-range in contralateral S1, in L2/3 and L5a that targeted L1 (0.3%, 15 out of 6,173 neurons, 3 brains).

Local input to L1 and L5 neurons assessed with fb and rabies virus
The experiments so far revealed that subtypes of L5 neurons project in a specific manner to L1.These observations raise a question: Are neurons that project to L1 presynaptic to the L5 pyramidal neurons?To address this question, we combined fb application in L1 with rabies-based retrograde monosynaptic tracing for the IT and PT neurons.Monosynaptic rabies tracing labels neurons connected by one synapse to the source (starter) neurons (Wickersham et al. 2007;Luo et al. 2008;Kim et al. 2016).

Distribution of input along the somato-dendritic axis
To first examine whether the rabies virus targets synapses along the entire soma-dendritic arbor of the starter neuroni.e. to examine whether boutons of presynaptic neurons could be distributed in L1 to L6-we used a modified rabies virus linked to synaptophysin (Wickersham et al. 2013).With this approach, presynaptic cell bodies and their axons/boutons were labeled with synaptophysin linked to RFP (presynaptic) and the starter population was additionally labeled with cytosolic GFP (postsynaptic) and nuclear cerulean f luorescent protein (CFP; Supplementary Fig. 7A).We examined the distribution of boutons and presynaptic (RFP) label in a circuit where it was known that a class of presynaptic neurons target specific dendritic loci.For this assay, we used the Gpr26-Cre line, which expresses Cre in a class of CA1 pyramidal neurons (Supplementary Fig. 7B).Entorhinal neurons target the apical tuft dendrites of CA1 pyramidal neurons (Masurkar et al. 2017), so finding presynaptic, RFP-labeled neurons in entorhinal cortex would indicate that the rabies strategy is effective even for inputs that target the tuft dendrites.When we examined entorhinal cortex, we found presynaptic, RFP-labeled neurons, suggesting that rabies retrograde approach can work for synapses at the tuft dendrites (Supplementary Fig. 7C).
Next, we used the same approach in S1 cortex of Sim1-Cre and Tlx3-Cre mice (Supplementary Fig. 7D to F). Starter cells expressing CFP were found in L5, and GFP expressing cells and RFP-positive boutons were detected in all layers, including L1, across the cortical column (Tlx3-Cre, n = 369 synapses, Sim1-Cre, n = 158 synapses, bouton counts from 10 boxes for each layer with a size of 20 × 10 μm).Presynaptic neurons were also detected in motor cortex, S2, and thalamus (Supplementary Fig. 7F).These experiments suggest that the rabies approach targets the entire somato-dendritic compartment-including synapses on apical tuft dendrites in L1-of the starter pyramidal neurons.

L1 projecting neurons, presynaptic to PT and IT neurons
With these results in hand, we performed injections in Tlx3-Cre and Sim1-Cre mice with the standard AAV-rabies virus approach in combination with fb application to see whether neurons that were presynaptic to L5 neurons also project to L1 (Fig. 5A and Supplementary Fig. 8A and B).In Tlx3-Cre mice, fb application in L1-labeled 9,596 neurons, mCherry was found in 6,138 cortical neurons presynaptic to the Tlx3 neurons.There were 449 neurons that were double labeled, i.e. were presynaptic to Tlx3 neurons and projected to L1.In Sim1-Cre mice, there were 13,131 fb labeled, 6,495 presynaptic neurons, and 823 double-labeled neurons, i.e. those that were presynaptic to Sim1-Cre neurons and projected to S1 L1.The percentage of fb+/presynaptic was 7.5% in Tlx3 and 9.0% in Sim1; the percentage of fb+/starter was 0.9% in Tlx3 and 5.2% in Sim1 (Fig. 5B and Supplementary Table 3a  and b).

Local intracortical pattern
Total presynaptic input was 6,138 neurons in Tlx3 mice and 6,495 neurons in Sim1 mice.Example images for the injection site in Tlx3-Cre and Sim1-Cre brains show starter cells in yellow and mCherry + fb double-labeled cells in cyan.Locally, in Tlx3-Cre brains, the presynaptic-labeled neurons were mainly seen in L2/3 and L5 (1-way ANOVA, * * * P > 0.001).Neurons that were double labeled, i.e. were presynaptic to Tlx3 neurons and labeled with fb, were most common in L2/3 and L5a and L5b (Fig. 5C).In Sim1-Cre brains, most of the local presynaptic input was from L2/3 and L5b (Fig. 5D) and double-labeled neurons were detected mainly in L2/3 and L5b.

A control for fb, when combined with rabies approach
We obtained similar results for fb application on L1 in the transgenic lines as we do with fb application in wild-type mice (Supplementary Fig. 8C to F).To confirm whether fb application on L1 interferes with rabies, we applied rabies virus to Tlx3-Cre and Sim1-Cre brains without fb application (Supplementary Fig. 9A to D and Supplementary Table 3c).The pattern of fb distribution was similar in both approaches.In Tlx3-Cre brains, presynaptic neurons were mainly found in L2/3 and L5a; in Sim1-Cre brains, presynaptic neurons were mainly found in L2/3 and L5b.

Long-range input to L1 and L5 neurons
Earlier work has shown that in addition to targeting L1, longrange input from cortical areas like M1 targets infragranular layers (Zagha et al. 2013;Kinnischtzke et al. 2016;Yamawaki et al. 2021;Geng et al. 2022).This raises the question whether neurons that project to L1 from other brain areas are presynaptic to the L5 pyramidal neurons?Long-range presynaptic input to Tlx3-Cre and Sim1-Cre neurons in S1 was primarily from motor cortices, S2 cortex, and thalamus (Fig. 6 and Supplementary Table 3a and  b).In both, Tlx3 and Sim1 brains, long-range input arose from neurons in L2/3 and L5 and double-labeled neurons-i.e.neurons that were presynaptic to L5 cells and projected to L1 of S1were most commonly found in L2/3 (1-way ANOVA, * * * * P > 0.0001).We focused our count of presynaptic neurons on the areas that provide input to L1.In S2 of Tlx3-Cre brains, presynaptic input was from L2/3 and L5a neurons (with a smaller input from L6a and L6b), and 49 neurons were double labeled for fb and mCherry (Fig. 6A).None of the L6 neurons were double labeled.
L2/3 and L5a neurons in motor cortex provided most of the presynaptic input to S1 Tlx3 neurons, and few were double labeled (Fig. 6B).Some presynaptic neurons were found in L6a and none were double labeled.
In S2 of Sim1-Cre brains, L2/3 and L5b neurons generated the bulk of the presynaptic input and 31 neurons were double labeled (Fig. 6D).Presynaptic neurons in motor cortex were found in L2/3, L5a, and L5b, but only 13 were double labeled (Fig. 6E).Presynaptic neurons to both L5 lines were found in thalamus with slightly more input to Tlx3 neurons than to Sim1 neurons.Double-labeled neurons accounted for a small percentage of labeled thalamic neurons in both lines (Tlx3: 1.8%, Sim1: 4.1%, Fig. 6C and F).

Local versus long-range input to L1
Taken together these results suggest that the bulk of the presynaptic input to Tlx3 and Sim1 neurons was from local neurons (Tlx3: 80%, Sim1: 63%).A large portion of the input targeting L1 that was presynaptic to L5 neurons arose from local neurons.A fraction of the neurons presynaptic to L5 pyramidal neurons also projected to L1. Long-range cortical input to these pyramidal cells was from S2, motor, and visual cortices (Fig. 6G and H and Supplementary Table 3d).
When we extended the approach of counting fb-labeled neurons to the whole brain, and to assess the input from the different brain areas, we found that the local input to L1 of S1 arises predominantly from neurons underlying the application site (68%), with long-range input constituting 32% of the total input (n = 4, 11,956 ± 2,469 neurons per brain; Fig. 7A and Supplementary Table 4).The long-range input to L1 was from higher-order thalamus (7%), ipsilateral motor cortices (9%), S2 cortex (2%), visual cortices (3%), and perirhinal cortex (0.1% ipsi and 0.4% contra).Contralateral S1, motor, and perirhinal cortex also contributed input to L1 of somatosensory cortex.We did not quantify input from brainstem to L1.We obtained a similar input pattern when injecting retrobeads in L1: 67% was local, and 33% was long-range input (n = 3, 3,487 ± 868 neurons per brain; Fig. 7A and Supplementary Table 4).

Discussion
Layer 1 has been dubbed an enigma, a "crowning mystery" (Hubel 1982;Rudy et al. 2011) with its anatomy, connectivity, and function still barely charted.By combining traditional retrograde tracing, we quantitatively characterized aspects of organization and function of this layer in primary somatosensory cortex.Our work reveals the extent to which multiple classes of local and longrange neurons connect to L1, and it shows the functional consequences of L5 input to L1.

Input to L1 is primarily local and can modulate activity of L1 interneurons
It is widely accepted that the majority of synapses in any patch of cortex arise from cortical neurons with thalamic inputs, contributing ∼15% of all synapses in cortex, even in L4 (White and Keller 1989;Binzegger et al. 2004;Douglas and Martin 2004;Shepherd 2004).An earlier, detailed analysis of axons and boutons of 39 filled excitatory and inhibitory neurons of all layers of the cat visual cortex suggested that non-local sources could provide a majority of the input-"dark synapses"-to L1 (Binzegger et al. 2004;Douglas andMartin 2004, 2007;Voges et al. 2010;Boucsein et al. 2011).We do not directly address the origin of these mysterious boutons, but along with earlier studies that suggest that the bulk of input to L1 interneurons is from a local origin, our work suggests that local inputs contribute most of the input to L1 (Manns et al. 2004;Brown and Hestrin 2009;Wozny and Williams 2011;Narayanan et al. 2015;Sakmann 2017;Abs et al. 2018;Peng et al. 2021).An additional feature of cortical organization is a substantial recurrent local excitatory and inhibitory connectivity (Sachdev et al. 2012;Sherman and Usrey 2021).L1 neurons could, via inhibitory interactions in L1 disinhibit pyramidal neurons and interneurons in L2/3 (Anastasiades et al. 2021;Cohen-Kashi Malina et al. 2021).
Our work suggests that local recurrent inputs are important in L1 as well.Neurons with dendrites in L1-i.e.L2/3 and L5 pyramidal neurons or L2/3 VIP/SST neurons-and neurons with axons in L1 contribute local recurrent excitation and inhibition in this layer (Larkum et al. 2018;Wozny and Williams 2011).Using optogenetic approaches, we show that input to L1 from L5 IT pyramidal neurons has specific functional consequences.L5 IT neurons effectively activate L1 interneurons but do not have a measurable effect on the apical tuft dendrites of L5 PT pyramidal neurons.The negligible effect of L5 IT pyramidal neurons on apical dendrites of L5 PT neurons could ref lect weak connectivity between these cell types in L1, or a specific lack of connectivity between the 2 classes of neurons.L5 pyramidal neurons do not seem to generate local recurrent excitation to other L5 pyramidal neurons via their inputs to L1.It has been shown earlier that uncaging in dendrites can elicit measurable input to the soma (Dodt et al. 1998;Pettit and Augustine 2000;Frick et al. 2001;Shoham et al. 2005;Harris and Pettit 2007), results that are consistent with ours.Taken together, our data suggest that circuit mechanisms exist for driving local recurrent excitation and inhibition through synapses in L1 (Wozny and Williams 2011;

Long-range input to L1 is primarily from L2/3 and L5 and is cell type specific
Our work shows that the principal source of long-range input to L1 is from L2/3 neurons, with a smaller but significant input from L5 cortical pyramidal cells.While neurons in contralateral S1, M1, and perirhinal cortices, and ipsilateral S2, M1, M2, perirhinal and visual cortices provide long-range input to L1, the bulk of the long-range input to L1 is from S2, M1, M2, visual cortices, and thalamus (Fig. 7B).Our work also shows that 10% to 20% of the PT and IT neurons in S2 and M1 and contralateral S1 project to L1.While L6 cortico-cortical neurons provide some intracortical feedback, their input to L1 is negligible.In contrast to the results with fb applied on the cortical surface, when tracer was injected into S1 cortex, an additional class of cortico-cortical projection neurons in L6 was evident, data that relate to previous studies on hierarchical models with retrograde tracers (Markov et al. 2014;Vezoli et al. 2021).4.

Few neurons that are presynaptic to L5 pyramidal neurons project to L1
Previous work has shown that motor cortices, POm thalamus, and perirhinal cortex project to L1 and have substantial input to L5 (Aronoff et al. 2010;Avermann et al. 2012;Chen et al. 2013Chen et al. , 2015Chen et al. , 2016;;Zagha et al. 2013;van der Bourg et al. 2017;Yamashita et al. 2018;Doron et al. 2020;Geng et al. 2022).Our experiments with rabies virus show that input from higher-order somatosensory cortex, S2, to L5 PT and IT neurons constitutes more than 40% of the total long-range input to these neurons.Additionally, the rabies approach shows that some L6a and L6b neurons in S2 and M1 are presynaptic to the IT neurons.Nevertheless, even for S2, only a fraction of neurons that were presynaptic to PT and IT neurons also projected to L1.Note that Tlx3 and Sim1 neurons are only subclasses of all L5 neurons; consequently these connectivity numbers are likely to be underestimates.

Input to L1 in other cortices
The approach we have taken here to study input to L1 in S1 cortex has also been successfully applied to other cortices such as motor/prefrontal cortices (Herkenham 1979;Arbuthnott et al. 1990;Geng et al. 2022).Tracing of axons of filled thalamic neurons confirms that just as in S1 cortex where higher-order POm axons target L1 (Wimmer et al. 2010;Sermet et al. 2019), higher-order thalamic VM neurons target cortical L1 (Kuramoto et al. 2009(Kuramoto et al. , 2015(Kuramoto et al. , 2017) ) and can affect movement initiation (Takahashi et al. 2021).

Methodological considerations
Why use fb applied on the cortical surface?
Traditional anatomical approaches for assessing connectivity between cortical areas use retrograde tracers injected into the brain.More recently AAV-retro and pseudo-rabies approaches have been developed targeting Cre-expressing neurons.But input to L1 cannot be assessed with beads, cholera toxin, fb, or AAVretro when injected into cortex.The rabies approach is specific and can target specific classes of neurons, i.e.L1 interneurons or pyramidal cells with dendrites in L1, but there is no approach that targets pyramidal cell dendrites exclusively.Approaches targeting L1 interneurons can be specific (Abs et al. 2018), but only reveal input to a tiny fraction of all interneurons in cortex and even in L1, only to 60% of interneurons in L1 (Rudy et al. 2011;Schuman et al. 2019).
The use of fb and related tracers applied on the cortical surface was established ∼30 years ago by Cauller and colleagues (Cauller et al. 1998;Clancy and Cauller 1999;Keizer et al. 1983;Kuypers et al. 1980).This work showed long-range connectivity to L1 and revealed that local L6b and L5 connect to L1 of somatosensory cortex in rats.A concern of retrograde tracing with fb applied on the cortical surface is that fb could spread into L2/3 or even L4.To mitigate against this, we discarded brains in which fb label is found in VPm or VPl.VPm neurons project to L4 and can target the apical tufts specifically in L2/3 (Killackey and Ebner 1973;Wimmer et al. 2010;Feldmeyer 2012;Oberlaender et al. 2012;Sermet et al. 2019;Guest et al. 2021).This approach mitigates against but does not completely eliminate connections that might target L2/3.
Another concern with fb application is the potential for nonspecific labeling of neurons directly at the center of the application site, i.e. too much background or potential uptake of fb by dendrites of pyramidal neurons.To address this issue, we used four approaches: (i) We injected a different class of retrograde tracers-retrobeads-and compared the results we obtained with retrobeads to those with fb.When we used retrobeads, we observe sparser label locally and in distant cortical and subcortical structures, but the overall pattern of label with both retrobeads and fb is similar; (ii) We considered the functional consequences of local input to L1.Earlier physiology with single-cell fills of pyramidal neurons indicated that L2/3 axons were pervasive in L1 and generated substantial input to L1 (Wozny and Williams 2011).As with L2/3, single-cell fills revealed that L5a neurons have dense projections to L1, but while L5b axons are present in L1, the local axon is sparser (Sakmann 2017).In our study, we tested the physiological effect of these local axons on both L1 interneurons and L5 pyramidal cells and show a specific functional effect of L5a input to L1 interneurons.The interaction between L5 cells, their input to L5 dendrites, in L1 was negligible; (ii) We counted both at center of the application site where there was dense fb label, where it was possible to overestimate labeled neurons, and 140 microns off center of the application site.The overall pattern of label was similar in both cases; (iv) We injected fb into cortex and found a different laminar pattern of both local and long-range label.
One additional concern with all tracing approaches is the uptake of tracer by fibers of passage.While this is possible, it should be noted that all tracing methods seem to show a similar pattern of labeling, with neurons in the same cortical and subcortical areas labeled.Our labeling with fb is consistent from mice to mice, generating a similar laminar profile for inputs.Uptake by fibers of passage might have generated a less consistent pattern of label.

Definition of cortical layers
While L1 and L6b are easily defined, the borders that separate the other layers from each other are not easily identifiable.We defined layers using 100/200-micron bins.This approach has the benefit of consistency and reproducibility, but there can be errors arising from brain shrinkage, brain size, or selection of sections used for counting.To mitigate against this, we examined the proportion of fb-labeled neurons in a 30-micron strip at the center of the 100/200-micron bin defining each layer.The proportion of fb-labeled neurons in L4 changed with this approach, suggesting that some of the L4 label could be attributed to L2/3 or L5.

Rabies virus tracing
The incubation time we used for rabies tracing (7 days) might have been too short to obtain expression in all long-range neurons (Kim et al. 2016).Furthermore, even though we have ruled out leakage of rabies virus and found little nonspecific expression (Zolnik et al. 2020), the rabies virus approach could be marginally nonspecific.

Implication
The nub of the mystery concerning L1 can be reduced to 2 questions: (i) How is the input to this layer organized and (ii) What does the input do?Our work shows that L1 has some features that are common to neocortical organization.Most input to this layer is local, arising from excitatory and inhibitory neurons in the patch of cortex underlying L1.This implies that when feedforward, feedback, or both inputs drive somatic action potentials in L2/3 or L5, this activity is likely to directly generate synaptic activity in L1, thus likely to modulate the effect of a backpropagating action potentials that interact with contextual long-range input.
Our characterization of anatomical and functional input to L1 in S1 cortex adds to the growing knowledge and reveals key organizational principles for this circuit.Our work also points to a need to develop tools to further understanding cortical organization and function.Overall, our work suggests that local input to L1 may play a critical role in the information f lowing through the underlying cortex.

Fig. 1 .
Fig. 1.Local input to L1. A) Schematic of experiment.Retrograde tracer fb was applied on L1 of S1 cortex, and fb-labeled neurons in S1 cortex were counted 7 days after incubation.B) Uptake of fb was evident across all cortical layers, with prominent uptake in L2/3 and L5.The tabulation of fblabeled neurons in different lamina under the application site is shown in the bar graph.C) Percentage of L1 projecting fb-labeled neurons in each layer.Nissl-stained (green) and double-labeled (fb and nissl) neurons were counted to obtain a percentage of labeled neurons in each lamina (30-micron bin, red dotted line).L2/3 and L5 had the highest percentage of double-labeled neurons.D, E) Injection of retrobeads in L1 shows a similar pattern of local and long-range input to L1 as with fb application.Total number of neurons counts are shown in each panel.Statistical analysis with 1-way ANOVA, Bonferroni post hoc test, * * * * P < 0.0001.Total number of neurons counted and mice used (in brackets) are shown in each panel, analysis details in Supplementary Table 1a.Scale bars in A, e 500 μm, in B to D 100 μm.

Fig. 2 .
Fig. 2. Long-range input to S1 L1 from other cortical areas.A to E) Percentage of total fb-labeled neurons (left), and percentage of labeled neurons in relation to nissl in each layer (right).A) S1 contralateral cortex, B) M2 cortex, C) M1 cortex, D) S2 cortex, E) visual (V1, V2L) cortices, F) perirhinal cortices ipsilateral and contralateral.Counts of fb-labeled neurons ipsilateral to application site (left) and contralateral (right).Quantification of the laminar pattern of long-range input to L1 shows the highest percentage of fb-labeled neurons in cortical L2/3 and L5a.G) Fast blue label (total numbers) in POm thalamus, zona incerta, and midline thalamus shows that these 3 subcortical nuclei have projections to L1. Abbreviations in g: MD, midline; POm, posteromedial complex of thalamus; VPl, ventrolateral nucleus of thalamus; VPm, ventromedial nucleus of thalamus; Zi, zona incerta.Each dot in the graphs represents one brain section.Data are from 4 brains.Statistical analysis with 1-way ANOVA, Bonferroni post hoc test, * * * * P < 0.0001.Analysis details in Supplementary Table1b.Scale bars in A to F 100 μm, in G 500 μm, in zoom-ins 50 μm.
Fig. 2. Long-range input to S1 L1 from other cortical areas.A to E) Percentage of total fb-labeled neurons (left), and percentage of labeled neurons in relation to nissl in each layer (right).A) S1 contralateral cortex, B) M2 cortex, C) M1 cortex, D) S2 cortex, E) visual (V1, V2L) cortices, F) perirhinal cortices ipsilateral and contralateral.Counts of fb-labeled neurons ipsilateral to application site (left) and contralateral (right).Quantification of the laminar pattern of long-range input to L1 shows the highest percentage of fb-labeled neurons in cortical L2/3 and L5a.G) Fast blue label (total numbers) in POm thalamus, zona incerta, and midline thalamus shows that these 3 subcortical nuclei have projections to L1. Abbreviations in g: MD, midline; POm, posteromedial complex of thalamus; VPl, ventrolateral nucleus of thalamus; VPm, ventromedial nucleus of thalamus; Zi, zona incerta.Each dot in the graphs represents one brain section.Data are from 4 brains.Statistical analysis with 1-way ANOVA, Bonferroni post hoc test, * * * * P < 0.0001.Analysis details in Supplementary Table1b.Scale bars in A to F 100 μm, in G 500 μm, in zoom-ins 50 μm.
Fig. 2. Long-range input to S1 L1 from other cortical areas.A to E) Percentage of total fb-labeled neurons (left), and percentage of labeled neurons in relation to nissl in each layer (right).A) S1 contralateral cortex, B) M2 cortex, C) M1 cortex, D) S2 cortex, E) visual (V1, V2L) cortices, F) perirhinal cortices ipsilateral and contralateral.Counts of fb-labeled neurons ipsilateral to application site (left) and contralateral (right).Quantification of the laminar pattern of long-range input to L1 shows the highest percentage of fb-labeled neurons in cortical L2/3 and L5a.G) Fast blue label (total numbers) in POm thalamus, zona incerta, and midline thalamus shows that these 3 subcortical nuclei have projections to L1. Abbreviations in g: MD, midline; POm, posteromedial complex of thalamus; VPl, ventrolateral nucleus of thalamus; VPm, ventromedial nucleus of thalamus; Zi, zona incerta.Each dot in the graphs represents one brain section.Data are from 4 brains.Statistical analysis with 1-way ANOVA, Bonferroni post hoc test, * * * * P < 0.0001.Analysis details in Supplementary Table1b.Scale bars in A to F 100 μm, in G 500 μm, in zoom-ins 50 μm.

Fig. 4 .
Fig. 4. Local and long-range input from classes of L5 and L6b, and inhibitory cells to L1.A to G) Example images of fb uptake in Tlx3-Cre, Sim1-Cre, Ctgf-Cre, Drd1a-Cre, Sst-Cre, Vip-Cre, and Pv-Cre brains.Ipsilaterally, 29% of Tlx3 S1 project to L1.Twenty percent of Sim1 neurons project to L1. Fiftytwo percent of L6b Ctgf-positive neurons project to L1.Only 13% of L6b Drd1a neurons project to L1. Sst and Vip neurons in all layers take up fb on the ipsilateral side.More than 80% of L2/3 Vip neurons take up fb.Pv neurons take up fb in layers 2/3 and 5. H) Example images of long-range fb uptake in Tlx3-Cre brain.Contralaterally, 8% of Tlx3 neurons in S1 project to L1. Tlx3 neurons in motor and S2 cortices also project to L1.Quantification shown in graphs.I) Example images of long-range fb uptake in Sim1-Cre mice.Contralaterally, few Sim1 neurons in S1 project to L1. Sim1 neurons in motor cortices and S2 cortices also project to L1.Quantification shown in graphs.K) No Ctgf or Drd1 neurons contralaterally of S1 project to L1.Each dot in the graphs represents one brain section.Total number of neurons counted in each mouse line and mice used (in brackets) are shown in each panel, analysis details in Supplementary Table2cand d.Fast blue pseudo colored in cyan.Scale bars 50 μm.
Fig. 4. Local and long-range input from classes of L5 and L6b, and inhibitory cells to L1.A to G) Example images of fb uptake in Tlx3-Cre, Sim1-Cre, Ctgf-Cre, Drd1a-Cre, Sst-Cre, Vip-Cre, and Pv-Cre brains.Ipsilaterally, 29% of Tlx3 S1 project to L1.Twenty percent of Sim1 neurons project to L1. Fiftytwo percent of L6b Ctgf-positive neurons project to L1.Only 13% of L6b Drd1a neurons project to L1. Sst and Vip neurons in all layers take up fb on the ipsilateral side.More than 80% of L2/3 Vip neurons take up fb.Pv neurons take up fb in layers 2/3 and 5. H) Example images of long-range fb uptake in Tlx3-Cre brain.Contralaterally, 8% of Tlx3 neurons in S1 project to L1. Tlx3 neurons in motor and S2 cortices also project to L1.Quantification shown in graphs.I) Example images of long-range fb uptake in Sim1-Cre mice.Contralaterally, few Sim1 neurons in S1 project to L1. Sim1 neurons in motor cortices and S2 cortices also project to L1.Quantification shown in graphs.K) No Ctgf or Drd1 neurons contralaterally of S1 project to L1.Each dot in the graphs represents one brain section.Total number of neurons counted in each mouse line and mice used (in brackets) are shown in each panel, analysis details in Supplementary Table2cand d.Fast blue pseudo colored in cyan.Scale bars 50 μm.
Fig. 4. Local and long-range input from classes of L5 and L6b, and inhibitory cells to L1.A to G) Example images of fb uptake in Tlx3-Cre, Sim1-Cre, Ctgf-Cre, Drd1a-Cre, Sst-Cre, Vip-Cre, and Pv-Cre brains.Ipsilaterally, 29% of Tlx3 S1 project to L1.Twenty percent of Sim1 neurons project to L1. Fiftytwo percent of L6b Ctgf-positive neurons project to L1.Only 13% of L6b Drd1a neurons project to L1. Sst and Vip neurons in all layers take up fb on the ipsilateral side.More than 80% of L2/3 Vip neurons take up fb.Pv neurons take up fb in layers 2/3 and 5. H) Example images of long-range fb uptake in Tlx3-Cre brain.Contralaterally, 8% of Tlx3 neurons in S1 project to L1. Tlx3 neurons in motor and S2 cortices also project to L1.Quantification shown in graphs.I) Example images of long-range fb uptake in Sim1-Cre mice.Contralaterally, few Sim1 neurons in S1 project to L1. Sim1 neurons in motor cortices and S2 cortices also project to L1.Quantification shown in graphs.K) No Ctgf or Drd1 neurons contralaterally of S1 project to L1.Each dot in the graphs represents one brain section.Total number of neurons counted in each mouse line and mice used (in brackets) are shown in each panel, analysis details in Supplementary Table2cand d.Fast blue pseudo colored in cyan.Scale bars 50 μm.

Fig. 5 .
Fig. 5. Local synaptic input to S1 L1 in Tlx3 and Sim1 neurons in combination with fb.A) Schematic of injection scheme.Fourteen days after the first injection with AAV-Syn-f lex-nGToG-EGFP (step 1), fb was applied on L1 in parallel to the injection of rabies virus (steps 2 and 3).Both, rabies virus and fb were allowed to express for 7 days.B) Number of starter, mCherry (presynaptic), fb-labeled and double-labeled neurons (mCherry + fb), and percentages of double labeled/presynaptic, and double labeled/starter in single mice.In Tlx3-Cre mice, the presynaptic neurons derived from 479 starter neurons (out of 706 GFP expressing neurons); in Sim1-Cre mice, the presynaptic neurons derived from 292 starter neurons (out of 643 GFP expressing neurons).C, D) Example images of virus expression and fb uptake in S1 cortex of C) Tlx3-Cre and D) Sim1-Cre brains.Insets showing starter neurons, presynaptic neurons, and fb double-labeled neurons.Percentages of presynaptic neurons in S1 cortex and double-labeled fb and rabies, neurons in Tlx3-Cre, and Sim1-Cre brains in graphs.Each dot in the graphs represents one brain section.Total number of neurons counted are shown in each panel, data from 3 mice each genotype, data shown as mean ± SEM.Statistical analysis with 1-way ANOVA, Bonferroni post hoc test, * * * * P > 0.0001.Analysis details in Supplementary Table3aand b.Scale bars in C, D 100 μm, in zoom-ins 50 μm.

Fig. 6 .
Fig. 6.Long-range synaptic input to S1 L1 in Tlx3 and Sim1 neurons in combination with fb.A to C) Images and quantification of presynaptic neurons (rabies and fb double labeled) for Tlx3 brains in A) S2 cortex ipsilateral, in B) motor cortices, and in C) thalamus.D to F) Images and quantification of presynaptic neurons (rabies and fb double labeled) for Sim1 brains, in D) S2 cortex ipsilateral, in E) motor cortices ipsilateral, in F) thalamus.G) Pie charts showing presynaptic input to S1 for both IT neurons (left) and L1 projecting neurons that are presynaptic to L5 neurons (right).H) Pie charts showing presynaptic input to S1 for both PT neurons (left) and L1 projecting neurons that are presynaptic to L5 neurons (right).Most input to IT and PT neurons was local.The long-range input to S1 L1 was from S2, visual (V1 and V2L), M1, and thalamus.For IT neurons, the bulk of the neurons presynaptic to these cells that also projected to L1 were local neurons (89%).Presynaptic input to PT neurons in S1 was divided into 63% local and 37% long-range.The long-range input was from S2, visual (V1 and V2L), and M1.Presynaptic input to PT neurons that also targeted L1 arose from local neurons (93% of the total input).Total number of neurons counted in each mouse line, data from 3 brains each genotype, shown as mean ± SEM.Each dot in the graphs (in A to F) represents one brain section.Abbreviations in a: I, ipsilateral; c, contralateral.Statistical analysis with 1-way ANOVA, Bonferroni post hoc test, * * * * P > 0.0001.Analysis details in Supplementary Table 3a, b, and d.Scale bar in A, B, D, E, 100 μm, in C, F, 500 μm.

Fig. 7 .
Fig. 7. Comparison of anatomical input to L1. A) Schematic of input to L1 assessed with fb (left) and retrobeads (right).The bulk of input to L1 was from local neurons and a fraction of the input arising from long-range sources.Thalamus (including POm, zona incerta), ipsilateral and contralateral motor cortices, contralateral S1, ipsilateral S2, bilateral perirhinal ), and ipsilateral visual cortices all provide input to L1. B) Schematic of input to L1. Proportion of local and long-range input and classification of long-range input to L1. Local input arose mainly from L2/3 and L5.Long-range input arose from S2, motor, and visual cortices ipsilateral to the application site and was less pronounced from the contralateral S1 and M1.Abbreviations in a: I, ipsilateral; c, contralateral.Data from 4 mice for fb, 3 mice for retrobeads, analysis details in Supplementary Table4.