Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleResearch Article: Methods/New Tools, Novel Tools and Methods

Regulatory Elements Inserted into AAVs Confer Preferential Activity in Cortical Interneurons

Anna N. Rubin, Ruchi Malik, Kathleen K. A. Cho, Kenneth J. Lim, Susan Lindtner, Sarah E. Robinson Schwartz, Daniel Vogt, Vikaas S. Sohal and John L. R. Rubenstein
eNeuro 16 November 2020, 7 (6) ENEURO.0211-20.2020; DOI: https://doi.org/10.1523/ENEURO.0211-20.2020
Anna N. Rubin
1Nina Ireland Laboratory of Developmental Neurobiology, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Anna N. Rubin
Ruchi Malik
2Center for Integrative Neuroscience, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ruchi Malik
Kathleen K. A. Cho
2Center for Integrative Neuroscience, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenneth J. Lim
1Nina Ireland Laboratory of Developmental Neurobiology, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Susan Lindtner
1Nina Ireland Laboratory of Developmental Neurobiology, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sarah E. Robinson Schwartz
2Center for Integrative Neuroscience, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel Vogt
1Nina Ireland Laboratory of Developmental Neurobiology, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vikaas S. Sohal
2Center for Integrative Neuroscience, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Vikaas S. Sohal
John L. R. Rubenstein
1Nina Ireland Laboratory of Developmental Neurobiology, Department of Psychiatry, University of California San Francisco Weill Institute for Neurosciences, University of California San Francisco, San Francisco, CA 94158
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for John L. R. Rubenstein
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Cortical interneuron (CIN) dysfunction is thought to play a major role in neuropsychiatric conditions like epilepsy, schizophrenia and autism. It is therefore essential to understand how the development, physiology, and functions of CINs influence cortical circuit activity and behavior in model organisms such as mice and primates. While transgenic driver lines are powerful tools for studying CINs in mice, this technology is limited in other species. An alternative approach is to use viral vectors such as AAV, which can be used in multiple species including primates and also have potential for therapeutic use in humans. Thus, we sought to discover gene regulatory enhancer elements (REs) that can be used in viral vectors to drive expression in specific cell types. The present study describes the systematic genome-wide identification of putative REs (pREs) that are preferentially active in immature CINs by histone modification chromatin immunoprecipitation and sequencing (ChIP-seq). We evaluated two novel pREs in AAV vectors, alongside the well-established Dlx I12b enhancer, and found that they drove CIN-specific reporter expression in adult mice. We also showed that the identified Arl4d pRE could drive sufficient expression of channelrhodopsin for optogenetic rescue of behavioral deficits in the Dlx5/6+/− mouse model of fast-spiking CIN dysfunction.

  • AAV
  • cortical interneurons
  • enhancers
  • fast spiking
  • regular spiking

Significance Statement

Uncovering the contribution of cortical interneuron (CIN) dysfunction to neuropsychiatric conditions is pivotal to generating therapies for these debilitating disorders. To this end, it is important to study the development, physiology, and function of CINs in model organisms. While transgenic driver lines are powerful tools to study CINs in mice, this technology is limited in other species. As an alternative, gene regulatory enhancer elements (REs) can be used in viral vectors to drive expression in multiple species including primates. The present study describes the genome-wide identification of REs preferentially active in CINs and the use of these RE AAVs for CIN-specific targeting in adult mice. The methodology and novel REs described here provide new tools for studying and targeting CINs.

Introduction

GABAergic cortical interneurons (CINs) constitute a diverse population of biochemically, physiologically and morphologically distinct neuronal subtypes that have specialized functions within cortical circuitry (Tremblay et al., 2016; Fishell and Kepecs, 2020). Based on their neurochemical and cellular properties, CINs can be broadly divided into three major subtypes: (1) parvalbumin (PV)-expressing CINs, which have fast-spiking physiological properties and innervate the perisomatic or axonal regions of excitatory pyramidal neurons (PNs); (2) somatostatin (SST)-expressing CINs, which typically have regular-spiking properties and target distal dendritic regions of PNs; and (3) vasoactive peptide (VIP)-expressing CINs, which have irregular-spiking properties and cause disinhibition of PNs by preferentially targeting other CINs. Multiple lines of evidence highlight that these CIN subtypes play dissociable roles in cortical circuits. In particular, fast-spiking PV+ CINs promote cortical γ oscillations, which are necessary for cognitive functions such as attention, working memory, and cognitive flexibility (Cardin et al., 2009; Sohal et al., 2009; Cho et al., 2015), and SST+ CINs play important roles in cortical β oscillations, input integration, and sensory processing (Yavorska and Wehr, 2016; Chen et al., 2017).

Our understanding of the development and functional heterogeneity of CINs owes heavily to the use of transgenic mouse lines expressing Cre recombinase in molecularly specified CIN subtypes (Taniguchi et al., 2011). Although these Cre-driver lines have been extremely valuable to the field, the current reliance on them has some limitations. First, while Cre-driver lines can be used in combination with other recombinases (such as Flp) to label different subtypes of neurons, the ability to manipulate and record from multiple classes of CINs simultaneously or study interactions of CINs with specific types of PNs has been lacking. Second, the Cre-driver lines are not able to distinguish between additional subclasses within the PV, SST, and VIP CIN subtypes, which have distinct physiological properties and molecular markers (Nigro et al., 2018). Third, while Cre-driver lines are widely available for mice, this technology is currently unsuitable for research in other species such as non-human primates.

A large number of studies implicate loss or alteration of the activity of CIN subtypes, particularly PV+ fast-spiking CINs, in neurologic diseases such as schizophrenia, epilepsy and autism (Lewis et al., 2012; Cho and Sohal, 2014; Vogt et al., 2015, 2017; Jiang et al., 2016; Hashemi et al., 2017; Malik et al., 2019). To develop targeted treatments for these disorders, it is imperative to develop technologies that will allow manipulation of the activity of CIN subtypes in humans.

The use of recombinant AAV vectors has become a widespread technique for labeling neurons as well as monitoring and modifying their electrophysiological activity (Betley and Sternson, 2011). For instance, AAV1, AAV5, and AAV9 disperse throughout the brain readily and can infect most neurons (Burger et al., 2004; Taymans et al., 2007; Nathanson et al., 2009b; Watakabe et al., 2015). AAVs have the added benefit of being translatable across species and are becoming a promising avenue for human therapeutic intervention (Haggerty et al., 2020).

Despite these advantages, limited viral genome capacity has hampered the generation of AAVs to target specific CIN subclasses. Because of their small size (often <1 kb) and specific patterns of activity, gene regulatory enhancer elements (REs) such as enhancers have proven to be a good tool for driving CIN-specific expression by viral methods. In particular, Dlx1/2 RE I12b, which labels CINs in transgenic mice (Potter et al., 2009), has been used in lentiviral (Arguello et al., 2013; Vogt et al., 2014) and AAV vectors to drive the expression of fluorophores and channelrhodopsin (Lee et al., 2014a,b; Cho et al., 2015). Subsequent work showed that another RE located at the Dlx5/6 locus has similar properties and demonstrated additional applications of this approach (Dimidschstein et al., 2016). However, viral tools targeting subclasses of CINs are in their infancy, and many still require the use of intersectional recombinase-based strategies (Mehta et al., 2019).

Here, we describe a systematic genome-wide identification of REs putatively active in immature CINs. From these, we selected two candidates and tested their utility in AAV vectors. We found that these RE AAVs drove robust, CIN-specific expression when injected into the medial prefrontal cortex (mPFC). Notably, we observed that the RE AAVs differentially target distinct CIN populations; specifically, they have a different preference for CINs with fast-spiking versus regular-spiking properties. We also demonstrate that the AAV with higher specificity for fast-spiking CINs can drive sufficient expression for in vivo optogenetic rescue of cognitive flexibility deficits in Dlx5/6+/− heterozygous mutant mice.

Materials and Methods

Mice

All animal care, procedures, and experiments were conducted in accordance with the NIH guidelines and approved by the University of California San Francisco animal care committee’s regulations. Mice were group housed (two to five siblings) in a temperature-controlled environment (22–24°C), had ad libitum access to food and water, and were reared in normal lighting conditions (12/12 h light/dark cycle). Mice of either sex were used, except in the case of behavioral experiments, which included only male mice.

Gad67-GFP (MGI ID: 3590301; Tamamaki et al., 2003), I12b-Cre (MGI ID: 3833422; Potter et al., 2009), and Ai14 (MGI ID: 3809524; Madisen et al., 2010) mice were maintained on a mixed background outcrossed to CD-1. Wild-type CD-1 mice were used for electrophysiology experiments. Dlx5/6+/− mice (MGI ID: 5583955; Robledo et al., 2002) were backcrossed to C57Bl/6 mice for at least six generations for behavioral experiments.

Epigenomic experiments

Histone modification chromatin immunoprecipitation and sequencing (ChIP-seq)

Immature CINs and non-CINs were isolated from P2 Gad67-GFP heterozygotes for native ChIP for histone posttranslational modifications. The pups were anesthetized on ice and decapitated. Neocortical tissue was dissected in ice-cold Earle’s balanced saline solution (EBSS) and cut into small pieces. The tissue was dissociated using the Papain Dissociation System (Worthington) according to the manufacturer’s instructions. After addition of the inhibitor solution, the tissue was gently aspirated 10–15 times with a P1000 pipette and filtered through a 40-µm filter to achieve a single-cell suspension. The cells were spun down and resuspended in EBSS for fluorescence-activated cell sorting (FACS) on a FACSAria II machine (BD Biosciences). Sorted GFP-positive CINs and GFP-negative non-CINs were collected separately into DMEM with 20% fetal bovine serum (FBS).

The cells were washed once then resuspended in buffer 1 (0.3 m sucrose, 60 mm KCl, 15 mm NaCl, 5 mm MgCl2, 0.1 mm EGTA, 15 mm Tris-HCl pH7.5, 1 mm DTT, 0.1 mm PMSF, 10 mm sodium butyrate, and 1× EDTA-free protease inhibitors). The cells were lysed by adding an equal volume of buffer 2 (buffer 1 + 0.4% v/v NP-40) and incubating on ice for 7 min. Nuclei were spun down and resuspended in 0.32 m sucrose, 50 mm Tris-HCl, pH7.5, 4 mm MgCl2, 1 mm CaCl2, 0.1 mm PMSF, and 10 mm sodium butyrate. Micrococcal nuclease was added for 6 min at 37°C to cleave chromatin into mononucleosome and dinucleosome fragments, and the reaction was stopped by addition of 20 mm EDTA. Nuclei were spun down, and the supernatant containing the soluble S1 chromatin fraction was saved. Nuclei were resuspended in dialysis buffer (1 mm Tris-HCl, pH 7.5, 0.2 mm EDTA, 0.2 mm PMSF, 10 mm Na butyrate, and 1× EDTA-free protease inhibitors) and rotated overnight at 4°C. After spinning down the supernatant containing the S2 chromatin fraction was combined with the S1 fraction for ChIP.

