GIRK2 splice variants and neuronal G protein-gated K+ channels: implications for channel function and behavior

Many neurotransmitters directly inhibit neurons by activating G protein-gated inwardly rectifying K+ (GIRK) channels, thereby moderating the influence of excitatory input on neuronal excitability. While most neuronal GIRK channels are formed by GIRK1 and GIRK2 subunits, distinct GIRK2 isoforms generated by alternative splicing have been identified. Here, we compared the trafficking and function of two isoforms (GIRK2a and GIRK2c) expressed individually in hippocampal pyramidal neurons lacking GIRK2. GIRK2a and GIRK2c supported comparable somato-dendritic GIRK currents in Girk2 −/− pyramidal neurons, although GIRK2c achieved a more uniform subcellular distribution in pyramidal neurons and supported inhibitory postsynaptic currents in distal dendrites better than GIRK2a. While over-expression of either isoform in dorsal CA1 pyramidal neurons restored contextual fear learning in a conditional Girk2 −/− mouse line, GIRK2a also enhanced cue fear learning. Collectively, these data indicate that GIRK2 isoform balance within a neuron can impact the processing of afferent inhibitory input and associated behavior.

and electrophysiological studies involving wild-type and Girk2 −/− mice have confirmed the presence of GIRK2 in mouse hippocampal neurons 10,15 . Using RNA-Seq and laser-captured tissue samples from the adult mouse, we probed for Girk2a and Girk2c mRNAs in the CA1 region of the hippocampus. Because the 3′-UTR regions of Girk2a and Girk2c mRNAs are distinct and may influence their subcellular trafficking 29 , we compared transcript levels in cell body (stratum pyramidale) and neuropil (stratum radiatum) samples. Both transcripts were detected in neuropil, albeit at significantly lower levels than seen in CA1 cell body samples (q = 0.0012 for Girk2a and 0.0125 for Girk2c), suggesting that an active RNA trafficking mechanism exists for both isoforms ( Fig. 2A). Higher levels of Girk2c mRNA relative to Girk2a mRNA were detected, with the difference reaching statistical significance in neuropil (P = 0.04) but not cell bodies (P = 0.07). Interestingly, we detected considerable read coverage past the annotated terminal exon for both transcripts, suggesting transcription occurs past the annotated 3′ ends in hippocampal pyramidal neurons.
We next assessed whether GIRK2a and GIRK2c proteins are present in mouse hippocampus and cultured hippocampal neurons. A commercially-available GIRK2 antibody (pan-GIRK2) targeting a domain in both isoforms revealed a doublet in intact hippocampus and cultured hippocampal neurons from wild-type but not Girk2 −/− mice (Fig. 2B). We also used a custom GIRK2c-specific antibody, raised against the unique 11-amino acid domain of GIRK2c, to probe for the GIRK2c isoform in hippocampal samples. The GIRK2c antibody recognized GIRK2c, but not other GIRK subunits, when expressed in HEK cells ( Supplementary Fig. S1). In hippocampal samples, this antibody revealed a single band that aligned with the upper band in the doublet seen with the pan-GIRK2 1, P < 0.001) were observed, but there was no interaction between isoform and compartment (F 1,8 = 0.07, P = 0.80). Symbols: *P < 0.05 vs. GIRK2a (within compartment); +,++ P < 0.05 and 0.01, respectively, vs. cell body (within isoform). (B) Representative GIRK2 immunoblots (from 3 independent experiments) of adult hippocampus and neonatal hippocampal culture (10-12 DIV) from wildtype (WT) and Girk2 −/− mice, probed with pan-GIRK2 or GIRK2c antibodies (upper blots), as well as a β-actin antibody (lower blots). (C,D) Representative images showing GIRK2 immunolabeling (as revealed with the pan-GIRK2 antibody) in cultured hippocampal neurons from wild-type (C) and Girk2 −/− (D) mice. Neurons were co-stained with a MAP2 antibody to highlight neuronal morphology and dendritic processes. Scale bar: 50 microns.

