Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Research Articles, Cellular/Molecular

CaMKIIα-Positive Interneurons Identified via a microRNA-Based Viral Gene Targeting Strategy

Marianna K. Keaveney, Bahar Rahsepar, Hua-an Tseng, Fernando R. Fernandez, Rebecca A. Mount, Tina Ta, John A. White, Jim Berg and Xue Han
Journal of Neuroscience 9 December 2020, 40 (50) 9576-9588; https://doi.org/10.1523/JNEUROSCI.2570-19.2020
Marianna K. Keaveney
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bahar Rahsepar
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Bahar Rahsepar
Hua-an Tseng
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fernando R. Fernandez
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Fernando R. Fernandez
Rebecca A. Mount
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tina Ta
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John A. White
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for John A. White
Jim Berg
2Allen Institute for Brain Science, Seattle, Washington 98109
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xue Han
1Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Single-cell analysis is revealing increasing diversity in gene expression profiles among brain cells. Traditional promotor-based viral gene expression techniques, however, cannot capture the growing variety among single cells. We demonstrate a novel viral gene expression strategy to target cells with specific miRNA expression using miRNA-guided neuron tags (mAGNET). We designed mAGNET viral vectors containing a CaMKIIα promoter and microRNA-128 (miR-128) binding sites, and labeled CaMKIIα+ cells with naturally low expression of miR-128 (Lm128C cells) in male and female mice. Although CaMKIIα has traditionally been considered as an excitatory neuron marker, our single-cell sequencing results reveal that Lm128C cells are CaMKIIα+ inhibitory neurons of parvalbumin or somatostatin subtypes. Further evaluation of the physiological properties of Lm128C cell in brain slices showed that Lm128C cells exhibit elevated membrane excitability, with biophysical properties closely resembling those of fast-spiking interneurons, consistent with previous transcriptomic findings of miR-128 in regulating gene networks that govern membrane excitability. To further demonstrate the utility of this new viral expression strategy, we expressed GCaMP6f in Lm128C cells in the superficial layers of the motor cortex and performed in vivo calcium imaging in mice during locomotion. We found that Lm128C cells exhibit elevated calcium event rates and greater intrapopulation correlation than the overall CaMKIIα+ cells during movement. In summary, the miRNA-based viral gene targeting strategy described here allows us to label a sparse population of CaMKIIα+ interneurons for functional studies, providing new capabilities to investigate the relationship between gene expression and physiological properties in the brain.

SIGNIFICANCE STATEMENT We report the discovery of a class of CaMKIIα+ cortical interneurons, labeled via a novel miRNA-based viral gene targeting strategy, combinatorial to traditional promoter-based strategies. The fact that we found a small, yet distinct, population of cortical inhibitory neurons that express CaMKIIα demonstrates that CaMKIIα is not as specific for excitatory neurons as commonly believed. As single-cell sequencing tools are providing increasing insights into the gene expression diversity of neurons, including miRNA profile data, we expect that the miRNA-based gene targeting strategy presented here can help delineate many neuron populations whose physiological properties can be readily related to the miRNA gene regulatory networks.

  • calcium/calmodulin kinase type II
  • cell type classification
  • GCaMP6
  • optogenetics

Introduction

The brain is composed of many specialized cell subtypes, which together give rise to complex circuit dynamics and functions. Brain cells have traditionally been classified according to many criteria, including location, morphology, neurotransmitter contents, and electrophysiological properties (Zeng and Sanes, 2017). With advances in genetic tools, neuron types have been increasingly defined based on the expression of molecular markers, which allows for detailed characterization of functionality when associated with physiological phenotypes. However, neurons that express a common marker gene can still exhibit great diversity (DeFelipe et al., 2013; Tremblay et al., 2016), and many marker genes are only loosely related to the electrophysiology of the neuron subtype that they label. For example, CaMKIIα, the α-subunit of the serine/threonine kinase CaMKII, has conventionally been used as a marker gene for excitatory neurons (Benson et al., 1992; Jones et al., 1994; Han et al., 2009; Nathanson et al., 2009; Han, 2012; Liu and Murray, 2012). However, the function of CaMKIIα protein has not been specifically linked to the excitatory nature of these neurons.

Transcriptomic profiling has started to reveal that the gene expression diversity of neurons is too great to be encompassed by single marker genes (Roy et al., 2018). Early efforts focused on profiling gene expression from small cell populations taken from precise dissections (Belgard et al., 2011; Hawrylycz et al., 2012) or pools of sorted cells expressing common markers (Sugino et al., 2006; Cahoy et al., 2008; Zhang et al., 2014). However, population-level analyses have since been shown to mask the heterogeneities among individual cells. More recently, single-cell transcriptomic profiling has enabled the characterization of individual brain cells, revealing an increasing number of molecularly distinct classes of neurons (Pollen et al., 2014; Macosko et al., 2015; Usoskin et al., 2015; Zeisel et al., 2015; Tasic et al., 2016, 2018; Economo et al., 2018; Hodge et al., 2019). For example, two recent single-cell sequencing (seq) studies analyzed ∼20,000 single brain cell nuclei, and identified 133 and 75 transcriptomically distinct neuron subtypes in mouse primary visual cortex (Tasic et al., 2018) and in human brain, respectively (Hodge et al., 2019). Combining single-cell transcriptional profiling and axonal reconstructions, Economo et al. (2018) subdivided excitatory pyramidal tract neurons into two subtypes with distinct projection patterns and motor functions. These single-cell studies underscore the genetic and functional diversity among brain cells and highlight the importance of developing new cell targeting strategies that can link transcriptomic profiles and functional properties.

Traditional promotor-based viral labeling techniques, however, cannot capture the increasing gene expression diversity among single cells. Recent advances in sequencing technologies have enabled comprehensive analysis of the relative levels and combinational patterns of microRNA (miRNA) expression with unprecedented scope and accuracy. miRNAs are small, 20- to 24-nt-long, noncoding RNAs that suppress gene expression by hybridizing to their complementary recognition sites within target mRNA transcripts (Bartel, 2009). A given miRNA tends to regulate sets of related genes, and its expression is usually correlated with that of its targets and can thus be closely linked to gene expression profiles. Accumulating evidence suggests that many miRNAs play important roles in regulating synaptic plasticity and neural network dynamics, and as well are implicated in neurologic and psychiatric disorders (Harraz et al., 2011; Xu et al., 2013). For example, miRNA expression has been used to discriminate between neurons and glia cells (Colin et al., 2009; Karali et al., 2011), and has been known to vary significantly across brain regions (Bak et al., 2008; Minami et al., 2014) and between cell types (He et al., 2012). Recently, miRNAs have been used as the basis for genetically targeting known neural subsets(Sayeg et al., 2015; Keaveney et al., 2018), as well as for the diagnosis of Huntington's disease (Reed et al., 2018) and therapy (Pfister et al., 2018). Thus, viral vector designs integrating miRNA regulation could provide a promising new strategy to study the functions of cells with a specific gene expression profile. Here, we demonstrate the use of the miRNA-guided neuron tag (mAGNET; Sayeg et al., 2015; Keaveney et al., 2018), combinatorial to traditional promotor-based gene targeting strategies for targeted gene expression in neurons on the basis of miRNA expression.

Materials and Methods

Viral vector cloning