Histone modification ChIP was performed with antibodies against H3K27ac (Millipore catalog #05-1334, RRID:AB_1977244) or H3K27me3 (Millipore catalog #07-449, RRID:AB_310624). Chromatin from ∼80,000 nuclei was used for each condition. Antibody/chromatin complexes were purified using Dynabeads (Invitrogen) and washed extensively in Wash buffer 1 (50 mm Tris-HCl, pH 7.5, 10 mm EDTA, 125 mm NaCl, 5 mm sodium butyrate, and 0.1% v/v Tween 20), Wash buffer 2 (50 mm Tris-HCl, pH 7.5, 10 mm EDTA, 125 mm NaCl, 5 mm sodium butyrate, and 0.1% v/v NP-40), and TE. Complexes were eluted in 1% SDS, 10 mm sodium bicarbonate buffer at 65°C for 10 min. Eluted chromatin was treated with RNase (10 mg/200 ml reaction, 15 min at 37°C) and Proteinase K (10 mg/200 ml reaction, 60 min at 55°C) and cleaned using a ChIP DNA Clean & Concentrator kit (Zymo Research catalog #D5205). DNA libraries were prepared using the Ovation Ultralow DR Multiplex System (NuGEN Technologies) according to the manufacturer’s instructions. The libraries were size selected (180- to 300-bp range) using the BluePippin cassette system (Sage Science) and analyzed on a Bioanalyzer HS DNA chip (Agilent). The libraries were sequenced as single-ended 50-nucleotide reads on a HiSeq 4000 (Illumina).

ChIP-qPCR

H3K27ac enrichment was validated by ChIP-qPCR. H3K27ac ChIP was performed on sorted CIN and non-CIN cells as described above, and the purified chromatin was used for qPCR. We performed triplicate reactions in 10 µl volumes using the PerfeCTa SYBR Green 2× FastMix (Quanta Biosciences) on a CFX384 Real Time PCR system (Bio-Rad). The primers used are shown in Table 1. Triplicate Ct values were averaged, and fold-enrichment over matched input controls was calculated using the 2-ΔΔCt method (Livak and Schmittgen, 2001): Ct values were first normalized to a Dlx2 upstream non-conserved control region (Dlx2-non), i.e., ΔCt(sample) = Ct(target) – Ct(Dlx2-non); ChIP DNA samples were then normalized to the corresponding input DNA, i.e., ΔΔCt = ΔCt(ChIP) – ΔCt(input). Finally, the fold enrichment for each CIN sample was divided by fold enrichment for the matched non-CIN sample and log2 transformed, resulting in positive fold change values for CIN-enriched loci and negative values for non-CIN enriched loci.

View this table:
  • View inline
  • View popup
Table 1

Primer sequences used for ChIP-qPCR

View this table:
  • View inline
  • View popup
Table 2

Primer sequences used for viral vector cloning

Transcription factor (TF) ChIP

Embryonic day (E)13.5 basal ganglia TF ChIP-seq data were downloaded from the NCBI GEO database: DLX2 (accession: GSE124936; Lindtner et al., 2019), LHX6 (accession: GSE85704), and NKX2-1 (accession: GSE85704; Sandberg et al., 2016).

Sequence alignment

Histone modification ChIP and TF ChIP sequencing reads were aligned to the mouse genome (mm10) using BWA (v0.7.15) mem (Li, 2013). Histone modification ChIP replicates were aligned separately for further processing (see Histone modification ChIP peak calling and postprocessing). We also generated replicate-merged histone modification ChIP-seq alignments for visualization in the main figures. The two DLX2 replicates had already been merged (Lindtner et al., 2019); we merged the three LHX6 and NKX2-1 replicates during alignment for consistency with the DLX2 TF ChIP data. We used samtools (v1.10) to convert the resulting SAM files to BAM files as well as handle the MAPQ filtering and sorting (Li et al., 2009). We used a quality filter of 30.

To visualize our reads on the UCSC Genome Browser we generated bigwig files using deepTools (v3.4.3) bamCoverage (Ramírez et al., 2016). We used a bin size of 10 and a fragment extension length of 200. We also normalized to RPKM and ignored duplicate reads.

Histone modification ChIP peak calling and postprocessing

For our histone samples, treating replicates as separate samples, we identified significant peaks against matched input controls using MACS2 (v2.1.1) callpeak (Zhang et al., 2008). We disabled model-based peak calling and local significance testing and used a fixed fragment extension length of 200 bp. We used broad peak calling on both H3K27ac and H3K27me3 ChIP-seq samples with a broad-cutoff of 0.01 and a q value cutoff of 0.01.

We then tested the consistency of our replicates using the irreproducible discovery rate (IDR) R package (v1.2; Li et al., 2011), using the settings suggested in their package description (μ = 2.6, σ = 1.3, ρ = 0.8, p = 0.7; https://cran.r-project.org/web/packages/idr/idr.pdf). We used log10 fold-change scores from MACS2 as our input to calculate the IDR scores.

After confirming good distributions of signal-to-noise in our replicates, we then identified differential peaks between the CIN and non-CIN peak sets for each replicate using the MACS2 differential peak calling (bdgdiff) method with a minimum peak length of 150 bp (Zhang et al., 2008). This identified three sets of peaks for each replicate: CIN-enriched, non-CIN-enriched and common peaks.

Within each single bed file, peaks separated by 500 bp or less were merged to reduce variation in peak definitions because of coverage differences across the complex regions of enrichment often seen in histone modification ChIP-seq (unlike TF peaks, which tend to be discrete peaks) using bedtools merge (Quinlan and Hall, 2010). Finally, we defined high-confidence (HC) peak sets for each condition by finding the intersection between peak sets for the two replicates using bedtools intersect (Quinlan and Hall, 2010).

TF ChIP-seq postprocessing

For our TF ChIP-seq samples, we also used MACS2 (v2.1.1) callpeak (Zhang et al., 2008). We disabled model-based peak calling and local significance testing. We used a fixed fragment extension length of 200 bp and a q value cutoff of 0.01.

We merged our DLX2, LHX6, and NKX2-1 TF ChIP-seq peaks to form a set of the union of the three TFs and a set of the intersect of all three TFs. We then used the bedtools intersect method to locate the overlapping regions between our set of differentially called H3K27ac ChIP-seq peaks and our TF ChIP-seq sets.

ATAC-seq and conservation data

We obtained previously published chromatin accessibility data from ATAC-seq on purified CINs and PNs from GEO (accession: GSE63137; Mo et al., 2015). We used the provided bigwigs and called peaks without further processing. Bigwigs from replicate 1 are shown in the figures Cconserved elements predicted by phastCons (60-way vertebrate conservation) were obtained through the UCSC browser (https://genome.ucsc.edu).

Genomic annotation and motif analysis

We used HOMER (Heinz et al., 2010) to call motifs in our peak sets, using a size window of 200. We used the findMotifsGenome.pl and annotatePeaks.pl methods to acquire sets of enriched TF binding motifs for all our experiments and genomic annotations for our peaks. We also used the -genomeOntology setting to acquire genome ontology tables. TF binding motif fold-enrichment was calculated as the % of targets/% of background.

Gene-region associations and Gene Ontology (GO) analysis were performed on the HC peak sets with Genomic Regions Enrichment of Annotations Tool (GREAT) v4.0 (http://great.stanford.edu/public/html/index.php) using the default basal+extension association rule (McLean et al., 2010). GO term significance was based on default settings of binomial and hypergeometric false discovery rate (FDR) < 0.05 and region fold enrichment > 2.

Data availability

The histone modification ChIP-seq data generated in this study have been deposited in the NCBI GEO database and are accessible through GEO Series accession number GSE158428.

AAV methods

AAV production

The plasmid (p)AAV-I12b-Bg-ChR2-EYFP was reported previously (Cho et al., 2015). pAAV-I12b-Bg-ChR2-EYFP was modified, with the Arl4d and Dlgap1 enhancers replacing I12b, using Gibson assembly cloning (Gibson et al., 2009). Briefly, Arl4d_RE and Dlgap1_RE were amplified from mouse genomic DNA with the addition of overlap sequences using the primers in Table 2. The Beta-globin (Bg) minimal promoter region was amplified from pAAV-I12b-Bg-ChR2-EYFP. AgeI and SpeI cut sites were also introduced 3′ to the enhancers and 5′ to the Bg minimal promoter. Finally, pAAV-I12b-Bg-ChR2-EYFP was cut with MluI and BamHI, and the enhancer, Bg promoter and backbone were assembled for each plasmid using the In-Fusion HD Cloning kit (Takara Bio USA) according to the manufacturer’s instructions. Recombinant AAVs were packaged with serotype AAV5 by Virovek, Inc. AAV5-hSyn-ChR2-EYFP was obtained from the UNC Vector Core.

Virus injections

Adult mice were anesthetized using isoflurane (2.5% induction, 1.2–1.5% maintenance, in 95% oxygen) and placed in a stereotaxic frame (David Kopf Instruments). Body temperature was maintained using a heating pad. An incision was made to expose the skull for stereotaxic alignment using bregma and lambda as vertical references. For the Dlx5/6+/− mice for behavioral experiments, the scalp and periosteum were removed from the dorsal surface of the skull and scored with a scalpel to improve implant adhesion.

The mice were injected in the mPFC near the border between the prelimbic and infralimbic cortices (1.7 anterior-posterior, ±0.3 mediolateral, and −2.6 dorsoventral mm relative to bregma) with AAV-Arl4d-ChR2-EYFP, AAV-Dlgap1-ChR2-EYFP, AAV-I12b-ChR2-EYFP or AAV-hSyn-ChR2-EYFP (600–1000 nl per hemisphere). Viruses were infused at 100–150 nL/min through a 35-gauge, beveled injection needle (World Precision Instruments) using a microsyringe pump (World Precision Instruments, UMP3 UltraMicroPump). After infusion, the needle was kept at the injection site for 5–10 min and then slowly withdrawn. After surgery, mice were allowed to recover until ambulatory on a heated pad, then returned to their homecage. Approximately eight weeks after injection, wild-type mice were used for electrophysiological characterization, I12b-Cre;Ai14 mice were used for immunohistochemistry and cell counts, and Dlx5/6+/− mice were used for behavioral experiments.

Electrophysiology methods

Acute cortical slice preparation

Mice were anesthetized with an intraperitoneal injection of euthasol and transcardially perfused with an ice-cold cutting solution containing the following: 210 mm sucrose, 2.5 mm KCl, 1.25 mm NaH2PO4, 25 mm NaHCO3, 0.5 mm CaCl2, 7 mm MgCl2, and 7 mm dextrose (bubbled with 95% O2-5% CO2, pH ∼7.4). Approximately, three coronal slices (250 µm thick) of the brain containing the mPFC were obtained using a vibrating blade microtome (VT1200S, Leica Microsystems Inc.). Slices were allowed to recover at 34°C for 30 min followed by 30-min recovery at room temperature in a holding solution containing the following: 125 mm NaCl, 2.5 mm KCl, 1.25 mm NaH2PO4, 25 mm NaHCO3, 2 mm CaCl2, 2 mm MgCl2, 12.5 mm dextrose, 1.3 mm ascorbic acid, and 3 mm sodium pyruvate.

Whole-cell patch clamp recordings

Somatic whole-cell current clamp and voltage clamp recordings were obtained from submerged slices perfused in heated (32–34°C) artificial CSF (aCSF) containing the following: 125 mm NaCl, 3 mm KCl, 1.25 mm NaH2 PO4, 25 mm NaHCO3, 2 mm CaCl2, 1 mm MgCl2, and 12.5 mm dextrose (bubbled with 95% O2/5% CO2, pH ∼7.4). Neurons were visualized using DIC optics fitted with a 40× water-immersion objective (BX51WI, Olympus microscope). PNs and EYFP-expressing CINs located in layer 5 were targeted for patching. Patch electrodes (2–4 MΩ) were pulled from borosilicate capillary glass of external diameter 1 mm (Sutter Instruments) using a Flaming/Brown micropipette puller (model P-2000, Sutter Instruments). For current clamp recordings, electrodes were filled with an internal solution containing the following the following: 120 mm K-gluconate, 20 mm KCl, 10 mm HEPES, 4 mm NaCl, 7 mm K2-phosphocreatine, 0.3 mm Na-GTP, and 4 mm Mg-ATP (pH ∼7.3 adjusted with KOH). Biocytin (Vector Laboratories) was included (0.1–0.2%) for subsequent histologic processing. For voltage clamp recordings, the internal solution contained the following the following: 130 mm Cs-methanesulfonate, 10 mm CsCl, 10 mm HEPES, 4 mm NaCl, 7 mm phosphocreatine, 0.3 mm Na-GTP, 4 mm Mg-ATP, and 2 mm QX314-Br (pH ∼7.3 adjusted with CsOH).

Electrophysiology data were recorded using Multiclamp 700B amplifier (Molecular Devices). Voltages have not been corrected for measured liquid junction potential (∼8 mV). Upon successful transition to the whole-cell configuration, the neuron was given at least 5 min to stabilize before data were collected. During current clamp recordings, series resistance and pipette capacitance were appropriately compensated. Series resistance was usually 10–20 MΩ, and experiments were terminated if series resistances exceeded 25 MΩ.

Electrophysiology protocols and data analysis

All data analyses were performed using custom routines written in IGOR Pro (Wavemetrics).

Resting membrane potential (RMP) was measured as the membrane voltage measured in current clamp mode immediately after reaching the whole-cell configuration. Input resistance (Rin) was calculated as the slope of the linear fit of the voltage-current plot generated from a family of hyperpolarizing and depolarizing current injections (−50 to +20 pA, steps of 10 pA). Firing output was calculated as the number of action potentials (APs) fired in response to 800-ms-long depolarizing current injections (25–500 pA). Firing frequency was calculating as the number of APs fired per second. Rheobase was measured as the minimum current injection that elicited spiking. Firing traces in response to 50-pA current above the rheobase were used for analysis of single AP properties: AP threshold, maximum dV/dt (rate of rise of AP), AP amplitude, AP half-width and fast afterhyperpolarization (fAHP) amplitude. Threshold was defined as the voltage at which the value of third derivative of voltage (with respect to time) is maximum. AP amplitude was measured from threshold to peak, with the half-width measured at half this distance. fAHP was measured from the threshold to the negative voltage peak after the AP. Index of spike-frequency accommodation (SFA) was calculated as the ratio of the last inter-spike interval to the first inter-spike interval. Recorded CINs were classified as fast-spiking or regular-spiking based on electrophysiological properties. Specifically, CINs were classified as fast spiking if they satisfied at least three of the following four criteria: AP half-width < 0.5 ms, max firing frequency > 50 Hz, fAHP amplitude > 15 mV, or SFA < 2.

To measure optogenetically evoked spiking in EYFP+ CINs, and optogenetically evoked IPSCs in layer 5 PNs, we stimulated ChR2 using 5 ms long single light pulses (light power: 4 mW/mm2) generated by a lambda DG-4 high-speed optical switch with a 300 W Xenon lamp (Sutter Instruments) and an excitation filter centered around 470 nm, delivered to the slice through a 40× objective (Olympus).

Histology

Post hoc immunohistochemistry (IHC) and biocytin labeling

Patched slices were fixed overnight in 4% PFA then transferred to PBS. All incubations were conducted on a rotating platform at room temperature unless stated otherwise. The sections were rinsed twice in PBS, then blocked and permeabilized for 3 h in PBS with 10% FBS, 0.5% Triton X-100, and 0.05% sodium azide. Sections were immunostained overnight with primary antibodies: rabbit anti-PV (1:1000, Swant catalog #PV-27, RRID:AB_2631173) and rat anti-GFP (1:1000, Nacalai Tesque catalog #04404-84, RRID:AB_10013361) in PBS with 0.1% Triton X-100, 10% FBS and 0.025% sodium azide. Sections were washed 2 × 30 min in PBS with 0.25% Triton X-100, and 2 × 30 min in PBS. Goat anti-rabbit Alexa Fluor 488 secondary antibody (1:750, Invitrogen, catalog #A-11034) and Streptavidin-647 (1:500, Invitrogen, catalog #S-32357) were added with Hoescht 33 342 nuclear counterstain (1:2000, Invitrogen, catalog #H3570) for 4–6 h at room temperature, then overnight at 4°C. After washing 2 × 30 min in PBS with 0.25% Triton X-100, and 2 × 30 min in PBS, sections were mounted on Superfrost Plus slides and coverslipped with Fluorescence Mounting Medium (DAKO catalog #S3023).

EYFP and CIN marker colocalization

AAV-injected I12b-Cre; Ai14 mice were anesthetized with an intraperitonial injection of ketamine/xylazine and transcardially perfused with ice-cold PBS followed by fixative solution (4% PFA in 0.24 m sodium phosphate buffer, pH 7.2). Brains were dissected out and postfixed in fixative solution overnight at 4°C, followed by cryopreservation for 24 h in 15% sucrose in PBS and 24 h in 30% sucrose in PBS. Coronal sections of 40-µm thickness were obtained using a freezing microtome (SM2000 R, Leica Microsystems Inc.) and stored in sectioning solution (40% PBS, 30% glycerol, 30% ethylene glycol) at −20°C until use.

Sections were immunostained using a standard floating section IHC protocol with gentle agitation on a rotary shaking platform. Briefly, sections were rinsed in PBS, blocked for 1 h at room temperature in PBST (PBS + 0.25% Triton X-100) with 10% FBS, and incubated with primary antibody diluted in PBST + 10% FBS overnight at 4°C. Primary antibodies used were: rabbit anti-PV (1:1000, Swant catalog #PV-27, RRID:AB_2631173), rat anti-SST (1:400, Millipore catalog #MAB354, RRID:AB_2255365), rabbit anti-VIP (1:250, ImmunoStar catalog #20077, RRID:AB_572270), rat anti-GFP (1:1000, Nacalai Tesque catalog #04404-84, RRID:AB_10013361), and rabbit anti-GFP (1:5000, Abcam catalog #ab6556, RRID:AB_305564). The following day, sections were washed 2 × 15 min in PBST and 2 × 15 min in PBS, incubated for 2 h at room temperature with Alexa Fluor-conjugated secondary antibodies (1:750, Invitrogen). Finally, sections were washed 2 × 15 min in PBST and 2 × 15 min in PBS, mounted on Superfrost Plus slides and coverslipped with Fluorescence Mounting Medium (DAKO catalog #S3023).

Imaging

Low-magnification epifluorescent images were taken using a Coolsnap camera (Photometrics) mounted on a Nikon Eclipse 80i microscope using NIS Elements acquisition software (Nikon). Confocal images were taken with 20× air and 40× oil objectives on an Andor Borealis CSU-W1 spinning disk confocal mounted on a Nikon Ti Microscope and captured with an Andor Zyla sCMOS camera and Micro-Manager software (Open Imaging). The raw images were preprocessed with ImageJ software (v2.0.0) to adjust brightness/contrast and convert to eight-bit RGB. Confocal z-stacks were stitched laterally to create composite stacks using the Grid/Collection stitching ImageJ plugin with linear blending (Preibisch et al., 2009).

Cell counts

Cell counting was performed in Adobe Photoshop CC on a single stitched 20× confocal image plane. EYFP+ and marker+ cells were counted in the mPFC at the rostrocaudal levels shown in Extended Data Figure 4-1. For each virus, at least three rostrocaudal sections were counted for each of at least two animals. EYFP or marker-positive cells were counted using only the corresponding single-color channel before determining colocalization of markers among these cells. Counts are presented as % (EYFP+ marker+ cells/total EYFP+ cells). Counts normalized to CIN subtype abundance are also presented, i.e., (EYFP+ marker+ cells/total EYFP+ cells) ÷ (tdTomato+ marker+ cells/total tdTomato+ cells). The normalized counts for AAV-hSyn-ChR2-EYFP injected animals include EYFP+ CINs only (co-expressing tdTomato), rather than the total EYFP+ population.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

ChIP-seq reveals differential enrichment of histone modifications in P2 CINs and non-CINs. A, Schematic showing the histone modification ChIP-seq workflow using CINs and non-CINs dissected and FACS-purified from P2 Gad67-GFP neocortex. B–E, Enrichment of H3K27ac (active mark, cyan and teal) and H3K27me3 (repressive mark, red and dark red) around genes expressed in CINs (Gad2), PNs (Emx2) oligodendrocyte-lineage cells (Olig1, Olig2), and astrocytes (Aldh1l1). Gad2En1 is a validated enhancer with activity in adult CINs (yellow highlight; Pla et al., 2017). F, G, H3K27ac enrichment is similar in CINs and non-CINs across the gene body for Id4 and Nr2f1 but differentially enriched at enhancers with validated activity in the basal ganglia (green highlights) or cortex (red highlights). Whole-mount transient transgenics from the VISTA browser show the pattern of enhancer activity at E11.5 (https://enhancer.lbl.gov; Visel et al., 2013). ChIP-qPCR validation of H3K27ac enrichment at highlighted enhancers is shown in Extended Data Figure 1-1. HC called peaks are shown as solid bars under the corresponding bigwig tracks.

Extended Data Figure 1-1

Validation of histone modification ChIP-seq. A, IDR analysis showing replicate reproducibility for histone modification ChIP-seq. B, Table detailing the number of called peaks from histone modification ChIP-seq experiments by biological replicate. HC peaks were those called for both replicates. C, Validation of H3K27ac enrichment at REs and pREs by ChIP-qPCR, expressed as log2 fold change of CIN enrichment over non-CIN enrichment. Fold changes were not significantly different to zero by one-sample t test (p > 0.05). Download Figure 1-1, EPS file.

Rule shift methods

Surgery

Male Dlx5/6+/− adult mice (15 weeks old at time of experiment) received bilateral injections of AAV-Arl4d-ChR2-EYFP into the mPFC as described above. After injection of virus, a 200/240 mm (core/outer) diameter, NA = 0.22, dual fiber-optic cannula (Doric Lenses, DFC_200/240–0.22_2.3 mm_GS0.7_FLT) was slowly inserted into mPFC until the tip of the fiber reached a dorsoventral depth of −2.25 mm relative to bregma. Implants were affixed onto the skull using Metabond Quick Adhesive Cement (Parkell). We waited at least eight weeks after injection before behavioral experiments to allow for virus expression.

Rule shift task

For the cognitive flexibility task, mice were singly-housed and habituated to a reverse light/dark cycle and food intake was restricted until the mice reached 80–85% of the ad libitum feeding weight. At the start of each trial, the mouse was placed in its home cage to explore two bowls, each containing one odor and one digging medium, until it dug in one bowl, signifying a choice. As soon as a mouse began to dig in one bowl, the other bowl was removed, so there was no opportunity for “bowl switching.” The bait was a piece of a peanut butter chip (∼5–10 mg in weight) and the cues, either olfactory (odor) or somatosensory and visual (texture of the digging medium which hides the bait), were altered and counterbalanced. All cues were presented in small animal food bowls (All Living Things Nibble bowls, PetSmart) that were identical in color and size. Digging media were mixed with the odor (0.01% by volume) and peanut butter chip powder (0.1% by volume). All odors were ground dried spices (McCormick garlic and coriander), and unscented digging media (Mosser Lee White Sand Soil Cover, Natural Integrity Clumping Clay cat litter).

After mice reached their target weight, they underwent 1 d of habituation. On this day, mice were given 10 consecutive trials with the baited food bowl to ascertain that they could reliably dig and that only one bowl contained food reward. All mice were able to dig for the reward. Then, on days 1 and 2, mice performed the task (this was the testing done for experiments). After the task was done for the day, the bowls were filled with different odor-medium combinations and food was evenly distributed among these bowls and given to the mouse so that the mouse would disregard any associations made earlier in the day.

Mice were tested through a series of trials. The determination of which odor and medium to pair and which side (left or right) contained the baited bowl was randomized (subject to the requirement that the same combination of pairing and side did not repeat on more than three consecutive trials) using http://random.org. On each trial, while the particular odor-medium combination present in each of the two bowls may have changed, the particular stimulus (e.g., a particular odor or medium) that signaled the presence of food reward remained constant over each portion of the task (initial association and rule shift). If the initial association paired a specific odor with food reward, then the digging medium would be considered the irrelevant dimension. The mouse is considered to have learned the initial association between stimulus and reward if it makes eight correct choices during 10 consecutive trials. Each portion of the task ended when the mouse met this criterion. Following the initial association, the rule shift portion of the task began, and the particular stimulus associated with reward underwent an extradimensional shift. For example, if an odor had been associated with reward during the initial association, then a digging medium was associated with reward during the rule shift portion of the task. The mouse is considered to have learned this extradimensional rule shift if it makes eight correct choices during 10 consecutive trials. When a mouse makes a correct choice on a trial, it is allowed to consume the food reward before the next trial. Following correct trials, the mouse is transferred from the home cage to a holding cage for ∼10 s while the new bowls were set up (intertrial interval). After making an error on a trial, a mouse was transferred to the holding cage for ∼2 min (intertrial interval). All animals performed the initial association in a similar number of trials (average: 10–15 trials).

In vivo optogenetic stimulation

A 473 nm blue laser (OEM Laser Systems) was coupled to the dual fiber-optic cannula (Doric Lenses) through a 200 mm diameter dual fiber-optic patchcord with guiding socket (Doric Lenses) and 1 × 2 intensity division fiber-optic rotary joint (Doric Lenses), and adjusted such that the final light power was ∼0.5 mW total, summed across both fibers and averaged over light pulses and the intervening periods. A function generator (Agilent 33500B Series waveform Generator) connected to the laser generated a 40 Hz train of 5 ms pulses.

Light stimulation began once mice reached the 80% criterion during the initial association portion of the task. Mice then performed three additional initial association trials before the rule shift portion of the task began.

Statistical analysis

Statistical analyses were performed using Prism 7 (GraphPad) and detailed in the corresponding figure legends. Student’s t tests were used to make single-variable comparisons. For multivariate comparisons, we used one-way ANOVA with Tukey’s corrected p values for multiple comparisons. Chi-square or Fisher’s exact tests were used to test for significance among contingency-type data; *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. Comparisons with no asterisk had p > 0.05 and were considered not significant.

Results

Identification of candidate CIN gene REs

To identify putative REs (pREs) that could be used to drive expression in immature CINs, and potentially in mature CINs, we used genome-wide epigenetic profiling to identify regions of the genome that have signatures of an active transcriptional state (Fig. 1A). We used FACS to separate GFP+ immature CINs and GFP- non-CIN cells (enriched for glutamatergic PNs and glia) from acutely dissected and dissociated cortical tissue from postnatal day (P)2 Gad67-GFP heterozygote pups. Chromatin was extracted from the purified cell populations and processed for native histone modification ChIP with antibodies raised against H3K27ac, a mark of active REs, and H3K27me3, a mark of poised and repressive REs, followed by next-generation DNA sequencing (Creyghton et al., 2010; Rada-Iglesias et al., 2011). Genomic regions with enrichment (“peaks”) of H3K27ac and H3K27me3 were called by comparing the ChIP samples with an input control. IDR analysis confirmed high concordance between all pairs of replicates (Extended Data Fig. 1-1A; Li et al., 2011).

We thus defined high-confidence (HC) peaks as those present in both biological replicates. A total of 35,768 HC H3K27ac peaks were called for the CIN samples and 34,953 for the non-CINs. These peaks represent the pREs that could be driving gene expression in these cells, including proximal REs (promoters) and distal REs (enhancers). We then performed differential peak calling to identify regions with significantly more H3K27ac in CINs versus non-CINs, and vice versa. By this method, we identified 2937 HC CIN-specific pREs and 718 HC non-CIN pREs; 25,556 HC peaks were common to both populations. The HC peak sets were used for all downstream analyses. Peak numbers by replicate can be found in Extended Data Figure 1-1B.

A total of 14,915 HC H3K27me3 peaks were called for CIN samples and 17,550 were called for non-CINs. These peaks represent putatively poised and/or repressed REs, which would be expected to be associated with lower expression of their target genes. By differential peak calling, we identified 168 HC regions with more H3K27me3 in CINs than non-CINs and 285 HC regions with higher H3K27me3 in non-CINs. Peak numbers by replicate can be found in Extended Data Figure 1-1B.

We validated our histone ChIP-seq data first by confirming patterns of histone modification enrichment at known marker gene loci. As expected, CIN-relevant genes (e.g., Dlx1, Dlx2, Gad1, and Gad2) had significant enrichment of H3K27ac (greater activation) and depletion of H3K27me3 (less repression) in CINs compared with non-CINs (Figs. 1B, 2A; Extended Data Fig. 2-1B). Conversely, markers for PNs (Emx2), oligodendrocytes (Olig1 and Olig2) and astrocytes (Aldh1l1) showed greater enrichment of H3K27ac over the gene body and proximal sequences in the non-CIN samples (Fig. 1C–E). Statistically significant differential enrichment of H3K27ac was also found at validated enhancers (distal REs) active in CINs or their progenitors in the basal ganglia. For example, H3K27ac was enriched specifically in CINs at Gad2En1, hs1175 near Id4, and hs1060 near Nr2f1 (Fig. 1B,F,G; Visel et al., 2013; Silberberg et al., 2016; Pla et al., 2017). Conversely, H3K27ac was enriched in one of the two non-CIN samples at hs1533 near Id4, and hs1172 near Nr2f1 (Fig. 1F,G), which have cortical activity at E11.5 (Visel et al., 2013; Pattabiraman et al., 2014).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Motif analysis and genomic annotation of regions enriched for histone H3K27ac. A, De novo motif analysis of CIN-enriched H3K27ac peaks by HOMER, showing the top five enriched motifs, their best-match known motifs, fold-enrichment (Fold-en) over background, and significance of enrichment. De novo motif analyses of other peak sets are shown in Extended Data Figure 2–1A–C. B, Heatmap showing the enrichment of genomic annotations (curated list) generated by HOMER for the different sets of H3K27ac peaks compared with background. The full list of annotations and enrichment values are shown in Extended Data Figure 2-1D. CIN-enriched_unionTFs peaks are CIN-enriched peaks that overlap with at least one DLX2, LHX6, or NKX2-1 ChIP binding site. CIN-enriched_intersectTFs peaks are CIN-enriched peaks that overlap with DLX2, LHX6, and NKX2-1 binding sites. C, Proportion of peaks in each set associated with genomic annotations generated by HOMER. GO analyses of genes associated with the CIN-enriched, non-CIN-enriched, and common peaks are shown in Extended Data Figures 2-1E,F, 2-2. TTS: transcription termination site, TSS: transcription start site, UTR: untranslated region.

Extended Data Figure 2-1

Motif analysis and genomic annotation of regions enriched for histone H3K27ac. A–C, De novo motif analysis by HOMER, showing the top five enriched motifs, their best-match known motifs, fold-enrichment over background, and significance of enrichment. D, Genomic ontology generated by HOMER, showing log ratio enrichment over background. Significantly overrepresented (or underrepresented) annotations (p < 0.01) are shown in bold. E, F, GO analysis of gene-associated H3K27ac-enriched regions by GREAT for CIN-specific and non-CIN peaks. The top ten significantly enriched biological process terms are shown for each peak set ranked by binomial p value. Relevant terms are highlighted in cyan. The full list of enriched GO biological process terms can be found in Extended Data Figure 2-2. Download Figure 2-1, EPS file.

Extended Data Figure 2-2

Full lists of significantly enriched GO terms from GREAT analysis for the HC common, CIN-specific, and non-CIN H3K27ac peaks. Download Figure 2-2, XLSX file.

We then confirmed differential enrichment of H3K27ac at these REs by ChIP-qPCR on sorted CIN and non-CIN samples from P2 Gad67-GFP cortices. Gad2En1, hs1175, and hs1060 showed a trend for enrichment in the CIN samples, while hs1533 and hs1172 were more enriched in non-CINs as expected (Extended Data Fig. 1-1C). However, these enrichments were not significant by one-sample t test (p > 0.05).

We hypothesized that REs specifically active in CINs would be bound by TFs with known roles in CIN development. Therefore, we performed de novo motif discovery with HOMER to identify TF binding motifs enriched in the CIN-enriched H3K27ac peak set and compared these to motifs enriched in the non-CIN-enriched and common peaks (Fig. 2A; Extended Data Fig. 2-1A,B).

The top CIN-enriched motifs included SOX motifs, a DLX motif and a MEIS motif and were highly significant (p ≤ 1e-21; Fig. 2A). DLX TFs are central to GABAergic cell development (Hu et al., 2017; Lindtner et al., 2019). MEIS1 is transcriptionally downstream of DLX1/2 and ARX (another important CIN TF) and is expressed in CINs migrating to the cortex (Long et al., 2009; Friocourt and Parnavelas, 2011). The enrichment of these binding motifs suggests that the CIN-enriched peaks represent functional REs that may regulate CIN gene expression. The non-CIN-enriched de novo motifs had low p values (≤1e-11) and were flagged as possible false positives by HOMER (Extended Data Fig. 2-1A). They were also found in <5% of the target sequences, compared with 17–43% of target sequences for the top five CIN-enriched motifs. This suggests a lack of any strong consensus TF binding motifs among the non-CIN-enriched peak set, which may be a function of the cell type heterogeneity of the non-CIN sample.

To determine whether CIN-enriched and non-CIN-enriched peaks are associated with particular functional regions of the genome, we also performed genomic annotation with HOMER. All the peak sets (CIN, non-CIN, CIN-enriched, non-CIN-enriched and common peaks) showed enrichment over background at promoters, 5′ untranslated regions (5′UTRs) and introns, as expected for H3K27ac ChIP (Fig. 2B; Extended Data Fig. 2-1D).

Compared with common (not differentially enriched) peaks, both CIN-enriched and non-CIN-enriched peaks were less often associated with promoter regions/transcription start sites (TSSs), exons and 5′UTRs, and more often associated with intergenic regions (Fig. 2C). This supports the notion that chromatin activation state is more dynamically regulated at enhancers than promoters, and, therefore, enhancers may show more cell type-specific activity than promoters.

We then used the GREAT to assign H3K27ac peaks to nearby genes they may regulate (McLean et al., 2010). GO analysis of gene-associated H3K27ac-enriched regions with GREAT showed significant enrichment of relevant biological process GO terms for each population. The top ten GO terms for CIN-specific pREs included terms related to forebrain neurogenesis and specifically GABAergic CIN differentiation, while terms further down the list were also associated with neuronal migration and axonogenesis (Extended Data Figs. 2-1E, 2-2). HC non-CIN pREs were strongly associated with glial development (Extended Data Fig. 2-1F). The significantly enriched GO terms for the common peak set were overwhelmed by terms related to generic cellular processes such as transcription and translation (Extended Data Fig. 2-2). However, many genes related to synaptogenesis and pan-neuronal development and function were associated with multiple common H3K27ac peaks (e.g., Tubb2a/b, Nlgn1/2/3, Ctnnb1, Rbfox2, and Camk2d).

To refine our list of candidate pREs, we leveraged previously published TF ChIP-seq data from embryonic basal ganglia (Sandberg et al., 2016; Lindtner et al., 2019). DLX2, LHX6, and NKX2-1 are TFs important for CIN development, whose expression within the telencephalon is primarily limited to the GABAergic lineages generated in the basal ganglia anlage (Kessaris et al., 2014; Hu et al., 2017). We hypothesized that DLX2, LHX6, and NKX2-1 binding at pREs increases the probability that the RE will be active in maturing GABAergic lineage neurons including CINs.

A total of 40% of CIN-specific pREs (1163 out of 2937) overlapped with at least one called DLX2, LHX6, or NKX2-1 peak (CIN-enriched_unionTF peak set). The most significantly enriched de novo motif found at CIN-enriched_unionTF pREs (TAATTACVVS) was a close match for known DLX1/2/5/6 motifs, and also very similar to published LHX6 consensus motifs (Extended Data Fig. 2-1C; Sandberg et al., 2016).

Sequence conservation across species is another predictor of functional REs (Nord et al., 2015). We found that 72% (2115 out of 2937) of our HC CIN-specific pREs overlapped with predicted vertebrate conserved elements identified by phastCons (obtained through the UCSC browser). Furthermore, chromatin accessibility is a prerequisite for RE activation and gene expression. We therefore integrated published chromatin accessibility data from ATAC-seq on sorted adult cortical cells to further inform our selection of REs to test (Mo et al., 2015).

As a proof of principle, we looked at the I12b enhancer, which is located between Dlx1 and Dlx2. I12b has well-characterized activity in the ganglionic eminences and their derivates, including CINs, in transgenic mice and when placed in viral vectors (Potter et al., 2009; Arguello et al., 2013; Lee et al., 2014b; Cho et al., 2015; Vogt et al., 2015). We found that the entire Dlx1, Dlx2 and I12b genomic region was highly enriched for H3K27ac in the purified P2 CINs, but not in the non-CIN population (Fig. 3A). I12b was also bound by DLX2, LHX6, and NKX2-1 in E13.5 basal ganglia cells based on TF ChIP (Fig. 3A). Additionally, ATAC-seq demonstrated that the region was accessible in adult PV+ and VIP+ CINs, but not in purified excitatory CAMK2A+ PNs (Fig. 3A; Mo et al., 2015).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Integrated epigenetic signatures of candidate CIN pREs. A, The I12b RE region (yellow highlight) is highly enriched for H3K27ac in CINs compared with the non-CINs. It also has Dlx2, LHX6, and NKX2-1 TF binding peaks in E13.5 basal ganglia (BG) and shows enhanced accessibility by ATAC-seq in P30 PV+ and VIP+ CINs, but not in CAMK2A+ PNs. B, The Arl4d pRE (yellow highlight) is highly enriched for H3K27ac in CINs compared with non-CINs, has DLX2, LHX6, and NKX2-1 TF binding peaks and shows enhanced accessibility by ATAC-seq in P30 PV+ CINs but not VIP+ CINs or CAMK2A+ PNs. It also overlaps a vertebrate conserved element. C, The Dlgap1 intronic pRE (yellow highlight) shows greater enrichment of H3K27ac in CINs than in non-CINs and modest binding peaks for DLX2, LHX6, and NKX2-1. The Dlgap1 pRE overlaps an ATAC peak in P30 PV+ CINs, but not in VIP+ CINs or PNs, and a vertebrate conserved element. D, The Syn1 promoter region (red highlight) has comparable H3K27ac enrichment in CINs and non-CINs and is highly accessible in P30 PV+ CINs, VIP+ CINs and CAMK2A+ PNs. The highly conserved homologous human region drives AAV expression in many PNs and some CINs (see Fig. 5). HC called peaks are shown as solid bars under the corresponding bigwig track. Additional pREs are shown in Extended Data Figure 3-1. Coverage and called peaks separated by replicate for the Arl4d and Dlgap1 loci are shown in Extended Data Figure 3-2. E, In situ hybridization in adult cortex shows CIN-specific staining of Dlx1. Arl4d and Dlgap1 mRNA appear to be expressed in CINs [especially obvious in neocortical layer 1 (red arrows) and hippocampal stratum radiatum (blue arrows)] and many PNs. Images are taken from the Allen Brain Atlas (https://www.brain-map.org; Lein et al., 2007). Scale bars: 500 µm.

Extended Data Figure 3-1

Epigenetic signatures of additional candidate CIN pREs. A, The Arx locus contains several validated enhancers. Enhancers with activity in the basal ganglia (green highlights) have HC H3K27ac CIN-specific peaks, DLX2, LHX6, and NKX2-1 TF binding peaks, and enhanced accessibility by ATAC-seq in P30 PV+ and/or VIP+ CINs but not in CAMK2A+ PNs. Enhancers active in the cortex (red highlights) or outside the telencephalon (blue highlights) lack this particular signature: hs123 is also accessibly by ATAC in CAMK2A+ PNs; hs122 has H3K27ac enrichment in both CINs and non-CINs. Sections showing transient transgenic enhancer activity at E11.5 are taken from the VISTA Enhancer Browser (https://enhancer.lbl.gov; Visel et al., 2013). B–F, Epigenetic signatures of CIN-specific pREs (yellow highlights) near the CIN genes Gad1, Slc32a1 (Vgat), Htr3a, Maf, and Sst. Download Figure 3-1, TIF file.

Extended Data Figure 3-2

H3K27ac enrichment at pREs separated by replicate. Browser tracks showing merged coverage and HC called peaks (white background), and coverage and called peaks by replicate (grey background) for CINs and non-CINs. Differential called peaks (CIN-enriched or non-CIN-enriched) are shown at the top. A, The Arl4d pRE (yellow highlight) shows strong H3K27ac enrichment in CINs in both replicates and is called as CIN-enriched in both replicates. B, The Dlgap1 pRE (yellow highlight) shows modest H3K27ac enrichment, with a called CIN peak and differential CIN-enriched peak in replicate 1. Download Figure 3-2, EPS file.

Based on the positive results with the I12b locus, we screened other loci for similar epigenetic signatures. The Arx locus contains several REs with validated enhancer activity in the developing telencephalon at E11.5: hs119 and hs121 are active in the basal ganglia, while hs122 and hs123 are active in cortical progenitors (VISTA REs, https://enhancer.lbl.gov; Visel et al., 2013). hs121 (also known as UAS3) also remains active in immature and adult CINs (Colasante et al., 2008). All four enhancers had enrichment of DLX2 TF binding in E13.5 basal ganglia, but only the REs active in the basal ganglia and immature CINs had both HC CIN-specific H3K27ac peaks at P2 and chromatin accessibility in adult PV+ or VIP+ CINs but not in CAMK2A+ PNs (Extended Data Fig. 3-1A). This highlights the utility of integrating multiple epigenetic marks to identify strong CIN-specific pRE candidates.

Among the many possibilities we focused on pREs located near two genes: Arl4d and Dlgap1. Both Arl4d and Dlgap1 mRNA are expressed in the embryonic basal ganglia (the ganglionic eminences), and their expression is downregulated on loss of Nkx2-1 expression in the medial ganglionic eminence (Sandberg et al., 2016). Arl4d expression is also reduced in the basal ganglia of Dlx1/2 knock-outs (Lindtner et al., 2019). In the adult, Arl4d is expressed in many cells throughout the cortex, including CINs most visible in neocortical layer 1 and the hippocampus (Fig. 3E). DLGAP1 is expressed at excitatory postsynaptic densities in neurons throughout the forebrain (Fig. 3E; Rasmussen et al., 2017). Although both of these genes show widespread expression in forebrain neurons, enhancers are often active in a restricted subset of their target gene’s expression domain (Visel et al., 2013).

The Arl4d pRE is located ∼30 kb downstream of the Arl4d TSS. H3K27ac enrichment was significantly higher in P2 CINs compared with non-CINs (Fig. 3B). The locus also contains TF binding peaks for DLX2, LHX6, and NKX2-1 and had ATAC-seq peaks in adult PV CINs but not in VIP CINs or PNs (Fig. 3B). The Dlgap1 pRE is located within its first intron. This locus showed modest enrichment of H3K27ac in P2 CINs (called in one replicate, see Extended Data Fig. 3-2B), as well as binding peaks for DLX2, LHX6, and NKX2-1 in embryonic basal ganglia (Fig. 3C). Both the Arl4d and Dlgap1 pREs also overlapped with vertebrate conserved elements predicted by phastCons, suggesting their activity could translate across species.

H3K27ac enrichment was further validated by ChIP-qPCR. I12b, the Arl4d pRE, and the Dlgap1 pRE all showed relative enrichment in CINs versus non-CINs (Extended Data Fig. 1-1C). The magnitude of enrichment was similar for the three REs, although it was not significant by one-sample t test (p > 0.05).

Many further candidate CIN pREs were identified based on these integrated epigenetic data. Several examples found near the CIN-expressed genes Gad1, Slc32a1 (Vgat), Htr3a, Maf, and Sst are shown in Extended Data Figure 3-1B–F. The epigenetic signature at the mouse Syn1 promoter, homologous to the human region used as a control RE active in PNs (see CIN RE AAVs preferentially label PV-positive neurons in the mPFC), did not contain any HC H3K27ac peaks or called TF binding peaks but was highly accessible in PV and VIP CINs as well as in PNs by ATAC-seq (Fig. 3D).

Physiologic properties of neurons targeted by the RE AAVs

CINs can be classified into diverse subclasses based on molecular markers, morphology, circuit function and electrophysiological properties (Kepecs and Fishell, 2014; Tremblay et al., 2016). To test whether our candidate pREs (Arl4d_RE and Dlgap1_RE) can be used to generate viral tools for the observation and manipulation of physiologically distinct subpopulations of CINs, we made AAV constructs using these pREs to drive expression of a channelrhodopsin-enhanced YFP fusion protein (ChR2-EYFP). The resulting AAVs (AAV-Arl4d-ChR2-EYFP and AAV-Dlgap1-ChR2-EYFP) and a previously characterized RE AAV (AAV-I12b-ChR2-EYFP; Cho et al., 2015) were injected into the deep layers of the mPFC.

We obtained current clamp recordings from transduced layer 5 neurons, identified by EYFP expression (Extended Data Fig. 4-1A). Reliable short latency (<10 ms) firing in response to flashes of blue light (470 nm) was used to validate AAV-mediated ChR2 expression in the recorded EYFP+ CINs (Extended Data Fig. 4-1E). Notably, the RE AAVs specifically targeted CINs and not PNs in layer 5 of the mPFC, as optogenetic stimulation of RE AAV-labeled neurons elicited robust IPSCs in EYFP negative layer 5 PNs, which were completely blocked by bath application of 10 μm gabazine (GABAA receptor antagonist; Extended Data Fig. 4-1F).

To determine whether different RE AAVs differentially target CINs with distinct physiological properties, recorded cells were classified as either fast spiking or regular spiking based on their firing output and AP properties (Extended Data Fig. 4-1B–D). This analysis showed that a large proportion of CINs targeted by AAV-Arl4d-ChR2-EYFP had fast-spiking properties (11/13 neurons; Fig. 4A). Specifically, the fast-spiking CINs targeted by AAV-Arl4d-ChR2-EYFP had shorter AP half-width, large fAHP amplitude, low SFA and high peak firing frequency in comparison to the smaller proportion of EYFP+ regular-spiking CINs (Extended Data Fig. 4-1E–I). Similar to the transduction specificity of AAV-Arl4d-ChR2-EYFP, a large proportion of AAV-I12b-ChR2-EYFP labeled CINs also had fast-spiking physiological properties (9/12 neurons; Fig. 4). In contrast, AAV-Dlgap1-ChR2-EYFP labeled a larger proportion of regular-spiking CINs (7/12 neurons) than fast-spiking CINs (5/12 neurons; Fig. 4). Overall, the proportions of regular versus fast-spiking CINs transduced by AAV-Arl4d-ChR2-EYFP versus AAV-Dlgap1-ChR2-EYFP were significantly different (p = 0.0414, Fisher’s exact test).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Electrophysiological properties of CINs targeted by AAV-Arl4d, AAV-Dlgap1, and AAV-I12b. A, Pie charts showing the percentages of fast-spiking (FS) and regular-spiking (RS) CINs targeted by the different pREs. AAV-Arl4d-ChR2-EYFP: 11/13 FS, 2/13 RS; AAV-Dlgap1-ChR2-EYFP: 5/12 FS, 7/12 RS; AAV-I12b-ChR2-EYFP: 9/12 FS, 3/12 RS. B–I, Subthreshold, AP, and firing properties of the layer 5 CINs targeted by the enhancer AAVs. Black dots denote FS CINs, and red dots denote RS CINs. Data are shown as mean ± SEM. Additional methodology and post hoc immunolabeling of recorded cells is shown in Extended Data Figure 4-1.

Extended Data Figure 4-1

Electrophysiological characterization of layer 5 CINs targeted by RE-AAVs. A, DIC and fluorescent images showing that EYFP+ layer 5 CINs in mPFC were targeted during patch clamp recordings. Scale bar: 25 µm. B, Example voltage traces in response to subthreshold current injections (–50 to +50 pA, top) and to suprathreshold current injections (250 and 375 pA) in an EYFP+ fast-spiking (FS) CIN. C, Same as B for an EYFP+ regular-spiking (RS) CIN. D, Example voltage traces showing single APs recorded from EYFP+ FS and RS CINs. E, left, Recording configuration showing that current clamp recordings were obtained from layer 5 EYFP+ CINs. Right, APs were reliably elicited in EYFP+ CINs in response to blue light (470 nm) flashes during patch clamp recordings. F, left, Recording configuration showing that voltage clamp recordings (+10 mV) were obtained from EYFP-negative layer 5 PNs. IPSCs (black trace) were recorded in PNs in response to blue light flashes. Currents were completely blocked by application of gabazine (10 µm). G, Representative images of post hoc labeling for biocytin, PV, and EYFP in recorded cells for each enhancer AAV injection. Arrows indicate PV+ recorded cells and open arrowhead indicates a PV– recorded cell. Scale bar: 25 µm. H, Pie charts showing the percentages of biocytin-labeled FS and RS CINs expressing PV by post hoc labeling. Some cells were not recovered after patching. Number of cells of each type is indicated. Download Figure 4-1, tif file.

Since fast-spiking physiological properties are associated with CINs that express PV (Kawaguchi et al., 1987; Rudy et al., 2011), we immunostained for PV expression in a subset of recorded EYFP+ neurons (Extended Data Fig. 4-1G,H). In line with the increased targeting of fast-spiking CINs by AAV-Arl4d-ChR2-EYFP and AAV-I12b-ChR2-EYFP, a large proportion of EYFP+ CINs transduced with these enhancer AAVs expressed PV (Arl4d: 5/5 labeled neurons; I12b: 8/12 labeled neurons). Surprisingly, many AAV-Dlgap1-ChR2-EYFP-labeled neurons also expressed PV (7/8 neurons; Fig. 4A; Extended Data Fig. 4-1G,H).

CIN RE AAVs preferentially label PV-positive neurons in the mPFC

To further characterize the molecular identity of the RE AAV-labeled neurons, we injected the AAV-RE-ChR2-EYFP viruses into the mPFC of adult I12b-Cre;Ai14 mice, in which the majority of CINs are labeled by tdTomato reporter expression (Extended Data Fig. 5-1B; Potter et al., 2009). As an additional control virus, we used AAV-hSyn-ChR2-EYFP, which labels both PNs and CINs. After seven to eight weeks, we performed immunohistochemistry (IHC) for EYFP and CIN subtype markers SST, PV, and VIP (Fig. 5A). Marker and EYFP colocalization was assessed in the mPFC of 2–3 bilaterally injected animals, compiled from three to four rostrocaudal sections from each animal (Extended Data Fig. 5-1A).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

IHC characterization of CINs targeted by RE-AAVs injected into adult I12b-Cre; Ai14 mPFC. A, Single confocal planes showing colocalization of EYFP from the RE-AAVs with endogenous tdTomato and CIN subtype markers PV, SST, and VIP within the mPFC. tdTomato expression from the I12b-Cre; Ai14 reporter line marks most CINs. Dotted white lines indicate border between cortical layer 6 and corpus callosum, and approximate layer markings are shown in the low-magnification tdTomato panels (top). Yellow boxes indicate areas shown at high magnification. Arrows indicate EYFP+, marker+ cells; open arrowheads indicate EYFP+ only cells. Scale bars: 50 µm. (tdTomato and PV images are from the same section for each animal.) B, Quantification of the % of EYFP+ cells targeted by the RE AAVs that express tdTomato and the CIN subtype markers PV, SST, and VIP. One-way ANOVA (tdTomato: I12b vs hSyn, Arl4d vs hSyn, Dlgap1 vs hSyn; ****p < 0.0001); (PV: I12b vs hSyn, **p = 0.0018; Arl4d vs hSyn, ***p = 0.0007; Dlgap1 vs hSyn, **p = 0.0010); (SST: I12b vs hSyn, **p = 0.0025; Arl4d vs hSyn, **p = 0.0074; Dlgap1 vs hSyn, **p = 0.0097); (VIP: I12b vs hSyn, *p = 0.0488). C, % of EYFP+ cells expressing CIN subtype markers normalized to the abundance of each subtype. Values used for normalization are summarized in Extended Data Figure 4-1C. A one-sample t test was used to determine significant difference of normalized values from 1 (SST: AAV-Dlgap1, **p = 0.0063; AAV-hSyn *p = 0.0323; VIP: AAV-I12b, *p = 0.0339; AAV-Arl4d, **p = 0.0028; AAV-Dlgap1, **p = 0.0032; AAV-hSyn, *p = 0.0110). Data are shown as mean ± SEM, as well as individual values for each animal. Non-significant p values not shown. Quantification of EYFP and marker colocalization broken down by rostro-caudal level is shown in Extended Data Figure 5-1D–F.

Extended Data Figure 5-1

IHC characterization of CINs targeted by RE-AAVs by rostrocaudal level of the mPFC of I12b-Cre; Ai14 mice. A, Coronal schematics of the sections A–D used for histology. Counts were taken from the areas highlighted in yellow. Cg, cingulate cortex; PrL, prelimbic cortex; IL, infralimbic cortex; VO, ventral orbital cortex; MO medial orbital cortex; DP, dorsal peduncular cortex; CPu, caudate putamen (striatum); cc, corpus callosum. Adapted from Paxinos and Franklin (2001). B, Quantification of % coverage of the three main CIN subtypes by endogenous tdTomato expression in the mPFC of I12b-Cre; Ai14 mice. Similar coverage of CIN populations by GFP in I12b-Cre; Z/EG somatosensory cortex was reported previously (Potter et al., 2009). C, Quantification of the abundance of PV, SST, and VIP CINs as % of tdTomato+ cells in the mPFC of I12b-Cre; Ai14 mice. D, Quantification of the % of EYFP-positive cells targeted by the RE AAVs that are CINs, i.e., label for endogenous tdTomato. E, Quantification of the % of EYFP-positive cells targeted by the RE AAVs that express the CIN subtype markers PV, SST, and VIP. Counts are shown separated by rostro-caudal levels A–D. F, % of EYFP-positive cells expressing CIN subtype markers normalized to the abundance of each subtype (% marker/EYFP divided by % marker/tdTomato), separated by rostro-caudal levels A–D. Data are shown as mean ± SEM, as well as individual values for each animal. Download Figure 5-1, EPS file.

AAV-I12b-ChR2-EYFP reporter expression was restricted to GABAergic CINs as previously described (Cho et al., 2015); 90% of EYFP+ cells co-expressed tdTomato (Fig. 5B), despite CINs accounting for only ∼20% of the neurons in the cortex. In line with our electrophysiological characterization, AAV-Arl4d-ChR2-EYFP and AAV-Dlgap1-ChR2-EYFP also drove EYFP expression predominantly in CINs; EYFP cells were 78% and 75% tdTomato+, respectively (Fig. 5B). As expected, AAV-hSyn-ChR2-EYFP reporter expression was seen in many more PNs than CINs, and only 6% of EYFP+ cells were positive for tdTomato. The pattern of EYFP labeling produced by AAV-hSyn-ChR2-EYFP was significantly different from the CIN-restricted expression of the other three viruses (I12b vs hSyn, Arl4d vs hSyn, and Dlgap1 vs hSyn: ****p < 0.0001; one-way ANOVA with Tukey’s correction for multiple comparisons).

We then immunostained with the CIN markers PV, SST, and VIP to determine whether our RE AAVs target specific subtypes of CINs. EYFP+ cells labeled by AAV-Arl4d-ChR2-EYFP and AAV-Dlgap1-ChR2-EYFP were 57% and 53% PV+, respectively, while AAV-I12b-ChR2-EYFP-labeled cells were 52% PV+ (Fig. 5B). SST-expressing neurons accounted for 11–16% of the EYFP+ cells and VIP-expressing neurons made up <3% of EYFP+ cells for all three CIN viruses (Fig. 5B). The majority of AAV-Syn-ChR2-EYFP-labeled cells were PNs and thus did not colocalize with PV, SST, or VIP. Counts separated by rostro-caudal position are shown in Extended Data Figure 5-1D–F.

To determine whether any of the RE AAVs preferentially label PV+ CINs, EYFP-CIN marker colocalization was normalized to the percent of marker+ CINs out of the total tdTomato+ CIN count in the same section (Fig. 5C; Extended Data Fig. 5-1C). AAV-I12b-ChR2-EYFP, AAV-Arl4d-ChR2-EYFP and AAV-Dlgap1-ChR2-EYFP-labeled populations were 1.8–twofold enriched for PV+ CINs compared with the overall CIN population, although this enrichment was not significant. SST CINs were significantly underrepresented in the EYFP+ population labeled by AAV-Dlgap1-ChR2-EYFP (p = 0.0063, one-sample t test), and VIP CINs were underrepresented in the populations labeled by all three CIN viruses (I12b, p = 0.0339; Arl4d, p = 0.0028; Dlgap1, p = 0.0032). This suggests some degree of preferential labeling of PV+ CINs by all three viruses.

RE AAV-driven expression in a CIN-dependent behavioral task: rule shift

We next investigated the function of CINs targeted by the Arl4d pRE using a cognitive flexibility assay in which mice learn a “rule shift.” This task requires mice to choose between bowls using odor-based or texture-based cues to find hidden food rewards and then shift from an odor-based rule to a texture-based one, or vice versa. These rule shifts depend on prefrontal γ oscillations, which are believed to be generated by fast-spiking PV+ CINs (Cardin et al., 2009; Sohal et al., 2009; Cho and Sohal, 2014; Cho et al., 2015).

In this study we took advantage of Dlx5/6+/− mice, which have abnormal PV CIN physiology, decreased task-evoked γ oscillations, and fail to suppress perseverative responses during rule shifts (Cho et al., 2015). It was previously shown that optogenetic stimulation of mPFC CINs using AAV-I12b-ChR2-EYFP can completely normalize rule-shift performance in Dlx5/6+/− mice. This effect occurs when stimulation is delivered at γ frequencies (40 or 60 Hz; Cho et al., 2015). Therefore, we tested whether another RE AAV that preferentially targets fast-spiking PV CINs, AAV-Arl4d-ChR2-EYFP, can also produce ChR2 expression that rescues rule-shift performance in Dlx5/6+/− mice (Fig. 6A). On day 1, we assayed rule-shift performance in the absence of optogenetic stimulation to quantify baseline performance of these mutant mice. As expected, the mice were specifically impaired during the rule-shift portion of the task, making a large number of perseverative errors (Fig. 6B,C). The number of trials to criterion and perseverative errors were very similar to previously published results for adult Dlx5/6+/− mice (Cho et al., 2015). Strikingly, delivery of optical stimulation at 40 Hz on day two dramatically improved rule-shift performance compared with day 1 in Dlx5/6+/− mice and significantly reduced the number of perseverative errors (Fig. 6B,C). The improvements were comparable to those seen with stimulation after injection of AAV-I12b-ChR2-EYFP in the Dlx5/6+/− mutant mice, whereas stimulation following injection of mock AAV (AAV-I12b-mCherry) did not rescue this behavioral deficit (Cho et al., 2015).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Optogenetic stimulation of AAV-Arl4d-ChR2-EYFP-transduced neurons improves cognitive flexibility in Dlx5/6−/− mice. A, Dlx5/6+/− mice (n = 3) had bilateral AAV-Arl4d-ChR2-EYFP injections and fiber-optic implants into the mPFC. Low-magnification image of EYFP expression after bilateral injection of AAV-Arl4d-ChR2-EYFP. Experimental design: day 1: no stimulation; day 2: 40-Hz stimulation during the rule-shift portion of the task. B, 40-Hz stimulation on day 2 rescues rule-shift behavior in mutant mice (t(2) = 12.85, **p = 0.006). C, 40-Hz stimulation on day 2 decreased the number of perseverative errors compared with no stimulation on day 1 (t(2) = 5.196, *p = 0.035). There was no change in random errors (t(2)= 1.512, p = 0.2697). Two-tailed, paired t tests were used. Individual values for each animal are shown.

Discussion

Elegant combinatorial transgenic approaches have made it possible to study the functions of defined cell populations in vivo in mice, including for subtypes of CINs (Dymecki et al., 2010; Taniguchi et al., 2011; He et al., 2016). However, these approaches cannot be easily adopted across species. The number of cell populations that can be independently accessed in the same animal is also limited with these approaches.

AAV targeting strategies have been lauded as an alternative to overcome these obstacles and to facilitate translation to human genetic therapies, but the lack of AAVs that can target specific cell types remains a significant impediment to realizing these goals within the heterogeneous cellular environment of the brain (Betley and Sternson, 2011). Here we have reported two novel REs (Arl4d_RE and Dlgap1_RE) that can be used to drive CIN expression of reporter and effector genes using the AAV delivery system. We also identified nearly three thousand CIN-specific pREs genome wide, providing a starting point for the development of further tools for CIN functional studies.

AAV vectors are restricted in their genomic carrying capacity, which limits the lengths of sequences that can be used to direct the expression of reporters and effectors. Several attempts have been made to use short promoter sequences from neuronal marker genes. This strategy has been successful for pan-neuronal or pan-excitatory neuron targeting with a variety of viral vectors, but less successful for labeling subtypes of neurons, including CINs (Kügler et al., 2003a,b; Dittgen et al., 2004; Nathanson et al., 2009a; Sohal et al., 2009). Distal REs (enhancers) can also be short enough to be used in AAV vectors. REs can drive cell-type-specific expression within the developing forebrain (Visel et al., 2013; Pattabiraman et al., 2014; Silberberg et al., 2016), and REs of the Dlx genes have been used in AAV to successfully target CINs in adult mice (Lee et al., 2014a,b; Cho et al., 2015; Dimidschstein et al., 2016; Mehta et al., 2019). Thus, we sought to identify REs that could show specific activity within CIN subtypes.

Epigenetic screen for candidate CIN enhancers

We first generated a genome-wide list of pREs likely to be active in immature CINs using an epigenomic approach. We leveraged genetic GFP labeling of immature CINs to purify these neurons before performing ChIP-seq for active (H3K27ac) and repressive (H3K27me3) histone modifications. By performing histone modification ChIP-seq on the non-CIN population in parallel, we were able to identify pREs likely to be not only active in CINs, but also specific to CINs compared with PNs and other cortical cells such as astrocytes. This served as a starting point for identifying cell type-specific pREs. We then combined this with a directed “candidate gene” based approach looking for pREs near genes expressed in CINs. Finally, we integrated our H3K27ac histone ChIP-seq data with published epigenomic datasets from embryonic (TF ChIP-seq) and adult (ATAC-seq) CINs to increase the likelihood that these candidate pREs would be active across CIN developmental stages.

Other groups have taken a similar approach using mainly chromatin accessibility data from ATAC-seq experiments to identify REs for use in AAV targeting CINs (Hrvatin et al., 2019; Mich et al., 2020; Vormstein-Schneider et al., 2020). Candidate regions with enrichment of H3K27ac are likely to give a better success rate than accessible regions determined by ATAC-seq alone because H3K27ac is a more specific mark of active enhancers. This makes our dataset a valuable resource for others seeking to identify additional CIN-specific REs.

One of the main drawbacks of using a genome-wide screen to identify pREs is that it may identify many false positives, as individual epigenetic markers are not sufficient to prove that a locus will function in the cell type under investigation (Nord et al., 2015; Hrvatin et al., 2019; Mehta et al., 2019; Mich et al., 2020; Vormstein-Schneider et al., 2020). However, we postulated that by integrating multiple sources of data we could increase our confidence in the candidate pREs. Thus, we combined our H3K27Ac screen with previously published DLX2, LHX6, and NKX2-1 TF ChIP-seq (Sandberg et al., 2016; Lindtner et al., 2019) to help identify enhancers active in developing CINs. This, however, still does not ensure activity in adult CINs in the absence of further testing.

In addition, H3K27ac enrichment in the CIN population compared with the non-CIN population does not guarantee complete specificity of RE activity to CINs. The non-CIN population is highly heterogeneous, thus a pRE active in CINs and a minority non-CIN cell type, e.g., a subtype of PN, may not show significant enrichment in the non-CIN bulk sample. Any screen will inevitably produce false negatives as well. Sample quality and sequencing depth may affect how many false negatives are produced.

Another consideration is the fact that genomic context can affect expression; even if endogenously active REs are accurately identified, their activity could be altered in heterologous systems such as extrachromosomal AAV vectors. A recent study compared REs in AAVs and transgenic mouse lines to drive expression in glutamatergic cells of the entorhinal cortex and found that only half of the tested AAVs targeted cells in a similar pattern to the corresponding transgenic line (Nair et al., 2020). Despite improvements in computational RE prediction methods, this is likely to mean that individual candidates will continue to need empirical testing. An elegant, recently published method addresses the testing bottleneck by combining massively parallel reporter assays with single-nucleus RNA-seq to test hundreds of candidate enhancer AAVs in a single animal (Hrvatin et al., 2019). They found that a small fraction of tested sequences drove SST+ CIN-specific expression as predicted from chromatin accessibility signatures.

Finally, because we used the entire neonatal CIN population (purified from Gad67-GFP mice) for H3K27ac ChIP-seq, our screen may not efficiently identify REs whose activity is specific to individual adult CIN subtypes. Our approach can be extended to improve on this by using, for example, sorted CINs from SST-Cre, PV-Cre, and VIP-Cre transgenic animals as has been done for ATAC-Seq (Mo et al., 2015). The major hurdle to this is obtaining sufficient numbers of cells for histone modification ChIP-Seq, which generally requires at least ∼100,000 cells. However, new methods such as CUT&RUN and CUT&TAG are improving the ability to perform epigenetic analyses with cell numbers in the 100–1000 s (Skene and Henikoff, 2017; Skene et al., 2018; Kaya-Okur et al., 2019).

Discovery of two new CIN-specific REs

We have shown that our two selected candidate pREs, Arl4d and Dlgap1, can drive robust expression when placed in AAV vectors. The strength of reporter expression and specificity to the GABAergic population are comparable to AAV using the well-characterized I12b RE.

Our slice physiology experiments demonstrated that Arl4d and I12b pRE-driven EYFP is predominantly expressed in fast-spiking CINs. In addition, we have shown that AAV-Arl4d-ChR2-EYFP can rescue behavioral deficits linked to reduced fast-spiking activity in Dlx5/6+/− mutants through optogenetic stimulation. Thus, these RE AAVs can drive channelrhodopsin expression at sufficient levels for functional/behavioral studies, which can sometimes be a problem with AAVs targeting specific cell types (Sohal et al., 2009). This demonstrates the wide-ranging utility of our RE AAVs for the neuroscience community.

We also showed that AAV-Dlgap1-ChR2-EYFP-targeted cells were more likely to exhibit regular-spiking than fast-spiking properties in slice physiology experiments. This suggests that specificity for different physiologically defined CIN types might be achievable by selecting different REs. In particular, Dlgap1_RE and Arl4d_RE seem to differentially label fast-spiking versus regular-spiking neurons. However, when we assessed this further by immunohistochemistry, we found that all three enhancers preferentially labeled PV-expressing neurons. Thus, additional work will be necessary to parse out the differences in specificity reported using electrophysiological and immunohistological analyses.

One possible explanation for the difference between the electrophysiology results and immunostaining for AAV-Dlgap1-ChR2-EYFP is that some cells could express EYFP more strongly than others. These would likely be the cells chosen for patching. In contrast, immunohistochemistry used antibody labeling for EYFP, which could amplify labeling in weakly expressing cells. Thus, immunohistochemistry may sample a broader population of transduced neurons than electrophysiology experiments that depend on live fluorescence. Along these lines, Mehta et al. (2019) found that the h12R enhancer drove expression at high and low levels with a bimodal distribution when they assessed reporter expression by in situ hybridization. In addition, fast-spiking properties and PV expression, although highly correlated, do not show a 1:1 relationship (Moore and Wehr, 2013). For example, PV-negative fast-spiking Chandelier cells have been reported in the mPFC (Taniguchi et al., 2013), and PV+ CINs with non-fast-spiking properties have also been reported (Blatow et al., 2003). One possibility is that AAV-Dlgap1-ChR2-EYFP may preferentially label regular-spiking PV+ cells. Of course, our electrophysiology dataset has a smaller sample size than our immunohistochemistry dataset. As such, the discordance between these two methods may simply represent statistical anomalies in our electrophysiological dataset.

Future utility of CIN-specific AAVs

AAVs with CIN-specific conserved REs offer a significant advantage over the existing mouse transgenic lines because of the ease of use of these tools across species, including in humans and non-human primates (Dimidschstein et al., 2016; Mehta et al., 2019; Mich et al., 2020; Vormstein-Schneider et al., 2020). In addition, RE AAVs injected in small volumes would facilitate effector protein expression and/or gene knock-out in restricted brain regions, and therefore could be used to uncover the physiology and function of CINs in different cortical regions.

REs could potentially be further engineered by mutation or addition of TF binding sites to increase activity or specificity. The REs described in the present study could also be tested in the context of different AAV serotypes, as the specificity of AAV-driven expression can be affected by the serotype and titer of the virus used (Nathanson et al., 2009b). Overall, we have identified novel CIN specific REs and our study describes methodology for identifying potentially hundreds of additional enhancers that could drive CIN subtype-specific expression in AAV vectors.

Acknowledgments

Acknowledgements: Next-generation sequencing was carried out by the Center for Advanced Technology at University of California San Francisco.

Footnotes

  • J.L.R.R. is a co-founder, stockholder, and currently on the scientific board of Neurona, a company studying the potential therapeutic use of interneuron transplantation. All other authors declare no competing financial interests.

  • This work was supported by National Institute of Mental Health Grants U01MH105948 (to J.L.R.R. and V.S.S.) and R01 MH081880 (to J.L.R.R.) and by a National Alliance for Research on Schizophrenia and Depression Young Investigator Grant, Leichtung family Investigator, Brain and Behavior Research Foundation (R.M.). Confocal microscopy was performed at the University of California San Francisco (UCSF) Nikon Imaging Center, supported by the National Institutes of Health (NIH) S10 Shared Instrumentation Grant 1S10OD017993-01A1. Cell sorting was performed at the UCSF Helen Diller Family Comprehensive Cancer Center Laboratory for Cell Analysis, using instrumentation supported by NIH Grant P30CA082103.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Arguello A, Yang X, Vogt D, Stanco A, Rubenstein JLR, Cheyette BNR (2013) Dapper antagonist of catenin-1 cooperates with Dishevelled-1 during postsynaptic development in mouse forebrain GABAergic interneurons. PLoS One 8:e67679. doi:10.1371/journal.pone.0067679 pmid:23826333
    OpenUrlCrossRefPubMed
  2. ↵
    Betley JN, Sternson SM (2011) Adeno-associated viral vectors for mapping, monitoring, and manipulating neural circuits. Hum Gene Ther 22:669–677. doi:10.1089/hum.2010.204 pmid:21319997
    OpenUrlCrossRefPubMed
  3. ↵
    Blatow M, Rozov A, Katona I, Hormuzdi SG, Meyer AH, Whittington MA, Caputi A, Monyer H (2003) A novel network of multipolar bursting interneurons generates theta frequency oscillations in neocortex. Neuron 38:805–817. doi:10.1016/s0896-6273(03)00300-3 pmid:12797964
    OpenUrlCrossRefPubMed
  4. ↵
    Burger C, Gorbatyuk OS, Velardo MJ, Peden CS, Williams P, Zolotukhin S, Reier PJ, Mandel RJ, Muzyczka N (2004) Recombinant AAV viral vectors pseudotyped with viral capsids from serotypes 1, 2, and 5 display differential efficiency and cell tropism after delivery to different regions of the central nervous system. Mol Ther 10:302–317. doi:10.1016/j.ymthe.2004.05.024 pmid:15294177
    OpenUrlCrossRefPubMed
  5. ↵
    Cardin JA, Carlén M, Meletis K, Knoblich U, Zhang F, Deisseroth K, Tsai L-H, Moore CI (2009) Driving fast-spiking cells induces gamma rhythm and controls sensory responses. Nature 459:663–667. doi:10.1038/nature08002 pmid:19396156
    OpenUrlCrossRefPubMed
  6. ↵
    Chen G, Zhang Y, Li X, Zhao X, Ye Q, Lin Y, Tao HW, Rasch MJ, Zhang X (2017) Distinct inhibitory circuits orchestrate cortical beta and gamma band oscillations. Neuron 96:1403–1418.e6. doi:10.1016/j.neuron.2017.11.033
    OpenUrlCrossRefPubMed
  7. ↵
    Cho KKA, Sohal VS (2014) Optogenetic approaches for investigating neural pathways implicated in schizophrenia and related disorders. Hum Mol Genet 23:R64–R68. doi:10.1093/hmg/ddu225 pmid:24824218
    OpenUrlCrossRefPubMed
  8. ↵
    Cho KKA, Hoch R, Lee AT, Patel T, Rubenstein JLR, Sohal VS (2015) Gamma rhythms link prefrontal interneuron dysfunction with cognitive inflexibility in Dlx5/6(+/-) mice. Neuron 85:1332–1343. doi:10.1016/j.neuron.2015.02.019 pmid:25754826
    OpenUrlCrossRefPubMed
  9. ↵
    Colasante G, Collombat P, Raimondi V, Bonanomi D, Ferrai C, Maira M, Yoshikawa K, Mansouri A, Valtorta F, Rubenstein JLR, Broccoli V (2008) Arx is a direct target of Dlx2 and thereby contributes to the tangential migration of GABAergic interneurons. J Neurosci 28:10674–10686. doi:10.1523/JNEUROSCI.1283-08.2008 pmid:18923043
    OpenUrlAbstract/FREE Full Text
  10. ↵
    Creyghton MP, Cheng AW, Welstead GG, Kooistra T, Carey BW, Steine EJ, Hanna J, Lodato MA, Frampton GM, Sharp PA, Boyer LA, Young RA, Jaenisch R (2010) Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc Natl Acad Sci USA 107:21931–21936. doi:10.1073/pnas.1016071107 pmid:21106759
    OpenUrlAbstract/FREE Full Text
  11. ↵
    Dimidschstein J, Chen Q, Tremblay R, Rogers SL, Saldi GA, Guo L, Xu Q, Liu R, Lu C, Chu J, Grimley JS, Krostag AR, Kaykas A, Avery MC, Rashid MS, Baek M, Jacob AL, Smith GB, Wilson DE, Kosche G, et al. (2016) A viral strategy for targeting and manipulating interneurons across vertebrate species. Nat Neurosci 19:1743–1749. doi:10.1038/nn.4430 pmid:27798629
    OpenUrlCrossRefPubMed
  12. ↵
    Dittgen T, Nimmerjahn A, Komai S, Licznerski P, Waters J, Margrie TW, Helmchen F, Denk W, Brecht M, Osten P (2004) Lentivirus-based genetic manipulations of cortical neurons and their optical and electrophysiological monitoring in vivo. Proc Natl Acad Sci USA 101:18206–18211. doi:10.1073/pnas.0407976101 pmid:15608064
    OpenUrlAbstract/FREE Full Text
  13. ↵
    Dymecki SM, Ray RS, Kim JC (2010) Mapping cell fate and function using recombinase-based intersectional strategies. Meth Enzymol 477:183–213.
    OpenUrlCrossRefPubMed
  14. ↵
    Fishell G, Kepecs A (2020) Interneuron types as attractors and controllers. Annu Rev Neurosci 43:1–30. doi:10.1146/annurev-neuro-070918-050421 pmid:31299170
    OpenUrlCrossRefPubMed
  15. ↵
    Friocourt G, Parnavelas JG (2011) Identification of Arx targets unveils new candidates for controlling cortical interneuron migration and differentiation. Front Cell Neurosci 5:28. doi:10.3389/fncel.2011.00028 pmid:22355284
    OpenUrlCrossRefPubMed
  16. ↵
    Gibson DG, Young L, Chuang R-Y, Venter JC, Hutchison CA, Smith HO (2009) Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods 6:343–345. doi:10.1038/nmeth.1318 pmid:19363495
    OpenUrlCrossRefPubMed
  17. ↵
    Haggerty DL, Grecco GG, Reeves KC, Atwood B (2020) Adeno-associated viral vectors in neuroscience research. Mol Ther Methods Clin Dev 17:69–82. doi:10.1016/j.omtm.2019.11.012 pmid:31890742
    OpenUrlCrossRefPubMed
  18. ↵
    Hashemi E, Ariza J, Rogers H, Noctor SC, Martínez-Cerdeño V (2017) The number of parvalbumin-expressing interneurons is decreased in the prefrontal cortex in autism. Cereb Cortex 27:1931–1943. doi:10.1093/cercor/bhw021 pmid:26922658
    OpenUrlCrossRefPubMed
  19. ↵
    He M, Tucciarone J, Lee S, Nigro MJ, Kim Y, Levine JM, Kelly SM, Krugikov I, Wu P, Chen Y, Gong L, Hou Y, Osten P, Rudy B, Huang ZJ (2016) Strategies and tools for combinatorial targeting of GABAergic neurons in mouse cerebral cortex. Neuron 91:1228–1243. doi:10.1016/j.neuron.2016.08.021 pmid:27618674
    OpenUrlCrossRefPubMed
  20. ↵
    Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK (2010) Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell 38:576–589. doi:10.1016/j.molcel.2010.05.004 pmid:20513432
    OpenUrlCrossRefPubMed
  21. ↵
    Hrvatin S, Tzeng CP, Nagy MA, Stroud H, Koutsioumpa C, Wilcox OF, Assad EG, Green J, Harvey CD, Griffith EC, Greenberg ME (2019) A scalable platform for the development of cell-type-specific viral drivers. Elife 8:e48089. doi:10.7554/eLife.48089
    OpenUrlCrossRefPubMed
  22. ↵
    Hu JS, Vogt D, Sandberg M, Rubenstein JL (2017) Cortical interneuron development: a tale of time and space. Development 144:3867–3878. doi:10.1242/dev.132852 pmid:29089360
    OpenUrlAbstract/FREE Full Text
  23. ↵
    Jiang X, Lachance M, Rossignol E (2016) Involvement of cortical fast-spiking parvalbumin-positive basket cells in epilepsy. Prog Brain Res 226:81–126. doi:10.1016/bs.pbr.2016.04.012 pmid:27323940
    OpenUrlCrossRefPubMed
  24. ↵
    Kawaguchi Y, Katsumaru H, Kosaka T, Heizmann CW, Hama K (1987) Fast spiking cells in rat hippocampus (CA1 region) contain the calcium-binding protein parvalbumin. Brain Res 416:369–374. doi:10.1016/0006-8993(87)90921-8 pmid:3304536
    OpenUrlCrossRefPubMed
  25. ↵
    Kaya-Okur HS, Wu SJ, Codomo CA, Pledger ES, Bryson TD, Henikoff JG, Ahmad K, Henikoff S (2019) CUT&Tag for efficient epigenomic profiling of small samples and single cells. Nat Commun 10:1930–1910. doi:10.1038/s41467-019-09982-5 pmid:31036827
    OpenUrlCrossRefPubMed
  26. ↵
    Kepecs A, Fishell G (2014) Interneuron cell types are fit to function. Nature 505:318–326. doi:10.1038/nature12983 pmid:24429630
    OpenUrlCrossRefPubMed
  27. ↵
    Kessaris N, Magno L, Rubin AN, Oliveira MG (2014) Genetic programs controlling cortical interneuron fate. Curr Opin Neurobiol 26:79–87. doi:10.1016/j.conb.2013.12.012 pmid:24440413
    OpenUrlCrossRefPubMed
  28. ↵
    Kügler S, Kilic E, Bähr M (2003a) Human synapsin 1 gene promoter confers highly neuron-specific long-term transgene expression from an adenoviral vector in the adult rat brain depending on the transduced area. Gene Ther 10:337–347. doi:10.1038/sj.gt.3301905 pmid:12595892
    OpenUrlCrossRefPubMed
  29. ↵
    Kügler S, Lingor P, Schöll U, Zolotukhin S, Bähr M (2003b) Differential transgene expression in brain cells in vivo and in vitro from AAV-2 vectors with small transcriptional control units. Virology 311:89–95. doi:10.1016/S0042-6822(03)00162-4 pmid:12832206
    OpenUrlCrossRefPubMed
  30. ↵
    Lee AT, Gee SM, Vogt D, Patel T, Rubenstein JL, Sohal VS (2014a) Pyramidal neurons in prefrontal cortex receive subtype-specific forms of excitation and inhibition. Neuron 81:61–68. doi:10.1016/j.neuron.2013.10.031 pmid:24361076
    OpenUrlCrossRefPubMed
  31. ↵
    Lee AT, Vogt D, Rubenstein JL, Sohal VS (2014b) A class of GABAergic neurons in the prefrontal cortex sends long-range projections to the nucleus accumbens and elicits acute avoidance behavior. J Neurosci 34:11519–11525. doi:10.1523/JNEUROSCI.1157-14.2014 pmid:25164650
    OpenUrlAbstract/FREE Full Text
  32. ↵
    Lein ES, Hawrylycz MJ, Ao N, Ayres M, Bensinger A, Bernard A, Boe AF, Boguski MS, Brockway KS, Byrnes EJ, Chen L, Chen L, Chen TM, Chin MC, Chong J, Crook BE, Czaplinska A, Dang CN, Datta S, Dee NR, et al. (2007) Genome-wide atlas of gene expression in the adult mouse brain. Nature 445:168–176. doi:10.1038/nature05453 pmid:17151600
    OpenUrlCrossRefPubMed
  33. ↵
    Lewis DA, Curley AA, Glausier JR, Volk DW (2012) Cortical parvalbumin interneurons and cognitive dysfunction in schizophrenia. Trends Neurosci 35:57–67. doi:10.1016/j.tins.2011.10.004 pmid:22154068
    OpenUrlCrossRefPubMed
  34. ↵
    Li H (2013) Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997v2.
  35. ↵
    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R; 1000 Genome Project Data Processing Subgroup (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079. doi:10.1093/bioinformatics/btp352 pmid:19505943
    OpenUrlCrossRefPubMed
  36. ↵
    Li Q, Brown JB, Huang H, Bickel PJ (2011) Measuring reproducibility of high-throughput experiments. Ann Appl Stat 5:1752–1779. doi:10.1214/11-AOAS466
    OpenUrlCrossRefPubMed
  37. ↵
    Lindtner S, Catta-Preta R, Tian H, Su-Feher L, Price JD, Dickel DE, Greiner V, Silberberg SN, Mckinsey GL, McManus MT, Pennacchio LA, Visel A, Nord AS, Rubenstein JLR (2019) Genomic resolution of DLX-orchestrated transcriptional circuits driving development of forebrain GABAergic neurons. Cell Rep 28:2048–2063.e8. doi:10.1016/j.celrep.2019.07.022 pmid:31433982
    OpenUrlCrossRefPubMed
  38. ↵
    Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25:402–408. doi:10.1006/meth.2001.1262 pmid:11846609
    OpenUrlCrossRefPubMed
  39. ↵
    Long JE, Cobos I, Potter GB, Rubenstein JLR (2009) Dlx1&2 and Mash1 transcription factors control MGE and CGE patterning and differentiation through parallel and overlapping pathways. Cereb Cortex 19 [Suppl 1]:i96–i106. doi:10.1093/cercor/bhp045 pmid:19386638
    OpenUrlCrossRefPubMed
  40. ↵
    Madisen L, Zwingman TA, Sunkin SM, Oh SW, Zariwala HA, Gu H, Ng LL, Palmiter RD, Hawrylycz MJ, Jones AR, Lein ES, Zeng H (2010) A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. Nat Neurosci 13:133–140. doi:10.1038/nn.2467 pmid:20023653
    OpenUrlCrossRefPubMed
  41. ↵
    Malik R, Pai EL-L, Rubin AN, Stafford AM, Angara K, Minasi P, Rubenstein JL, Sohal VS, Vogt D (2019) Tsc1 represses parvalbumin expression and fast-spiking properties in somatostatin lineage cortical interneurons. Nat Commun 10:4994–4916. doi:10.1038/s41467-019-12962-4 pmid:31676823
    OpenUrlCrossRefPubMed
  42. ↵
    McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G (2010) GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol 28:495–501. doi:10.1038/nbt.1630 pmid:20436461
    OpenUrlCrossRefPubMed
  43. ↵
    Mehta P, Kreeger L, Wylie DC, Pattadkal JJ, Lusignan T, Davis MJ, Turi GF, Li W-K, Whitmire MP, Chen Y, Kajs BL, Seidemann E, Priebe NJ, Losonczy A, Zemelman BV (2019) Functional access to neuron subclasses in rodent and primate forebrain. Cell Rep 26:2818–2832.e8. doi:10.1016/j.celrep.2019.02.011 pmid:30840900
    OpenUrlCrossRefPubMed
  44. ↵
    Mich JK, Graybuck LT, Hess EE, Mahoney JT, Kojima Y, Ding Y, Somasundaram S, Miller JA, Weed N, Omstead V, Bishaw Y, Shapovalova NV, Martinez RA, Fong O, Yao S, Mortrud M, Chong P, Loftus L, Bertagnolli D, Goldy J, et al. (2020) Functional enhancer elements drive subclass-selective expression from mouse to primate neocortex. bioRxiv 29:555318.
    OpenUrl
  45. ↵
    Mo A, Mukamel EA, Davis FP, Luo C, Henry GL, Picard S, Urich MA, Nery JR, Sejnowski TJ, Lister R, Eddy SR, Ecker JR, Nathans J (2015) Epigenomic signatures of neuronal diversity in the mammalian brain. Neuron 86:1369–1384. doi:10.1016/j.neuron.2015.05.018 pmid:26087164
    OpenUrlCrossRefPubMed
  46. ↵
    Moore AK, Wehr M (2013) Parvalbumin-expressing inhibitory interneurons in auditory cortex are well-tuned for frequency. J Neurosci 33:13713–13723. doi:10.1523/JNEUROSCI.0663-13.2013 pmid:23966693
    OpenUrlAbstract/FREE Full Text
  47. ↵
    Nair RR, Blankvoort S, Lagartos MJ, Kentros C (2020) Enhancer-driven gene expression (EDGE) enables the generation of viral vectors specific to neuronal subtypes. iScience 23:100888. doi:10.1016/j.isci.2020.100888 pmid:32087575
    OpenUrlCrossRefPubMed
  48. ↵
    Nathanson JL, Jappelli R, Scheeff ED, Manning G, Obata K, Brenner S, Callaway EM (2009a) Short promoters in viral vectors drive selective expression in mammalian inhibitory neurons, but do not restrict activity to specific inhibitory cell-types. Front Neural Circuits 3:19. doi:10.3389/neuro.04.019.2009 pmid:19949461
    OpenUrlCrossRefPubMed
  49. ↵
    Nathanson JL, Yanagawa Y, Obata K, Callaway EM (2009b) Preferential labeling of inhibitory and excitatory cortical neurons by endogenous tropism of adeno-associated virus and lentivirus vectors. Neuroscience 161:441–450. doi:10.1016/j.neuroscience.2009.03.032 pmid:19318117
    OpenUrlCrossRefPubMed
  50. ↵
    Nigro MJ, Hashikawa-Yamasaki Y, Rudy B (2018) Diversity and Connectivity of Layer 5 Somatostatin-Expressing Interneurons in the Mouse Barrel Cortex. J Neurosci 38:1622–1633.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    Nord AS, Pattabiraman K, Visel A, Rubenstein JLR (2015) Genomic perspectives of transcriptional regulation in forebrain development. Neuron 85:27–47. doi:10.1016/j.neuron.2014.11.011 pmid:25569346
    OpenUrlCrossRefPubMed
  52. ↵
    Pattabiraman K, Golonzhka O, Lindtner S, Nord AS, Taher L, Hoch R, Silberberg SN, Zhang D, Chen B, Zeng H, Pennacchio LA, Puelles L, Visel A, Rubenstein JLR (2014) Transcriptional regulation of enhancers active in protodomains of the developing cerebral cortex. Neuron 82:989–1003. doi:10.1016/j.neuron.2014.04.014 pmid:24814534
    OpenUrlCrossRefPubMed
  53. ↵
    Pla R, Stanco A, Howard MA, Rubin AN, Vogt D, Mortimer N, Cobos I, Potter GB, Lindtner S, Price JD, Nord AS, Visel A, Schreiner CE, Baraban SC, Rowitch DH, Rubenstein JLR (2017) Dlx1 and Dlx2 promote interneuron GABA synthesis, synaptogenesis, and dendritogenesis. Cereb Cortex 29:1–19.
    OpenUrl
  54. ↵
    Potter GB, Petryniak MA, Shevchenko E, Mckinsey GL, Ekker M, Rubenstein JLR (2009) Generation of Cre-transgenic mice using Dlx1/Dlx2 enhancers and their characterization in GABAergic interneurons. Mol Cell Neurosci 40:167–186. doi:10.1016/j.mcn.2008.10.003 pmid:19026749
    OpenUrlCrossRefPubMed
  55. ↵
    Preibisch S, Saalfeld S, Tomancak P (2009) Globally optimal stitching of tiled 3D microscopic image acquisitions. Bioinformatics 25:1463–1465. doi:10.1093/bioinformatics/btp184 pmid:19346324
    OpenUrlCrossRefPubMed
  56. ↵
    Quinlan AR, Hall IM (2010) BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26:841–842. doi:10.1093/bioinformatics/btq033 pmid:20110278
    OpenUrlCrossRefPubMed
  57. ↵
    Rada-Iglesias A, Bajpai R, Swigut T, Brugmann SA, Flynn RA, Wysocka J (2011) A unique chromatin signature uncovers early developmental enhancers in humans. Nature 470:279–283. doi:10.1038/nature09692 pmid:21160473
    OpenUrlCrossRefPubMed
  58. ↵
    Ramírez F, Ryan DP, Grüning B, Bhardwaj V, Kilpert F, Richter AS, Heyne S, Dündar F, Manke T (2016) deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res 44:W160–W165. doi:10.1093/nar/gkw257 pmid:27079975
    OpenUrlCrossRefPubMed
  59. ↵
    Rasmussen AH, Rasmussen HB, Silahtaroglu A (2017) The DLGAP family: neuronal expression, function and role in brain disorders. Mol Brain 10:43. doi:10.1186/s13041-017-0324-9 pmid:28870203
    OpenUrlCrossRefPubMed
  60. ↵
    Robledo RF, Rajan L, Li X, Lufkin T (2002) The Dlx5 and Dlx6 homeobox genes are essential for craniofacial, axial, and appendicular skeletal development. Genes & Development 16:1089–1101.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    Rudy B, Fishell G, Lee S, Hjerling-Leffler J (2011) Three groups of interneurons account for nearly 100% of neocortical GABAergic neurons. Dev Neurobiol 71:45–61. doi:10.1002/dneu.20853 pmid:21154909
    OpenUrlCrossRefPubMed
  62. ↵
    Sandberg M, Flandin P, Silberberg S, Su-Feher L, Price JD, Hu JS, Kim C, Visel A, Nord AS, Rubenstein JLR (2016) Transcriptional networks controlled by NKX2-1 in the development of forebrain GABAergic neurons. Neuron 91:1260–1275. doi:10.1016/j.neuron.2016.08.020 pmid:27657450
    OpenUrlCrossRefPubMed
  63. ↵
    Silberberg SN, Taher L, Lindtner S, Sandberg M, Nord AS, Vogt D, Mckinsey GL, Hoch R, Pattabiraman K, Zhang D, Ferran JL, Rajkovic A, Golonzhka O, Kim C, Zeng H, Puelles L, Visel A, Rubenstein JLR (2016) Subpallial enhancer transgenic lines: a data and tool resource to study transcriptional regulation of GABAergic cell fate. Neuron 92:59–74. doi:10.1016/j.neuron.2016.09.027 pmid:27710791
    OpenUrlCrossRefPubMed
  64. ↵
    Skene PJ, Henikoff S (2017) An efficient targeted nuclease strategy for high-resolution mapping of DNA binding sites. Elife 6:e21856. doi:10.7554/eLife.21856
    OpenUrlCrossRefPubMed
  65. ↵
    Skene PJ, Henikoff JG, Henikoff S (2018) Targeted in situ genome-wide profiling with high efficiency for low cell numbers. Nat Protoc 13:1006–1019. doi:10.1038/nprot.2018.015 pmid:29651053
    OpenUrlCrossRefPubMed
  66. ↵
    Sohal VS, Zhang F, Yizhar O, Deisseroth K (2009) Parvalbumin neurons and gamma rhythms enhance cortical circuit performance. Nature 459:698–702. doi:10.1038/nature07991 pmid:19396159
    OpenUrlCrossRefPubMed
  67. ↵
    Tamamaki N, Yanagawa Y, Tomioka R, Miyazaki J-I, Obata K, Kaneko T (2003) Green fluorescent protein expression and colocalization with calretinin, parvalbumin, and somatostatin in the GAD67-GFP knock-in mouse. J Comp Neurol 467:60–79.
    OpenUrlCrossRefPubMed
  68. ↵
    Taniguchi H, He M, Wu P, Kim S, Paik R, Sugino K, Kvitsiani D, Kvitsani D, Fu Y, Lu J, Lin Y, Miyoshi G, Shima Y, Fishell G, Nelson SB, Huang ZJ (2011) A resource of Cre driver lines for genetic targeting of GABAergic neurons in cerebral cortex. Neuron 71:995–1013. doi:10.1016/j.neuron.2011.07.026 pmid:21943598
    OpenUrlCrossRefPubMed
  69. ↵
    Taniguchi H, Lu J, Huang ZJ (2013) The spatial and temporal origin of chandelier cells in mouse neocortex. Science 339:70–74. doi:10.1126/science.1227622 pmid:23180771
    OpenUrlAbstract/FREE Full Text
  70. ↵
    Taymans JM, Vandenberghe LH, Haute CVD, Thiry I, Deroose CM, Mortelmans L, Wilson JM, Debyser Z, Baekelandt V (2007) Comparative analysis of adeno-associated viral vector serotypes 1, 2, 5, 7, and 8 in mouse brain. Hum Gene Ther 18:195–206. doi:10.1089/hum.2006.178 pmid:17343566
    OpenUrlCrossRefPubMed
  71. ↵
    Tremblay R, Lee S, Rudy B (2016) GABAergic interneurons in the neocortex: from cellular properties to circuits. Neuron 91:260–292. doi:10.1016/j.neuron.2016.06.033 pmid:27477017
    OpenUrlCrossRefPubMed
  72. ↵
    Visel A, Taher L, Girgis H, May D, Golonzhka O, Hoch RV, McKinsey GL, Pattabiraman K, Silberberg SN, Blow MJ, Hansen DV, Nord AS, Akiyama JA, Holt A, Hosseini R, Phouanenavong S, Plajzer-Frick I, Shoukry M, Afzal V, Kaplan T, et al. (2013) A high-resolution enhancer atlas of the developing telencephalon. Cell 152:895–908. doi:10.1016/j.cell.2012.12.041 pmid:23375746]
    OpenUrlCrossRefPubMed
  73. ↵
    Vogt D, Hunt RF, Mandal S, Sandberg M, Silberberg SN, Nagasawa T, Yang Z, Baraban SC, Rubenstein JLR (2014) Lhx6 directly regulates Arx and CXCR7 to determine cortical interneuron fate and laminar position. Neuron 82:350–364. doi:10.1016/j.neuron.2014.02.030 pmid:24742460
    OpenUrlCrossRefPubMed
  74. ↵
    Vogt D, Cho KKA, Lee AT, Sohal VS, Rubenstein JLR (2015) The parvalbumin/somatostatin ratio is increased in Pten mutant mice and by human PTEN ASD alleles. Cell Rep 11:944–956. doi:10.1016/j.celrep.2015.04.019 pmid:25937288
    OpenUrlCrossRefPubMed
  75. ↵
    Vogt D, Cho KKA, Shelton SM, Paul A, Huang ZJ, Sohal VS, Rubenstein JLR (2017) Mouse Cntnap2 and human CNTNAP2 ASD alleles cell autonomously regulate PV+ cortical interneurons. Cereb Cortex 51:1–12.
    OpenUrl
  76. ↵
    Vormstein-Schneider D, Lin JD, Pelkey KA, Chittajallu R, Guo B, Arias-Garcia MA, Allaway K, Sakopoulos S, Schneider G, Stevenson O, Vergara J, Sharma J, Zhang Q, Franken TP, Smith J, Ibrahim LA, M Astro KJ, Sabri E, Huang S, Favuzzi E, et al. (2020) Viral manipulation of functionally distinct interneurons in mice, non-human primates and humans. Nat Neurosci 50:825–828.
    OpenUrl
  77. ↵
    Watakabe A, Ohtsuka M, Kinoshita M, Takaji M, Isa K, Mizukami H, Ozawa K, Isa T, Yamamori T (2015) Comparative analyses of adeno-associated viral vector serotypes 1, 2, 5, 8 and 9 in marmoset, mouse and macaque cerebral cortex. Neurosci Res 93:144–157. doi:10.1016/j.neures.2014.09.002 pmid:25240284
    OpenUrlCrossRefPubMed
  78. ↵
    Yavorska I, Wehr M (2016) Somatostatin-expressing inhibitory interneurons in cortical circuits. Front Neural Circuits 10:76. doi:10.3389/fncir.2016.00076 pmid:27746722
    OpenUrlCrossRefPubMed
  79. ↵
    Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS (2008) Model-based analysis of ChIP-Seq (MACS). Genome Biol 9:R137–R139. doi:10.1186/gb-2008-9-9-r137 pmid:18798982
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Karl Herrup, University of Pittsburgh School of Medicine

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. The following reviewer(s) agreed to reveal their identity: Danny Leung.

The manuscript addresses the question of the role of cortical interneuron (CIN) dysfunction in neuropsychiatric conditions such as epilepsy, schizophrenia and autism. In approaching this topic the authors set as their goal the development of a new and more powerful tool that uses viral vectors to discover gene regulatory enhancer elements (REs) that are operative in CIN and other specific cell types. The work uncovered two potential REs, which are compared to the better established Dlx I12b enhancer. Both REs drove CIN-specific reporter expression in adult mice and proved useful in specific experimental settings.

The reviewers of the work both felt that this was a well-constructed study that applies basic epigenomic analysis to improve the technology used to drive CIN specific gene expression. They found the emphasis on the genome-wide screen for putative REs to be interesting and potentially valuable. Minor concerns were expressed that the REs defined here were not all that different in function from published constructs. Technical issues were also raised, in particular the adequacy of the descriptions of techniques offered in the Methods as well as the specific approaches to the bioinformatics. On balance, however, there was considerable enthusiasm for working with the authors to improve the work and see it published in eNeuro.

With regard to the bioinformatics, more details must be provided in the Methods to inform the reader as to how the material was generated. For example, in Figure 2, it is implicit but not spelled out that the authors have re-processed the TF and ATAC data with methods similar to the ChIP data that were generated. Further, it is not evident that there is ATAC-seq signal at the highlighted region as mentioned. Beyond the need for more methodological details, concerns were raised about the peak calling. Two replicates were performed, which would usually mean an irreproducible discovery rate (favored by groups such as the ENCODE consortia) would be determined. With more replicates, it would have been possible to use a package such as diffBind and statistical testing provided from DESeq2/edgeR to establish peaks of low variability, Instead, the authors used MACS2 for differential peak calling. This is acceptable given the q<0.01 used for final peaks. However, there needs to be considerably more information provided in the Methods to allow the reader to assess the results. For example, it is unclear if they have called bdgdiff for each replicate using the 4 bedgraphs: t1 (CIN) t2 (nonCIN) c1 input for CIN and c2 input for nonCIN, essential in making the decision to call peaks as high confidence.

Also, it would be better if there were validation experiments performed for the ChIP-seq datasets. An independent ChIP-qPCR at some of the loci shown in the screenshots in Figure 1 would serve this purpose. Moreover, it would allow for quantitative comparison of enrichment like those mentioned in Figure 2B. The authors used GREAT to assign the associations between an RE and a gene. Although it is a standard for GO analysis of non-coding sequences, GREAT largely only considers linear distance between 2 loci. It would help if, in addition, the authors were to compare definitions with existing lists of enhancer-promoter pairs. For instance, the ENCODE project has an extensive list from various different tissues and cell types. Along the same lines, when asserting that an RE “could regulate Dhx8 in addition to or instead of, Arl4d", incorporation of existing definitions and referring to published higher-order chromatin interactions such as those defined by Hi-C data would help.

The authors aimed to define CIN-specific REs by doing ChIP-seq in sorted cells. This is a sound approach but the non-CIN cells are likely far more heterogenous, making it difficult to ensure that the signal is not simply diluted and that the RE is actually not CIN-specific. Flow cytometry analysis to determine the percentage of co-occurrence of EYFP+ and other cell type markers could address this point and help quantify the specificity of the RE in their in vivo assay. While they would strengthen the manuscript, if these experiments are not possible, it will be essential to point out this limitation in the discussion.

For the RE list that incorporates K27ac, TF ChIP-seq and ATAC-seq datasets, a motif enrichment analysis would be informative. It would also be useful to provide the genomic annotations for CIN-peaks compared to expected distributions (fold-elevated). The authors could analyze whether the identified peaks were more likely to be in intergenic regions or proximal genes in dense regions. Consideration should also be given to comparing the REs they have defined to those found in other mammals (primates/humans/etc) since their motivation is cross-species tools to target CINs specifically. They might also wish to consider running a de novo motif analysis (MEME/DREAM suite or something similar) on the CIN specific, non-CIN, and joint peaks to strengthen the link to specific TFs.

Minor Comments:

1. On line 11, histone ChIP-seq should be histone modification ChIP-seq

2. Line 135 mentions reverse crosslinked chromatin but the protocol was a native ChIP.

3. On Line 446, H3K27ac+ should be H3K27 enriched, as the signal is not binary.

4. Clarify for the reader the meaning of the data showing the enrichment of H3K27ac over gene body and proximal sequence. This is a less standard approach to the questions posed and has no widely accepted meaning in the field; it would help if this line of reasoning was expanded.

5. Line 471 should probably read RE and not GE.

6. Genome browser screenshots lack y-axis values.

7. It would be helpful to include a track on the genome browser screenshots that show the peak definitions.

8. Figure 1A. It is unclear in the cartoon as to what “conservation” applies to. Evolutionary? Sequences only for CINs?

9. Figure 4A. The subtype markers used were not explicitly specified. Please clarify this point. It would also be helpful to see the enrichment of K27ac at the promoters of marker genes.

10. Figure 5. An additional negative (mock AAV) or positive (AAV-I12b-ChR2-EYFP) control would go a long way towards supporting this experiment. The authors should address the question of whether the elements they identify in this paradigm perform similarly to Dlx-I12b (given levels of staining/infection)

Back to top

In this issue

eneuro: 7 (6)
eNeuro
Vol. 7, Issue 6
November/December 2020
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Regulatory Elements Inserted into AAVs Confer Preferential Activity in Cortical Interneurons
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Regulatory Elements Inserted into AAVs Confer Preferential Activity in Cortical Interneurons
Anna N. Rubin, Ruchi Malik, Kathleen K. A. Cho, Kenneth J. Lim, Susan Lindtner, Sarah E. Robinson Schwartz, Daniel Vogt, Vikaas S. Sohal, John L. R. Rubenstein
eNeuro 16 November 2020, 7 (6) ENEURO.0211-20.2020; DOI: 10.1523/ENEURO.0211-20.2020

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Regulatory Elements Inserted into AAVs Confer Preferential Activity in Cortical Interneurons
Anna N. Rubin, Ruchi Malik, Kathleen K. A. Cho, Kenneth J. Lim, Susan Lindtner, Sarah E. Robinson Schwartz, Daniel Vogt, Vikaas S. Sohal, John L. R. Rubenstein
eNeuro 16 November 2020, 7 (6) ENEURO.0211-20.2020; DOI: 10.1523/ENEURO.0211-20.2020
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • AAV
  • cortical interneurons
  • enhancers
  • fast spiking
  • regular spiking

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Article: Methods/New Tools

  • Automated Segmentation of the Mouse Body Language to Study Stimulus-Evoked Emotional Behaviors
  • The Role of Genetically Distinct Central Amygdala Neurons in Appetitive and Aversive Responding Assayed with a Novel Dual Valence Operant Conditioning Paradigm
  • In Utero Electroporated Neurons for Medium-Throughput Screening of Compounds Regulating Neuron Morphology
Show more Research Article: Methods/New Tools

Novel Tools and Methods

  • Automated Segmentation of the Mouse Body Language to Study Stimulus-Evoked Emotional Behaviors
  • The Role of Genetically Distinct Central Amygdala Neurons in Appetitive and Aversive Responding Assayed with a Novel Dual Valence Operant Conditioning Paradigm
  • In Utero Electroporated Neurons for Medium-Throughput Screening of Compounds Regulating Neuron Morphology
Show more Novel Tools and Methods

Subjects

  • Novel Tools and Methods

  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.