Subcellular distribution of GIRK2a and GIRK2c in Girk2 −/− pyramidal neurons. While GIRK2a
and GIRK2c exhibited comparable labeling patterns when expressed in cultured rat hippocampal neurons 22 , endogenous GIRK2 might mask isoform-dependent differences in the trafficking and/or function of exogenously expressed subunits. Thus, we sought to examine the subcellular distribution of each isoform when expressed individually in cultured hippocampal neurons from Girk2 −/− mice. AAV vectors harboring coding sequence for GIRK2a or GIRK2c, inserted downstream of the CaMKIIα promoter, were used to drive expression of untagged GIRK2a or GIRK2c (and EGFP, bi-cistronically) in Girk2 −/− neurons. Strong overlap between recombinant GIRK2 and endogenous CaMKIIα immunolabeling was observed ( Supplementary Fig. S2), confirming that the expression of GIRK2 was restricted primarily to pyramidal neurons.
Consistent with a previous report 35 , the pan-GIRK2 antibody yielded a punctate labeling pattern in wild-type cultures, with extensive labeling seen in dendritic processes (Fig. 2C). No such labeling was seen in hippocampal cultures from Girk2 −/− mice (Fig. 2D), confirming antibody specificity. The labeling patterns for recombinant GIRK2a and GIRK2c expressed in Girk2 −/− pyramidal neurons were overlapping but distinct (Fig. 3). While both subunits were detected in the soma and dendrites of infected neurons, GIRK2a labeling was more restricted to the cell body and proximal dendrites (Fig. 3A). In contrast, GIRK2c labeling was uniformly distributed throughout the neuron, including distal dendrites (Fig. 3B). Quantification of GIRK2 isoform labeling confirmed that GIRK2c exhibited a higher density of labeling than GIRK2a in secondary and tertiary dendrites of infected hippocampal pyramidal neurons (Fig. 3C).
GIRK3 has been shown to mediate interactions between GIRK channels and proteins that can regulate GIRK channel trafficking (sorting nexin 27, or SNX-27) 36 and receptor coupling efficiency (regulator of G protein signaling 2, or RGS2) 37 . As GIRK3 is also expressed in hippocampal neurons 7, 15, 38 , we examined the subcellular distribution of GIRK2a and GIRK2c expressed in cultured hippocampal pyramidal neurons from mice lacking both GIRK2 and GIRK3 (Girk2 −/− /Girk3 −/− mice). The labeling patterns for GIRK2a and GIRK2c in neurons from Girk2 −/− /Girk3 −/− mice were comparable to those seen in Girk2 −/− neurons ( Supplementary Fig. S3), indicating that the GIRK2 isoform-dependent difference in subcellular distribution is not dependent on GIRK3.
GIRK2-containing channels are enriched near excitatory synapses 15,39 , and some overlap between GIRK2 and PSD-95 was reported in cultured hippocampal neurons 35 . As the PDZ interaction motif found in GIRK2c may promote the association of GIRK channels with excitatory synapses, we next compared the overlap between The insets highlight different densities of GIRK2a and GIRK2c puncta along proximal/primary (Ai,Bi) and distal/ secondary (Aii,Bii) dendritic segments. Scale bars: 5 microns. (C) Quantification of GIRK2a and GIRK2c labeling in dendrites from infected Girk2 −/− pyramidal neurons. GIRK2 fluorescence intensity was measured in 2-3 primary (t 27 = 1.4, P = 0.17), secondary (t 29 = 2.1, *P < 0.05), and tertiary (t 27 = 5.2, ***P < 0.001) dendritic segments from 5 different neurons expressing each subunit. Fluorescence intensity is expressed as arbitrary units (AU) normalized to segment length. Only GIRK2 labeling that overlapped directly with the dendritic segment (identified by MAP2 labeling) was quantified in this analysis. endogenous PSD-95 and either GIRK2a or GIRK2c in Girk2 −/− pyramidal neurons (Fig. 4A,B). No significant difference in the extent of overlap between PSD-95 and either GIRK2a or GIRK2c was detected at baseline (Fig. 4C, control). We also explored the possibility that incubation with morphine, a treatment previously shown to increase the overlap between GIRK2 and PSD-95 in cultured hippocampal neurons 35 , might reveal an isoform-dependent difference. Consistent with published observations, we found that morphine treatment increased total GIRK2 protein level by approximately 2-fold relative to vehicle-treated controls ( Supplementary  Fig. S4A,B), while having no apparent effect on PSD-95 labeling intensity ( Supplementary Fig. S4C). Morphine treatment increased the overlap between PSD-95 and both GIRK2 isoforms relative to control conditions, with GIRK2c showing a small but significantly higher degree of overlap with PSD-95 than did GIRK2a following treatment (Fig. 4C, morphine; Supplementary Fig. S4D,E). GIRK currents measured in Girk2 −/− pyramidal neurons expressing GIRK2a or GIRK2c. We next asked whether GIRK2a and GIRK2c conferred unique functionality to GIRK channels in cultured hippocampal pyramidal neurons. Whole-cell currents evoked by GABA B R stimulation were measured in wild-type and Girk2 −/− pyramidal neurons expressing EGFP, and in Girk2 −/− pyramidal neurons expressing GIRK2a or GIRK2c. Baclofen evoked a prominent inward current in wild-type pyramidal neurons that was not seen in Girk2 −/− control neurons expressing EGFP alone (Fig. 5A,B). Viral expression of either GIRK2 isoform in Girk2 −/− pyramidal neurons yielded baclofen-induced currents that were larger than currents in wild-type controls. Activation (but not deactivation) kinetics of the baclofen-induced currents in Girk2 −/− pyramidal neurons expressing GIRK2a or GIRK2c also differed from those in wild-type control neurons, though no isoform-dependent differences were found (Fig. 5C,D). Similarly, while there was no isoform-dependent difference in channel sensitivity (EC 50 ) to GABA B R activation, GABA B R-GIRK coupling efficiency in both cases was diminished relative to wild-type controls (Fig. 5E).
Girk2 −/− pyramidal neurons expressing GIRK2a or GIRK2c, and in wild-type control neurons (Fig. 5J), indicating that viral expression of either GIRK2 isoform restores a normal level of GIRK1-containing heteromeric channel activity.
Synaptic GIRK currents in Girk2 −/− pyramidal neurons expressing GIRK2a or GIRK2c. The lack of natural connectivity and normal afferent inputs to cultured neurons complicates efforts to study the full range of physiological contributions made by GIRK channels to pyramidal neuron excitability. Thus, we next examined the functional properties of GIRK2a and GIRK2c in organotypic hippocampal slices. Organotypic hippocampal slices were cultured from young wild-type or Girk2 −/− mice (P8-9), and the CA1 region was targeted with control or GIRK2 isoform-specific expression viruses (Fig. 6A). Following a 1-wk incubation period, we measured somatodendritic and synaptic GIRK-dependent responses in infected dorsal CA1 pyramidal neurons. Consistent with our observations in infected Girk2 −/− pyramidal neurons in culture, we observed larger baclofen-induced somatodendritic currents in dorsal Girk2 −/− CA1 pyramidal neurons expressing either GIRK2a or GIRK2c, as compared to wild-type controls (Fig. 6B).
Slow inhibitory postsynaptic currents (slow IPSCs) in dorsal CA1 pyramidal neurons were evoked by electrical stimulation within proximal (Schaeffer collateral/SC) or distal (perforant path/PP) dendritic fields. SC and PP slow IPSCs were reliably observed in wild-type CA1 pyramidal neurons (Fig. 6C), and were blocked completely by the GABA B R receptor antagonist CGP54626 (2 μM, not shown). Slow IPSCs were not observed in Girk2 −/− CA1 neurons expressing EGFP. While viral expression of either GIRK2a or GIRK2c in Girk2 −/− slices restored both SC and PP slow IPSCs in CA1 pyramidal neurons, the relative amplitudes of slow IPSCs measured in SC and PP differed between the groups (Fig. 6C). Indeed, the ratio of slow IPSC amplitudes, measured in PP and SC for each neuron, was significantly smaller in Girk2 −/− CA1 pyramidal neurons expressing GIRK2a, as compared to both wild-type control neurons, and Girk2 −/− pyramidal neurons expressing GIRK2c (Fig. 6D). Recently, we reported that mice lacking GIRK2 in forebrain pyramidal neurons (CaMKIICre(+):Girk2 fl/fl mice) were deficient in contextual fear learning, a hippocampal-dependent behavior 6 . As CaMKIICre(+):Girk2 fl/fl mice exhibited a loss of GIRK-dependent signaling in dorsal but not ventral CA1 pyramidal neurons 6 , we asked whether viral restoration of GIRK2a or GIRK2c in dorsal CA1 pyramidal neurons of CaMKIICre(+):Girk2 fl/fl mice could rescue normal fear learning. Adult CaMKIICre(+):Girk2 fl/fl mice received bilateral infusions of Cre-dependent variants of the GIRK2a and GIRK2c expression viruses into the dorsal CA1 (Fig. 7A). After a suitable post-surgical recovery period (2 wk), we tested the efficacy of our viral treatment by measuring baclofen-induced somato-dendritic currents in CA1 pyramidal neurons in acutely isolated slices. Currents in dorsal CA1 pyramidal neurons from CaMKIICre(+):Girk2 fl/fl mice treated with control virus were small (Fig. 7B,C), and comparable in size to those reported previously for untreated CaMKIICre(+):Girk2 fl/fl mice 6 . Dorsal CA1 pyramidal neurons from CaMKIICre(+):Girk2 fl/fl mice expressing GIRK2a or GIRK2c exhibited significantly larger baclofen-induced currents (Fig. 7B,C). Indeed, currents measured in dorsal CA1 pyramidal neurons expressing GIRK2a or GIRK2c were notably larger than those reported previously in wild-type or CaMKIICre(−):Girk2 fl/fl controls 6,15 , suggesting that the GIRK2 isoforms were over-expressed in these neurons.
We utilized a trace fear conditioning protocol to assess the impact of viral treatment on fear learning behavior, as both contextual and cue fear learning in this model are supported by the dorsal hippocampus 44,45 . Untreated