Standard molecular cloning procedures (PCR and restriction enzyme digestion) were used to construct adeno-associated virus (AAV) vectors. To produce AAV-CaMKIIα-GCaMP6f, the 1.3 kb CaMKIIα promoter (plasmid #64545, Addgene) and the 1.35 kb GCaMP6f gene (plasmid #40755, Addgene; Chen et al., 2013) were amplified via PCR and inserted between the ITRs (inverted terminal repeats) of an AAV-2 genomic vector (plasmid #50954, Addgene; Lin et al., 2013) via EcoRI, AscI, and SpeI restriction sites. To produce AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET, four instances of the microRNA-128 (miR-128) recognition site (exact reverse complement of the mouse miRNA sequence) were added to the 3′ untranslated region (UTR) of the GCaMP6f transgene by annealing single-stranded DNA oligonucleotides synthesized by DNA (Integrated DNA Technologies) and cloning into NotI sites. To produce AAV-CaMKIIα-mRuby2-4 × 128-mAGNET, the mRuby2 reporter (plasmid #54771, Addgene; Lam et al., 2012) was processed by PCR and used to replace GCaMP6f in AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET via AscI and SpeI sites. Viral vector sequences are available at GenBank (accession numbers: MW139899, for AAV-CaMKII-mRuby2-4 × 128-mAGNET; MW139900, for AAV-CaMKII-GCaMP6f; and MW139901, for AAV-CaMKII-GCaMP6f-4 × 128-mAGNET).

Virus preparation

Recombinant AAV2.9 was packaged via triple transfections of the AAV-2 transfer vector (AAV-CaMKIIα-GCaMP6f, AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET, or AAV-CaMKIIα-mRuby2-4 × 128-mAGNET), cap-9 plasmid encoding the Rep and Cap proteins (serotype 9), and pXX680 helper packaging plasmid, into HEK293FT cells (R700-07, Thermo Fisher Scientific) using TransIT-293 transfection reagent (MIR 2700, Mirus Bio). Two days post-transfection, HEK cells were washed and lysed via multiple freeze–thaw cycles and sonication, and virus was purified from the lysate via polyethylene glycol precipitation (LV810A-1, System Biosciences). Viral titer was determined via real time PCR (RT-PCR) (AAVpro Titration Kit version 2, Takara) to be ∼2E12 viral genomes (vg)/ml. Viral transduction was first tested in primary neuron cultures before in vivo injections were performed.

Patch-seq experiments

Tissue processing

All procedures were done in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Allen Institute for Brain Sciences Institutional Animal Care and Use Committees. Mice (male and female) between the ages of postnatal day 45 (P45) and P70 were anesthetized with 5% isoflurane and intracardially perfused with 25 or 50 ml of ice-cold slicing artificial CSF [ACSF; 0.5 mm calcium chloride (dehydrate), 25 mm d-glucose, 20 mm HEPES buffer, 10 mm magnesium sulfate, 1.25 mm sodium phosphate monobasic monohydrate, 3 mm myo-inositol, 12 mm N-acetyl-l-cysteine, 96 mm N-methyl-d-glucamine chloride (NMDG-Cl), 2.5 mm potassium chloride, 25 mm sodium bicarbonate, 5 mm sodium l-ascorbate, 3 mm sodium pyruvate, 0.01 mm taurine, and 2 mm thiourea, pH 7.3, continuously bubbled with 95% O2/5% CO2]. Coronal slices (350 μm) were generated (model VT1200S Vibratome, Leica Biosystems). Slices were then transferred to an oxygenated and warmed (34°C) slicing ACSF for 10 min, then transferred to room temperature holding ACSF [2 mm calcium chloride (dehydrate), 25 mm d-glucose, 20 mm HEPES buffer, 2 mm magnesium sulfate, 1.25 mm sodium phosphate monobasic monohydrate, 3 mm myo-inositol, 12.3 mm N-acetyl-l-cysteine, 84 mm sodium chloride, 2.5 mm potassium chloride, 25 mm sodium bicarbonate, 5 mm sodium L-ascorbate, 3 mm sodium pyruvate, 0.01 mm taurine, and 2 mm thiourea, pH 7.3, continuously bubbled with 95% O2/5% CO2] for the remainder of the day until transferred for patch-clamp recordings.

Patch-clamp recording

Slices were bathed in warm (34°C) recording ACSF [2 mm calcium chloride (dehydrate), 12.5 mm d-glucose, 1 mm magnesium sulfate, 1.25 mm sodium phosphate monobasic monohydrate, 2.5 mm potassium chloride, 26 mm sodium bicarbonate, and 126 mm sodium chloride, pH 7.3, continuously bubbled with 95% O2/5% CO2). The bath solution contained blockers of fast glutamatergic (1 mm kynurenic acid) and GABAergic synaptic transmission (0.1 mm picrotoxin). Thick-walled borosilicate glass (catalog #G150F-3, Warner Instruments) electrodes were manufactured (catalog #PC-10, Narishige) with a resistance of 4–5 MΩ. Before recording, the electrodes were filled with ∼1.0–1.5 µl of internal solution with biocytin [110 mm potassium gluconate, 10.0 mm HEPES, 0.2 mm ethylene glycol-bis (2-aminoethylether)-N,N,N′,N′-tetraacetic acid, 4 mm potassium chloride, 0.3 mm guanosine 5′-triphosphate sodium salt hydrate, 10 mm phosphocreatine disodium salt hydrate, 1 mm adenosine 5′-triphosphate magnesium salt, 20 µg/ml glycogen, 0.5 U/µl RNase inhibitor (Takara, 2313A), and 0.5% biocytin (catalog #B4261, Sigma-Aldrich), pH 7.3]. The pipette was mounted on a Multiclamp 700B amplifier headstage (Molecular Devices) fixed to a micromanipulator (PatchStar, Scientifica).

Electrophysiology signals were recorded using an ITC-18 Data Acquisition Interface (HEKA). Commands were generated, signals processed, and amplifier metadata acquired using MIES (https://github.com/AllenInstitute/MIES/), written in Igor Pro (WaveMetrics). Data were filtered (Bessel) at 10 kHz and digitized at 50 kHz. Data were reported uncorrected for the measured –14 mV liquid junction potential between the electrode and bath solutions.

After formation of a stable seal and break-in, the resting membrane potential of the neuron was recorded (typically within the first minute). A bias current was injected, either manually or automatically using algorithms within the MIES data acquisition package, for the remainder of the experiment to maintain that initial resting membrane potential. Bias currents remained stable for a minimum of 1 s before each stimulus current injection. Upon completion of electrophysiology recording, the nucleus was slowly extracted from the neuron. The pipette containing internal solution, cytosol, and nucleus was removed from the pipette holder, and the contents were expelled into a PCR tube containing the lysis buffer (catalog #634894, Takara).

cDNA amplification and library construction

We used the SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (catalog #634894, Takara) to reverse transcribe poly(A) RNA and amplify full-length cDNA according to the manufacturer instructions. We performed reverse transcription and cDNA amplification for 21 PCR cycles in eight-well strips, in sets of 12–24 strips at a time. At least one control strip was used per amplification set, which contained four wells without cells and four wells with 10 pg control RNA. Control RNA was either Mouse Whole Brain Total RNA (catalog #MR-201, Zyagen) or control RNA provided in the SMART-Seq v4 kit. All samples proceeded through Nextera XT DNA Library Preparation (catalog # FC-131–1096, Illumina) using the Nextera XT Index Kit V2 Set A (FC-131–2001). Nextera XT DNA Library prep was performed according to manufacturer instructions except that the volumes of all reagents including cDNA input were decreased to 0.4× or 0.5× by volume. Subsampling was conducted to a read depth of 0.5 million per cell.

Sequencing data processing

Fifty bp paired-end reads were aligned to GRCm38 (mm10) using a RefSeq annotation gff file retrieved from National Center for Biotechnology Information on January 18, 2016 (https://www.ncbi.nlm.nih.gov/genome/annotation_euk/all/). Sequence alignment was performed using STAR version 2.5.3 (Dobin et al., 2013) in two-pass mode. PCR duplicates were masked and removed using STAR option “bamRemoveDuplicates.” Only uniquely aligned reads were used for gene quantification. Gene counts were computed using the R Genomic Alignments package summarize Overlaps function using “IntersectionNotEmpty” mode for exonic and intronic regions separately (Lawrence et al., 2013).

Mapping to neuron class or subclass

Neurons were mapped to class or subclass as described in the study by Gouwens et al. (2020). The Patch-seq transcriptomes were mapped to the mouse primary visual reference taxonomy tree (Tasic et al., 2018) in a top-down manner, trying to resolve broad classes first, followed by subclasses. For each single-cell transcriptome, starting from the root and at each branch point, we computed its correlation with the reference cell types using the markers associated with the given branch point (i.e., the genes that best distinguished groups at each branch of the tree) and then chose the most correlated branch. To determine the confidence of mapping, we applied 100 bootstrapped iterations at each branch point, and in each iteration, 70% of the reference cells and 70% of the markers were randomly sampled for mapping. This was done to evaluate the degree to which mappings were robust to the effects of individual genes or reference cells. The percentage of times a cell was mapped to a given leaf or branch point in the reference taxonomy was defined as the corresponding mapping probability. For each cell, the mapping probability was sorted, and then the cell was assigned to the class (for glutamatergic neurons) or the subclass (for GABAergic neurons) to which that cell mapped with highest probability. Each cell included in this dataset was mapped with a combined probability of >0.95 to the specific cell types within each assigned class or subclass.

Animal surgery for in vivo calcium imaging and slice physiology

All procedures were performed in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Boston University Institutional Animal Care and Use and Biosafety Committees. First, surgery was performed to inject a mixture of AAV viruses. Mice were adult female C57BL/6 wild-type mice (Taconic Biosciences) that were 8 weeks old at the time of virus injection. AAV viruses (AAV-CaMKIIα-GCaMP6f + AAV-CaMKIIα-mRuby2-4 × 128-mAGNET or AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET + AAV-CaMKIIα-mRuby2-4 × 128-mAGNET) were premixed at equal particle numbers. Next, 2.5E8 vg of each virus (∼600 nl total) was injected into motor cortex (mediolateral, +1.250; anteroposterior, +1.500; dorsoventral, −0.500). Seven days after the first surgery, mice underwent a second surgery for implantation of a sterilized custom imaging window that was previously assembled. The window consisted of a stainless steel imaging cannula [outer diameter (OD), 3.17 mm; inner diameter, 2.36 mm; height, 2 mm; catalog #B004TUE45E, AmazonSupply), fitted with a circular coverslip [size 0; OD, 3 mm; Deckgläser Cover Glasses, catalog #64–0726 (CS-3R-0), Warner Instruments] adhered using a UV-curable optical adhesive (Norland Products). To access the motor cortex, the skull directly above the injection site was carefully removed with a dental drill, then the imaging window was placed flush against the brain tissue and secured in place with Kwik-Sil and acrylic cement. During the same surgery, a custom aluminum head plate was attached to the skull with additional acrylic cement, posterior to the imaging window.

For the patch-clamp recording from parvalbumin (PV) and Lm128C neurons, another set of animals, either PV-tdTomato (n = 4) or PV-cre (n = 3) male and female mice, were injected with AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET or an equal particle mixture of AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET + AAV-flex-tdTomato, respectively. Injections were made at motor cortex (same coordinates as above).

Slice electrophysiology

Coronal sections of motor cortex were prepared from two sets of animals. One group was the CaMKIIα mice coinjected with two AAVs and subjected to in vivo calcium imaging recordings. The other group included PV-cre mice, coinjected with flex-tdTomato and AAV-CaMKIIα-GCaMP6f-4 × 128-mAGNET. Animals were anesthetized with isoflurane, decapitated, and brains were removed and placed in oxygen-bubbled sucrose-substituted ACSF: 185 mm NaCl, 2.5 mm KCl, 1.25 mm NaH2PO4, 10 mm MgCl2, 25 mm NaHCO3, 12.5 mm glucose, and 0.5 mm CaCl2. 400 µm slices were cut via microtome (catalog #VT1200, Leica Microsystems) and incubated in 30°C ACSF consisting of 125 mm NaCl, 25 mm NaHCO3, 25 mm d-glucose, 2 mm KCl, 2 mm CaCl2, 1.25 mm NaH2PO4, and 1 mm MgCl2 for 20 min and then cooled to room temperature.

To differentiate among the different cell types, slice recordings were performed using a two-photon imaging system (Thorlabs) with a mode-locked Ti:Sapphire laser (Chameleon Ultra II, Coherent) set to wavelengths between 920 and 950 nm, which was used to excite the Alexa Fluor 488, GCamP6f and tdTomato using a 20×, numerical aperture 1.0 (Olympus) objective lens. In one group, Lm128C (labeled with the red fluorescent marker mRuby2) and CaMKIIα cells (labeled with green fluorescent GCaMP6f), recordings were made from both green (overall CaMKIIa+ cell population) and yellow cells (CaMKIIα+ and mRuby2+). In the other experimental group, PV cells expressed red fluorescent tdTomato, Lm128C/PV− cells expressed green GCaMP6f, and Lm128C/PV+ cells coexpressed the two markers and appeared yellow.

Slice electrophysiology recording were done using a patch-clamp electrode with tips (0.6–1 μm) pulled via horizontal puller (Sutter Instrument), and pipette resistance between 4 and 6 MΩ. Pipettes were filled with intracellular fluid containing the following: 120 mm K-gluconate, 20 mm KCl, 10 mm HEPES, 7 mm diTrisPhCr, 4 mm Na2ATP, 2 mm MgCl2, 0.3 mm Tris-GTP, and 0.2 mm EGTA, buffered to pH 7.3 with KOH. Alexa Fluor 488 hydrazide (Thermo Fisher Scientific) was added to this solution (0.15% w/v) to visualize the pipette under the two-photon imaging system.

For each recording, pipette capacitance neutralization and bridge balance compensation were performed before breaking in. Data were acquired using Multiclamp 700B (Molecular Devices) and a Digidata 1550 (Molecular Devices) with a sampling rate of 10 kHz. The mean ± SE access resistance of all recordings was 31.6 ± 5.2 MΩ (n = 36 cells).

Electrophysiology data analysis

Brain slice electrophysiology data analysis was performed in python with custom-written scripts using the pyABF package. Spike onsets were identified by a threshold of 10 on their first derivative, with their peak as the maximum value. Spike half-widths are the times between the two halves of the peak voltage amplitude. The rate of rise was calculated by dividing the change in voltage by the time from spike onset to peak. Rate of decay corresponds to time for voltage change from peak to the offset of a spike, defined as the minimum value of the first derivative of the voltage. To calculate the firing frequency versus injected current (FI) gain, the change in firing frequency from lowest firing rate to the highest firing rate was divided by the amount of current injected. Sag was calculated by subtracting the average membrane potential during the first 50 ms of hyperpolarizing current (−50 pA) from the voltage during the last 100 ms of the step. The adaptation ratio was obtained by dividing the mean interspike intervals (ISIs) during the last 200 ms of spiking by the mean ISIs during the first 200 ms of spiking.

Calcium imaging data acquisition

Animal habituation

Following full recovery from the second surgery, animals were habituated to being head fixed and running on a 3D treadmill. The mice were habituated for 5, 10, 15, and 20 min/d, respectively, for 4 d.

Image acquisition with scientific complementary metal–oxide–semiconductor cameras

Animals were positioned underneath the microscope, and imaged while freely navigating the spherical treadmill. Image acquisition was performed via custom microscope equipped with a scientific complementary metal–oxide–semiconductor (sCMOS) camera (ORCA-Flash4.0 LT Digital CMOS camera C11440-42U, Hamamatsu). GCaMP6f fluorescence excitation was accomplished with a 5 W LED (460 nm; model LZ1-00B200, LED Engin). Custom microscope included a Leica N Plan 10 × 0.25 PH1 microscope objective lens, a dual-band excitation filter (catalog #FF01-468/553–25), a dichroic mirror (catalog #FF493/574-Di01-25 × 36), and a dual-band emission filter (catalog #FF01-512/630–25, Semrock), and a commercial SLR lens focused to infinity as the tube lens (80-200 mm f/4 AI-s Zoom NIKKOR lens, Nikon).

A custom MATLAB script was used to trigger frame capture and to synchronize the acquisition of image data with movement information. TTL (transistor–transistor logic) trigger pulses were delivered to the camera using a common input/output interface (catalog #USB-6259, National Instruments). Image acquisition was performed using HC Image Live (Hamamatsu). Because triggering was initiated using a MATLAB software loop, the exact sampling intervals varied based on demands of the Windows 7 operating system. In general, the sampling rate was ∼20 Hz. For each image frame, exposure time was fixed at 20 ms. Image data were stored as the multipage tagged image file format (mpTIFF). For a recording session of 20 min, ∼48 GB of image data were stored, spreading across 12 mpTIFF video files. The acquisition software was configured to buffer all frames in computer RAM during image acquisition to optimize acquisition speed. For each animal, full-session recordings of GCaMP6f fluorescence (20 min) were performed using the specifications noted above.

Movement data acquisition

The spherical treadmill was constructed following the design of Dombeck et al. (2007). Briefly, the treadmill consisted of a 3D printed plastic housing and a Styrofoam ball supported with air. Movement was monitored using two computer USB mouse motion sensors affixed to the plastic housing at the equator of the Styrofoam ball. Each mouse sensor was mounted 3–4 mm away from the surface of the ball to prevent the interference of ball movement. The LED sensors projected on the ball surface 78° apart. The x- and y-surface displacement measured by each mouse sensor was acquired using a separate computer running a Linux OS (minimal CentOS 6). A multithreaded python script that asynchronously reads and accumulates mouse motion events was used to send packaged <dx,dy> data at 100 Hz to the image acquisition computer via a RS232 serial link. Motion data were received on the imaging computer using a MATLAB script that stored the accumulated <dx,dy> between frame triggers synchronized to each acquired image frame.

Calcium imaging data analysis

Motion correction, manual ROI selection, and ΔF/F time series calculation

Motion correction was performed using a custom python script as described in previous studies (Mohammed et al., 2016; Gritton et al., 2019). We first performed a series of image processing procedures to enhance the contrast and the characters of the image. Specifically, for each imaging session, we first generated a reference image by projecting the mean values of every pixel in the first mpTIFF video file saved containing ∼2000 image frames. We then enhanced the contrast of the reference image and each frame as following. We removed the low-frequency component (potential nonuniform background) with a Gaussian filter (python scipy package, ndimage.Gaussian_filter, σ = 50). We then enhanced the boundaries of the high-intensity areas to sharpen the image. In brief, we low-pass filtered the image with a Gaussian filter twice at two levels (σ = 2 and then 1). The intensity differences in two filtered images, which represent the boundaries of high-intensity areas, were multiplied by 100 and added back to the first low pass-filtered image (σ = 2), resulting in a boundary-sharpened image. Finally, to prevent the potential bleaching that may affect the overall intensity of the whole image, we normalized the intensity of each image by shifting the mean intensity to zero and dividing by the SD of the intensity. We then obtained the displacement between the location of the maximum coefficient and the center of the image by calculating the cross-correlations between the enhanced reference image and each enhanced frame. For each frame, the shift that countered the displacement was then applied to the original unenhanced image to complete the motion correction.

We then manually selected ROIs, using the maximum-minus-mean image for each GCaMP6f recording session. The center of each neural soma was registered using a circle with a radius of 6 pixels, corresponding to 7.8 μm. Raw fluorescence traces for each ROI were then extracted from motion-corrected video using a custom MATLAB script. For each ROI, we also normalized the trace by shifting the mean value to zero and scaling so that the maximum intensity equals 100%.

Calcium event identification and waveform calculations

Because of the baseline shift of GCaMP6f fluorescence during each calcium imaging session, and the heterogeneity of signal quality among different ROIs within each field of view, it is difficult to detrend each fluorescence trace without affecting calcium event waveforms. Thus, we identified calcium waveforms using frequency domain features of the GCaMP6f fluorescence traces, where each calcium event is associated with an increase in low-frequency power (<2Hz), because of the time constant of GCaMP6f responses.

We first calculated the spectrogram from raw fluorescence traces (MATLAB chronux, http://chronux.org/; mtspecgramc with tapers = [2 3] and window=[1 0.05]), and averaged the oscillatory power <2 Hz. To detect any significant changes in power, we calculated the change in the power at each time point and identified all positive outliners corresponding to increases in power (MATLAB function “isoutliner”). For outliners that occurred at consecutive time points, we chose the first outliner as the onset of the power change. We then determined the offset of the power change as the first time point where the power decreases. We then mapped the onset and the offset time points of the power change back to the raw fluorescence traces for calcium event detection. We first identified the peak of each calcium event. Briefly, we examined the time points from the offset of the power change, to the second time point when the raw fluorescence trace decreased. Within this time window, we identified the peak fluorescence as the peak of a calcium event. The onset of a calcium event is determined as the time point when the raw fluorescence is at the minimum between the power change onset event and the calcium event peak. We excluded calcium events with an amplitude (the signal difference between the peak and the onset) <3 SDs of the fluorescence trace in a 10 s time window before the calcium event onset. In a few occasions, our event detection procedure may detect the same event more than once; therefore, we further refine the results by combining events that share overlapping time points.

To calculate the rise and fall time of the calcium event of each ROI, we first aligned all calcium events to their peak and averaged the traces within 5 s before and after the peak as the representative waveform (length, 10 s). The rise and fall times were then determined using the representative calcium waveform for each ROI. The rise time is defined as the time from the time point when the intensity increases above threshold, defined as 1 SD above the mean intensity during the 10 s window, to the peak of the waveform. The fall time is the time between the peak and the time point when the intensity dropped below the threshold. We excluded any ROIs that did not have well defined rise or fall times, such as the trace of the waveform did not fall below threshold after peak.

The size of each calcium event is calculated as the area of each calcium waveform above the threshold, calculated as described at each waveform of a given ROI and then averaged across all calcium events for a given ROI.

Pairwise correlation analysis between neuron pairs

We first binarized the calcium intensity trace by converting the time points within the rising phase of a calcium event (from event onset to peak, as identified above) to ones, and all other time points to zeros. The pairwise correlation coefficient from neuron A to neuron B is defined as the time points sharing ones in both A and B, divided by the number of total time points of ones in A. To obtain the representative pairwise correlation coefficient for each neuron, we averaged the pairwise correlation coefficient from the neuron to every other neuron.

Movement data analysis

Spherical treadmill calibration

To map the Styrofoam ball movement to actual mouse linear velocity, we first calibrated the ball. We pinned the two sides of the ball at the equator perpendicular to each of the two computer mouse sensors, and then rotated the ball over a fixed distance perpendicular to the sensor to obtain the translation factors for calculating physical distance. The two computer mouse sensors were calibrated separately. Each computer mouse sensor reading was composed of a vertical reading and a lateral reading.

Linear velocity calculation

We calculated the linear displacement of the animal in two perpendicular directions. One direction (y-axis of animal movement) is calculated from the vertical reading of one sensor. The other direction (x-axis of animal movement) is calculated from the vertical readings of both sensors, using the following equations: Embedded Image Embedded Image

Where L is the vertical reading from the left sensor, R is the vertical reading from the right sensor, and θ is the angle between the two sensors relative to the center of the ball that was fixed at 78°. The changes in linear distance were computed as follows: Embedded Image … D is the change in linear distance between two data sampling points, where each reading is triggered by TTL pulses as described above. To compute linear velocity, we divided the linear distance over time between sample points.

Identification of sustained periods of movement with high and low speed

We identified periods of high versus low linear speed similar to that described by Gritton et al. (2019). Briefly, we first smoothed each linear velocity trace using a low-pass Butterworth filter with a low-pass frequency of 1.5 Hz and then identified the time periods at which the smoothed linear velocity exceeded 5 cm/s for at least 1 s (MATLAB commands “butter” and then “filtfilt” for zero-phase filtering). These periods were considered to be periods of movement with high speed. Similarly, periods during which the smoothed linear velocity of the mouse stayed at <5 cm/s for at least 1 s were considered periods of movement with low speed. If a period overlapped with the first or the last data point of a recording session, it is excluded from further analysis, because of the uncertainty of the start or the end of the period.

Calcium event properties across cell types during different movement periods

Calcium events considered to be in the high-speed or low-speed movement periods were used for event rate and pairwise correlation calculations as long as the onset of a calcium event occurred during one of these time periods. Calcium event properties (event rate and pairwise correlation) were averaged for each ROI, and the differences between high- and low-speed time periods were assessed using a paired t test across all neurons.

Immunohistology

After the conclusion of in vivo calcium imaging recordings, animals were perfused with 4% PFA (except for animals used for slice electrophysiology). Brains were postfixed in 4% PFA at room temperature for 3 h, cryoprotected in 30% sucrose, then flash frozen in OCT (optimal cutting temperature) compound. Frozen brains were sectioned (coronal) at 50 μm on a cryostat, and slices were stored in PBS + 0.05% sodium azide. Brain slices were stained with antibodies against CaMKIIα (1:50; catalog #sc-13 082, Santa Cruz Biotechnology), followed by Alexa Fluor 633 goat anti-rabbit or goat-anti-guinea pig secondary antibody (1:1000; catalog #A21071 or #A21105, Thermo Fisher Scientific). All stained brained slices were imaged via confocal microscopy.

Confocal imaging and quantification

Confocal imaging was performed on an Olympus FV1000 scanning confocal microscope, using a 10× lens or a 60× water-immersion lens. Z-stacks were taken by imaging at 2 μm intervals throughout the 50 µm slices. Consistent laser settings were used for all imaging sessions, as follows: 488 nm laser 3.5%, 350 V; 543 nm laser 13%, 600 V; and 633 nm laser 5%, 600V. GCaMP6f+ and/or mRuby2+ cell bodies were first identified throughout the stack, then categorized as immunopositive or immunonegative for specific antibody makers. Immunopositive cell counts were pooled across slices stained for the same marker for each animal and normalized to the total Lm128C cell count across those slices. In this way, the proportion (percentage) of Lm128C cells that colocalize with each immunomarker was calculated independently. For each animal, we analyzed at least one to four nonoverlapping confocal stacks from two to six slices.

Data availability

The data that support the findings of this study are available from the corresponding author upon request. Viral vector sequences are available from GenBank (accession numbers: MW139899 for AAV-CaMKII-mRuby2-4X128-mAGNET; MW139900 for AAV-CaMKII-GCaMP6f; and MW139901 for AAV-CaMKII-GCaMP6f-4X128-mAGNET).

Results

Viral labeling of a sparse population of CaMKIIα+ and low in miR-128 expression brain cells (Lm128C cells)

CaMKIIα has conventionally been used as a cell type-specific marker gene for excitatory neurons in the cortex (Benson et al., 1992; Jones et al., 1994; Han et al., 2009; Nathanson et al., 2009; Han, 2012; Liu and Murray, 2012), and cortical CaMKIIα+ neurons have been shown to express miR-128 at a high level (He et al., 2012). However, when we incorporated miR-128 target sites in our earlier mAGNET designs, we consistently observed a small fraction of labeled CaMKIIα+ neurons, indicating a lack of miR-128 expression in these CaMKIIα+ neurons (Sayeg et al., 2015). These results suggest that CaMKIIα+ cells can be further divided into subtypes based on their miR-128 expression levels. To characterize this unique neuron subset that is CaMKIIα+ and low in miR-128 expression (low-miRNA-128 CaMKIIα cells, Lm128C cells), we designed an AAV-Lm128C-mAGNET viral construct. AAV-Lm128C-mAGNET uses a 1.3 kb mouse CaMKIIα promoter (Dittgen et al., 2004) and contains four repeats of the miR-128 target site (perfectly complimentary to the mouse miR-128 consensus sequence) within the 3′ UTR of the transgene (Fig. 1A,B). Using this Lm128C-mAGNET design, all infected CaMKIIα+ cells produce mRNA transcripts because of the chosen promoter. However, in cells with a naturally high expression of miR-128, gene transcripts are knocked down, whereas in cells with low levels of miR-128, gene transcripts survive to be translated.

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

Lm128C versus overall CaMKIIα neurons in mouse motor cortex. A, AAV-Lm128C-GCaMP6f-mAGNET and AAV-Lm128C-mRuby-mAGNET vector designs. Lm128C animals (n = 5) were coinjected with these two AAVs for both in vivo calcium activity recordings and static labeling of Lm128C cells. B, AAV-CaMKIIα-GCaMP6f-Control and AAV-Lm128C-mRuby-mAGNET vector designs. CaMKIIα animals (n = 5) were coinjected with these two AAVs for in vivo calcium activity recordings of overall CaMKIIα+ cells and static labeling of Lm128C cells. C, Schematic of stereotaxic virus injection site in mouse motor cortex. D, Representative confocal images of GCaMP6f fluorescence from Lm128C cells, immunofluorescence of a CaMKIIα antibody, and a merged overlay. E, Representative low-magnification (top) and high-magnification (bottom) confocal microscopy images of the virus injection site in a CaMKIIα mouse, showing GCaMP6f expression in overall CaMKIIα+ neurons, mRuby2 expression in Lm128C neurons, and a merged overlay. F, Same as E, for a Lm128C mouse. G, Quantification of labeled Lm128C cells. Percentage of GCaMP6f+ Lm128C cells that are also positive for CaMKIIα immunofluorescence, in Lm128C animals (left, gray); fraction of mRuby2+ Lm128C cells over the total GCaMP6f-expressing cells in CaMKIIα animals (middle, blue); percentage of GCaMP6f+ Lm128C cells that also express mRuby2 in Lm128C animals (right, red). Scatter plot shows individual animals (circles) and mean (horizontal bar) ± SD across animals (n = 5 mice each for each comparison, and at least 50 neurons were quantified per animal).

To statically label Lm128C neurons, we delivered a red fluorescent reporter mRuby2 in the AAV-Lm128C-mAGNET viral vectors (Fig. 1A,B). Additionally, to facilitate in vivo calcium activity recordings from either Lm128C neurons or the overall CaMKIIa+ population, we also codelivered another AAV bearing the calcium indicator GCaMP6f under the regulation of the CaMKIIα promoter, with or without the cassette of miR-128 binding sites (Fig. 1A,B). CaMKIIα mice were injected with AAV-CaMKIIα-GCaMP6f-Control (Fig. 1A) to label all CaMKIIα+ neurons with GCaMP6f. Lm128C mice were injected with AAV-Lm128C-GCaMP6f-mAGNET (Fig. 1B) to label only Lm128C neurons with GCaMP6f. Both groups were also coinjected with AAV-Lm128C-mRuby-mAGNET to label Lm128C cells with mRuby2. All injections consisted of ∼600 µl of a 1:1 virus mixture (1E9 vg/virus) in superficial motor cortex (Fig. 1C).

At the conclusion of in vivo experiments, we processed the brain tissue and characterized mRuby2-labeled Lm128C cells via confocal microscopy. We found that Lm128C neurons are sparsely distributed across all cortical layers in mouse motor cortex (Fig. 1E,F). Staining with an anti-CaMKIIα antibody showed that 99 ± 1% of Lm128C neurons expressed CaMKIIα (Fig. 1D,G), confirming that CaMKIIα promoter exclusively labeled CaMKIIα-expressing neurons. We then further quantified the percentage of GCaMP6f-expressing cells (overall CaMKIIα+ neurons) that also expressed mRuby2 (Lm128C marker) in CaMKIIα animals (n = 1334 cells from five mice) and found that Lm128C neurons constitute 8 ± 2% of the overall CaMKIIα+ neuron population (Fig. 1E,G). Similar colocalization analysis in Lm128C animals revealed that 98 ± 2% of GCaMP6f-expressing neurons (Lm128C marker) also expressed mRuby2 (Lm128C marker; n = 347 cells from five mice), confirming that coinjected AAVs transduced nearly all of the same cells (Fig. 1F,G).

Single-cell transcriptomic analysis reveals that Lm128C cells express CaMKIIα and GABA transporters, and can be mapped to mainly parvalbumin-expressing or somatostatin-expressing subtypes

To assess the gene expression patterns and the corresponding physiological properties of Lm128C cells, we performed single-cell mRNA sequencing of Lm128C cells. We injected mice with AAV-Lm128C-GCaMP6f-mAGNET in the superficial motor cortex region and performed patch-seq analysis of GCaMP6f+ Lm128C cells. We found that 16 of the 18 cells transcriptomically analyzed were GABAergic, expressing glutamate decarboxylase genes [glutamic acid decarboxylase-1 (GAD1) and GAD2], and the GABA transporter gene Slc6a1 (Fig. 2A). Two of the 18 cells analyzed exhibited low levels of expression of Slc17a7, a vesicular glutamate transporter gene, indicative of excitatory nature. Based on the gene clusters generated in our previous study using the single-cell RNA-seq dataset in the mouse visual cortex (Tasic et al., 2018), the GABAergic Lm128C neurons mapped to the somatostatin (SST)-expressing (11 cells) or PV-expressing (5 cells) branches of the transcriptomic tree. The electrophysiological features obtained in these neurons are in general consistent with the faster spiking properties of these interneuron subtypes (Fig. 2B,C).

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

Single-cell sequencing analysis of Lm128C cells. A, Heat map showing gene expression values for individual Lm128C cells in mouse motor cortex slices using the Patch-seq technique (n = 18). Color in top bar indicates the subclass that each cell maps to, and the corresponding dendrogram. Orange, SST; red, PV; green, glutamatergic (Glu). B, Representative electrophysiology traces from cells that map to SST, PV, and glutamatergic subclasses. One-second-long current injections to hyperpolarized potentials show passive properties, while spike train properties are shown with injections for each cell 60 pA above rheobase. C, Box plots of population electrophysiology responses. Top, FI gain, calculated as the change in firing frequency from lowest firing rate to highest firing rate, divided by the current injected. Middle, Rate of spike decay, calculated as the time from peak voltage to the minimum value of the first derivative of voltage. Bottom, Spike half-width, calculated as the time between the two midpoints of the spike waveforms. Individual data points for each cell are overlaid, with color corresponding to the same mapping as in A and B, with gray indicating cells without successful transcriptomic processing.

Consistent with our antibody staining and colocalization analysis with CaMKIIα (Fig. 1D,G), all but one cell analyzed contained significant levels of CaMKIIα mRNA. miR-128 has been shown to regulate the expression of several ion channels that govern membrane excitability (Tan et al., 2013); thus, we focused our analysis on these genes. We observed consistent expression of most ion channel genes that have been shown to be regulated by miR-128, and are known to regulate membrane excitability, such as Cacna1a/b/c/d and Cacna2d3 (families of voltage-gated calcium channels), Kcnq2/3 (potassium channels), and Scn1a/1b/2b (sodium channels; Fig. 2A). The presence of these proteins that are known to be negatively regulated by miR-128 confirmed that Lm128C cells are indeed low in miR-128 expression. Together, these results demonstrate that Lm128C cells exhibit electrophysiological and transcriptomic features consistent with a population of inhibitory interneurons, including both PV and SST subtypes, and with gene expression profiles consistent with low miR-128 expression.

Lm128C neurons exhibit biophysical properties that are characteristic of fast-spiking interneurons

To further characterize the biophysical properties of Lm128C cells, we performed two-photon guided patch-clamp electrophysiological recordings in motor cortex brain slices prepared from mice previously used for in vivo calcium imaging experiments. We identified CaMKIIα cells as those only positive for GCaMP6f, and Lm128C cells as those positive for both GCaMP6f and mRuby2 (as detailed in Fig. 1A). Electrophysiology analysis was performed on neurons in the superficial cortex region, primarily layer 2/3 cells that can also be analyzed using wide-field calcium imaging. We measured the voltage responses to step current injections in both CaMKIIα (Fig. 3A, blue) and Lm128C cells (Fig. 3A, orange). The slope of the spike frequency–current (F–I) relationship was significantly steeper for Lm128C cells (Fig. 3B), suggesting heightened membrane excitability. Additionally, Lm128C cells exhibit narrower spike half-width, and faster spike decay rate than the overall CaMKIIα cell population (Fig. 3Ci). Increased excitability, narrower spikes, and greater afterhyperpolarization are characteristics of fast-spiking interneurons. Together, these results further confirm that the electrophysiological properties of the overall CaMKIIα cells are consistent with those typically observed in pyramidal cells, whereas the properties of Lm128C neurons are consistent with those typically observed in fast-spiking interneurons (McCormick et al., 1985; Suter et al., 2013; Gonzalez-Sulser et al., 2014; Ferrante et al., 2017). miR-128 has been shown to regulate the expression of several ion channels and signaling components of the extracellular signal-regulated kinase 2 network that governs neuronal excitability, and deficiency in this miRNA impacts motor function (Tan et al., 2013). The observed elevated membrane excitability in Lm128C neurons, with naturally low levels of miR-128 expression, is consistent with the role of miR-128 in regulating membrane excitability.

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

Distinct electrophysiological properties of Lm128C and overall CaMKIIα neurons. A, Representative voltage responses (top) of a CaMKIIα neuron (blue) and an Lm128C neuron (orange) to stepped current injections (bottom). Bi, FI for each neuron recorded (n = 5 each). Bii, FI gain (***p = 0.001). Ci, Representative average spike shape of a CaMKIIα neuron (blue) and an Lm128C neuron (orange). Dots indicate spike threshold, and spike waveform midpoint used for calculating spike half-width. Cii, Rate of spike decay (***p = 0.0009). Ciii, Spike half-width (***p = 0.0006). All p values were calculated using two-tailed unpaired Student's t test.

Lm128C interneurons regardless of PV expression have fast-spiking biophysical profiles similar to that of PV+ interneurons

Single-cell sequencing analysis suggested that Lm128C neurons can be mapped to either the PV or the SST subtype based on their transcriptomic properties (Fig. 2). To further assess biophysical differences between the two subclasses of Lm128C cells, we performed patch-clamp recording in PV-tdTomato mice injected with AAV-Lm128C-GCaMP6f-mAGNET (n = 4), and PV-cre mice (n = 3) injected with both AAV-flex-tdTomato and AAV-Lm128C-GCaMP6f-mAGNET. In both cases, PV+ cells expressed tdTomato. Lm128C/PV− neurons, presumable SST-mapped Lm128C cells, are positive for GCaMP6f and negative for tdTomato, LM128c/PV+ neurons are positive for both GCaMP6f and tdTomato, and PV+ cells are positive for tdTomato and negative for GCaMP6f. Electrophysiology results revealed that all three populations exhibited steep F–I curves, narrow spike shapes, and large afterhyperpolarizations following the repolarization of spikes (Fig. 4B,C), consistent with the fast-spiking features, and that observed in Figure 3. Lm128C neurons, both Lm128C/PV− and LM128C/PV+ subtypes, appear to have slightly narrower spikes than PV+ neurons; however, the difference is not significant (Fig. 4C). Together, these results demonstrate that Lm128C cells exhibit fast-spiking electrophysiological properties that are traditionally associated with GABAergic PV+ interneurons.

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

Electrophysiological properties of Lm128C interneuron subtypes. A, Representative voltage responses (top) of a PV+ interneuron (navy), an Lm128C neuron positive for PV (Lm128C/PV+, orange), and an Lm128C neuron negative for PV (Lm128C/PV–, yellow), to stepped current injections (bottom). Bi, FI for each neuron recorded (n = 10 neurons for PV+, n = 9 for LM128C/PV+, and n = 7 for Lm128C/PV–). Bii, FI gain for the three neuron groups, calculated as the change in firing frequency from the lowest firing rate to the highest firing rate divided by the amount of current injected. C, Sag for the three neuron groups. D, Adaptation ratio for the three neuron groups (*p = 0.03). Ei, Representative average spike shapes of the three neuron populations. Eii, Spike half-width, calculated as the time between the two midpoints of the spike waveforms. Eiii, Rate of decay, the time from peak voltage to the minimum value of the first derivative of voltage. p Values were calculated using two-tailed unpaired Student's t test, NS, not significant.

Although both PV and SST cells can exhibit a fast-firing and narrow spike shape phenotype, previous studies have demonstrated that these two cell types can be differentiated by two additional electrophysiological features, the adaptation ratio (Large et al., 2016) and the membrane voltage sag (Yavorska and Wehr, 2016). We estimated the membrane voltage sag associated with a step hyperpolarization as the absolute difference between the steady-state decrease in voltage and the largest decrease in voltage following a hyperpolarizing current step. Specifically, sag is calculated as the average voltage during the first 50 ms in response to a −50 pA current injection, subtracted from the last 100 ms of the −50 pA current injection step. We found no differences among the three groups in terms of the size of the membrane voltage sag (Fig. 4C). In contrast, the adaptation ratio, calculated as the ratio of interspike intervals during the last 100 ms of spiking divided by the first 100 ms of spiking, was significantly greater in LM128C/PV− cells than in LM128C/PV+ neurons (Fig. 4D), consistent with previous observations in SST cells (Large et al., 2016). Together, these results further support the Patch-seq transcriptomic results that Lm128C cells are predominantly inhibitory neurons of the PV and SST subtypes.

Lm128C cells exhibit distinct calcium event dynamics, event rate, and pairwise correlation during locomotion compared with overall CaMKIIα neurons in mouse motor cortex

To illustrate the utility of the mAGNET labeling strategy in functional analysis of neurons during behavior, we performed calcium imaging from Lm128C versus overall CaMKIIa cell populations in mice during locomotion. Since Lm128C cells comprise a minor fraction (∼8%; Fig. 1G) of the neurons labeled in CaMKIIα animals, the effects observed in CaMKIIα animals are primarily dominated by neurons with high levels of miR-128. Thus, our labeling scheme allows us to differentiate the physiological differences between neurons with high and low levels of miR-128 expression.

The experimental timeline for in vivo calcium imaging recordings is summarized in Figure 5D. First, mice were coinjected with a combination of AAVs (Fig. 1A,B), including the genetically encoded calcium indicator GCaMP6f targeted either to the population of overall CaMKIIα neurons (“CaMKIIα” animals) or specifically to CaMKIIα neurons with low miR-128 (“Lm128C” animals) using the AAV-Lm128C-mAGNET vector (Fig. 1A,B). We later implanted a head bar and an imaging window centered over the injection site in motor cortex (Fig. 5A). During in vivo calcium imaging, mice were awake, head fixed, and engaged in voluntary locomotion on a spherical treadmill (Dombeck et al., 2007; Mohammed et al., 2016; Gritton et al., 2019) positioned under a wide-field fluorescence microscope (Fig. 5B,C).

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

In vivo calcium imaging of overall CaMKIIα and Lm128C neurons in motor cortex. A, Schematic of an imaging window and a head bar implanted on the mouse head. B, In vivo calcium imaging paradigm: animals were awake and head fixed under the imaging microscope during voluntary locomotion on a spherical treadmill. Two USB motion sensors were used for movement tracking. C, Custom GCaMP6f imaging setup with an sCMOS camera. D, Experimental timeline for both the CaMKIIα and Lm128C mice groups. E, Representative wide-field images from in vivo calcium imaging recordings (maximum–mean pixel intensity across all frames) for a CaMKIIα mouse (left) and an Lm128C mouse (right). Top, GCaMP6f fluorescence. Bottom, GCaMP6f fluorescence overlaid with manually identified ROIs corresponding to individual GCaMP6f-expressing neurons. F, Representative calcium traces (top) and waveforms (bottom) from a CaMKIIα neuron (left) and an Lm128C neuron (right). Black traces (bottom) represent the average event waveform, and the red line indicates the average activation threshold. G–J, Calcium event rate (G), event rise time (H), and fall time (I), and intrapopulation pairwise correlation for the overall CaMKIIα population and Lm128C neurons (J). Bar plots represent the mean ± SD across recording sessions (n = 10 recording sessions for each group, 2 recording sessions/mouse; n = 5 mice/group), p values were determined via a two-tailed unpaired Student's t test.

We motion corrected GCaMP6f fluorescence videos, manually selected ROIs from the maximum mean fluorescence intensity images (Fig. 5E), extracted the normalized fluorescence intensity traces for these ROIs, and identified individual calcium events. Consistent with the trend observed in the postfixed brain slices (Fig. 1), recordings from CaMKIIα mice contained as many as approximately six times the number of GCaMP6f-expressing neurons as recordings from Lm128C mice (Fig. 5E, compare right, left). We found that the calcium event rate was significantly elevated for Lm128C neurons, compared with the overall CaMKIIα population (Fig. 5G). In addition, the durations of calcium events were significantly shorter for Lm128C neurons compared with the CaMKIIα population, with faster rise and fall times (Fig. 5F,H,I), suggesting that these neurons exhibit faster intracellular calcium dynamics. Furthermore, Lm128C neurons show significantly elevated pairwise correlation with other Lm128C cells when compared with the correlation among the overall population of CaMKIIα neurons (Fig. 5J). Overall, the features exhibited by Lm128C neurons in behaving mice, including elevated event frequency, faster event dynamics, and greater synchrony, are typically associated with nonpyramidal interneurons in the cortex (Chen et al., 1996; Hu and Agmon, 2015), consistent with our Patch-seq transcriptomic analysis and electrophysiological characterizations in brain slices (Figs. 2–4).

Lm128C cells are distinctly modulated by movement compared with overall CaMKIIα neurons during locomotion

Because motor cortex neurons are known to be modulated by movement, we examined whether the overall CaMKIIα neurons and Lm128C neurons exhibit distinct relationships with movement parameters. We identified periods when mice were moving at high speed (>5 cm/s for at least 1 s) versus at low speed (<5 cm/s for at least 1 s) and compared calcium responses of the overall CaMKIIα cells and Lm128C neurons during these time periods (Fig. 6A). We found that Lm128C neurons exhibited significantly elevated event rate during periods of high-speed movement, whereas the overall CaMKIIα neuron population exhibited similar event rates during both periods (Fig. 6B). In addition, both the overall CaMKIIα population and Lm128C neurons exhibited larger event size, calculated as the area under the waveform, during high-speed movement compared with that during low-speed movement (Fig. 6C). These results are consistent with the general observation that many motor cortex neurons increase their activity when movement is more rigorous (Armstrong and Drew, 1984; Ebina et al., 2018). The fact that Lm128C neurons showed a significant increase in event rate, but not the overall CaMKIIα neuron population, may be related to the fact that Lm128C neurons have shorter event duration. Interestingly, we also found that the overall CaMKIIα population exhibited greater pairwise correlation during high-speed movement compared with that during low-speed movement, whereas Lm128C cells showed no change in pairwise correlation (Fig. 6D). Together, these results demonstrate that the overall CaMKIIα population and Lm128C neurons were both positively modulated by movement speed. The overall CaMKIIα population could further increase their intrapopulation correlation when movement speed increases, whereas Lm128C neurons showed no change in correlation when movement speed changes.

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

Differential modulation of calcium activity by movement speed in the overall CaMKIIα population and Lm128C neurons in the motor cortex. A1, A2, Representative examples of mouse movement speed and corresponding calcium traces during a recording session, for a CaMKIIα mouse (A1) and a Lm128C mouse (A2). Periods of high-speed movement are highlighted in red, and periods of low-speed movement are highlighted in green. B–D, Event rate (events per minute; B), event size (area under the curve; C), and intrapopulation pairwise correlation (D) during periods of high-speed and low-speed movement for the overall CaMKIIα population and Lm128C neurons. Bar plots present the mean ± SE (CaMKIIα: n = 1850 neurons for B and D, n = 1659 neurons for C, Lm128C: n = 330 neurons for B and D, n = 322 neurons for C; n = 5 mice for CaMKIIα group, and n = 5 mice for Lm128C group; two-tailed paired Student's t test), NS, not significant.

Discussion

We here demonstrate a novel viral gene expression strategy, using mAGNET designs, combinatorial with traditional promotor-based strategies, for targeted gene expression in neurons with specific endogenous miRNA expression profiles. This strategy enables in vivo labeling and functional characterization of unique neuron populations. Using this viral strategy, we labeled a subtype of CaMKIIα+ inhibitory interneurons on the basis of miRNA profile (low miR-128) and illustrated their unique electrophysiological and calcium activity patterns in mouse motor cortex. This study provides direct evidence that CaMKIIα is present in a set of inhibitory neurons, and identified a novel neural subtype characterized on the basis of miRNA profile in vivo. Importantly, this study opens the door for a new avenue of studying cells with distinct expression profiles through miRNA regulation. As single-cell analysis tools provide increasing insight into the gene expression diversity of neurons, including miRNA profile data, we expect that miRNAs will emerge as important cellular markers. Furthermore, as the regulation networks of individual miRNAs are further explored, this will in turn allow the design of mAGNETs to delineate unique neuron subtypes whose physiological properties can be readily related to the miRNA gene regulatory networks dictated by the classification method.

We designed mAGNET AAV viral vectors containing the CaMKIIα promoter to label a subset of CaMKIIα+ neurons with naturally low expression of miR-128 (low-miRNA-128 CaMKIIα cells, Lm128C cells). CaMKII is a serine/threonine kinase essential for synaptic plasticity and memory, and a major component of the postsynaptic density in many neurons (Herring and Nicoll, 2016). The α subunit (CaMKIIα) is generally regarded as a marker gene for pyramidal excitatory neurons in the mammalian cortex (Benson et al., 1992; Jones et al., 1994; Han et al., 2009; Nathanson et al., 2009; Han, 2012; Liu and Murray, 2012), although little is known about the gene expression diversity among cortical CaMKIIα+ neurons. Among the miRNAs expressed in mammalian neocortical CaMKIIα+ neurons, miR-128 has been shown to be the second most highly expressed (He et al., 2012) and, overall, is one of the most abundant miRNAs in the mammalian brain (Lagos-Quintana et al., 2002). This miRNA has been shown to regulate gene networks that govern neuronal excitability(Tan et al., 2013), and the fast-spiking feature observed in Lm128C cells is consistent with the role of miR-128 in regulating membrane excitability. Thus, miRNA-based viral vector designs integrating miRNA regulation and promoter sequences represent a novel strategy to study the functions of cells with more defined gene expression profiles.

It is interesting to observe the expression of CaMKIIα in GABAergic cortical interneurons, whose properties have not been specifically characterized previously. Over the years, in situ hybridization and immunolabelling of CaMKIIα and Gad1 (the rate-limiting enzyme for GABA production) in the mammalian cortex have suggested that these proteins are expressed in mutually exclusive cells (Benson et al., 1992; Jones et al., 1994). For these reasons, CaMKIIα has come to be associated with cortical excitatory cells. The fact that we found a small, yet distinct, population of cortical inhibitory neurons that express CaMKIIα demonstrates that CaMKIIα is not as specific for excitatory neurons as commonly believed. Although CaMKIIα has been associated with principal excitatory neurons in the mammalian cortex, its primary function in mediating postsynaptic long-term potentiation is not inherently specific to excitatory neurons (Herring and Nicoll, 2016). While extensive work has been directed toward the involvement of CaMKII at excitatory synapses, this does not preclude the expression of CaMKII in interneurons.

Similarly, PV has been widely adopted as a marker gene for a subset of GABAergic cells known as fast-spiking interneurons (Celio, 1986) and has therefore been assumed to not express CaMKIIα. However, CaMKIIα and GAD are known to colocalize in neurons within some brain regions, such as the stria terminalis and cerebellum (Benson et al., 1992; Liu and Murray, 2012). Several recent studies suggest that there may be greater heterogeneity and overlap between cortical neurons expressing traditional excitatory and inhibitory genetic markers than previously thought. For example, Nathanson et al. (2009) reported that some GAD67+ neurons were labeled in mouse somatosensory cortex with the CaMKIIα promoter delivered from an AAV2.1 vector, although they did not further characterize these cells. Furthermore, single-cell transcriptomic profiling of cortical neurons has recently been used to identify seven transcriptionally distinct subsets of PV+ neurons in mouse primary visual cortex and acknowledged the existence of cells that coexpress a combination of markers including the excitatory marker Slc7a7 as well as the inhibitory markers Gad1, Gad2, and PV (Tasic et al., 2016).

Footnotes

  • X.H. was supported by National Institutes of Health (NIH) Grants 1DP2-NS-082126 and 1R01-NS-087950-01), and National Science Foundation Grant CBET-1848029; and M.K.K. was supported by NIH Grant F31-NS-095465-02. Patch-seq experiments were funded by the Allen Institute for Brain Science. We thank Drs. Howard Gritton and Mike Romano for guidance regarding in vivo calcium imaging recordings and data analysis, and Dr. Kyle Hansen for sharing MATLAB code facilitating manual selection of ROIs from preprocessed calcium images. We also thank the LAS team, Tissue Processing, and Molecular Biology teams at the Allen Institute for Brain Science for experimental support; and Brian Lee, Kristen Hadley, Rusty Mann, Lindsay Ng, Ram Rajanbabu, Jessica Trinh, and Katherine Baker at the Allen Institute for Brain Science for contributing toward Patch-seq recordings.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Xue Han at xuehan{at}bu.edu

SfN exclusive license.

References

  1. ↵
    1. Armstrong DM,
    2. Drew T
    (1984) Discharges of pyramidal tract and other motor cortical neurones during locomotion in the cat. J Physiol 346:471–495. doi:10.1113/jphysiol.1984.sp015036 pmid:6699782
    OpenUrlCrossRefPubMed
  2. ↵
    1. Bak M,
    2. Silahtaroglu A,
    3. Møller M,
    4. Christensen M,
    5. Rath MF,
    6. Skryabin B,
    7. Tommerup N,
    8. Kauppinen S
    (2008) MicroRNA expression in the adult mouse central nervous system. RNA 14:432–444. doi:10.1261/rna.783108 pmid:18230762
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Bartel DP
    (2009) MicroRNAs: target recognition and regulatory functions. Cell 136:215–233. doi:10.1016/j.cell.2009.01.002 pmid:19167326
    OpenUrlCrossRefPubMed
  4. ↵
    1. Belgard TG,
    2. Marques AC,
    3. Oliver PL,
    4. Abaan HO,
    5. Sirey TM,
    6. Hoerder-Suabedissen A,
    7. García-Moreno F,
    8. Molnár Z,
    9. Margulies EH,
    10. Ponting CP
    (2011) A transcriptomic atlas of mouse neocortical layers. Neuron 71:605–616. doi:10.1016/j.neuron.2011.06.039 pmid:21867878
    OpenUrlCrossRefPubMed
  5. ↵
    1. Benson DL,
    2. Isackson PJ,
    3. Gall CM,
    4. Jones EG
    (1992) Contrasting patterns in the localization of glutamic acid decarboxylase and Ca2+/calmodulin protein kinase gene expression in the rat central nervous system. Neuroscience 46:825–849. doi:10.1016/0306-4522(92)90188-8 pmid:1311814
    OpenUrlCrossRefPubMed
  6. ↵
    1. Cahoy JD,
    2. Emery B,
    3. Kaushal A,
    4. Foo LC,
    5. Zamanian JL,
    6. Christopherson KS,
    7. Xing Y,
    8. Lubischer JL,
    9. Krieg PA,
    10. Krupenko SA,
    11. Thompson WJ,
    12. Barres BA
    (2008) A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J Neurosci 28:264–278. doi:10.1523/JNEUROSCI.4178-07.2008 pmid:18171944
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Celio MR
    (1986) Parvalbumin in most gamma-aminobutyric acid-containing neurons of the rat cerebral cortex. Science 231:995–997. doi:10.1126/science.3945815 pmid:3945815
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Chen TW,
    2. Wardill TJ,
    3. Sun Y,
    4. Pulver SR,
    5. Renninger SL,
    6. Baohan A,
    7. Schreiter ER,
    8. Kerr RA,
    9. Orger MB,
    10. Jayaraman V,
    11. Looger LL,
    12. Svoboda K,
    13. Kim DS
    (2013) Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499:295–300. doi:10.1038/nature12354 pmid:23868258
    OpenUrlCrossRefPubMed
  9. ↵
    1. Chen W,
    2. Zhang JJ,
    3. Hu GY,
    4. Wu CP
    (1996) Electrophysiological and morphological properties of pyramidal and nonpyramidal neurons in the cat motor cortex in vitro. Neuroscience 73:39–55. doi:10.1016/0306-4522(96)00009-7 pmid:8783228
    OpenUrlCrossRefPubMed
  10. ↵
    1. Colin A,
    2. Faideau M,
    3. Dufour N,
    4. Auregan G,
    5. Hassig R,
    6. Andrieu T,
    7. Brouillet E,
    8. Hantraye P,
    9. Bonvento G,
    10. Déglon N
    (2009) Engineered lentiviral vector targeting astrocytes in vivo. Glia 57:667–679. doi:10.1002/glia.20795 pmid:18942755
    OpenUrlCrossRefPubMed
  11. ↵
    1. DeFelipe J,
    2. López-Cruz PL,
    3. Benavides-Piccione R,
    4. Bielza C,
    5. Larrañaga P,
    6. Anderson S,
    7. Burkhalter A,
    8. Cauli B,
    9. Fairén A,
    10. Feldmeyer D,
    11. Fishell G,
    12. Fitzpatrick D,
    13. Freund TF,
    14. González-Burgos G,
    15. Hestrin S,
    16. Hill S,
    17. Hof PR,
    18. Huang J,
    19. Jones EG,
    20. Kawaguchi Y, et al
    . (2013) New insights into the classification and nomenclature of cortical GABAergic interneurons. Nat Rev Neurosci 14:202–216. doi:10.1038/nrn3444 pmid:23385869
    OpenUrlCrossRefPubMed
  12. ↵
    1. Dittgen T,
    2. Nimmerjahn A,
    3. Komai S,
    4. Licznerski P,
    5. Waters J,
    6. Margrie TW,
    7. Helmchen F,
    8. Denk W,
    9. Brecht M,
    10. Osten P
    (2004) Lentivirus-based genetic manipulations of cortical neurons and their optical and electrophysiological monitoring in vivo. Proc Natl Acad Sci U S A 101:18206–18211. doi:10.1073/pnas.0407976101 pmid:15608064
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Dobin A,
    2. Davis CA,
    3. Schlesinger F,
    4. Drenkow J,
    5. Zaleski C,
    6. Jha S,
    7. Batut P,
    8. Chaisson M,
    9. Gingeras TR
    (2013) STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29:15–21. doi:10.1093/bioinformatics/bts635 pmid:23104886
    OpenUrlCrossRefPubMed
  14. ↵
    1. Dombeck DA,
    2. Khabbaz AN,
    3. Collman F,
    4. Adelman TL,
    5. Tank DW
    (2007) Imaging large-scale neural activity with cellular resolution in awake, mobile mice. Neuron 56:43–57. doi:10.1016/j.neuron.2007.08.003
    OpenUrlCrossRefPubMed
  15. ↵
    1. Ebina T,
    2. Masamizu Y,
    3. Tanaka YR,
    4. Watakabe A,
    5. Hirakawa R,
    6. Hirayama Y,
    7. Hira R,
    8. Terada SI,
    9. Koketsu D,
    10. Hikosaka K,
    11. Mizukami H,
    12. Nambu A,
    13. Sasaki E,
    14. Yamamori T,
    15. Matsuzaki M
    (2018) Two-photon imaging of neuronal activity in motor cortex of marmosets during upper-limb movement tasks. Nat Commun 9:1879. doi:10.1038/s41467-018-04286-6 pmid:29760466
    OpenUrlCrossRefPubMed
  16. ↵
    1. Economo MN,
    2. Viswanathan S,
    3. Tasic B,
    4. Bas E,
    5. Winnubst J,
    6. Menon V,
    7. Graybuck LT,
    8. Nguyen TN,
    9. Smith KA,
    10. Yao Z,
    11. Wang L,
    12. Gerfen CR,
    13. Chandrashekar J,
    14. Zeng H,
    15. Looger LL,
    16. Svoboda K
    (2018) Distinct descending motor cortex pathways and their roles in movement. Nature 563:79–84. doi:10.1038/s41586-018-0642-9 pmid:30382200
    OpenUrlCrossRefPubMed
  17. ↵
    1. Ferrante M,
    2. Tahvildari B,
    3. Duque A,
    4. Hadzipasic M,
    5. Salkoff D,
    6. Zagha EW,
    7. Hasselmo ME,
    8. McCormick DA
    (2017) Distinct functional groups emerge from the intrinsic properties of molecularly identified entorhinal interneurons and principal cells. Cereb Cortex 27:3186–3207. doi:10.1093/cercor/bhw143 pmid:27269961
    OpenUrlCrossRefPubMed
  18. ↵
    1. Gonzalez-Sulser A,
    2. Parthier D,
    3. Candela A,
    4. McClure C,
    5. Pastoll H,
    6. Garden D,
    7. Surmeli G,
    8. Nolan MF
    (2014) GABAergic projections from the medial septum selectively inhibit interneurons in the medial entorhinal cortex. J Neurosci 34:16739–16743. doi:10.1523/JNEUROSCI.1612-14.2014
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Gouwens NW,
    2. Sorensen SA,
    3. Baftizadeh F,
    4. Budzillo A,
    5. Lee BR,
    6. Jarsky T,
    7. Alfiler L,
    8. Arkhipov A,
    9. Baker K,
    10. Barkan E,
    11. Berry K,
    12. Bertagnolli D,
    13. Bickley K,
    14. Bomben J,
    15. Braun T,
    16. Brouner K,
    17. Casper T,
    18. Crichton K,
    19. Daigle TL,
    20. Dalley R, et al
    . (2020) Toward an integrated classification of neuronal cell types: morphoelectric and transcriptomic characterization of individual GABAergic cortical neurons. bioRxiv. Advance online publication. Retrieved November 8, 2020. doi: 10.1101/2020.02.03.932244. doi:10.1101/2020.02.03.932244
    OpenUrlCrossRef
  20. ↵
    1. Gritton HJ,
    2. Howe WM,
    3. Romano MF,
    4. DiFeliceantonio AG,
    5. Kramer MA,
    6. Saligrama V,
    7. Bucklin ME,
    8. Zemel D,
    9. Han X
    (2019) Unique contributions of parvalbumin and cholinergic interneurons in organizing striatal networks during movement. Nat Neurosci 22:586–597. doi:10.1038/s41593-019-0341-3 pmid:30804530
    OpenUrlCrossRefPubMed
  21. ↵
    1. Han X
    (2012) Optogenetics in the nonhuman primate. Prog Brain Res 196:215–233. doi:10.1016/B978-0-444-59426-6.00011-2 pmid:22341328
    OpenUrlCrossRefPubMed
  22. ↵
    1. Han X,
    2. Qian X,
    3. Bernstein JG,
    4. Zhou HH,
    5. Franzesi GT,
    6. Stern P,
    7. Bronson RT,
    8. Graybiel AM,
    9. Desimone R,
    10. Boyden ES
    (2009) Millisecond-timescale optical control of neural dynamics in the nonhuman primate brain. Neuron 62:191–198. doi:10.1016/j.neuron.2009.03.011 pmid:19409264
    OpenUrlCrossRefPubMed
  23. ↵
    1. Harraz MM,
    2. Dawson TM,
    3. Dawson VL
    (2011) MicroRNAs in Parkinson's disease. J Chem Neuroanat 42:127–130. doi:10.1016/j.jchemneu.2011.01.005 pmid:21295133
    OpenUrlCrossRefPubMed
  24. ↵
    1. Hawrylycz MJ,
    2. Lein ES,
    3. Guillozet-Bongaarts AL,
    4. Shen EH,
    5. Ng L,
    6. Miller JA,
    7. van de Lagemaat LN,
    8. Smith KA,
    9. Ebbert A,
    10. Riley ZL,
    11. Abajian C,
    12. Beckmann CF,
    13. Bernard A,
    14. Bertagnolli D,
    15. Boe AF,
    16. Cartagena PM,
    17. Chakravarty MM,
    18. Chapin M,
    19. Chong J,
    20. Dalley RA, et al
    . (2012) An anatomically comprehensive atlas of the adult human brain transcriptome. Nature 489:391–399. doi:10.1038/nature11405 pmid:22996553
    OpenUrlCrossRefPubMed
  25. ↵
    1. He M,
    2. Liu Y,
    3. Wang X,
    4. Zhang MQ,
    5. Hannon GJ,
    6. Huang ZJ
    (2012) Cell-type-based analysis of microRNA profiles in the mouse brain. Neuron 73:35–48. doi:10.1016/j.neuron.2011.11.010 pmid:22243745
    OpenUrlCrossRefPubMed
  26. ↵
    1. Herring BE,
    2. Nicoll RA
    (2016) Long-term potentiation: from CaMKII to AMPA receptor trafficking. Annu Rev Physiol 78:351–365. doi:10.1146/annurev-physiol-021014-071753 pmid:26863325
    OpenUrlCrossRefPubMed
  27. ↵
    1. Hodge RD,
    2. Bakken TE,
    3. Miller JA,
    4. Smith KA,
    5. Barkan ER,
    6. Graybuck LT,
    7. Close JL,
    8. Long B,
    9. Johansen N,
    10. Penn O,
    11. Yao Z,
    12. Eggermont J,
    13. Hollt T,
    14. Levi BP,
    15. Shehata SI,
    16. Aevermann B,
    17. Beller A,
    18. Bertagnolli D,
    19. Brouner K,
    20. Casper T, et al
    . (2019) Conserved cell types with divergent features in human versus mouse cortex. Nature 573:61–68. doi:10.1038/s41586-019-1506-7 pmid:31435019
    OpenUrlCrossRefPubMed
  28. ↵
    1. Hu H,
    2. Agmon A
    (2015) Properties of precise firing synchrony between synaptically coupled cortical interneurons depend on their mode of coupling. J Neurophysiol 114:624–637. doi:10.1152/jn.00304.2015 pmid:25972585
    OpenUrlCrossRefPubMed
  29. ↵
    1. Jones EG,
    2. Huntley GW,
    3. Benson DL
    (1994) Alpha calcium/calmodulin-dependent protein kinase II selectively expressed in a subpopulation of excitatory neurons in monkey sensory-motor cortex: comparison with GAD-67 expression. J Neurosci 14:611–629. doi:10.1523/JNEUROSCI.14-02-00611.1994
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Karali M,
    2. Manfredi A,
    3. Puppo A,
    4. Marrocco E,
    5. Gargiulo A,
    6. Allocca M,
    7. Corte MD,
    8. Rossi S,
    9. Giunti M,
    10. Bacci ML,
    11. Simonelli F,
    12. Surace EM,
    13. Banfi S,
    14. Auricchio A
    (2011) MicroRNA-restricted transgene expression in the retina. PLoS One 6:e22166. doi:10.1371/journal.pone.0022166 pmid:21818300
    OpenUrlCrossRefPubMed
  31. ↵
    1. Keaveney MK,
    2. Tseng HA,
    3. Ta TL,
    4. Gritton HJ,
    5. Man HY,
    6. Han X
    (2018) A microRNA-based gene-targeting tool for virally labeling interneurons in the rodent cortex. Cell Rep 24:294–303. doi:10.1016/j.celrep.2018.06.049 pmid:29996091
    OpenUrlCrossRefPubMed
  32. ↵
    1. Lagos-Quintana M,
    2. Rauhut R,
    3. Yalcin A,
    4. Meyer J,
    5. Lendeckel W,
    6. Tuschl T
    (2002) Identification of tissue-specific microRNAs from mouse. Curr Biol 12:735–739. doi:10.1016/S0960-9822(02)00809-6 pmid:12007417
    OpenUrlCrossRefPubMed
  33. ↵
    1. Lam AJ,
    2. St-Pierre F,
    3. Gong Y,
    4. Marshall JD,
    5. Cranfill PJ,
    6. Baird MA,
    7. McKeown MR,
    8. Wiedenmann J,
    9. Davidson MW,
    10. Schnitzer MJ,
    11. Tsien RY,
    12. Lin MZ
    (2012) Improving FRET dynamic range with bright green and red fluorescent proteins. Nat Methods 9:1005–1012. doi:10.1038/nmeth.2171 pmid:22961245
    OpenUrlCrossRefPubMed
  34. ↵
    1. Large AM,
    2. Kunz NA,
    3. Mielo SL,
    4. Oswald AM
    (2016) Inhibition by somatostatin interneurons in olfactory cortex. Front Neural Circuits 10:62. doi:10.3389/fncir.2016.00062 pmid:27582691
    OpenUrlCrossRefPubMed
  35. ↵
    1. Lawrence M,
    2. Huber W,
    3. Pagès H,
    4. Aboyoun P,
    5. Carlson M,
    6. Gentleman R,
    7. Morgan MT,
    8. Carey VJ
    (2013) Software for computing and annotating genomic ranges. PLoS Comput Biol 9:e1003118. doi:10.1371/journal.pcbi.1003118 pmid:23950696
    OpenUrlCrossRefPubMed
  36. ↵
    1. Lin JY,
    2. Knutsen PM,
    3. Muller A,
    4. Kleinfeld D,
    5. Tsien RY
    (2013) ReaChR: a red-shifted variant of channelrhodopsin enables deep transcranial optogenetic excitation. Nat Neurosci 16:1499–1508. doi:10.1038/nn.3502 pmid:23995068
    OpenUrlCrossRefPubMed
  37. ↵
    1. Liu XB,
    2. Murray KD
    (2012) Neuronal excitability and calcium/calmodulin-dependent protein kinase type II: location, location, location. Epilepsia 53 [Suppl 1]:45–52. doi:10.1111/j.1528-1167.2012.03474.x pmid:22612808
    OpenUrlCrossRefPubMed
  38. ↵
    1. Macosko EZ,
    2. Basu A,
    3. Satija R,
    4. Nemesh J,
    5. Shekhar K,
    6. Goldman M,
    7. Tirosh I,
    8. Bialas AR,
    9. Kamitaki N,
    10. Martersteck EM,
    11. Trombetta JJ,
    12. Weitz DA,
    13. Sanes JR,
    14. Shalek AK,
    15. Regev A,
    16. McCarroll SA
    (2015) Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell 161:1202–1214. doi:10.1016/j.cell.2015.05.002 pmid:26000488
    OpenUrlCrossRefPubMed
  39. ↵
    1. McCormick DA,
    2. Connors BW,
    3. Lighthall JW,
    4. Prince DA
    (1985) Comparative electrophysiology of pyramidal and sparsely spiny stellate neurons of the neocortex. J Neurophysiol 54:782–806. doi:10.1152/jn.1985.54.4.782 pmid:2999347
    OpenUrlCrossRefPubMed
  40. ↵
    1. Minami K,
    2. Uehara T,
    3. Morikawa Y,
    4. Omura K,
    5. Kanki M,
    6. Horinouchi A,
    7. Ono A,
    8. Yamada H,
    9. Ohno Y,
    10. Urushidani T
    (2014) miRNA expression atlas in male rat. Sci Data 1:140005. doi:10.1038/sdata.2014.5 pmid:25977763
    OpenUrlCrossRefPubMed
  41. ↵
    1. Mohammed AI,
    2. Gritton HJ,
    3. Tseng HA,
    4. Bucklin ME,
    5. Yao Z,
    6. Han X
    (2016) An integrative approach for analyzing hundreds of neurons in task performing mice using wide-field calcium imaging. Sci Rep 6:20986. doi:10.1038/srep20986 pmid:26854041
    OpenUrlCrossRefPubMed
  42. ↵
    1. Nathanson JL,
    2. Yanagawa Y,
    3. Obata K,
    4. Callaway EM
    (2009) 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
  43. ↵
    1. Pfister EL,
    2. DiNardo N,
    3. Mondo E,
    4. Borel F,
    5. Conroy F,
    6. Fraser C,
    7. Gernoux G,
    8. Han X,
    9. Hu D,
    10. Johnson E,
    11. Kennington L,
    12. Liu P,
    13. Reid SJ,
    14. Sapp E,
    15. Vodicka P,
    16. Kuchel T,
    17. Morton AJ,
    18. Howland D,
    19. Moser R,
    20. Sena-Esteves M, et al
    . (2018) Artificial miRNAs reduce human mutant Huntingtin throughout the striatum in a transgenic sheep model of Huntington's disease. Hum Gene Ther 29:663–673. doi:10.1089/hum.2017.199 pmid:29207890
    OpenUrlCrossRefPubMed
  44. ↵
    1. Pollen AA,
    2. Nowakowski TJ,
    3. Shuga J,
    4. Wang X,
    5. Leyrat AA,
    6. Lui JH,
    7. Li N,
    8. Szpankowski L,
    9. Fowler B,
    10. Chen P,
    11. Ramalingam N,
    12. Sun G,
    13. Thu M,
    14. Norris M,
    15. Lebofsky R,
    16. Toppani D,
    17. Kemp DW 2nd.,
    18. Wong M,
    19. Clerkson B,
    20. Jones BN, et al
    . (2014) Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex. Nat Biotechnol 32:1053–1058. doi:10.1038/nbt.2967 pmid:25086649
    OpenUrlCrossRefPubMed
  45. ↵
    1. Reed ER,
    2. Latourelle JC,
    3. Bockholt JH,
    4. Bregu J,
    5. Smock J,
    6. Paulsen JS,
    7. Myers RH
    (2018) MicroRNAs in CSF as prodromal biomarkers for Huntington disease in the PREDICT-HD study. Neurology 90:e264–e272. doi:10.1212/WNL.0000000000004844
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Roy AL,
    2. Conroy R,
    3. Smith J,
    4. Yao Y,
    5. Beckel-Mitchener AC,
    6. Anderson JM,
    7. Wilder EL
    (2018) Accelerating a paradigm shift: the common fund single cell analysis program. Sci Adv 4:eaat8573. doi:10.1126/sciadv.aat8573 pmid:30083611
    OpenUrlFREE Full Text
  47. ↵
    1. Sayeg MK,
    2. Weinberg BH,
    3. Cha SS,
    4. Goodloe M,
    5. Wong WW,
    6. Han X
    (2015) Rationally designed microRNA-based genetic classifiers target specific neurons in the brain. ACS Synth Biol 4:788–795. doi:10.1021/acssynbio.5b00040 pmid:25848814
    OpenUrlCrossRefPubMed
  48. ↵
    1. Sugino K,
    2. Hempel CM,
    3. Miller MN,
    4. Hattox AM,
    5. Shapiro P,
    6. Wu C,
    7. Huang ZJ,
    8. Nelson SB
    (2006) Molecular taxonomy of major neuronal classes in the adult mouse forebrain. Nat Neurosci 9:99–107. doi:10.1038/nn1618 pmid:16369481
    OpenUrlCrossRefPubMed
  49. ↵
    1. Suter BA,
    2. Migliore M,
    3. Shepherd GM
    (2013) Intrinsic electrophysiology of mouse corticospinal neurons: a class-specific triad of spike-related properties. Cereb Cortex 23:1965–1977. doi:10.1093/cercor/bhs184 pmid:22761308
    OpenUrlCrossRefPubMed
  50. ↵
    1. Tan CL,
    2. Plotkin JL,
    3. Venø MT,
    4. von Schimmelmann M,
    5. Feinberg P,
    6. Mann S,
    7. Handler A,
    8. Kjems J,
    9. Surmeier DJ,
    10. O'Carroll D,
    11. Greengard P,
    12. Schaefer A
    (2013) MicroRNA-128 governs neuronal excitability and motor behavior in mice. Science 342:1254–1258. doi:10.1126/science.1244193 pmid:24311694
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Tasic B,
    2. Menon V,
    3. Nguyen TN,
    4. Kim TK,
    5. Jarsky T,
    6. Yao Z,
    7. Levi B,
    8. Gray LT,
    9. Sorensen SA,
    10. Dolbeare T,
    11. Bertagnolli D,
    12. Goldy J,
    13. Shapovalova N,
    14. Parry S,
    15. Lee C,
    16. Smith K,
    17. Bernard A,
    18. Madisen L,
    19. Sunkin SM,
    20. Hawrylycz M, et al
    . (2016) Adult mouse cortical cell taxonomy revealed by single cell transcriptomics. Nat Neurosci 19:335–346. doi:10.1038/nn.4216 pmid:26727548
    OpenUrlCrossRefPubMed
  52. ↵
    1. Tasic B,
    2. Yao Z,
    3. Graybuck LT,
    4. Smith KA,
    5. Nguyen TN,
    6. Bertagnolli D,
    7. Goldy J,
    8. Garren E,
    9. Economo MN,
    10. Viswanathan S,
    11. Penn O,
    12. Bakken T,
    13. Menon V,
    14. Miller J,
    15. Fong O,
    16. Hirokawa KE,
    17. Lathia K,
    18. Rimorin C,
    19. Tieu M,
    20. Larsen R, et al
    . (2018) Shared and distinct transcriptomic cell types across neocortical areas. Nature 563:72–78. doi:10.1038/s41586-018-0654-5 pmid:30382198
    OpenUrlCrossRefPubMed
  53. ↵
    1. Tremblay R,
    2. Lee S,
    3. 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
  54. ↵
    1. Usoskin D,
    2. Furlan A,
    3. Islam S,
    4. Abdo H,
    5. Lönnerberg P,
    6. Lou D,
    7. Hjerling-Leffler J,
    8. Haeggström J,
    9. Kharchenko O,
    10. Kharchenko PV,
    11. Linnarsson S,
    12. Ernfors P
    (2015) Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. Nat Neurosci 18:145–153. doi:10.1038/nn.3881 pmid:25420068
    OpenUrlCrossRefPubMed
  55. ↵
    1. Xu B,
    2. Hsu PK,
    3. Stark KL,
    4. Karayiorgou M,
    5. Gogos JA
    (2013) Derepression of a neuronal inhibitor due to miRNA dysregulation in a schizophrenia-related microdeletion. Cell 152:262–275. doi:10.1016/j.cell.2012.11.052 pmid:23332760
    OpenUrlCrossRefPubMed
  56. ↵
    1. Yavorska I,
    2. Wehr M
    (2016) Somatostatin-expressing inhibitory interneurons in cortical circuits. Front Neural Circuits 10:76. doi:10.3389/fncir.2016.00076 pmid:27746722
    OpenUrlCrossRefPubMed
  57. ↵
    1. Zeisel A,
    2. Muñoz-Manchado AB,
    3. Codeluppi S,
    4. Lönnerberg P,
    5. Manno GL,
    6. Juréus A,
    7. Marques S,
    8. Munguba H,
    9. He L,
    10. Betsholtz C,
    11. Rolny C,
    12. Castelo-Branco G,
    13. Hjerling-Leffler J,
    14. Linnarsson S
    (2015) Brain structure. Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. Science 347:1138–1142. doi:10.1126/science.aaa1934 pmid:25700174
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Zeng H,
    2. Sanes JR
    (2017) Neuronal cell-type classification: challenges, opportunities and the path forward. Nat Rev Neurosci 18:530–546. doi:10.1038/nrn.2017.85 pmid:28775344
    OpenUrlCrossRefPubMed
  59. ↵
    1. Zhang YK,
    2. Chen K,
    3. Sloan SA,
    4. Bennett ML,
    5. Scholze AR,
    6. O'Keeffe S,
    7. Phatnani HP,
    8. Guarnieri P,
    9. Caneda C,
    10. Ruderisch N,
    11. Deng S,
    12. Liddelow SA,
    13. Zhang C,
    14. Daneman R,
    15. Maniatis T,
    16. Barres BA,
    17. Wu JQ
    (2014) An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci 34:11929–11947. doi:10.1523/JNEUROSCI.1860-14.2014 pmid:25186741
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

The Journal of Neuroscience: 40 (50)
Journal of Neuroscience
Vol. 40, Issue 50
9 Dec 2020
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience 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.
CaMKIIα-Positive Interneurons Identified via a microRNA-Based Viral Gene Targeting Strategy
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
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
CaMKIIα-Positive Interneurons Identified via a microRNA-Based Viral Gene Targeting Strategy
Marianna K. Keaveney, Bahar Rahsepar, Hua-an Tseng, Fernando R. Fernandez, Rebecca A. Mount, Tina Ta, John A. White, Jim Berg, Xue Han
Journal of Neuroscience 9 December 2020, 40 (50) 9576-9588; DOI: 10.1523/JNEUROSCI.2570-19.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
Request Permissions
Share
CaMKIIα-Positive Interneurons Identified via a microRNA-Based Viral Gene Targeting Strategy
Marianna K. Keaveney, Bahar Rahsepar, Hua-an Tseng, Fernando R. Fernandez, Rebecca A. Mount, Tina Ta, John A. White, Jim Berg, Xue Han
Journal of Neuroscience 9 December 2020, 40 (50) 9576-9588; DOI: 10.1523/JNEUROSCI.2570-19.2020
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

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

Keywords

  • calcium/calmodulin kinase type II
  • cell type classification
  • GCaMP6
  • optogenetics

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 Articles

  • Sex differences in histamine regulation of striatal dopamine
  • Gestational Chlorpyrifos Exposure Imparts Lasting Alterations to the Rat Somatosensory Cortex
  • The Neurobiology of Cognitive Fatigue and Its Influence on Effort-Based Choice
Show more Research Articles

Cellular/Molecular

  • Sex differences in histamine regulation of striatal dopamine
  • A Critical Role of Neuroligin 2 C-Terminus in OCD and Social Behavior
  • Time-Dependent Actions of Corticosterone on Infralimbic Cortex Pyramidal Neurons of Adult Male Rats
Show more Cellular/Molecular
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • 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
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci 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 JNeurosci.