|Home | About | Journals | Submit | Contact Us | Français|
Systematic genetic access to GABAergic cell types will facilitate studying the function and development of inhibitory circuitry. However, single gene-driven recombinase lines mark relatively broad and heterogeneous cell populations. Although intersectional approaches improve precision, it remains unclear whether they can capture cell types defined by multiple features. Here we demonstrate that combinatorial genetic and viral approaches target restricted GABAergic subpopulations and cell types characterized by distinct laminar location, morphology, axonal projection, and electrophysiological properties. Intersectional embryonic transcription factor drivers allow finer fate mapping of progenitor pools that give rise to distinct GABAergic populations, including laminar cohorts. Conversion of progenitor fate restriction signals to constitutive recombinase expression enables viral targeting of cell types based on their lineage and birth time. Properly designed intersection, subtraction, conversion, and multi-color reporters enhance the precision and versatility of drivers and viral vectors. These strategies and tools will facilitate studying GABAergic neurons throughout the mouse brain.
He et al. present strategies and tools for the combinatorial targeting of GABAergic neurons in the mouse brain. Specific cell type targeting can be achieved by combining 2 to 3 driver-reporter alleles with viral vectors that engage multiple cell-defining features.
Understanding the functional organization of neural circuits requires comprehensive knowledge of the basic cellular components. Whereas glutamatergic pyramidal neurons make up ~80% of neurons in rodent neocortex and constitute the major cortical processing streams and output channels, GABAergic neurons regulate the delicate balance and dynamic operations of local circuitry and route information flow in cortical networks (Harris and Shepherd, 2015; Huang, 2014; Kepecs and Fishell, 2014). Although a minority, GABAergic neurons consist of diverse cell types that reflect an elaborate division of labor of an intricate inhibitory control system (Jiang et al., 2015; Klausberger and Somogyi, 2008; Markram et al., 2015; Markram et al., 2004; Roux and Buzsaki, 2015). However, comprehensive identification and characterization of GABAergic cell types have remained an unfulfilled goal, largely due to the lack of experimental access to specific cell types in a systematic way.
Although GABAergic neurons have been classified according to their morphology, connectivity, physiological properties, gene expression and developmental origin (Petilla Interneuron Nomenclature et al., 2008), a combination of multiple correlated features is likely more informative in cell type definition (Kepecs and Fishell, 2014). The generation of mouse recombinases driver lines begins to provide reliable access to multiple GABAergic cell populations and lineages (Taniguchi et al., 2011), and has facilitated progress in studying the development, connectivity, and function of cortical inhibitory interneurons (e.g. Cardin et al., 2009; Fu et al., 2014; Lee et al., 2013; Pfeffer et al., 2013; Pi et al., 2013; Taniguchi et al., 2013). However, as there is no simple correlation between single genes and specific cell types, most individual driver lines often include multiple cell types (Jiang et al., 2015; Munoz et al., 2014; Pronneke et al., 2015; Xu et al., 2013). Thus a key challenge of broad significance is to increase the specificity of cell type targeting. Enhanced specificity will facilitate multi-faceted studies toward elucidating the biological basis for defining GABAergic cell types, a fundamental issue in understanding cortical circuitry, and will provide a ground truth starting point for comprehensive identification of cortical cell types.
Intersectional methods based on 2–3 genes or other cell features improve the precision of cell targeting (Dymecki et al., 2010; Fenno et al., 2014; Madisen et al., 2015). However, strategically designed combinatorial driver and reporter lines are currently very limited. More importantly, few strategies have been demonstrated to incorporate multiple cell-defining phenotypes toward more specific cell targeting. Furthermore, there are practical limits to the number of driver/reporter allele combinations suitable for mouse engineering and breeding. Thus it remains unclear to what extent genetic methods will achieve sufficient specificity in targeting GABAergic subpopulations and cell types.
Here we have designed multiple combinatorial strategies, implemented a suite of genetic and viral tools, and demonstrated more precise targeting of GABAergic neurons and progenitors. We designed intersectional transcription factor drivers that allow more precise fate mapping of progenitor pools giving rise to more restricted GABAergic populations, including laminar cohorts. We showed that intersection of Cre and Flp lines driven by two marker genes targeted more restricted GABAergic subpopulations, and that incorporation of viral methods allowed capturing specific cell types. We further implemented a method that converted a transient fate restriction signal in progenitors to constitutive recombinase expression, which enabled viral targeting of cell types in the mature brain based on cell lineage and birth time. Therefore, highly specific targeting is feasible by combining 2 to 3 driver-reporter alleles with viral vectors to engage a spectrum of cell-defining features that include lineage, birth time, marker genes, and anatomy. Together these strategies and tools will accelerate progress in studying GABAergic circuits in the cerebral cortex and in other brain areas.
Our strategy is to combine a small set of driver-reporter alleles and viral vectors that incorporate multiple cell features to increase targeting specificity (Figure 1, Table 1, Table S1–2). Our newly generated drivers include: 1) Nkx2.1-ires-FlpO (Nkx2.1-Flp) for intersectional fate mapping of medial ganglionic eminence and preoptic area progenitors, 2) vasointestinal peptide-ires-FlpO (VIP-Flp) and somatostatin-ires-FlpO (SST-Flp) for dissecting these two broad cell populations, 3) a highly efficient inducible PV-2A-CreER driver which is far superior to the previous version (Taniguchi et al., 2011). The new reporter lines include: 1) an intersectional reporter expressing nuclear-targeted GFP (HG), 2) a high performance intersection-subtraction reporter (IS), 3) a conversion reporter that converts Cre expression to Flp expression, 4) a RGBow reporter that expresses one of three membrane-tethered fluorescent proteins in a random manner upon Cre-mediated recombination. We have also designed several viral vectors tailored for combinatorial targeting with driver and reporter lines (Table 1, Table S1–2). All mouse lines are deposited to the Jackson Laboratory for distribution.
The embryonic subpallium contains a developmental plan embedded in progenitors along the lateral-, medial- and caudal- ganglionic eminence (LGE, MGE, CGE) and the preoptic area (POA) that orchestrates the reliable generation of all GABAergic neurons in rodent neocortex (Gelman and Marin, 2010; Rudy et al., 2011) (Figure 1C). Beyond these major domains, 18 subdomains have been identified according to combinatorial TF expression patterns (Flames et al., 2007). Furthermore, the subpallium contains different types of progenitors, including radial glial progenitors (RGs), which reside in the ventricular zone, (VZ) and intermediate progenitors (IPs), which reside in the subventricular zone (SVZ) (Brown et al., 2011; Hansen et al., 2013; Harwell et al., 2015). Whether these progenitor subdomains and cell types contribute to different aspects of GABAergic neuron production and diversity remains unclear.
The homeodomain transcription factor Nkx2.1 is expressed in MGE and POA throughout embryonic neurogenesis in both RGs and IPs (Marin and Rubenstein, 2001; Sussel et al., 1999). We tested our Nkx2.1-Flp driver with multiple reporters and demonstrated that it captured all major groups of MGE and POA derived cortical GABAergic neurons (Figure S1).
The expression patterns of Nkx6.2 and Er81 at E12.5-E13.5 mark the dorsal and ventral MGE, respectively (among other LGE, CGE and POA domains, Figure 1C) (Flames et al., 2007; Sousa et al., 2009). By combining the Nkx6.2-CreER or Er81-CreER driver with Nkx2.1-Flp and the Ai65 intersectional reporter, we fate mapped dorsal vs. ventral MGE progenitors at E13.5. Herein we refer to Nkx2.1-Flp;Nkx6.2-CreER;Ai65 allele combination as Nkx2.1/Nkx6.2 and Nkx2.1-Flp;Er81-CreER;Ai65 as Nkx2.1/Er81 (Figure 2A–B); the same naming convention will be followed for other allele combinations throughout the text. In adult progeny, typically several hundred neurons were sparsely labeled in each cortex across multiple cortical areas and layers (Figure 2C–D, Figure S3). While neurons derived from Nkx2.1/Nkx6.2 vs. Nkx2.1/Er81 progenitors show laminar bias in their cortical settlement, the most notable difference between them is the morphological characteristics of their cell composition. Significantly more putative layer2/3 (L2/3) Martinotti cells were labeled in Nkx2.1/Nkx6.2 (~20%, Figure 2Ci, Figure S3A) compared to Nkx2.1/Er81(<6%, Figure 2Di). Neurons residing in L1/2 with dense horizontally extending axons (Figure 2Cii) and a rarer type of L6 neurons that extended axons along the corpus callosum (Figure 2Ciii) were only observed in Nkx2.1/Nkx6.2 but not found in Nkx2.1/Er81. Most of the labeled cells in Nkx2.1/Er81 were multipolar cells that resided in deep layers (especially L5) with radiating axons (Figure 2Dii–iii, Figure S3C, D). While Nkx2.1/Nkx6.2 progenitors produced only a small number of chandelier cells (ChCs) in both superficial and deep layers (Figure 2Civ, S3F-H), Nkx2.1/Er81progenitors produced more that were almost exclusively superficial layer ChCs (Figure 2Div, Figure S3E). Together, these results suggested that progenitors in dorsal vs. ventral MGE are fated to generate distinct sets of GABAergic neurons. Nkx2.1-Flp thus allows intersectional fate mapping of finer grain progenitor pools in MGE and POA defined by additional transcription factors.
The bHLH transcription factor Ascl1 is primarily expressed by IPs throughout the embryonic subpallium. Inducible fate mapping from IPs should allow more precise labeling of temporal cohorts of neurons born at the time of CreER induction than from RGs that undergo multiple rounds of neurogenesis. By combining the Ascl1-CreER knockin driver (Kim et al., 2011) with Nkx2.1-Flp, we demonstrated the feasibility of fate mapping MGE IPs at multiple embryonic time points. Tamoxifen induction at E12.5, E13.5 and E14.5 labeled GABAergic neurons that occupied cortical layers in a clear “inside-out” pattern: neurons born at E12.5 mainly occupy L5 to L6, at E13.5 mainly L4, and at E14.5 L2/3 (Figure 2 E–H). These laminar patterns of cell distribution were significantly sharper than those fate mapped from the Olig2-CreER driver, which included both RGs and IPs (Miyoshi et al., 2007). Therefore, the Nkx2.1-Flp driver allows differential fate mapping of IPs in MGE and POA defined by appropriate TF drivers. This provides an improved tool for analyzing the relationship between neuronal birthdate and mature phenotypes that contribute to cell identity.
Among GABAergic neurons derived from CGE, the VIP expressing cells make up to 15% of all cortical interneurons (Miyoshi et al., 2010; Rudy et al., 2011) and display multiple morphological and physiological features (Cauli et al., 2004; Pronneke et al., 2015). The generation of a VIP-ires-Cre driver (Taniguchi et al., 2011) allowed a direct demonstration that certain VIP interneurons preferential inhibit other GABAergic neurons (Lee et al., 2013; Pfeffer et al., 2013; Pi et al., 2013), and enabled studies that revealed the function of this dis-inhibitory circuit in the context of behavior (Fu et al., 2014; Pi et al., 2013).
On the other hand, the VIP population consists of several subpopulations that show distinct morphology (e.g. bipolar type, small basket type (Bayraktar et al., 2000)), intrinsic electrophysiological properties (Lee et al., 2013; Pronneke et al., 2015) and differential marker expression (e.g. calretinin, cholecystokinin, corticotropin-releasing hormone (Kubota, 2014; Kubota et al., 2011)). To demonstrate the morphological heterogeneity within the VIP population, we used the VIP-ires-Cre driver to activate a triple-color reporter RGBow (Figure 1D) that we have generated to achieve multi-color labeling of VIP cells in the same animal (Figure 3B).
To further define the VIP cell subpopulations, we generated a VIP -Flp line (Figure 3A) which showed a near identical recombination specificity and efficiency as VIP-ires-Cre (Figure S4 A, B). We combined this line with calretinin(CR) - and cholecystokinin(CCK)-ires-Cre (Taniguchi et al., 2011), to carry out intersectional targeting (Figure 3, Figure S4, Table S2–4). The VIP/CR intersection resulted in a significant enrichment of cells located in lower L2/3, with smaller fractions in deep layers (Figure 3A, C). These cells showed prominent vertically oriented bipolar morphology (18 out of 24 biocytin filled single cells recovered from patch clamp recording in brain slice), consistent with previous studies (Cauli et al., 2014; Jiang et al., 2015); they extended vertically oriented dendrites that span much of the cortical thickness and projected vertically restricted axons (Figure 3Ei–ii). In rare cases, neurons with multipolar morphology can also be observed (Figure 3Eiii). Previous studies reported that VIP interneurons as a group show a range of different firing patterns including irregular spiking (IS) which has also been frequently reported for CR expressing interneurons (Cauli et al., 1997; Cauli et al., 2000), bursting non-adapting (b-NA), and fast adapting (fAD) (Figure 3G) (Jiang et al., 2015; Lee et al., 2010; Miyoshi et al., 2010; Pronneke et al., 2015). The VIP/CR intersection significantly enriched for cells with IS firing pattern, with a subpopulation also showing b-NA or fAD (Figure 3G). Both the prominently bipolar axon morphology and the irregular-spiking intrinsic features are characteristic of interneuron-selective cells (Lee et al., 2013), suggesting that the VIP/CR intersection significantly enriched for this subpopulation.
In contrast, most of the VIP/CCK cell somata were enriched in upper L2/3 near the L1 and L2 boundary, with a smaller fraction in deep layers (Figure 3A, D). 15 out of 17 VIP/CCK neurons recovered from cell filling show multipolar morphology. Many of these cells extend highly exuberant local axon arbors within L2/3, but their multi-polar dendrites extended to L1 (Figure 3Fi–iii). We also recovered rare cells that projected axons to L1 (Figure 3Fiv). The intrinsic properties of a substantial proportion of VIP/CCK cells were either bursting or fast-adapting, which were less frequently observed in VIP/CR cells (Figure 3G, Table S3). These morphological and physiological features suggest that VIP/CCK cells likely correspond to a set of VIP− and CCK-expressing small basket cells (Kawaguchi and Kubota, 1998; Kubota, 2014; Ma et al., 2011). It is notable that this population appears to be found in a large scale recording in juvenile rat (Markram et al., 2015) but was undetected by similar studies in the mouse even using the VIP-ires-Cre driver (Jiang et al., 2015; Pronneke et al., 2015). Together, these results demonstrate the presence of multiple distinct subpopulations within the VIP population and the feasibility of intersectional targeting of these subpopulations.
SST-expressing neurons constitute approximately half of MGE-POA derived GABAergic neurons and up to ~30% of all GABAergic neurons in many cortical areas (Rudy et al., 2011). SST cells are present throughout L2 to L6, with notable enrichment in deep layers (Pfeffer et al., 2013), and consist of multiple heterogeneous subpopulations (Ma et al., 2006; Munoz et al., 2014; Tomioka et al., 2005; Xu et al., 2013). Although the SST-ires-Cre driver confers experimental access to SST cells, the diverse cell types within this broad population are not well understood.
There is evidence that some SST and CR double positive cells may correspond to dendrite-targeting Martinotti cells (Wang et al., 2004; Xu et al., 2006) and SST and nNOS double positive cells to long projection cells (Tomioka et al., 2005). To examine these correlations, we generated a SST-Flp driver that allows intersectional dissection of the SST population. The SST-Flp driver gave rise to a nearly identical recombination pattern to SST-ires-Cre and labeled SST cells across L2 to L6 (Figure 4A). In sharp contrast, SST/CR intersection revealed a subpopulation with striking laminar distribution patterns: SST/CR cells were highly enriched in L2/3 and L5a, with less abundance in L5b and L6 but were virtually absent in L4 (Figure 4A, B, Figure 5B). These cells often had prominent laminar-targeted axon projection and arborization, with L2/3 and L5a cells projecting axons to L1 (Figure 4C), and L6 cells likely projecting axons to L4/5a (Markram et al., 2015). Single cell reconstruction confirmed that L2/3 and L5 SST/CR cells have Martinotti cell morphology with an ascending axon that arborizes in L1 with additional branching in L2/3 (Figure 4C), consistent with the population pattern.
We compared the intrinsic properties of cells sampled in SST/CR with those sampled from the whole SST population in SST-Cre;Ai14 mice. In the SST population (n= 71) we observed: 1) cells with a strongly adapting firing pattern (~44%), 2) cells with narrower spikes that can sustain much higher firing frequencies resembling the L4-targeting SST cells in the X-94 mouse line (~47%) (Ma et al., 2006; Xu et al., 2013), and 3) cells with very high input resistance that produced rebound spikes following a hyperpolarization, resembling the SST cells reported in the X-98 mouse line (~9%) (Ma et al., 2006) (Figure 4D). In the SST/CR population, intrinsic electrophysiological properties are much more homogeneous, with 90% strongly characterized as adapting firing pattern (AD) (Figure 4D). Therefore, the SST/CR intersection significantly enriched, although did not exclusively purify, Martinotti cells in L2/3 and L5, and it essentially excluded L4-targeting and long range projection SST neurons.
Nitric oxide (NO) is a versatile signaling molecule in the brain synthesized by the neuronal isoform of nitric oxide synthase (nNOS) (Tricoire et al., 2013). In mature cerebral cortex, all nNOS cells are GABAergic and consist of two major populations (Kubota et al., 2011; Perrenoud et al., 2012). Type II neurons weakly stain for nNOS and have small somata distributed in multiple layers; they include several subpopulations (e.g. neurogliaform cells) (Magno et al., 2012; Olah et al., 2009; Perrenoud et al., 2012). In contrast, type I neurons have heavy nNOS stain and large somata distributed predominantly in L5 and L6 and, to a lesser extent, L2 (Kubota et al., 1994; Tricoire et al., 2013). These highly unique GABAergic neurons are not interneurons but extend axons to distant cortical and subcortical regions(Tomioka et al., 2005). They are selectively activated during slow wave sleep (Kilduff et al., 2011). The distribution, connectivity and function of type I nNOS neurons are poorly understood due to the lack of specific experimental access.
In SST/nNOS mice, a single dose tamoxifen induction at P16 labeled ~85% of type I nNOS neurons characterized by strong nNOS immunostaining and large cell bodies (Figure 5B, Figure S5B–C). SST/nNOS cells were prominently localized to L5 and L6, and to a much lesser extent at the L1/2 border, with a small number of cells also scattered in L3 and the white matter (Figure 5B, Figure S6A). The abundance and laminar distribution of type I cells varied according to cortical areas (Figure 5B, Table S5). At the population level, the axon arbors of type I cells were dense and extensive, forming a matrix that appeared to cover the entire cortical volume (Figure 5B, Figure S6 A–F). These axons were dense in L1 and also extended into the white matter (Figure S6B–E). This particularly rare cell type has not been described by previous physiological studies (Jiang et al., 2015; Markram et al., 2015). We found that the intrinsic properties of type I cells are quite homogeneous, characterized by an irregular accommodating spiking profile with a delayed first spike; their action potentials have relatively narrow spike widths (0.51 ms) and always showed a biphasic after potential consisting of a fast afterdepolarization(ADP) followed by an afterhyperpolarization (AHP) (Figure 5C). Single cell reconstruction of L6 SST/nNOS cells revealed that while their multipolar dendrites were largely restricted to the same layer, their axons arborized extensively with branches entering the white matter that could not be fully traced in brain slice (Figure 5D).
We next designed a genetic and viral intersection strategy to examine the distribution pattern and axon trajectory of SST/nNOS cells that make distant projections to a specific cortical area. In a SST-Flp;Ai65 mouse, we injected a retrograde-transported canine virus CAV2-Cre (Hnasko et al., 2006; Kremer, 2005) in the primary motor cortex (M1) (Figure 5E–H). Only Flp expressing SST cells that were infected with CAV2-Cre would activate the Ai65 reporter allele and express tdTomato. At the injection site, many SST cells expressed tdTomato as their local axons were infected by CAV2-Cre (Figure 5E square boxed area). In more distant cortical areas from M1, only SST cells that projected across a long distance to the injection site could acquire CAV2-Cre and activate tdTomato expression. We found that M1-projecting SST neurons were distributed in widespread cortical areas ranging from the rostral orbital frontal cortex to the primary and secondary visual cortex in the caudal end (Figure 5E–H, Movie S1). Their rostro-caudal axon projection distance exceeded 3 millimeters across several cortical areas, and cells were labeled more than 2.5 mm away laterally in the somatosensory area. We detected only 1 cell in the contralateral hemisphere in 3 mice. Approximately half of the SST projection neurons were strongly nNOS positive (Figure 5G), while the other half were either negative or very weakly stained, consistent with an earlier report (Tamamaki and Tomioka, 2010). As expected, these M1 long projection cells were highly enriched in L5 and L6, with another minor group of cells in L2/3 (Movie S1). In areas distant from the injection site, sparse cell labeling revealed their extensive axon arbors, which often projected branches into L1 (Figure S6F, G). Thus SST/nNOS neurons likely project to distant cortical areas through L1 as well as the white matter. The SST-Flp and CAV2-Cre intersection strategy thus provides a reliable and efficient method to systematically map the areal organization and projection pattern of this unique and intriguing class of cortical GABAergic projection neurons.
We further developed a combinatorial genetic-AAV strategy, termed “Triple-Trigger”, to establish viral-based experimental access that allows integration of multiple tools for studying the connectivity and function of SST/nNOS neurons (Figure 1D, Figure 5I). SST-Flp;nNOS-CreER mice were first bred with a Rosa26-loxPNeoSTOPloxP-tTA (LNL-tTA) (Wang et al., 2008) line to generate SST-Flp;nNOS-CreER;LNL-tTA mice, in which tamoxifen induction can convert CreER into constitutive tTA expression in nNOS neurons. We further generated a new AAV vector in which TVA-mcherry expression is dependent on both tTA and Flp activities (Figure 1D, Table 1, Figure 5I). We injected this AAV-tetO-fDIO-TVA-mcherry into the motor cortex of SST-Flp;nNOS-CreER;LNL-tTA mice, followed by tamoxifen induction of Cre recombination to activate tTA expression. Three weeks later, a pseudo-typed glycoprotein deficient rabies virus expressing GFP (RV-EnvA-ΔG-GFP) was injected in the same area to infect TVA expressing neurons. This method yielded intense GFP labeling of a small number (~10) of SST/nNOS neurons in L5/6 of motor cortex in a restricted region (<1mm diameter) (Figure 5J), which extended dense and widespread axon arbors towards L1. Descending axon branches from L1 were observed in sections far away from their parent cell bodies. The longest projecting axons were detected 2.4 mm away from the center of the injection site (Figure 5K, L).
In summary, nNOS-CreER and SST-Flp intersection captures a specific and unique type of GABAergic projection neurons that likely contributes to the regulation of neuronal populations in widespread cortical regions. Further incorporation of intersectional viral vectors allows versatile experimental access enabling investigators to map the projection and connectivity patterns, monitor the activity and manipulate the function of this enigmatic and likely important cell type.
In contrast to glutamatergic pyramidal neurons, which are organized with laminar patterns characteristic to the neocortex, GABAergic neurons appear more evenly distributed across the cortical thickness. However, it is unclear to what extent more specific GABAergic subpopulations or cell types are deployed to more restricted cortical layers.
Using serial two-photon tomography (Ragan et al., 2012), we quantified the laminar distribution of six GABAergic populations in somatosensory, motor and medial prefrontal cortices (Figure 6). We found that, while the broad VIP and SST populations showed laminar biases, each of the four subpopulations revealed significantly more restricted laminar distributions, with some variations in different cortical areas. Whereas VIP/CCK cells were highly enriched in upper L2/3, most VIP/CR cells were located in deeper L2/3 (Figure 3, ,6).6). SST/CR showed striking enrichment in L2/3, L5a, L6, and were nearly absent in L4 in barrel cortex (Figure 6, note the dip at 50% of cortical depth in barrel cortex, which was not present in the other two cortical regions; Table S5). Lastly, SST/nNOS neurons were concentrated in deep L5/6 (~74% combined in bin 6–8, Table S5), with another much smaller cohort in upper L2/3 (~8% in bin 2, Table S5). Together, these results indicate that GABAergic subpopulations and cell types are in fact deployed to restricted cortical layers with characteristic distribution patterns. It is possible that further purification of cell types might reveal even more specific laminar distribution patterns (e.g. ChCs (Taniguchi et al., 2013)). As the location of cell somata anchors the spatial distribution of their mostly local dendritic and axonal arbors that mediate connectivity (with a notable exception for the long projection axons of type I nNOS neurons), the laminar restriction of GABAergic neurons suggests that they are likely embedded in spatially segregated subnetworks.
ChCs are one of the most distinctive types of cortical GABAergic neurons that specifically innervate pyramidal neurons at their axon initial segment (Somogyi et al., 1982), the site of action potential generation, and might exert powerful control over the firing of a pyramidal neuron ensemble (Woodruff et al., 2010). Genetic fate mapping using the Nkx2.1-CreER driver revealed the developmental origin of ChCs and allowed reliable morphological labeling (Taniguchi et al., 2013). Since Nkx2.1 is mainly expressed in progenitors and no specific markers have been identified for postmitotic ChCs, it remains difficult to monitor and manipulate ChCs in mature cortex.
We have designed a strategy that converts the developmentally transient Nkx2.1-CreER induction to a constitutive Flp expression, which allows specific and reliable viral targeting of ChCs in mature animals. We have generated a mouse line, Rosa26-loxPSTOPloxP-Flpo (LSL-Flp), in which Flp expression is conditional upon Cre-mediated excision of a transcription STOP cassette (Figure 1C). In Nkx2.1-CreER;LSL-Flp mice, late embryonic tamoxifen induction results in constitutive Flp expression in postmitotic ChCs throughout the mouse’s life span (Figure 7A). Flp-dependent AAV-fDIO-GFP injection in the mature primary somatosensory cortex of E17.5 tamoxifen-induced Nkx2.1-CreER; LSL-Flp mice resulted in specific and intense labeling of ChCs (Figure 7B). Morphological reconstruction revealed that, interestingly, L2 ChCs consisted of at least two subsets defined by their axon morphology. Most L2 ChCs extended dense axon arbors within L2/3, characterized by highly distinct vertically oriented cartridges of synaptic boutons (Figure 7C). Another more rare subset projected additional axonal branches to L5 and L6, with elaborate local arbors studded with cartridges of boutons (Figure 7D), consistent with a previous report using in utero electroporation labeling of ChCs (Tai et al., 2014). Such dual laminar-projecting ChCs likely coordinate the firing of different sets of pyramidal neuron ensembles spatially segregated in different cortical layers. Combined with other viral tools that express functional effectors such as ChR2, DREADDs or pseudo-typed rabies virus receptor TVA, the conversion method allows manipulation of ChCs’ activity in vivo and tracing their local and long range inputs (J.T and J.L et al, in preparation).
Our previous study revealed that ChCs consist of both PV immunopositive (PV+) and PV immunonegative (PV−) subsets (Taniguchi et al., 2013)(Figure 7E). As PV levels may reflect different functional subsets of hippocampal and neocortical basket cells (Dehorter et al., 2015; Donato et al., 2013), it is possible that PV+ and PV− ChCs have different physiological and functional properties. We designed a novel intersection and subtraction reporter (IS reporter): Rosa26-loxPSTOPloxP-Frt-Tdtomato-Frt-3XGFP. When combined with Nkx2.1-CreER; PV-Flp, the IS reporter differentially labeled PV+ ChCs with tdTomato and PV− ChCs with GFP (Figure 7F–H). This simultaneous and differential labeling will facilitate the physiological characterization of PV+ vs PV− ChCs. More generally, the high performance IS reporter will allow differential labeling of cell types using appropriate Cre and Flp driver lines. For example, using the Nkx2.1-CreER; SST –Flp;IS mouse, we were able to differentially label SST+ vs. SST− cells born at similar time from the Nkx2.1 progenitors (Figure S7).
Efforts toward achieving a comprehensive catalogue of neuronal cell types in brain circuits face two fundamental challenges: 1) specificity – defining and working with pure cell types rather than mixtures of types, and 2) comprehensiveness – identifying many, if not all, cell types in the circuit such that their connectivity and interactions can be examined systematically toward understanding circuit operations (Seung and Sumbul, 2014). While researchers working in certain brain areas (e.g. the retina) have achieved consensus in solving these two issues and a cell catalogue appears within reach (Sanes and Masland, 2015; Seung and Sumbul, 2014), a cell catalogue of the neocortex is a longer term goal, and even the definition of cortical cell types remains contentious (DeFelipe et al., 2013). Framing these challenges in the context of understanding cortical GABAergic neurons, achieving sufficient specificity in cell targeting and cell type definition is a prerequisite for progress toward comprehensive cell type discovery. Given the scope of GABAergic neuron diversity, a major technical difficulty is simply the lack of reliable experimental access to a set of well-parsed cell populations such that different investigators can obtain and discuss results concerning the same cell population. As demonstrated by studies in the retina (Sanes and Masland, 2015) and invertebrate (Mauss et al., 2015) circuits, the more specific cell populations are studied, the more clear answers concerning cell types are likely to emerge.
Given the practical limitations in mouse engineering that severely constrain the number of gene combinations, it has remained unclear to what extent genetic methods will achieve sufficient specificity in targeting cell types, which are better defined by multiple features beyond 1 or 2 molecular markers. The key strength of our approach is not only to intersect gene expression patterns but, more importantly, to combine 2 or 3 strategically designed driver-reporter alleles with viral vectors to achieve multi-feature targeting of cell types. Further, in our scheme, gene expression patterns at different developmental stages of a cell’s life span, from transcription factors in progenitors to mature cell markers, are recruited for combinatorial targeting. These strategies and tools can be applied and modified for specific cell targeting in many brain circuits.
Cortical GABAergic neurons comprise three largely non-overlapping populations marked by PV, SST and 5HT3aR (~40% of the 5HT3aR population are marked by VIP) (Rudy et al., 2011). The generation of PV− (Madisen et al., 2015), SST− and VIP-Flp drivers allows intersectional targeting of more restricted subpopulations within these broad classes as well as simultaneous targeting of non-overlapping populations by combining with proper Cre lines. These knock-in Flp lines faithfully recapitulated the endogenous gene expression pattern and have identical recombination pattern as the corresponding knockin Cre lines (Taniguchi et al., 2011). This reliability of gene targeting-based approach is critical for the rational design of combinatorial cell targeting. Here we demonstrate that intersectional targeting with VIP− and SST-Flp substantially increases specificity. On the other hand, most subpopulations captured by two-gene intersections (e.g. VIP/CR, VIP/CCK, SST/CR) may still contain more than one cell type defined by axonal morphology and/or firing patterns (Figure 3, ,4).4). SST/nNOS may represent an exception and appears to capture a rare type of categorically distinct GABAergic projection neurons (Figure 5). Recent large-scale recordings in the mouse brain slice have yielded impressive progress in characterizing neocortical GABAergic cell types (Jiang et al., 2015), yet appear to have missed rare types such as SST/nNOS and VIP/CCK. A coordinated effort that combines large scale recording and intersectional genetic targeting will likely facilitate further advances. An increased repertoire of conditional viral vectors that confer tool-gene expression (e.g.ChR2) upon Cre, Flp, and/or tTA activation (Fenno et al., 2014; Madisen et al., 2015) will further enhance the versatility and experimental power of cell targeting and manipulation.
Previous studies using BrdU-based birth dating have provided evidence that cortical GABAergic neurons tend to populate the cortex in an inside-out sequence (Fairen et al., 1986; Miller, 1985; Valcanis and Tan, 2003). More recent experiments using the Olig2-CreER driver to fate map MGE progenitors (including both RGs and IPs) showed that early- and late-born interneurons do follow an inside-out layering trend and further display distinct, though diverse physiological properties (Miyoshi et al., 2007). It remains unclear to what extent cell birth date correlates and contributes to their laminar positioning and other cell phenotypes. Because IPs marked by Ascl1 are likely more fate-restricted than RGs and often only divide once to produce two neurons (Brown et al., 2011), Ascl1-CreER and Nkx2.1-Flp intersection allows more precise fate mapping of MGE IPs and birth dating of their progenies. Indeed our results revealed tight and striking laminar patterns and an inside-out sequence. This method thus enables studying the relationship among cell lineage, birth time/order, laminar location and cell phenotypes at a much higher resolution. By intersecting with other IP drivers (e.g. Dlx1-CreER), this will address whether IPs defined by different transcription factors generate distinct set of interneurons. This strategy further allows comprehensive discovery of MGE derived GABAergic neurons through sparse labeling and phenotyping.
Genetic strategies that engage the molecular mechanisms of developmental programs such as cell lineage and birth order are particularly effective in fate mapping and targeting cell types (Harris et al., 2015; Taniguchi et al., 2013). However, transcription factors that act in progenitors often cease expression or change expression patterns in postmitotic neurons and thus are not directly suited for gaining experimental access to mature cell types. The essence of our conversion strategy is to convert a transient fate restriction signal in progenitors to a permanent “genetic handle” (e.g. Flp, tTA) in progeny neurons that allows viral access to engage anatomical features. In addition, we demonstrate that TF expression in progenitors (e.g. Nkx2.1) can be combined with phenotype-defining mature markers (e.g. PV) to target highly specific subtypes (e.g. PV+ vs. PV− ChCs) using appropriate reporters (e.g. IS). Together, these strategies combine 2 or 3 driver-reporter alleles with viral vectors to engage a spectrum of cell-defining features, thereby significantly expanding the versatility and enhancing the specificity of cell type targeting.
The MGE gives rise to >60% of cortical GABAergic neurons that are categorically distinct from those generated in the CGE (Kessaris et al., 2014; Rudy et al., 2011). Because all MGE progenitors express Nkx2.1 and most MGE-derived interneurons can be parsed by PV or SST into largely non-overlapping populations (Kessaris et al., 2014; Rudy et al., 2011), these boundary-defining molecular markers provide the basis and scheme for a systematic genetic dissection. First, by combining PV− or SST-Flp drivers with proper Cre drivers, intersection/subtraction reporters, and viral vectors, it is possible to target increasingly restricted subpopulations and cell types based on mature cell markers. Second, by combining Nkx2.1-Flp driver with proper TF Cre drivers, it is possible to target restricted progenitor pools (e.g. spatial domains, RG vs. IP) and capture more specific cell populations, including laminar cohorts. Third, by combining a progenitor driver (e.g. Ascl1-CreER), a postmitotic cell driver (e.g. PV-Flp, SST-Flp), and proper reporters and viral vectors, it is possible to target cells jointly defined by lineage, birth order, marker expression, and anatomy. With increasing knowledge of cellular gene expression revealed by transcriptome analysis (Tasic et al., 2016; Zeisel et al., 2015), even a modest number of additional driver lines that expand combinatorial targeting will enable decisive progress.
In summary, although the diversity of cortical GABAergic neurons has been difficult to untangle, this diversity likely results from the manifestation of developmental programs that reliably generate a large but not infinitely large set of neuronal cardinal types, which subsequently further differentiate into diverse and distinct cell types. With genetic and viral strategies that engage these developmental and molecular programs, as well as anatomic features, it is increasingly feasible to enhance the specificity of cell targeting and move toward a more comprehensive identification of MGE-derived cell types - a major fraction of cortical GABAergic neurons.
New driver and reporter mouse lines were generated as described before (Taniguchi et al., 2011). All experimental procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of CSHL or NYU School of Medicine in accordance with NIH guidelines, or in accordance with institutional guidelines of Institutes of Brain Science, Fudan University, China.
Mice were anaesthetized (using avertin or chloral hydrate) and intracardially perfused with saline followed by 4% paraformaldehyde (PFA) in 0.1 M PB. Following 24hrs of post fixation at 4°C, brain slices at 50μm or 75μm thickness were sectioned with a Leica 1000s vibratome, stained and imaged as described before (Taniguchi et al., 2011).
Adult mice were anesthetized using a ketamine/xylazine mixture and stereotactic injections were performed via rodent stereotax as described (Huang et al., 2014).
Neurons in coronal slices (300 μm) containing the barrel field of the somatosensory cortex from P21–30 mice were recorded, filled with biocytin and reconstructed as described before (Xu et al., 2013).
We are grateful to Hongkui Zeng for providing Ai3 and Ai65 vectors and the Ai65 mouse line, Dawen Cai, Jeff Lichtman and Joshua Sanes for sharing the Brainbow vector and antibodies, Sang Yong Kim for help with generation of knock in mice, Adam Kepecs and Hiroki Taniguchi, Dhananjay Huilgol for comments on the manuscript. This work was supported in part by NIH 5U01 MH078844-05, CSHL Robertson Neuroscience Fund to Z.J.H., National Natural Science Foundation of China 31471037, 91432106, 31421091 and Shanghai Yangfan 14YF1400500 to M.H., NARSAD Young Investigator Grant to Y.K., NIH U01MH105971 to P.O., and NIH R01NS30989 and P01NS074972 to B.R., J.T, S.M.K. and J.M.L. are supported by NRSA F30 Medical Scientist Predoctoral Fellowships.
AUTHOR CONTRIBUTIONSZ.J.H conceived and organized the study. M.H. and Z.J.H. designed the experiments. M.H., P.W and Z.J.H designed and generated all new knockin mouse lines. J.T., M.H. and P.W. designed and generated all new AAVs. M.H., S.M.K, Y.C., Y.H., L.G., and J.M.L. conducted mouse breeding, tamoxifen induction, immunostaining, imaging and quantification. J.T., J.M.L., and L.G. conducted virus injection experiments. Y. K. P.O, M.H., J.M.L. performed STP experiments. S.L, M.N. I.K. and B.R performed electrophysiology recording and single cell reconstruction. M.H. and Z.J.H wrote the manuscript, with input from all other co-authors.