Discussion
In this study, we probed for differences between GIRK2a and GIRK2c in an array of expression systems. Our findings suggest that these two isoforms differ primarily in terms of subcellular trafficking, with GIRK2c achieving a more uniform distribution throughout neurons, including distal dendrites. This difference in subcellular trafficking likely underlies the differing abilities of GIRK2a and GIRK2c to support slow IPSCs evoked via the electrical stimulation of proximal and distal dendritic fields. Indeed, our findings suggest that GIRK2c is critical for ensuring robust synaptic inhibitory GIRK-dependent responses in the distal dendritic fields of CA1 pyramidal neurons.
We observed some interesting differences between GIRK-dependent signaling in transfected HEK cells and infected Girk2 −/− pyramidal neurons. Notably, GIRK1/GIRK2 channels expressed in HEK cells were more sensitive to GABA B R activation (i.e., they exhibited a lower EC 50 ) than the GIRK channels found in wild-type pyramidal neurons or Girk2 −/− neurons expressing GIRK2a or GIRK2c. This difference could be explained by a higher density of receptor and channel expressed in HEK cells, and/or by different receptor-channel ratios achieved in the different cell systems. The difference in coupling efficiency could also reflect the presence of negative regulatory elements in neurons that are not found in HEK cells. Indeed, we reported previously that a complex consisting of RGS7 and Gβ5 decreases the sensitivity of GIRK channels to GABA B R activation in hippocampal neurons [46][47][48] , while also prominently accelerating GABA B R-GIRK current deactivation kinetics. In support of this contention, the deactivation kinetics of GABA B R-GIRK currents were significantly faster in neurons than in HEK cells.
Another difference between HEK cell and neuron models is the presence of GIRK3, which has been reported to decrease the activity and surface trafficking of GIRK1/GIRK2 heteromers and GIRK2 homomers in expression systems 22,49 . We show here that while the relative dendritic labeling patterns for GIRK2a and GIRK2c were comparable in Girk2 −/− and Girk2 −/− /Girk3 −/− neurons, GPCR-GIRK currents measured in Girk2 −/− /Girk3 −/− neurons expressing GIRK2a or GIRK2c were larger than those measured in Girk2 −/− neurons. These results are consistent with the possible role of GIRK3 as a negative regulator of GIRK channel surface trafficking. In midbrain DA neurons, GIRK3 mediates an interaction between the GIRK channel, a GIRK2c/GIRK3 heteromer 13 , and SNX-27 36, 50, 51 , a protein implicated in the surface trafficking and internalization of several channels and receptors [52][53][54][55][56] . The interaction between SNX-27 and GIRK3 is conferred by a Class I PDZ interaction domain at the GIRK3 C-terminus, which is identical to that found in GIRK2c 33 . At present, the relevance of SNX-27 and GIRK3 to GIRK channel trafficking in hippocampal neurons is unclear.
Previous immunoelectron microscopy studies have shown that GIRK1 and GIRK2 labeling is concentrated at the peri-and extra-synaptic positions relative to the postsynaptic density (PSD) in hippocampal CA1. In addition, GIRK2 (but not GIRK1) labeling was detected within the PSD of CA1 pyramidal neurons 15 . The presence of a PDZ domain on the C-terminus of GIRK2c suggested that this isoform might promote the distribution of GIRK channels within or around the PSD, by virtue of its interaction with a scaffolding protein such as PSD-95. As noted above, GIRK2c was reported to bind PSD-95 31 , as well as SAP-97 34 . Moreover, SAP97 conferred both basal activity and G protein-dependent gating to GIRK2c-containing channels expressed in Xenopus oocytes 34 . Our data here, however, show that GIRK2a and GIRK2c assume comparable subcellular distributions relative to PSD-95 in mouse hippocampal pyramidal neurons. Thus, alternative splicing of the coding sequence of the Girk2 gene does not seem to influence the position of GIRK channels relative to the PSD. We did, however, detect significantly more Girk2c than Girk2a mRNA in the neuropil of CA1 in vivo, suggesting that RNA localization may play a role in setting the amount of GIRK2 isoforms in CA1 dendrites. Moreover, our RNA-Seq data suggest that the 3′ terminal exons of both isoforms may be extended in vivo, which could further regulate the targeting, translation, and turnover of these transcripts. Follow up studies looking at the significance of this extended transcription are needed.
The level of GIRK1 protein is significantly reduced in the hippocampus of constitutive Girk2 −/− mice 15,24 . Our experiments with the GIRK1-selective agonist ML297, however, suggest that viral expression of GIRK2 isoforms in Girk2 −/− pyramidal neurons restored a normal level of functional GIRK1/GIRK2 heteromeric channels. Baclofen-induced currents in Girk2 −/− pyramidal neurons expressing either GIRK2a or GIRK2c were significantly larger than those measured in wild-type controls, however, indicating that some of the GABA B R-GIRK current was carried by GIRK2 homomeric or GIRK2/GIRK3 heteromeric channels. Interestingly, expression of GIRK1 in VTA dopamine neurons, which normally express a GIRK2c/GIRK3 heteromeric channel, increased GIRK channel sensitivity to GABA B R activation 37 . Moreover, GIRK2/GIRK3 heteromers are less sensitive to G protein activation than GIRK2 homomers 37 . Thus, an altered proportion of GIRK1-containing and GIRK1-lacking channels in Girk2 −/− pyramidal neurons expressing GIRK2a or GIRK2c could explain why some features of the baclofen-induced current, including the EC 50 , differed from those in wild-type neurons. Notably, adenosine and 5HT 1A R-dependent currents in Girk2 −/− neurons expressing GIRK2a or GIRK2c did not differ in terms of amplitude from wild-type controls, suggesting that adenosine and 5HT 1A receptors may couple preferentially to GIRK1-containining channels, while GABA B R may be able to couple effectively to GIRK channels irrespective of channel subtype.
Multiple forms of plasticity involving GIRK-dependent signaling in hippocampal neurons have been reported. For example, Jan and colleagues showed that stimuli capable of inducing long-term potentiation (LTP) also enhanced GABA B R-GIRK signaling in CA1 pyramidal neurons 57 . Subsequently, this group showed that NMDA receptor activation triggered a rapid increase in surface trafficking of GIRK channels in hippocampal cultures 58,59 . This adaptation involved the dephosphorylation of an N-terminal GIRK2 residue (Ser-9) found in both GIRK2a and GIRK2c, and yielded enhanced basal GIRK channel activity and GIRK channel responses to adenosine receptor activation, but had no effect on GABA B R-GIRK signaling. Slesinger and colleagues showed that prolonged morphine treatment in hippocampal cultures shifted the subcellular distribution of GIRK2 from dendritic shafts to spines, leading to enhanced basal and 5HT-induced GIRK currents, but no change in GABA B R-GIRK signaling 35 . Our data suggest that this morphine-induced adaptation does not require elements within the Girk2 promoter or untranslated domains in the Girk2 mRNAs, as only coding sequence was included in the viral reconstitution vectors. The unique GIRK2c C-terminus is also not required, as both GIRK2a and GIRK2c were similarly impacted by morphine treatment. Interestingly, the strengthening of GIRK-dependent signaling in VTA dopamine neurons evoked by burst firing was blocked by inclusion of a peptide mimic of the GIRK2c (and GIRK3) C-terminal domain 60 . These data indicate that the GIRK2c C-terminal domain may be a critical determinant of plasticity involving specific GIRK channel subtypes, including the GIRK1-lacking channels expressed in midbrain dopamine neurons 13,31 .
Constitutive Girk2 −/− mice exhibit an array of neurological phenotypes, including impaired fear learning 6,61 . The global nature of the genetic manipulation, however, complicates attempts to link the phenotype to a defined neuron population. Accordingly, we developed conditional Girk2 −/− (Girk2 fl/fl ) mice and have crossed them with transgenic Cre driver lines to permit GIRK channel ablation in a neuron-specific manner 6,61 . In CaMKIICre(+):Girk2 fl/fl mice, Girk2 ablation is driven by Cre recombinase, expressed under the control of the CaMKIIα promoter. While the CaMKIICre driver line is generally considered to promote gene ablation in postnatal forebrain pyramidal neurons, the scope of recombination achieved with this driver line depends on the "floxed" gene in question, with some studies showing highly restricted target loss in dorsal CA1 pyramidal neurons and others showing recombination in multiple forebrain pyramidal neuron populations [62][63][64][65][66][67] . We found that GIRK2 expression was preserved in most forebrain structures in CaMKIICre(+):Girk2 fl/fl mice, but that GIRK-dependent signaling was blunted in dorsal (but not ventral) CA1 pyramidal neurons 6 .
To further probe the contribution of GIRK-dependent signaling in the hippocampus to fear learning, we assessed the performance of CaMKIICre(+):Girk2 fl/fl mice in trace fear conditioning test. In this paradigm, the trace interval inserted between the presentation of the conditioned (cue/tone) and unconditioned (footshock) stimulus activates temporal memory mechanisms within the hippocampus, thus rendering both contextual and cue fear learning dependent on the hippocampus 44,68 . Viral-mediated over-expression of either GIRK2a or GIRK2c in dorsal CA1 pyramidal neurons rescued contextual fear learning in CaMKIICre(+):Girk2 fl/fl mice, suggesting that the fear learning deficit in CaMKIICre(+):Girk2 fl/fl mice is linked to GIRK2 loss in these neurons. Although we did not see a difference in cue fear learning between untreated CaMKIICre(−):Girk2 fl/fl and CaMKIICre(+):Girk2 fl/fl mice, over-expression of GIRK2a, but not GIRK2c, in dorsal CA1 pyramidal neurons significantly enhanced cue learning. The enhanced cue learning seen in CaMKIICre(+):Girk2 fl/fl mice over-expressing GIRK2a in dorsal CA1 pyramidal neurons, together with our evidence that the GIRK2 isoforms contribute in distinct fashion to the processing of inhibitory input to proximal and distal dendritic fields, suggests that trace cue learning may be differentially modulated by inputs from the PP and SC.
In summary, we demonstrate in this study that alternative splicing of the mouse Girk2 gene generates at least two variants, and that these splice variants can give rise to channel subtypes that differ in terms of their subcellular trafficking and the neurophysiological processes in which they participate. As GIRK-dependent signaling plays an important role in moderating glutamatergic input to neurons, our findings raise the intriguing possibility that altering the balance of GIRK2 splice variants could represent a means by which neurons fine-tune the processing of excitatory synaptic input, in a subcellular compartment-dependent manner.

