Abstract
Genetic studies have shown an association between smoking and variation at the CHRNA5/A3/B4 gene locus encoding the α5, α3, and β4 nicotinic receptor subunits. The α5 receptor has been specifically implicated because smoking-associated haplotypes contain a coding variant in the CHRNA5 gene. The Chrna5/a3/b4 locus is conserved in rodents and the restricted expression of these subunits suggests neural pathways through which the reinforcing and aversive properties of nicotine may be mediated. Here, we show that, in the interpeduncular nucleus (IP), the site of the highest Chrna5 mRNA expression in rodents, electrophysiological responses to nicotinic acetylcholine receptor stimulation are markedly reduced in α5-null mice. IP neurons differ markedly from their upstream ventral medial habenula cholinergic partners, which appear unaltered by loss of α5. To probe the functional role of α5-containing IP neurons, we used BAC recombineering to generate transgenic mice expressing Cre-recombinase from the Chrna5 locus. Reporter expression driven by Chrna5Cre demonstrates that transcription of Chrna5 is regulated independently from the Chrna3/b4 genes transcribed on the opposite strand. Chrna5-expressing IP neurons are GABAergic and project to distant targets in the mesopontine raphe and tegmentum rather than forming local circuits. Optogenetic stimulation of Chrna5-expressing IP neurons failed to elicit physical manifestations of withdrawal. However, after recent prior stimulation or exposure to nicotine, IP stimulation becomes aversive. These results using mice of both sexes support the idea that the risk allele of CHRNA5 may increase the drive to smoke via loss of IP-mediated nicotine aversion.
SIGNIFICANCE STATEMENT Understanding the receptors and neural pathways underlying the reinforcing and aversive effects of nicotine may suggest new treatments for tobacco addiction. Part of the individual variability in smoking is associated with specific forms of the α5 nicotinic receptor subunit gene. Here, we show that deletion of the α5 subunit in mice markedly reduces the cellular response to nicotine and acetylcholine in the interpeduncular nucleus (IP). Stimulation of α5-expressing IP neurons using optogenetics is aversive, but this effect requires priming by recent prior stimulation or exposure to nicotine. These results support the idea that the smoking-associated variant of the α5 gene may increase the drive to smoke via loss of IP-mediated nicotine aversion.
Introduction
Nicotine is a powerfully addictive drug, but under specific conditions of dosage and prior exposure, nicotine consumption can be either reinforcing or aversive. Understanding the variety of nicotinic receptors (nAChRs) that mediate these effects and the neural pathways that express these receptors is of key importance for understanding the mechanisms and genetics of nicotine addiction and for the design of new medications to treat this disorder. The strongest known genetic risk for increased tobacco consumption is associated with certain haplotypes at the Chrna5/a3/b4 gene locus that occur at high frequency in European populations (Berrettini and Doyle, 2012; Lassi et al., 2016). These high-risk haplotypes contain a nonsynonymous Chrna5 polymorphism, Chrna5(D398N), which appears to reduce the function of α4β2-containing (α4β2*) and α3β4* nicotinic receptors into which the risk-associated 398N variant is incorporated (Kuryatov et al., 2011; George et al., 2012; Tammimäki et al., 2012; Morel et al., 2014; Sciaccaluga et al., 2015; Koukouli et al., 2017).
The Chrna5/a3/b4 gene locus is conserved between humans and rodents and the tissue-specific expression of the receptors encoded by these genes suggests neural pathways of importance in nicotine addiction (Picciotto and Mineur, 2014). The channel-forming α3 and β4 receptor subunits are prominently expressed in the ventral medial habenula (MHbV) and the accessory α5 receptor subunit is highly expressed in the interpeduncular nucleus (IP) that receives MHbV input (Marks et al., 1992; Salas et al., 2003; Hsu et al., 2013). The α5 subunit is also detected in the ventral tegmental area (VTA), where it regulates striatal nicotine effects (Exley et al., 2012), and in deep layers of the cerebral cortex, where loss of α5 may impair performance in attentional tests (Bailey et al., 2010). Genetic studies have shown that mice lacking the α5 subunit (α5KO mice) have attenuated somatic signs of nicotine withdrawal (Jackson et al., 2008; Salas et al., 2009) and show self-administration and place preference with high doses of nicotine that are usually aversive (Jackson et al., 2010; Fowler et al., 2011; Morel et al., 2014).
Pharmacological studies also support a role for the MHbV and/or IP as mediating behavioral responses to nicotine, including reinforcement, aversion, and withdrawal (Salas et al., 2009; Zhao-Shea et al., 2013, 2015; Antolin-Fontes et al., 2015; Harrington et al., 2016), but are not by themselves sufficient to assess the role of the α5 subunit, which does not participate in the usual ligand binding site of nAChRs and thus has no specific agonist or antagonist. Here, we show that the IP neurons of α5KO mice have a markedly attenuated physiological response to acetylcholine and nicotine, making the IP a strong candidate for mediating some of the behavioral effects of this nicotinic receptor subunit. We did not observe an attenuation in the nicotine response of MHbV neurons in α5KO mice despite the prevailing idea that α5 expression is relevant to MHb function (Fowler et al., 2011; Frahm et al., 2011; Fowler and Kenny, 2014; Antolin-Fontes et al., 2015).
Prior studies have addressed α5 function in mice using BAC transgenic methods to manipulate the Chrna5/a3/b4 locus (Frahm et al., 2011; Hsu et al., 2013). However, this approach may misexpress the other genes in the α5 cluster, which can independently affect physiology and behavior. Here, we have used a BAC recombineering strategy to produce a new Chrna5Cre transgenic line to manipulate α5-expressing neurons without misexpression artifacts. Reporter expression driven by Chrna5Cre shows that Chrna5 transcription is independently regulated from the Chrna3/b4 genes transcribed on the opposite strand. Chrna5Cre-expressing neurons in the rostral IP (IPR) are GABAergic, but project to specific areas of the mesopontine raphe rather than forming local circuits within the IP. Optogenetic stimulation of α5-expressing IP neurons is aversive, but only after recent prior stimulation or exposure to a single dose of nicotine. Therefore, the activation of neurons expressing α5* receptors in different brain regions may govern the overall rewarding or aversive response to nicotine.
Materials and Methods
Animals.
Mice bearing a deletion of the Chrna5 gene have been described previously (Salas et al., 2003) and have been backcrossed to C57BL/6 (Charles River Laboratories) for >10 generations. Mice null for α5 (α5KO mice) were generated from crosses of heterozygous (α5+/−) mice. To express channelrhodopsin (ChR2-YFP) in the MHbV of α5KO and control mice, ChAT-ChR2 mice (ChAT-ChR2-YFP BAC; Jax #014546; Zhao et al., 2011) expressing channelrhodopsin in ChAT-positive cholinergic neurons were first interbred with α5KO mice to produce compound heterozygotes. Mice with the α5+/−/ChAT-ChR2 genotype were then bred to α5+/− mice to produce mice heterozygous for ChAT-ChR2 and either α5 WT (α5+/+, control) or α5KO for electrophysiological experiments. To confirm that the expression of ChAT-ChR2 did not alter the basal firing properties of MHbV neurons, we replicated the findings from ChR2 mice in tonically firing MHbV neurons of mice arising from the original α5+/− cross that did not express ChR2.
Three independent lines of Chrna5BACCre mice, hereafter referred to as Chrna5Cre mice, were generated by oocyte injection of a modified Chrna5/a3/b4 BAC, as described in the Results. The parent BAC containing Cre recombinase targeted to the Chrna5 locus was generated as part of the GENSAT project (Gong et al., 2003; Gerfen et al., 2013) and obtained from the BACPAC Resources Center, Children's Hospital Oakland Research Institute. The transgene encompassed ∼50 kb of genomic sequence spanning chr9:54,809,614–54,859,149 of the mouse genome (NCBI version 37) and extending from ∼19 kb upstream from the Chrna5 translation start site to ∼3.5 kb beyond the 3′ terminus of the Chrna5 transcribed region adjacent to the 3′ terminus of the Chrna3 transcript, which is transcribed on the opposite strand. The Cre expression cassette was inserted immediately after the initial ATG of the Chrna5 open reading frame at chr9:54,828,904. Recombineering of the parent Gensat BAC-Cre construct and oocyte injections were performed by Cyagen Biosciences. We resequenced all the relevant parts of the BAC transgenic construct to confirm the complete removal of the Chrna3 and Chrnb4 open reading frames and also to confirm that the original GENSAT modification of the Chrna5 locus with the insertion of the Cre cassette were all correct. Because the Chrna5Cre transgene contains no transcribed sequences from the Chrna3 or Chrnb4 loci, there is no direct mechanism by it can affect the expression levels of the endogenous Chrna5/a3/b4 genes. In addition, because the Chrna5Cre transgene is not integrated at the Chrna5/a3/b4 gene locus, cis-acting regulatory sequences contained within the BAC transgene will only affect the expression of the Cre transgene, not the endogenous genes. Three male founder transgenic mice were characterized by breeding with the reporter strain Ai6 (below). After the initial characterization of reporter gene expression, a single Chrna5Cre transgenic line (“line 3”) was used for all subsequent experiments.
Somatostatin (SST)-expressing neurons were identified using the SST-IRES-Cre targeted transgenic strain, Ssttm2.1(cre)Zjh/J (Jackson Laboratories, catalog #013044), which has been described previously (Taniguchi et al., 2011). Three Cre-inducible reporter strains, all targeted to the Gt(ROSA)26Sor locus, were used: Ai6 and Ai14 were used as Cre-dependent fluorescent reporters for the genetic identification of cell phenotypes; Ai6 (B6.Cg-Gt(ROSA)26Sortm6(CAG-ZsGreen1)Hze/J, Jax 007906) expressing a ZsGreen reporter and Ai14 (B6.Cg-Gt(ROSA)26Sortm14(CAG-tdTomato)Hze/J, Jax 007914) expressing a cytoplasmic tdTomato reporter, have been described previously (Madisen et al., 2010). Mice for Cre-driven optogenetic studies were generated using the mouse line Ai32 (B6.Cg-Gt(ROSA)26Sortm32(CAG-COP4*H134R/EYFP)Hze/J, Jax 024109), which conditionally expresses the Channelrhodopsin-2 variant ChR2(H134R)-EYFP (Madisen et al., 2012). Mice bearing a Chrna5Cre allele and an Ai32 allele were used all in vivo optogenetic experiments and are referred to as IPChR2 mice. All transgenic lines used in behavioral experiments were maintained on a C57BL/6 genetic background (Charles River Laboratories). Chrna5Cre mice were backcrossed at least three generations and Ai32 mice at least six generations with commercially supplied C57BL/6 mice before use in behavioral experiments. Experiments using α5KO or α5KO/ChAT-ChR2 mice and matched controls were performed under the guidelines of the Canadian Council on Animal Care and all experimental procedures were approved by the Faculty of Medicine Animal Care Committee at the University of Toronto. Experimental procedures in Chrna5Cre mice, inducible reporter lines, and controls were approved by the Institutional Animal Care and Use Committee of Seattle Children's Research Institute.
Electrophysiology methods.
Acute brain slice electrophysiology experiments in littermate WT (α5+/+) and homozygous knock-out (α5KO) mice were performed to assess the significance of α5 subunit expression on neurophysiological activity and cholinergic modulation of neurons in the IPR and the MHbV. Experiments included recordings from 12 WT mice and 13 α5KO mice. Coronal slices of brain tissue including the IPR or the MHbV were obtained by investigators blinded to genotype from male and female adult mice (mean ± SEM, postnatal day 121 ± 9).
Mice were anesthetized with chloral hydrate (400 mg/kg) and decapitated and their brains quickly removed and chilled in 4°C sucrose ACSF containing the following (in mm): 254 sucrose, 10 d-glucose, 24 NaHCO3, 2 CaCl2, 2 MgSO4, 3 KCl, and 1.25 NaH2PO4, pH 7.4. Coronal brain slices (∼bregma −3.5 for the IPR, ∼bregma 1.0–2.0 for the MHbV), 400 μm thick, were obtained using a Dosaka linear slicer (Sci Media) and recovered for ∼2 h in regular ACSF containing the following (in mm): 128 NaCl, 10 d-glucose, 26 NaHCO3, 2 CaCl2, 2 MgSO44, 3 KCl, and 1.25 NaH2PO4, pH 7.4. All solutions were oxygenated with 95% O2/5% CO2.
For whole-cell electrophysiological recordings, brain slices containing the IPR or MHbV were placed in a perfusion chamber on the stage of a BX50W1 microscope (Olympus) and perfused with oxygenated ACSF at a rate of 3–4 ml/min. Recording electrodes (2.5–4.5 MΩ) were filled with solution containing the following (in mm): 120 potassium gluconate, 5 KCl, 2 MgCl2, 4 K2-ATP, 0.4 Na2-GTP, 10 Na2-phosphocreatine, and 10 HEPES buffer, with pH adjusted to 7.3 using KOH. Data were acquired at 20 kHz, low-pass filtered at 3 kHz using pClamp software (Molecular Devices), and corrected for the liquid junction potential (14 mV). ChAT-ChR2-positive MHbV cells, which send cholinergic projections to GABAergic cells in the IP (Ren et al., 2011; Hsu et al., 2013), were identified visually using YFP fluorescence and patched in either hemisphere. Confirmation of neuronal identity was obtained both by observing the characteristic tonic firing at rest and confirming the presence of ChR2 using 473 nm blue light stimulation (5 ms pulse delivered while holding the neuron at −75 mV in voltage clamp). No differences were found in these light-evoked current responses between genotypes (WT, 604 ± 40 pA, N = 14, a5KO, 528 ± 81 pA, N = 13, t(25) = 0.86, p = 0.4).
To examine the electrophysiological impact of Chrna5 subunit expression in the IPR, we assessed the effects of bath application of ACh (acetylcholine chloride, Sigma-Aldrich; 100 μm to 1 mm, 15 s) and nicotine (nicotine hydrogen tartrate salt, Sigma-Aldrich; 1 μm, 60 s) on the amplitude of current responses using voltage-clamp recordings with neurons held at −75 mV. Rapid bath application (3–4 ml/min) was selected to survey receptors in a location-unbiased manner that would mimic the arrival of nicotine in the brain, which happens on a time course of seconds (Rose et al., 2010). However, we note that this approach may underestimate a rapidly desensitizing response. The role of synaptic transmission in mediating cholinergic response amplitude was assessed using the AMPA receptor antagonist CNQX (Tocris Bioscience; 20 μm), the NMDA receptor antagonist APV (D-2-amino-5-phosphonovaleric acid, Sigma-Aldrich; 50 μm), and the GABAA receptor antagonist bicuculline (bicuculline methiodide, Tocris Bioscience; 10 μm). We also assessed the role of the α4β2* nicotinic ACh receptor in cholinergic current responses using dihydro-β-erythroidine hydrobromide (DHβE, Sigma-Aldrich; 10 μm). To assess the role of Chrna5 subunit expression on nicotinic responses in the MHbV and IPR, bath application of nicotine was delivered during current-clamp recordings from resting membrane potential (MHbV) or from −70 mV (MHbV and IP).
All electrophysiological data are expressed as the mean ±SEM. Parametric or nonparametric statistical tests were used according to the Shapiro–Wilk test of normality. Intrinsic electrophysiological properties of WT and α5KO neurons, as well as current responses to light stimulation in MHbV ChAT-ChR2 neurons, were compared using independent t tests. Concentration response data were assessed using two-way ANOVA. Differences in responses to cholinergic and nicotinic stimulation between WT and α5KO neurons were assessed using Mann–Whitney nonparametric t tests. Differences in cholinergic responses before and after application of synaptic blockers or DHβE were assessed using Wilcoxon matched-pairs signed-rank nonparametric tests.
Methods used for measuring light-evoked action potentials in acute brain slices from adult mice have been described previously (Hsu et al., 2013). Neurons expressing ChR2-EYFP were identified visually in brain slices of IPChR2 mice aged 2–4 months by their endogenous fluorescence. Recordings were performed on a total of 8 adult IPChR2 mice. Glass pipettes (3–7 MΩ) filled with intracellular pipette solution containing the following (in mm): 140 K-gluconate acid, 1 CaCl2, 2 MgSO4, 10 EGTA, 4 Na2-ATP, 0.3 Na-GTP, and 10 HEPES, pH 7.2, were used for all recordings. The recording bath solution consisted of aCSF containing the following (in mm): 118 NaCl, 3 KCl, 25 NaHCO3, 1 NaH2PO4, 1 MgCl2, 1.5 CaCl2, and 30 glucose, pH 7.4, with added GABA blockers gabazine (10 μm) and 2-OH saclofen (50 μm) in all recordings. Light pulses (10 ms) were delivered from a 473 nm blue laser coupled to an optical fiber positioned in the recording bath 2 mm form the recording pipette tip, producing 1–2 mW/mm2 irradiance at the recording pipette tip. When possible, we measured the baseline firing of the cell and the light response of the cell to 5, 10, 20, and 50 Hz light pulses using cell-attached recording, followed by recording the response of the cell in current-clamp mode.
Tract tracing, immunofluorescence, and FISH.
For Cre-dependent labeling of cell bodies and axons, we used AAV1-CAG.Flex.tdTomato.WPRE.bGH (Addgene Plasmid #51503) or AAV1-CAG.Flex.EGFP.WPRE.bGH (Addgene Plasmid #51502, Oh et al., 2014), Viral stocks were prepared at the University of Pennsylvania Gene Therapy Program Vector Core (http://www.med.upenn.edu/gtp/vectorcore/). Enhanced labeling of the terminal fields of neuronal projections was performed by Cre-dependent viral expression of a synaptophysin-EGFP fusion protein (sypGFP). The plasmid pCAG.Flex.SypEGFP.WPRE (AAV-sypGFP) was constructed by replacing the EGFP moiety of pCAG-FLEX-EGFP-WPRE with the sypGFP construct from phSyn1(S)-FLEX-tdTomato-T2A-SypEGFP-WPRE (Addgene Plasmid #51509) by Julie Harris, Karla Hirokawa, and Hong Gu of the Allen Institute for Brain Science. In most experiments, the tdTomato axonal tracer and the sypGFP synaptic tracer viruses were coinjected.
The detailed methods used here for anterograde tract tracing with iontophoretic injection of AAV have recently been published in conjunction with the Allen Mouse Brain Connectivity Atlas (Harris et al., 2012; Oh et al., 2014). The injection coordinates used were as follows: Case IP-a1: AP 3.5 (from bregma), ML 0.1, DV 4.6; Case MnR-a1: AP 4.3, ML 0.2, DV 4.2; Case MnR-a2: AP 4.1, ML 0.3, DV 4.5; based on a standard atlas (Paxinos and Franklin, 2001). Animals were fixed by transcardial perfusion with 4% paraformaldehyde at 14–21 d after injection, brains were removed and equilibrated in graded sucrose solutions, frozen at −80°C in optimal cutting temperature solution, and cryosectioned at 25 μm for fluorescence/immunofluorescence imaging (Quina et al., 2017). Antisera used and references for the validation of their specificity include: Gpr151 C-terminus, rabbit polyclonal, Sigma-Aldrich SAB4500418, RRID:AB_10743815 (Broms et al., 2015); green fluorescent protein (GFP), goat polyclonal, Abcam ab6673, RRID:AB_305643; choline acetyltransferase (ChAT), goat polyclonal, Millipore AB144P, RRID:AB_2079751; substance P, rat monoclonal IgG2a, EMD Millipore MAB356, RRID:AB_94639; Somatostatin (S-14), rabbit polyclonal, Peninsula Laboratories T-4103, RRID:AB_518614; tryptophan hydroxylase 2 (Tph2), rabbit polyclonal, Millipore ABN60, rabbit polyclonal, RRID:AB_10806898; and cFos, rabbit polyclonal, Santa Cruz Biotechnology Sc-52, RRID:AB_2106783.
Double-label FISH (DFISH) was performed as part of the Allen Institute for Brain Science transgenic characterization project (Oh et al., 2014). Complete datasets used for expression Chrna5 + TH mRNA, used here to show Chrna5 expression alone (experiment 304864931), and Chrna5 + Gad1 mRNA (experiment 127833597) are available online at: http://connectivity.brain-map.org/transgenic.
Behavioral tests.
Three cohorts of mice were used for behavioral studies. All behavioral cohorts consisted of Chrna5Cre, Gt(ROSA)26SorAi32/+ experimental mice (IPChR2), and Gt(ROSA)26SorAi32/+ controls lacking the Cre driver from the same or simultaneously generated litters. Cohort 1 consisted of 6 male and 6 female experimental mice, and 4 male and 4 female control mice. Cohort 2 consisted of 3 male and 4 female experimental mice, and 3 male and 7 female control mice. Cohort 3 consisted of 6 male and 6 female experimental mice, and 8 male and 7 female control mice. All mice were between 2 and 3 months old at the time of surgery and between 3 and 6 months old at the time of the behavioral experiments. Cohort 1 underwent open-field locomotion, stimulus-primed real-time place preference (RTPP), and intracranial self-stimulation (ICSS). Cohort 2 underwent the light/dark box test, RTPP, RTPP controls with no light stimulus on day 2 of the protocol and priming of the response by light stimulus in the home cage, and testing for somatic signs of nicotine withdrawal. Cohort 3 underwent RTPP, RTPP control for priming of the response by light stimulus in the home cage, and priming of the RTPP response by acute administration of nicotine.
Mice were implanted with a 100 μm (cohorts 1 and 3) or 200 μm (cohort 2) fiber-optic cannula targeting the IP using a 20° angle and the target stereotaxic coordinates AP −3.3 to −3.5 (mm caudal to bregma), ML 0.00 (midline), and DV 4.60 (Paxinos and Franklin, 2001). After the behavioral experiments, mice were perfused with 4% paraformaldehyde and the brains of all subject animals were examined for ChR2-EYFP fluorescence in the IP to confirm the genotype and the consistent expression of the transgene and for cannula placement. Although some cannulas differed slightly from the intended coordinates, only one animal from cohort 1 and one animal from cohort 3 were excluded from the behavioral analysis due to a misplaced cannula; in these cases, the cannula had obstructed the cerebral aqueduct, causing hydrocephalus. Behavioral data are presented graphically as mean ± SEM.
Open-field locomotion.
Open-field locomotion was assessed in IPChR2 and control mice over three 10 min intervals: baseline, stimulation, and poststimulation. The stimulation period consisted of 30 s periods of 20 Hz light pulses alternating with 30 s periods of light off. The open-field test enclosure consisted of a 27.3 cm × 27.3 cm arena (ENV-510, Med Associates) housed in a sound attenuating chamber (ENV-018MD, Med Associates). For determination of time spent in the center and periphery of the enclosure, a center area of 17.5 × 17.5 cm was defined in the video image. The total distance traveled was analyzed using the video-tracking software EthoVision XT 10.0 (Noldus Information Technology).
Light/dark box.
Studies were conducted in a three-chamber place-preference box (ENV-3013, Med Associates). Incident light on the dark side of the chamber was blocked by covering it with black felt. Light intensity on the light side of the chamber was ∼750 lux and the black side was ∼2 lux. The 473 nm laser delivered 10 ms pulses at 20 Hz and 8 mW power to all subjects for the duration of the 15 min trial. Time spent in the dark chamber was measured for all subjects and compared between genotypes using unpaired t tests.
Video recording for the open-field test, RTPP test, and light/dark box test was performed with an ICD-49 B/W video camera with heliopan E35.5 RG 850 IR lens filter (Ikegami Tsushinki) under IR Illumination. In all experiments recorded with Noldus Ethovision, an Arduino UNO microprocessor was used to control the frequency and pulse duration of a 473 nm modulated laser diode (MLD, model 0473-06-01-0100-100, Cobolt) to deliver light pulses via a fiber-optic cannula.
Somatic signs associated with nicotine withdrawal.
To reproduce the methods used in prior work, somatic signs typically associated with nicotine withdrawal were assessed using a previously described optogenetic stimulation protocol (Zhao-Shea et al., 2013). Behaviors were assessed in a standard housing cage and animals were acclimated to the test cage for 45 min/d on the 2 d before the test. Experiments were performed in nicotine-naive mice. Optogenetic stimulation at 8 mW was delivered intermittently to the IP for 10 min, followed by 10 min of recording without stimulation. The stimulation consisted of repeated 50 Hz, 0.2 s trains of 5 ms pulses (10 pulses) delivered at 1 s intervals for 5 s, followed by lights-off periods of 5 s. Mice were video recorded throughout the 20 min trial and the entire period was scored manually by two independent raters who were blinded to genotype. Six somatic signs associated with withdrawal were scored: scratching, head or body shaking (“like a wet dog”), head nodding (at least two closely consecutive nods), backing (both hind limbs must step backward), rearing (on side or in middle of cage, both forepaws must leave bedding), and chewing. The average scores of two independent scorers blinded to genotype were used for each behavior and subject and the effect of genotype was assessed using unpaired t tests, with subsequent correction for multiple hypothesis testing.
ICSS.
To test ICSS reinforcement by IP stimulation, mice were placed into an operant chamber (ENV-307W, Med Associates) with two response wheels that record one event for every 90° of rotation (ENV- 113AMW, Med Associates). One response wheel was randomly assigned as the active wheel in the initial four 45 min trials and reversed for the second four trials. Two turning events recorded on the active wheel within a 2 s interval resulted in the delivery of 2 s of 473 nm light stimulation consisting of 10 ms pulses delivered at 20 Hz and 8mW power.
RTPP.
RTPP studies were conducted in a two-chamber place-preference box (ENV-010MC, Med Associates) in which mice received light stimulation on one side and could move freely between compartments. The 15 min sessions were initiated by placing the mouse in the center of the apparatus. For each subject, the active side was randomly assigned in the initial trial and reversed in each subsequent trial. Recording and laser stimulation were controlled with EthoVision XT using center point tracking. The active chamber was paired with a 473 nm laser that was activated, delivering 10 ms light pulses at a 20 Hz and 8 mW of total power, until the mouse crossed over to the inactive chamber, whereas no simulation was delivered in the inactive chamber. Active and inactive chamber occupancy and total distance traveled were then calculated for each 5 min interval of the 15 min session and for the entire session. In the standard RTPP experiments, two trials were performed on consecutive days, differing only by the reversal of the stimulated side. Two variants of the RTPP protocol were performed as controls. In one variant, one side of the chamber was stimulated on day 1 and no stimulus was present on day 2. In another variant, intermittent stimulation was provided for 15 min in a standard housing cage on day 1 and day 2 was conducted in the two-compartment chamber. To administer light stimulation in a standard housing cage, the stimulation intervals were pseudorandomized based on the typical distribution of stimulation periods delivered to control mice in the two compartment box, with a total stimulation time of 7.5 min (50% of trial). An identical 15 min intermittent stimulation protocol was used to induce the expression of cFos in the IP of IPChR2 mice (n = 3 IPChR2 mice, 3 cannulated controls). In the cFos induction trials, mice were acclimated to the test cage for 5 d before testing. Mice were perfused with paraformaldehyde 2 h after the end of the stimulation period.
Acute administration of nicotine and nicotine inhibition of home cage locomotion.
To assess the effects of nicotine on home cage locomotion and RTPP, mice were administered 0.3 mg/kg nicotine ditartrate subcutaneously and returned to a standard housing enclosure. Nicotine dosage was calculated based on the mass of the free base and administered as a 0.1 mg/ml solution (Bailey et al., 2010). Distance traveled was then recorded for 1–16 min after injection using Noldus Ethovision center point tracking. Two hours after injection, mice were placed in the two-compartment place-preference enclosure for a 15 min RTPP trial with optogenetic stimulation of the IP, as described above. The test was repeated using the same methods 10 d later. Due to the short half-life of nicotine in mouse plasma (<10 min; Siu and Tyndale, 2007), the timing of the locomotor test represents the acute effect of nicotine and the RTPP trial tests the effect of recent exposure to nicotine. Reversal trials in the nicotine-primed RTPP study were conducted on day 10 of the experiment to avoid any effect of activity-primed aversion.
Results
Electrophysiology of the IP response to ACh and nicotine
To address directly the significance of the expression of the Chrna5 subunit in the rostral interpeduncular nucleus, we performed acute slice recordings from sibling WT (α5+/+) and α5KO mice. The intrinsic properties of IPR neurons were recorded and compared by genotype (Table 1). There were no significant differences on measures of input resistance, resting membrane potential, capacitance, spike amplitude, or spike threshold.
We first investigated how the absence of Chrna5 subunit expression affected current response to bath application of ACh (1 mm). Data for both the WT and α5KO groups failed the Shapiro–Wilk test for normality (WT, n = 27, W = 0.75, p < 0.0001; α5KO, n = 27, W = 0.72, p < 0.0001, Fig. 1-1 relating to Fig. 1), necessitating the use of nonparametric analyses. The amplitude of current responses was markedly reduced in a5KO IP neurons (WT, 227 ± 50 pA, n = 27; α5KO, 37 ± 7 pA, n = 27, Mann–Whitney U = 70.0, p < 0.0001; Fig. 1A,B) and this genotype difference was observed across a range of concentrations (100 μm to 1 mm; WT, n = 8; α5KO n = 6; F1, 34 = 16.6, p = 0.0003; Fig. 1C,D).
Figure 1-1
IP neurons are predominantly GABAergic and receive synaptic inputs from the MHbV that are both glutamatergic and cholinergic (Ren et al., 2011; Hsu et al., 2013). The MHbV terminals are rich in nAChRs that could act presynaptically to facilitate glutamate release. In the presence of synaptic blockers (APV, CNQX, and bicuculline), there was a significant reduction in the current amplitude of responses to ACh in within-cell experiments (WT baseline, 252 ± 75, synaptic blockers, 147 ± 52 pA, n = 9; α5KO baseline, 41 ± 11, synaptic blockers, 19 ± 5 pA, n = 6; Wilcoxon W = −110.0, p = 0.006, data not shown). However, in the presence of synaptic blockers, there remains a strong and significant genotype difference in cholinergic currents (WT: 147 ± 52 pA, n = 9; α5KO: 18 ± 5 pA, n = 6; Mann–Whitney U = 4.0, p = 0.005; Fig. 1E,F). This demonstrates that glutamatergic synaptic transmission contributes to the ACh response of IP neurons, but their intrinsic response to ACh is strongly affected by the presence of the α5 subunit.
To determine the principal subtype of channel-forming nAChRs mediating the cholinergic response in the IPR and thus potentially associated with α5 subunits in this nucleus, we blocked ACh-induced currents with the β2* nAChR antagonist DHβE (10 μm, 10 min). Cholinergic currents were substantially reduced by DHβE in both genotypes (WT: baseline, 360 ± 143; DHβE, 75 ± 25 pA; n = 5; α5KO: baseline, 44 ± 14; DHβE, 8 ± 3 pA; n = 4; Wilcoxon W = −45.0 p = 0.004; Fig. 1G,H).
Because the current responses to ACh application were largely mediated by nAChRs, we next applied nicotine (1 μm) in the bath to assess differences in response between genotypes. Inward currents elicited by nicotine were significantly greater in WT neurons (WT, 122 ± 35 pA, n = 14; α5KO, 16 ± 4 pA, n = 12; Mann–Whitney U = 31.0, p = 0.005; Fig. 1I,J). These results indicate that the Chrna5 gene is critical for the normal response of IPR neurons to nicotinic stimulation.
Contrasting nicotine-induced excitability of WT and α5KO neurons in the IP and MHbV
To investigate the question of where α5 exerts its maximal effects within the habenula–IP pathway, we examined the ability of nicotine to excite IP neurons and MHbV cholinergic neurons in adult WT and α5KO mice (Fig. 2). Nicotine (1 μm) was bath applied to WT and α5KO neurons in the IPR in current-clamp conditions (Fig. 2A). The peak nicotine-elicited spike frequency was greater in neurons from WT (13.8 ± 5 Hz, n = 15) than α5KO mice (1.6 ± 1 Hz, n = 12, Mann–Whitney U = 46.0, p = 0.02; Fig. 2B–D; data failed the Shapiro-Wilk test for normality, necessitating the use of non-parametric analyses, see Fig. 2-1). Voltage-clamp recordings confirmed the significant difference in nicotinic currents observed above (WT, 144.9 ± 50 pA, n = 7, α5KO, 14.5 ± 9 pA, n = 12, Mann–Whitney U = 3.0, p < 0.001). Membrane and action potential properties did not differ significantly between genotypes except for greater input resistance in the α5KO neurons in this sample (WT, 348 ± 90 MΩ, n = 8, α5KO, 812 ± 164 MΩ, Mann–Whitney U = 18, p = 0.03), which would suggest that they were IP cells of a type that would otherwise have the most prominent responses to nicotine (Shih et al., 2014).
Figure 2-1
Because cholinergic neurons in the MHbV are the vital part of the habenula–IP pathway, we also tested the response of these neurons to nicotine. We recorded from the tonically firing, YFP-positive ChAT-ChR2 neurons in the MHbV of either hemisphere (Fig. 2E,F) and confirmed the presence of ChR2 (see Materials and Methods) in both WT and a5KO transgenic mice. Intrinsic properties of MHbV neurons were measured and compared by genotype and no significant differences between genotypes were found on any measures (Table 2). In the first set of experiments, nicotine (1 μm) was applied while neurons were firing tonically at rest in current clamp, increasing their firing frequency (Fig. 2F). In contrast to the IP recordings, however, peak nicotine-elicited firing frequency did not differ significantly between MHbV neurons from WT (9.3 ± 1.4 Hz, n = 12) and α5KO (7.9 ± 1.4 Hz, n = 13) mice (Mann–Whitney U = 58.0, p = 0.29; Fig. 2F–H). To avoid a potential ceiling effect, we next assessed nicotine-induced firing in current clamp from an initial membrane potential of −70 mV, similar to the starting membrane potential in the IP experiments. Nicotine depolarized MHbV cholinergic neurons, as illustrated in Figure 2I, yet no difference was found in nicotine-elicited peak firing frequency between genotypes (WT, 4.3 ± 0.8 Hz, n = 12, α5KO, 3.4 ± 1.1 Hz, n = 12, Mann–Whitney U = 49.0, p = 0.19; Fig. 2I–K). Additional voltage-clamp experiments with neurons held at −75 mV did not detect any significant difference in the response to nicotine between genotypes (WT, 20.6 ± 5 pA, n = 12, α5KO, 15.0 ± 3.8 pA, n = 11, Mann–Whitney U = 52.5, p = 0.4). In a final series of experiments, we recorded from nonlabeled, tonically firing neurons in the MHbV of WT and a5KO mice from the original a5+/− cross. In these MHbV neurons recorded at rest, nicotine-elicited peak firing frequencies also did not differ significantly between genotype (WT, 9.2 ± 1.6 Hz, n = 12, α5KO, 11.9 ± 1.2 Hz, n = 5, Mann–Whitney U = 20.5, p = 0.3), nor were there significant differences in MHbV membrane properties between genotypes (data not shown).
Generation and characterization of Chrna5Cre transgenic mice
The cell-specific expression of the α5/α3/β4 nAChRs determines the possible repertoire of pentameric receptors in a given neuron and thus the function of these receptors in specific brain regions. The Chrna3 and Chrnb4 genes are transcribed on the same DNA strand and, in rodents, they are highly expressed in the MHbV. Chrna5, which is transcribed on the opposite strand, is highly expressed in the IP and moderately expressed in the VTA and is also found in other CNS areas not necessarily correlated with transcription of Chrna3/b4, including layer 6 of the cerebral cortex, CA1 of the hippocampus, and the median/paramedian raphe (Wada et al., 1990; Marks et al., 1992; Winzer-Serhan and Leslie, 1997; Salas et al., 2003; Hsu et al., 2013). Mouse genomic BACs modified to express GFP or Cre-recombinase have been used to mark and manipulate gene expression in neurons expressing one or more of the genes in the Chrna5/a3/b4 cluster (Frahm et al., 2011; Hsu et al., 2013). Because the three genes at this locus are tightly linked, BAC transgenic models that target one of the three genes at this locus have the disadvantage that they may misexpress the other genes (Ables et al., 2017). To address these limitations and to provide a better model for studies of α5-expressing IP neurons, we constructed a new BAC-derived Chrna5Cre transgene that completely eliminates the Chrna3 and Chrnb4 transcribed sequences (Fig. 3) and thus cannot overexpress these subunits. The parent BAC clone RP24–318N16 was modified as part of the GENSAT project by the insertion of a Cre recombinase expression cassette immediately after the Chrna5 translation start site, blocking the native Chrna5 open reading frame (Gerfen et al., 2013). This BAC clone contains most of the intergenic region between Chrna5 and the upstream locus Psma4, the entire Chrna5 expressed region, and the entire Chrna3 and Chrnb4 loci, including several regions of high interspecies conservation (Fig. 3C,D). In humans, sequence variants in this genomic region appear to affect CHRNA5 expression levels (Smith et al., 2011). We also examined database information for CTCF binding sites within the locus determined by chromatin immunoprecipitation as part of the ENCODE project (Fig. 3E). Such sites are associated with chromosomal insulators that isolate enhancer function (Ali et al., 2016). We note that a cluster of CTCF-binding sites are found in the short intergenic region between Chrna5 and Chrna3, suggesting that these genes can be independently regulated and the 3′ end of the transgenic construct was chosen to coincide with these sites. We then derived a ∼50 kb subclone of the parent BAC that included all of the available Chrna5 upstream flanking sequence plus the inserted Cre-transcribed region and the Chrna5 exonic, intronic, and 3′ untranslated regions, ending between the closely opposed Chrna5 and Chrna3 loci (Fig. 3F,G, and Materials and Methods).
The BAC-derived transgenic construct Chrna5Cre was injected into mouse oocytes and three founder lines were derived. Founders were crossed with a Cre-inducible reporter strain, Ai6, which gives ZsGreen expression in the presence of Cre-recombinase (Fig. 4). Overall, the expression pattern of the Cre-induced reporter expression in Chrna5Cre/Ai6 mice resembled the previously described expression of a Chrna5GFP transgene derived from the same parent BAC (Hsu et al., 2013). In the IP, expression driven by all three founder lines closely resembled expression of the endogenous mRNA (Fig. 4A–F), including expression in the IPR, the dorsolateral IP subnucleus (IPDL), and the intermediate and central IP subnuclei (IPI/IPC) and lack of expression in lateral IP subnucleus (IPL). Chrna5-expressing neurons in IPR and IPDL consist predominantly of GABAergic projection neurons (Hsu et al., 2013), whereas IPI contains a preponderance of Chrna5-expressing GABAergic neurons that are local interneurons (Ables et al., 2017). In the VTA, expression of the transgene was very sparse (lines 1 and 2) or absent (line 3; Fig. 4G–I) and clearly did not reproduce the native pattern of Chrna5 mRNA expression in dopaminergic VTA neurons (Fig. 4J). A Chrna5GFP BAC transgene derived from the same genomic region is also not expressed in the VTA (Hsu et al., 2013), so it appears that enhancers controlling VTA expression reside outside the region included in the parent BAC.
Two Chrna5Cre transgenic lines showed some expression in the MHb and LHb, but in different sets of neurons (Fig. 4K,L). The third transgenic line showed negligible habenula expression (Fig. 4M). In the adult habenula, Chrna5 mRNA is weakly expressed and is not detectable by nonisotopic ISH or FISH (Fig. 4N and data not shown), but has been detected by 35S-ISH (Franceschini et al., 2002; Salas et al., 2003) and qPCR (Hsu et al., 2013). Because of the weak ISH signal and the lack of available antibodies for immunohistochemistry, the expression of Chrna5 has not been assigned to any specific habenula subnucleus. Marker expression in the two lines that showed habenula signal did not correlate with mRNA expression for the MHbV cholinergic markers vesicular acetylcholine transporter (VAChT) and Chrna3 (Fig. 4O,P). Because of this inconsistency, the habenula expression in these lines can be considered ectopic. Chrna5Cre line 3 was used for subsequent functional studies of the IP because of its lack of potentially confounding expression in the habenula and VTA.
Chrna5 mRNA is also expressed in layer 6 of the cerebral cortex and cortical neurons in α5KO mice show altered responses to acetylcholine (Bailey et al., 2010). Crosses of the three Chrna5Cre founder mice with the Ai6 reporter line revealed transgene expression in cortical layer 6 in one line (Fig. 5A–C), but with this was associated with significant ectopic expression in a radial pattern consistent with clonal expression (Fig. 5B). To determine the neurotransmitter phenotype of the Chrna5Cre-expressing neurons, Chrna5Cre-mice were generated expressing a Cre-dependent tdTomato reporter, Ai14 (Madisen et al., 2010), and a Cre-independent marker of GABAergic neurons, Gad67-GFP (Tamamaki et al., 2003). This is of interest because α5* receptors have been hypothesized to mediate nicotinic responses in VIP-expressing GABAergic interneurons in the mouse prelimbic cortex (Koukouli et al., 2017). Chrna5Cre-expression was not detected in GABAergic interneurons (Fig. 5D), implying that it is restricted to glutamatergic pyramidal neurons. This could be because the regulatory sequences contained in the Chrna5Cre transgene do not drive correct expression in GABAergic neurons or because the Chrna5 mRNA is not significantly expressed in these cells. For this reason, we examined the expression of mRNA for Gad1, a GABAergic marker, and Chrna5 in the prelimbic region of WT mice by DFISH (Fig. 5E). Coexpression of these markers was very rarely observed either within pyramidal cell layers 5 and 6, where most of the Chrna5-expressing neurons reside, or in more superficial layers, where Chrna5 neurons are scattered. However, occasional Gad1-expressing neurons appeared to express low levels of Chrna5 mRNA. This suggests that Chrna5Cre-driven reporter expression correctly reports that the α5 subunit is expressed in layer 6 cortical pyramidal neurons, but only rarely and/or at very low levels in the interneurons of any layer.
Figure 5-1
Remarkably, although reporter crosses of Chrna5Cre line 3 founder mice showed little cortical expression, these mice bred one generation into a C57BL/6 background (N1) showed extensive, appropriate expression in the cortex (Fig. 5-1). Analysis of the transmission frequency of the Chrna5Cre allele in line 3 was consistent with Mendelian transmission of a single insertion site (∼50% transmission), ruling out the presence of multiple insertion sites driving different expression patterns. We infer that the variable expression in the cortex is due to epigenetic effects on cortical expression of the Cre transgene. No such variability of Cre-dependent transgene expression was observed in the IP, where the induction of ROSA26-targeted ZsGreen (Ai6) and ChR2-YFP (Ai32) transgenes was 100% penetrant in all generations. The Chrna5Cre transgenic lines derived here contain only the Chrna5 gene and upstream regulatory sequences and do not contain coding or regulatory information from the Chrna3/b4 region of the combined Chrna5/a3/b4 locus. Therefore, these results demonstrate the autonomous transcriptional regulation of the Chrna5 gene in specific brain regions (IP, cortex) independent of enhancers that might reside in the Chrna3/b4 region of the gene locus. In other cell types, such as the VTA, information from outside the transgenic construct used here appears to be required for correct expression. These findings may help to provide a framework for the interpretation of the noncoding polymorphisms found in smoking-related haplotypes (Barrie et al., 2017).
Specific projections of Chrna5-expressing neurons in the IP and the mesopontine raphe
The IP subnuclei have distinct afferents and efferents. Chrna5Cre is expressed in IPR and IPC/IPI, which receive input from MHbV, and excluded from IPL, which receives input from the dorsal MHb (Quina et al., 2017). To map the specific efferents of the Chrna5Cre neurons, we used a Cre-dependent AAV tract-tracing strategy. Two viruses were coinjected, one expressing synaptophysin-EGFP (sypGFP), which prominently labels the presynaptic area of efferent fibers, and one expressing tdTomato, which labels cell bodies and fibers but is not specifically transported to synaptic areas. In each injected area, the majority of labeled neurons coexpressed both markers. In case IP-a1, AAV was injected into the IPR of a Chrna5Cre mouse (Fig. 6A,B), with a small number of labeled neurons also noted in the adjacent median/paramedian raphe (MnR/PMnR; Fig. 6C,D). Prior work has suggested that SST-expressing GABAergic neurons in IPR make presynaptic connections onto afferents from the MHbV, which may modulate habenula input (Zhao-Shea et al., 2013). These SST neurons in IPR have been shown to be a subpopulation of the Chrna5-expressing neurons (Hsu et al., 2013). However, a synaptic connection by these neurons within the IP has not been directly demonstrated. In case IP-a1, labeled cell bodies in the IP were restricted to IPR, but labeled synaptic puncta were observed in IPC and in symmetrical bands running between IPI and IPL (Fig. 6A). The synaptic labeling was sometimes associated with ChAT labeling in MHbV afferents (Fig. 6E–G), although it is not possible to determine with certainty whether these synapses are associated with these afferents or with local IP neurons. Paradoxically, no tdTomato-labeled fibers of passage were observed between IPR and the ventral synaptic puncta, suggesting that these synapses may be derived from the small population of labeled PMnR neurons. To test this hypothesis, in case MnR-a1, we injected sypGFP+tdTomato expressing viruses into the PMnR of a Chrna5Cre mouse caudal to the IP (Fig. 6H). Although fewer neurons were labeled compared with case IP-a1, dense synaptic labeling was observed in the ventral IP (Fig. 6I–K). To further characterize the raphe–IP projection, in case MnR-a2, we injected a GFP tract-tracing virus into the PMnR of an SSTCre mouse (Fig. 6L). This also produced strongly labeled afferents to IPC and the border between IPI and IPL (Fig. 6M). DFISH for Chrna5 and SST mRNA showed that the SST-expressing neurons in IPR are a subset of the Chrna5-expressing neurons found there (Fig. 6N). However, in the MnR/PMnR, expression of these transcripts does not overlap (Fig. 6O). We conclude that IPR neurons do not project within the IP as previously hypothesized, but rather, the IP receives input from distinct populations of Chrna5- and SST-expressing neurons in the raphe.
Prior work has shown that IP efferents project predominantly to the raphe and mesopontine tegmentum, with distinct projections for the IP subnuclei (Quina et al., 2017). Viral tract tracing of the efferents of Chrna5Cre-expressing neurons in IPR shows synaptic labeling in a subset of these areas (Fig. 7A–I). In the median raphe, synaptic puncta from IPR projections generally appear in PMnR, surrounding but usually not overlapping 5HT neurons in MnR proper (Fig. 7A,B). Labeling in the rostral part of dorsal raphe was sparse (DR; Fig. 7C,D). Intense synaptic labeling was observed in the rhabdoid nucleus (Rbd; Fig. 7C,E), which can be identified by immunofluorescence for the orphan receptor Gpr151. The expression of this orphan receptor in Rbd can probably be attributed to inputs from Gpr151-expressing excitatory neurons in the LHb (Broms et al., 2015), showing an interesting confluence of inhibitory IP and excitatory LHb inputs to this nucleus. Projections from Chrna5Cre IP neurons are sparse in the laterodorsal tegmental nucleus (LDTg) and are not closely associated with the cholinergic neurons there (Fig. 7F,G). Projections are absent from the central dorsal tegmental nucleus (DTg; Fig. 7H,I), although there is strong synaptic labeling in the surrounding the central pontine gray (CGPn) and near the midline, in the nucleus incertus (NI; Fig. 7I). Viral tract tracing of Chrna5Cre neurons in MnR (case MnR-a1; Fig. 7J–L) shows that these neurons also form synaptic connections to the mesopontine raphe and dorsal tegmentum in a pattern that overlaps that of the IPR efferents but spares Rbd. The caudal projections of Chrna5Cre IP neurons in IPR appear similar to IP projections described in the rat (Lima et al., 2017).
Activation of IP Chrna5 neurons is aversive
To test the behavioral function of the IP Chrna5 neurons, we generated an optogenetic model using Cre-inducible expression of ChR2. Chrna5Cre mice were interbred with mice homozygous for a floxed-stop ChR2-EYFP cassette targeted to the Gt(ROSA)26Sor gene locus (Madisen et al., 2012). ChR2 expression was observed in the IP, but not in adjacent structures such as the VTA, which might confound behavioral experiments (Fig. 8A). IPChR2 neurons and fibers were found in IPR, IPC, and IPI, areas that receive glutamatergic/cholinergic afferents from the MHbV (ChAT label; Fig. 8B), but not in IPL, which receives glutamatergic/peptidergic inputs from the MHbD (SP label; Fig. 8C). Therefore, IPChR2 mice allow specific functional testing of the MHbV-IP pathway most often implicated in the effects of nicotine.
To assess the light responsiveness of IPChR2 neurons, we recorded cellular responses in IPR neurons in acute brain slice preparations (see Materials and Methods) using whole-cell current-clamp (n = 3) and cell-attached (n = 5) recordings. Of eight recorded cells, only two were spontaneously firing at <1 Hz; the others were silent in the absence of light stimulus. Light responsiveness was tested using 2 s trains of 10 ms light pulses at stimulus frequencies of 5, 10, 20, and 50 Hz. Eight of 8 neurons tested were able to fire action potentials at 1:1 correspondence with the 5 Hz and 10 Hz stimuli and 6/8 neurons were able to follow a 20 Hz stimulus (Fig. 8C1–C3). However, in the presence of a 50 Hz stimulus, some neurons produced an action potential in response to the first light pulse and then entered a depolarization block (5/8 cells; Fig. 8C4), whereas others briefly fired action potentials at a 1:1 correspondence, but were unable to do so for 2 s (2/8 cells; Fig. 8D). Therefore, for all behavioral experiments, a 20 Hz stimulus frequency was used except when assessing somatic nicotine withdrawal signs, for which a specific published stimulus protocol was followed (see Materials and Methods).
To assess the effect of IP stimulation on behavior, IPChR2 mice and controls were implanted with fiber-optic cannulas just dorsal to or slightly within IPR (see Materials and Methods; Fig. 9A, Fig. 9-1). Light stimulation from cannulas in this dorsal position resulted in the induction of cFos expression in neurons throughout the IP (Fig. 9-2). Because visible light and drug administration can in themselves affect behavior, all behavioral experiments were structured as a comparison between genotypes, in which IPChR2 mice and controls received identical light stimulation and/or pharmacological treatment. To assess for any effect of IP activity on acute locomotion, IPChR2 mice were tested in an open-field arena (Fig. 9B–E). During the stimulus phase of the trial, a small-magnitude increase in distance traveled during the light on intervals was noted, but there was no effect of genotype [repeated-measures (RM)-ANOVA, genotype × light, effect of genotype: F(1,18) = 0.24, p = 0.63, effect of light: F(1,18) = 13.89, p = 0.0015; Fig. 9B,C], nor were significant differences noted between genotypes in the three phases of the experiment (multiple-hypothesis-corrected t test; baseline t(54) = 1.10, p = 0.83, stimulation t(54) = 0.49, p > 0.99, poststimulation t(54) = 0.24, p > 0.99; Fig. 9D). In the open-field test, avoidance of the center of the enclosure has been used as a model of anxiety (Bailey and Crawley, 2009). A small-magnitude decrease in center occupancy was observed during the lights-on intervals, but there was no effect of genotype (RM-ANOVA: genotype × light, effect of genotype: F(1,18) = 1.26, p = 0.28, effect of light: F(1,18) = 9.468, p = 0.0065; Fig. 9E). Similarly, no difference between genotypes was detected in another assay of anxiety, the light/dark box (unpaired t test between genotypes: t(28) = 0.97, p = 0.34; Fig. 9F).
Figure 9-1
Figure 9-2
Prior work has reported that optogenetic stimulation of GABAergic neurons in the IP can elicit somatic signs associated with nicotine withdrawal, which include scratching, shaking, nodding, backing, rearing, and chewing (Zhao-Shea et al., 2013). To exactly match the stimulation protocol used in the prior work, we used intermittent optogenetic stimulation at 50 Hz in a standard housing cage (see Materials and Methods). IPChR2 mice showed no difference from controls in any of the somatic signs tested (t test between genotypes for each behavior, p > 0.99 for all seven behaviors; Fig. 9G). In the prior study of somatic signs, the IP of Gad2Cre mice was injected with a conditional AAV vector expressing ChR2 in all GABAergic neurons in the injected area. The amount of injected virus (1 μl) was approximately twice the entire volume of the IP and is likely to have diffused into the surrounding tegmental nuclei, all of which contain GABAergic neurons, although the extent of viral spread was not discussed. Therefore, the prior result may represent off-target effects.
To determine whether the stimulation of α5-expressing IP neurons is reinforcing, aversive, or neutral, we used RTPP in a two-compartment shuttle box. To control for effects of preexisting preference or habituation, RTPP was performed with a 2 d reversal structure. On day 1, mice were assigned randomly to receive the light stimulation on the right or left side of the box, referred to here as the “A” side. On the subsequent day of testing, light stimulation was delivered on the opposite “B” side. No significant preference or avoidance was observed on day 1 of a 2 d RTPP experiment; however, on day 2, a significant avoidance of the stimulated side was observed in IPChR2 mice (RM-ANOVA: genotype × day, effect of genotype: F(1,18) = 14.71, p = 0.0012, effect of day: F(1,18) = 0.83, p = 0.38; Fig. 10A). Division of the test period into 5 min intervals for IPChR2 mice on the two test days revealed an incipient trend toward avoidance of the stimulated side on day 1, which became a significant difference between the start and end of the trial on day 2 (RM-ANOVA, day × epoch, effect of day: F(2,33) = 5.71, p = 0.0074, effect of day: F(1,33) = 8.67, p = 0.0059; Fig. 10B). In principle, the enhancement of avoidance seen on day 2 could either be a learning effect, in which a specific side preference is retained from day 1, or a “priming” effect, in which recent prior stimulation sensitizes the real-time avoidant response. However, given the reversal of the stimulated side between days, a choice to avoid the stimulated side on day 2 would result from a learned preference for the chamber stimulated on day 1. This seems unlikely based on the trend to develop avoidance of the stimulated side over the course of the test on both days. To test for the retention of a learned preference, we performed a 2 day experiment in which one side was stimulated on day 1 and the animals were returned to the shuttle box on day 2 without a light stimulus. No retained side preference was exhibited in the absence of the stimulus (RM-ANOVA, genotype × day, effect of genotype: F(1,24) = 2.13, p = 0.16, effect of day: F(1,24) = 2.00, p = 0.17; Fig. 10C). In a second cohort of mice, we tested the reversibility of the avoidance effect. A trend toward avoidance in IPChR2 mice on day 2 disappeared by day 10, demonstrating the transient nature of the priming effect (RM-ANOVA, day × genotype, effect of genotype: F(1,15) = 3.774, p = 0.07, effect of day: F(2,30) = 0.62, p = 0.55; Fig. 10D). Although the trend toward an effect of genotype did not reach statistical significance in this experiment, analysis of day 1 and day 2 data for the combined cohorts showed a clear aversive effect in IPChR2 mice on day 2 (RM-ANOVA, day × genotype, effect of genotype: F(1,35) = 24.86, p < 0.0001, effect of day: F(1,35) = 0.46, p = 0.50; Fig. 10E,F). Finally, to eliminate any positive or negative learned effect of side in the primed aversion effect, we provided the initial stimulation in a standard housing cage. Prior stimulation in a separate cage was sufficient to prime the aversive effect (t test between genotypes: t(41) = 2.47, p = 0.0176; Fig. 10G).
In prior studies, we used an optogenetic ICSS protocol to demonstrate that activation of the dorsal subnucleus of the medial habenula is reinforcing (Hsu et al., 2014, 2016). In this paradigm, mice are given access to paired response wheels, one of which triggers a laser stimulus to the area of interest. After 4 d of training, the active wheel is switched for four reversal trials. Increased wheel turns, preference for the active wheel after 1–2 d of training, and reversal of the preference with change of the active wheel are indicative of an ICSS effect. As expected given the aversive nature of IP stimulation in the shuttle box test, IPChR2 mice did not show a preference for the stimulus-associated wheel in this paradigm on day 4 or day 8, when a maximum response would be seen (Fig. 10H). The observed infrequent, sporadic wheel-turning behavior resulted in a very low frequency of light stimulation, which did not demonstrate reinforcement and could not be shown to be aversive because of the low rate of stimulation at baseline.
α5-expressing IP neurons receive input from cholinergic/glutamatergic neurons of the MHbV. The role of cholinergic transmission in this system has not been fully defined, but acetylcholine released at the MHbV-IP interface has the potential to both directly activate IP nicotinic receptors and to act presynaptically to increase glutamate release (Ren et al., 2011; Hsu et al., 2013; Frahm et al., 2015). This suggests that, like direct stimulation of the IP, exposure to nicotine could prime the aversive effect of IP activation in RTPP. Nicotine-primed RTPP was performed with a 2 day reversal design. The initial trial and reversal trial were separated by 10 day to avoid the effect of light stimulation priming, which is lost after this amount of time (Fig. 10D). On day 1, IPChR2 and control mice were given a subcutaneous injection of nicotine (0.3 mg/kg, see Materials and Methods) and placed a standard housing cage. The chosen dose of nicotine is by itself reinforcing or neutral, not aversive, in place preference assays (Neugebauer et al., 2011; Grieder et al., 2017). Acute locomotion, measured from 1 to 16 min after nicotine administration on day 1, was reduced compared with that observed in a prior trial in which nicotine was not delivered (RM-ANOVA, drug × genotype, effect of genotype: F(1,24) = 4.71, p = 0.04, effect of drug: F(1,24) = 46.01, p < 0.0001; Fig. 11A). RTPP was measured 2 h after nicotine administration. The rapid metabolism of nicotine in mice, with a half-life of <10 min (Matta et al., 2007), predicts that the residual blood level from this single dose of nicotine was low at the time of RTPP and is consistent with prior results (Bailey et al., 2010). The acute inhibition of locomotion after nicotine administration had resolved by the time RTPP was conducted. IPChR2 mice exposed subacutely to nicotine in this way showed a marked aversion to IP stimulation (RM-ANOVA, day × genotype, effect of day: F(1,21) = 0.70, p = 0.42; effect of genotype: F(1,21) = 7.856, p = 0.0107; Fig. 11B–D) with a strong reversal of side preference between day 1 and day 10 (Fig. 11C). In experiments with activity-primed RTPP, the aversive response increased over the 15 min trial of IP stimulation (RM-ANOVA, day × 5 min epoch, effect of day: F(2,54) = 5.184, p = 0.0087; effect of epoch: F(1,54) = 5.963, p = 0.0179; Fig. 10B), whereas in nicotine-primed RTPP, strong aversion was observed in IPChR2 mice in the first 5 min period of the trial on day 1 (Fig. 11E). We also noted that, after nicotine administration, control mice showed slight aversion to IP stimulation, spending ∼40% fractional time in the laser-paired compartment. This is consistent with the small effect of light seen in control animals in open-field locomotion and was not sufficient to mask the very strong genotype-dependent aversion observed in IPChR2 mice.
Discussion
Smokers self-titrate their nicotine dosage and the propensity to smoke and the intensity of smoking will result from the balance of dose-dependent reinforcement and aversion. Smokers with the α5 398N risk allele, a “mutant” allele in that it is not observed in other species, alter their smoking behavior to deliver more nicotine (Macqueen et al., 2014) and the importance of aversion in this equation is demonstrated by the lower nicotine aversion of these individuals (Jensen et al., 2015). Rodent models strongly implicate the mesolimbic and habenulopeduncular systems in determining the balance of nicotine reinforcement/aversion and resulting drive to consume the drug. Here, we have shown that inward currents elicited by ACh and nicotine are markedly decreased in α5KO mice in the IP, but not in the MHbV, and that specific stimulation of the IP neurons that express the α5 receptor can mediate aversion.
Role of the α5 subunit in cellular responses to ACh and nicotine
The α5 subunit is an accessory subunit that is not required for pentameric nAChR formation or ligand binding (Zoli et al., 2015; but see Jin et al., 2014 and Jain et al., 2016). In principle, the observed decrease in agonist-elicited currents in α5KO mice could be due to fewer available receptors, possibly through an effect on stability or trafficking, or because of reduced ion flux per receptor due to altered conductance or channel kinetics. Radioligand binding is unchanged in the IP of α5KO mice, as assayed by autoradiography (Baddick and Marks, 2011) and in synaptosomal membrane preparations (Beiranvand et al., 2014), suggesting that receptor levels are unchanged, although in these assays, the abundant nAChRs on MHbV afferents in the IP could obscure a loss of intrinsic binding. However, radioligand binding is also independent of α5 in other CNS regions that express this subunit, such as the VTA and cortical layer 6, and in synaptosome preparations from thalamus, hindbrain, and spinal cord (Jackson et al., 2010). One study using an α4-GFP chimeric subunit as a probe reported reduced α4* receptor levels in the VTA of α5KO mice (Chatterjee et al., 2013), but this is hard to reconcile with the ligand-binding studies. Therefore, it seems likely that most of the decrease in agonist-induced current observed in α5KO IP neurons results from reduced channel-specific ion flux rather than fewer available receptors. Direct studies of the effect of the α5 subunit on channel properties have not been performed in CNS neurons. However, incorporation of α5 in a4b2* channels in model cell systems results in increased Ca2+ permeability (Kuryatov et al., 1997; Gerzanich et al., 1998; Tapia et al., 2007) and inclusion of α5 in α3β4* channels in autonomic ganglia results in longer channel open times (Ciuraszkiewicz et al., 2013).
Together with other recent reports, these data add to a developing picture of where in the brain the α5 subunit is likely to mediate behavioral responses to nicotine and where it is not. Loss of α5 in mice results in marked decreases in agonist-elected currents in the rostral IP (shown here), prelimbic cortical layer 6 pyramidal neurons (Bailey et al., 2010, 2012; Verhoog et al., 2016), and the VTA (Chatterjee et al., 2013; Morel et al., 2014). We also show a sharp decline in nicotine-elicited action potential firing frequency in the IPR with the loss of α5. In contrast, we find that mice lacking the α5 subunit show no change in nicotine-evoked action potential firing frequency in adult MHbV neurons, consistent with earlier work on nicotine-elicited currents in juvenile mice (Dao et al., 2014). These results correlate strongly with the relative expression of the α5 mRNA in these areas (Hsu et al., 2013), pointing to the IP and not the MHb as the site where α5 subunit expression is relevant in the habenulopeduncular system. These results do not support a model in which the behavioral phenotype of α5-null mice (Fowler et al., 2011) or mice overexpressing the β4 receptor (Frahm et al., 2011) can be “rescued” or “balanced” by viral expression of the α5 subunit in the habenula.
Role of the IP in behavioral responses to nicotine
The IP is the interface between the MHb and its output to tegmental systems that may affect motivation and mood states, so it is important to understand the functional circuitry of the IP. Zhao-Shea et al. (2013) have proposed that IPR neurons expressing SST synapse on MHbV afferents, gating the glutamatergic/cholinergic input to the IP. Application of SST in acute slice preparations has been shown to inhibit EPSCs evoked in IP neurons by the stimulation of MHb afferents (Ables et al., 2017), but this does not identify the source of SST. The SST neurons in IPR are a subset of the α5-expressing neurons, which are in turn a subset of the GABAergic neurons that are the predominant cell type in the IP (Hsu et al., 2013; Quina et al., 2017). Here, using Cre-mediated tract tracing, we find no evidence that SST-expressing neurons form a local circuit within the IP to the areas of strongest MHbV input. Instead, the ventral IP subnuclei IPI and IPC receive synaptic input from distinct populations of neurons expressing α5 and SST in the raphe, which are also GABAergic (Hsu et al., 2013). These inputs terminate near MHbV afferents, but whether they synapse upon them or upon the local IP neurons is unknown.
Injection of the nicotinic antagonist mecamylamine into the IP precipitates somatic symptoms of nicotine withdrawal in nicotine-dependent mice (Salas et al., 2009). ACh is an excitatory neurotransmitter at the MHb–IP junction that can either act directly to excite postsynaptic IP neurons or act on MHbV autoreceptors to facilitate glutamate release (Ren et al., 2011; Hsu et al., 2013; Frahm et al., 2015). Therefore, mecamylamine administration probably acts to acutely decrease IP excitability in association with these withdrawal signs. Also consistent with this, inhibition of pacemaker currents in the MHb, which provides excitatory input to the IP, produces withdrawal symptoms (Görlich et al., 2013). Paradoxically, it has been reported that optogenetic stimulation of the IP results in somatic signs of nicotine withdrawal (Zhao-Shea et al., 2013). Here, however, no such signs were detected in IPChR2 mice using an identical stimulation protocol and, because of the nonspecific expression system used, the prior result is potentially explained by off-target effects.
IPR projects prominently to the mesopontine raphe (Quina et al., 2017), suggesting a role in modulating 5HT function. The reduced response of IP neurons to nicotine in α5KO mice correlates well with recent studies showing that, in WT mice, nicotine inhibits 5HT neuron firing, whereas in α5KO mice, nicotine increases 5HT neuron firing and hippocampal 5HT release (Besson et al., 2016). Loss of α5 function in the IP may thus result in a failure to appropriately brake 5HT release in response to nicotine. However, the α5* GABAergic neurons distributed throughout the MnR could also contribute to these effects (Hsu et al., 2013).
The role of the MHbV and IP in mediating the propensity to consume nicotine is clearly complex. Here, we have shown that direct, specific stimulation of the IP neurons that receive cholinergic MHbV input can produce aversion, which is greatly enhanced after nicotine administration. Using a self-stimulation paradigm, we could not find any level of IP stimulation that was reinforcing. Nicotinic activation of the MHb may contribute to reinforcement because blockade of α3β4* receptors there inhibits accumbal DA release in response to nicotine (McCallum et al., 2012; Eggan and McCallum, 2016), but this is unlikely to be affected by loss of the α5 subunit because it has no direct effect on the MHb response to nicotine. Although the MHbV provides excitatory input to the IP, activity of the IP appears to mediate nicotine aversion. Local infusion of nicotine into the IP is not reinforcing (Ikemoto et al., 2006) and silencing of the IP with local anesthetics (Fowler et al., 2011) leads to increased self-administration of nicotine. The mesolimbic DA system also clearly plays a role in the dose-dependent effects of nicotine on reinforcement and aversion (Exley et al., 2012). These effects are not necessarily paradoxical if viewed along a continuum of nicotine dose–response. In the absence of the α5 subunit, the VTA-mediated dose–response curve for reinforcement may be shifted to higher doses of nicotine. IP-mediated aversion may also be attenuated and fail to act as a brake on the consumption of normally aversive higher doses. The loss of α5 activity appears to affect the overall balance of reinforcement and aversion in a way that favors higher consumption in mouse self-administration (Fowler et al., 2011; Morel et al., 2014) and place preference tests (Jackson et al., 2010) and in humans (Macqueen et al., 2014; Jensen et al., 2015). It is this balance of reinforcement and aversion that is likely to be perturbed by the risk alleles of the α5 receptor subunit and contribute to the individual risk for tobacco use.
Footnotes
This work was supported by the National Institutes of Health (Grant R01-DA035838 to E.E.T. and Grant R01-MN093667 to E.E.T.), the Canadian Institutes of Health Research (Grant MOP-89825 to E.K.L.), a Canada Research Chair in Developmental Cortical Physiology (E.K.L.), and the Alzheimer's Society Research Program of Canada (Postdoctoral Fellowship to D.W.S.). The parent BAC construct used to generate Chrna5Cre transgenic mice was generated by the Gene Expression Nervous System Atlas (GENSAT) Project supported by National Institute of Neurological Disorders and Stroke Contracts N01NS02331 and HHSN271200723701C to The Rockefeller University. AAV-DIO-sypEGFP was generously provided by Julie Harris, Karla Hirokawa, and Hong Gu of the Allen Institute for Brain Science. We thank Mariella De Biasi for the original gift of α5KO mice, Julie Harris and Hongkui Zeng of the Allen Institute for Brain Science for access to unpublished data, and Lely Quina and Janice McNabb for expert technical assistance.
- Correspondence should be addressed to Eric E. Turner, Seattle Children's Research Institute, University of Washington, 1900 Ninth Avenue, Mail Stop C9S-10, Seattle, WA 98101. eric.turner{at}seattlechildrens.org