RNA-Seq and analysis.
Total RNA was extracted from laser-captured tissue microdissected from hippocampal CA1 cell body and neuropil layers from 3 adult (6-8 wk) male Amigo2-EGFP transgenic mice using the micro RNeasy kit (QIAGEN; Hilden, Germany), including on-column DNase-treatment. Total RNA samples were analysed for RNA integrity (RIN > 8.5) and concentration using the 2100 Bioanalyzer instrument (Agilent; Santa Clara, CA) and RNA 6000 Pico assay (Agilent). Stranded RNA-Seq libraries were prepared using 1-5 ng of total RNA per sample and region (1 ng for neuropil samples and 5 ng for cell body samples) with the Ovation RNA-Seq Systems 1-16 for mouse (NuGEN; San Carlos, CA), according to manufacturer instructions. Libraries were analyzed for size and concentration using the 2100 Agilent Bioanalyzer and the High Sensitivity DNA assay. Libraries were multiplexed and run on a NextSeq500 instrument (Illumina; San Diego, CA), acquiring an average of 66 ± 17 million 100 bp reads per sample. Reads were trimmed using Sickle software (https:// github.com/najoshi/sickle) and only paired reads with a quality score >20 and a minimum length of 20 bp were pseudo-aligned and quantified to merged GENCODE (M9) and RefSeq gene models using Kallisto 71 . Samples were normalized and pairwise comparisons (cell body to neuropil) were made using Sleuth software (https:// github.com/pachterlab/sleuth) 72 , with a false discovery rate (FDR, q value) of 0.05 (reported in Results). Scatter plots for Kncj6 transcripts corresponding to accession numbers NM_001025584.2 (Girk2a) and NM_010606.2 (Girk2c) were made using log 2 normalized counts. Primary cultures of hippocampal neurons were prepared as described 46 , with minor modifications. Briefly, hippocampi were extracted from neonatal (P1-4) pups and placed into an ice-cold modified Hank's Balanced Salt Solution (HBSS, Sigma-Aldrich) containing 1 mM kynurenic acid, 12 mM MgCl 2 , 3 mM HEPES, and 5.5 mM D-glucose. The tissue was digested for 20 min with papain and DNAse I in the modified HBSS solution at 37 °C with occasional inversion. Tissue was then washed twice with modified HBSS, twice with modified HBSS plus 20% FBS, and twice with growth medium: Neural-basal A medium containing 2% B27 supplement, 0.5 mM Glutamax, and antibiotic/antimycotic (Thermo Fisher Scientific). Hippocampi were mechanically-dissociated in growth medium using trituration in 1 mL pipettes. Neurons were pelleted by centrifugation (150 × g, for 5 min at 4 °C), and plated onto 8-mm glass coverslips pre-coated with 0.0005% poly-L-lysine (Sigma-Aldrich) in 48-well plates. Neurons were maintained in culture at 37 °C/5%CO 2 , and half of the medium was replaced with fresh growth medium every 3-4 d. After 10 d in culture, neurons were infected with control or GIRK expression viruses; electrophysiological analysis or immunofluorescence microscopy was conducted 4-5 d after infection.
Immunoblotting. Intact adult hippocampi and primary hippocampal cultures from wild-type and Girk2 −/− mice were homogenized in RIPA lysis buffer containing Halt phosphatase and protease inhibitor cocktail (Thermo Fisher Scientific) and centrifuged at 4 °C for 20 min at 16,000 × g. The supernatant was taken and protein concentrations were determined by BCA assay. Equal amounts of protein from each sample were mixed with SDS sample buffer and heated at 85 °C for 10 min. Samples were then separated by 12% SDS-PAGE and transferred to nitrocellulose membranes. Membranes were blocked in 5% milk/PBS, incubated overnight at 4 °C with rabbit polyclonal antibodies directed against GIRK1 (1:200, Alomone Labs; Jerusalem, Israel), GIRK2 (1:200, Alomone Labs), GIRK3 (1 μg/mL, Frontier Institute Co.; Ishikari, Japan), or β-actin (0.27 μg/mL, Abcam; Cambridge, MA), and diluted in 5% milk/PBS/0.1% Tween 20. The custom GIRK2c antibody was generated and affinity-purified using a peptide corresponding to the unique 11-amino acid C-terminal domain (DVANLENESKV, Pacific Immunology Corp.; Ramona, CA). Membranes were washed with PBS/0.1% Tween 20 and incubated with donkey anti-rabbit (0.2 μg/mL, LI-COR Biotechnology; Lincoln, NE) or donkey anti-mouse (0.5 μg/mL; LI-COR) secondary antibodies. Blots were developed using the Odyssey infrared imaging system (LI-COR), and the integrated density of each band was measured using NIH ImageJ software (Bethesda, MD).

Electrophysiology in cultured cells. Coverslips containing transfected HEK cells or infected hippocam-
pal neurons were transferred to a chamber containing a low-K + bath solution (in mM): 130 NaCl, 5.4 KCl, 1 CaCl 2 , 1 MgCl 2 , 5.5 D-glucose, 5 HEPES/NaOH (pH 7.4). Neurons were visualized using an Olympus IX-70 microscope. Fire-polished borosilicate patch pipettes (3)(4)(5) were filled with K-gluconate pipette solution (in mM): 140 K-gluconate, 2 MgCl 2 , 1.1 EGTA, 5 HEPES, 2 Na 2 -ATP, 0.3 Na-GTP, and 5 phosphocreatine (pH 7.2). Upon achieving whole-cell access, cells were held in voltage-clamp mode at −70 mV (with no correction for liquid junction potential). GIRK currents were measured in a high-K + bath solution containing (in mM): 120 NaCl, 25 KCl, 1 CaCl 2 , 1 MgCl 2 , 5.5 D-glucose, 5 HEPES/NaOH (pH 7.4); GPCR agonists were dissolved in the high-K + bath solution, which was then applied directly to the cell using a ValveLink 8.2 rapid perfusion system (AutoMate Scientific, Inc.; Berkeley, CA). Whole-cell currents were measured at room temperature with hardware (Axopatch-200B amplifier, Digidata 1322A) and software (pCLAMP v. 8.2) from Molecular Devices, LLC (Sunnyvale, CA). All currents were low-pass filtered at 2 kHz, digitized at 10 kHz, and stored on computer hard disk for subsequent analysis. Only those experiments in which the access resistances were stable and low (<18 MΩ) were included in the analysis. Peak and steady-state currents were measured for each experiment. Current activation and deactivation rates were extracted from 1-term Boltzmann/exponential fits of the recording segments corresponding to the onset and offset, respectively, of drug effect.
Electrophysiology in organotypic slices. Organotypic hippocampal slices were cut out of the culture insert and positioned in a recording chamber superfused (2-2.5 mL/min) with an external solution containing (in mM): 125 NaCl, 2.5 KCl, 25 NaHCO 3 , 1.3 MgCl 2 , 2 CaCl 2 , 0.4 ascorbic acid and 10 glucose (pH 7.4) and bubbled with 95% O 2 /5% CO 2 . Kynurenic acid (2 mM) and picrotoxin (100 μM) were added to the bath solution to block ionotropic glutamate and GABA A receptors, respectively. Bath and chamber temperatures were maintained at 29-30 °C using a TC344C Dual Automatic Temperature Controller (Warner Instruments; Hamden, CT). Whole-cell recordings (V hold = −60 mV) were made from EGFP-positive CA1 pyramidal neurons using borosilicate pipettes (3-5 MΩ) filled with a K-gluconate pipette solution. All measured and command potentials factored in a junction potential (−15 mV) predicted using JPCalc software (Molecular Devices, LLC). Recordings were amplified with an EPC10 HEKA amplifier, filtered at 2 kHz, digitized at 5-10 kHz using Patchmaster 2x73.2 software (HEKA Elektronik; Bellmore, NY) and stored on hard disk. Stimulation was delivered with a cluster-type stimulating electrode (FHC, Bowdoin, ME) and a DS3 stimulator (Digitimer North America, LLC; Ft. Lauderdale, FL). Stimulating electrodes were placed in the stratum radiatum or stratum lacunosum moleculare, the dendritic fields that receive Schaffer Collateral (SC) and Perforant Path (PP) inputs, respectively. Stimuli of increasing intensity (0.2-25 mA) were applied at a 0.05 Hz frequency to establish maximal slow IPSC amplitude. Slow IPSCs evoked by SC and PP were made in each neuron, and the order of stimulation was counterbalanced across experiments. Only experiments with stable (<20% variation) and low series resistances (<30 MΩ) were analyzed.
For quantification of GIRK2 fluorescence in dendritic processes, only infected neurons that were clearly separated from other infected cells were evaluated. MAP2 labeling was used to identify clearly the dendrites; 2-3 relatively straight segments (35-50 microns) of primary, secondary and tertiary dendrites exhibiting even MAP2 labeling from each neuron were selected for analysis. The primary dendrite fragments were located at least 20 microns from the soma. GIRK2 fluorescence intensity measurements within these processes were acquired following background subtraction, and data are presented as arbitrary units (AU). For GIRK2 isoform and PSD-95 co-localization analysis, the soma and surrounding area (20 microns from the edge of the soma) was removed from the image prior to analysis. Quantification was performed using the NIH ImageJ 1.29 plug-in, Puncta Analyzer. MAP2 labeling was used to delineate dendritic processes, and red (GIRK2) and green (PSD-95) channels were thresholded to highlight visible puncta without introduction of background noise. The percentage of co-localized GIRK puncta overlapping with total PSD-95 puncta was calculated for each analyzed neuron. Viral expression of GIRK2a and GIRK2c in dorsal CA1. CaMKIICre(+):Girk2 fl/fl mice (7-8 wk) were placed in a stereotaxic device (David Kopf Instruments; Tujunga, CA) under isoflurane anesthesia. Microinjectors were created by affixing a 33-gauge stainless steel hypodermic tube within a shorter 26-gauge stainless steel hypodermic tube, and were attached to polyethylene-20 tubing affixed to 10 μL Hamilton syringes. Microinjectors were lowered bilaterally through burr holes in the skull to the dorsal CA1 (bregma: −2.20 mm AP, ±1.60 mm ML, −1.40 mm DV) and 400 nL of AAV8-hSyn-DIO-GIRK2a-IRES-EGFP, AAV8-hSyn-DIO-GIRK2c-IRES-EGFP, or AAV8-hSyn-DIO-mCherry control virus (~1.3 × 10 13 viral particles/mL) was injected per side over 4 min using a Pump 11 Elite infusion pump (Harvard Apparatus; Holliston, MA). The syringe was left in place for 10 min following infusion to reduce solution backflow along the infusion track. Behavioral experiments were performed 2 wk after surgery to allow for full recovery and viral expression. The accuracy of viral targeting was assessed in all experiments using fluorescence confocal microscopy. Images were acquired with MetaMorph Advanced Acquisition v. 7.7.7.0 software. GIRK-dependent signaling in CA1 pyramidal neurons expressing GIRK2a or GIRK2c, indicated by EGFP expression, was assessed by measuring baclofen-induced somato-dendritic currents in acutely isolated slices, as described 6 . Fear conditioning. To most effectively isolate the contribution of the hippocampus to context fear learning, we utilized a trace fear conditioning protocol 45 . In brief, mice were acclimated to and handled within the testing room for 2 d prior to testing. For the training session (Day 1), mice were placed in fear conditioning chambers (Med Associates, Inc.; St Albans, VT) with metal bar floors and white wall insert that were cleaned with 75% ethanol. A 30-s white noise (65 dB/3 kHz) was used as the conditioned stimulus (CS, cue), and a 2-s footshock (0.50 mA) was used as the unconditioned stimulus (US). Training consisted of 2 CS-US pairings in which the US was delivered 30 s after the end of the CS; CS presentations were separated by 90 s (7.5 min total conditioning time). On Day 2, the level of contextual fear learning was assessed by returning mice to the training chambers and measuring freezing behavior for 5 min. On Day 3, the chambers were reconfigured (i.e., floor bars were covered with a plastic insert and white wall inserts were replaced with black tent inserts to alter both the color and shape of the chamber) and cleaned with 0.1% acetic acid. Mice were then assessed for freezing behavior in response to 2 × 3-min CS presentations. Freezing was assessed by Video Freeze software (Med Associates, Inc.) and data are presented as percent time spent freezing. Statistical analysis. Data are presented throughout as the mean ± SEM. Statistical analyses were performed using Prism 6 or 7 (GraphPad Software, Inc.; La Jolla, CA). Data were analyzed with Student's two-tailed t-test, 1-way ANOVA, 2-way ANOVA, or 2-way ANOVA with repeated measures, as appropriate. Tukey, Holm-Sidak, Dunnett's, and Bonferroni multiple comparison post hoc tests were used, if justified. Differences were considered significant if P < 0.05. Data availability. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.