Abstract
The habenular complex in the epithalamus consists of distinct regions with diverse neuronal populations. Past studies have suggested a role for the habenula in voluntary exercise motivation and reinforcement of intracranial self-stimulation but have not assigned these effects to specific habenula subnuclei. Here, we have developed a genetic model in which neurons of the dorsal medial habenula (dMHb) are developmentally eliminated, via tissue-specific deletion of the transcription factor Pou4f1 (Brn3a). Mice with dMHb lesions perform poorly in motivation-based locomotor behaviors, such as voluntary wheel running and the accelerating rotarod, but show only minor abnormalities in gait and balance and exhibit normal levels of basal locomotion. These mice also show deficits in sucrose preference, but not in the forced swim test, two measures of depression-related phenotypes in rodents. We have also used Cre recombinase-mediated expression of channelrhodopsin-2 and halorhodopsin to activate dMHb neurons or silence their output in freely moving mice, respectively. Optical activation of the dMHb in vivo supports intracranial self-stimulation, showing that dMHb activity is intrinsically reinforcing, whereas optical silencing of dMHb outputs is aversive. Together, our findings demonstrate that the dMHb is involved in exercise motivation and the regulation of hedonic state, and is part of an intrinsic reinforcement circuit.
Introduction
Changes in physical activity and the inability to enjoy rewarding or pleasurable experiences (anhedonia) are two of the hallmarks of major depression. Conversely, structured exercise is one of the best-documented nonpharmacological interventions for depressed mood (Trivedi et al., 2011). Consistent with this, rodents have a strong drive to engage in voluntary wheel running activity (WRA) (Novak et al., 2012), and access to a running wheel relieves behavioral states induced by acute and chronic stress that resemble human depression (Solberg et al., 1999; Greenwood et al., 2003; Duman et al., 2008). However, the neural pathways underlying exercise motivation and the maintenance of mood states are not well understood.
The medial habenula (MHb) is a paired midline structure located on the dorsal surface of the thalamus (epithalamus). Past studies using electrical stimulation, or using lesions of the habenula or its output tract, have assigned a variety of behavioral functions to this region, including primary reinforcement (Klemm, 2004; Lecourtier and Kelly, 2007). However, most of these studies have not distinguished the MHb from the adjacent lateral habenula, which is known to mediate negative motivational values, including reward prediction error and punishment (Jhou et al., 2009; Hong et al., 2011).
MHb neurons are glutamatergic, but two principal subnuclei can be identified by the use of acetylcholine as a cotransmitter in the ventral MHb (vMHb) and the expression of the neuropeptide substance P (SP) in the dorsal MHb (dMHb) (Contestabile et al., 1987; Lecourtier and Kelly, 2007; Ren et al., 2011). In addition, although the dMHb and vMHb project to the same region of the ventral midbrain, the interpeduncular nucleus (IP), they innervate distinct structures there (Quina et al., 2009; Hsu et al., 2013). Recent work has suggested a role for the cholinergic neurons of the vMHb in behavioral responses to nicotine (Grady et al., 2009; Salas et al., 2009; Görlich et al., 2013; Leslie et al., 2013). In contrast, no specific function has been assigned to the dMHb.
In the present study, we have investigated the behavioral functions of the dMHb using three mouse models: a developmental genetic lesion mediated by loss of function of the transcription factor Pou4f1 (Brn3a) that specifically ablates the dMHb; the Cre-mediated conditional expression of channelrhodopsin (ChR2H134R) in the dMHb, allowing the specific activation of this region with light; and conditional expression of halorhodopsin (eNpHR3.0) in the dMHb, allowing the silencing of the output of this nucleus to the IP. Ablation of the dMHb reduces voluntary WRA with relatively little effect on simple tests of gait and balance and no effect on basal locomotion, and also reduces sucrose preference, a measure of hedonic state. Consistent with a role for the dMHb in mediating positive mood states, optogenetic stimulation of the dMHb is strongly reinforcing in an intracranial self-stimulation (ICSS) paradigm. Conversely, optogenetic silencing of dMHb output to the IP produces acute place aversion. Collectively, our data suggest that the dMHb is part of a pathway required to motivate voluntary exercise and regulate hedonic state, and that dMHb activity can mediate primary reinforcement.
Materials and Methods
Generation of transgenic mice.
For Cre-mediated disruption of the Pou4f1 gene, mice were generated containing loxP sites flanking Pou4f1 exon 2 (“floxed” allele), which contains most of the Brn3a protein-coding sequence, including the DNA-binding domain. ES cell targeting was performed in 129 ES cells using a congenic targeting construct introducing the loxP sites and a TK-NEO selection cassette (see Fig. 1A). Correct targeting of the locus was demonstrated in the primary screen of recombinant ES cell clones by PCR across the 3′ targeting arm and confirmed by Southern blotting using an introduced Hind3 restriction site and a probe derived from the 5′ targeting arm. ES cells were injected into C57BL/6 blastocysts to generate chimeric founders, which were screened for germline transmission of the floxed allele. Transgenic services were provided by the University of California San Diego Transgenic Core Facility. Germ-line excision of the TK-NEO cassette in the chosen founder line was performed by crossing with mice bearing FLPeR allele, expressed from the Rosa26 locus (Farley et al., 2000), to generate the allele Pou4f1flox.
The Syt6-BAC Cre transgenic line (Syt6Cre) from the Gensat project, STOCK Tg(Syt6-Cre)KI148Gsat/Mmcd, was used for the manipulation of loxP-mediated transgenic events in the dMHb. The Syt6Cre transgenic line was obtained as cryopreserved sperm from the Mutant Mouse Regional Resource Center of the University of California, Davis, and recovered by in vitro fertilization of C57BL/6 mouse oocytes by the University of Washington Transgenic Resources Program. Both the Pou4f1flox and Syt6Cre mice were backcrossed >6 generations to C57BL/6NCrl mice (Charles River Laboratories) before use in behavioral experiments with the exception of FST, where the Syt6Cre mice were backcrossed for 4 generations to C57BL/6 mice. To generate mice with a tissue-specific null mutation of Pou4f1 in the dMHb, Syt6Cre mice were first interbred with a mouse strain containing a Pou4f1 targeted tauLacZ transgene, consisting of a fusion of bovine tau with Escherichia coli lacZ (Quina et al., 2005), to generate mice with the genotype Pou4f1tlacZ/+, Syt6Cre/+. Pou4f1tlacZ is also a null allele of Pou4f1. These mice were then backcrossed to Syt6Cre/+ animals to generate homozygous Cre driver mice with the genotype Pou4f1tlacZ/+, Syt6Cre/Cre. These mice were then interbred with Pou4f1flox/flox mice to yield the genotypes Pou4f1tlacZ/flox, Syt6Cre/+ (dMHbCKO) and Pou4f1flox/+, Syt6Cre/+ (control) mice in equal ratios. The initial characterization of the dMHb genetic lesion model was examined in control mice, with the genotype of Pou4f1flox/tlacZ, and dMHbCKO mice, with the genotype of Pou4f1flox/tlacZ, Syt6Cre/+ to allow immunofluorescent detection of the lacZ gene product, β-galactosidase, in both dMHbCKO mice and control mice.
The mouse strains used for Cre-inducible expression of Channelrhodopsin2 (Ai32, expressing ChR2(H134R)-YFP) and Halorhodopsin (Ai39, expressing eNpHR3.0-EYFP) have been described (Madisen et al., 2012). Both strains express their respective fusion proteins from a modified floxed-stop Gt(Rosa)26Sor locus. Ai32 mice were interbred with Syt6Cre mice to yield the double-heterozygous Ai32, Syt6Cre (dMHbChR2) optogenetic mice used for activation of the dMHb. Ai39 mice were interbred with Syt6Cre mice to yield double-heterozygous Ai39, Syt6Cre (dMHbNpHR) optogenetic mice used for silencing of dMHb efferents. Control mice for the optogenetic studies were derived from the same crosses but lacked either the opsin allele or the cre-driver allele. All mouse lines were maintained on a C57BL/6NCrl genetic background. Behavioral assays were performed on male mice 2–8 months of age.
Genotyping.
Genotyping was performed by real-time PCR using the following oligonucleotide pairs: Cre allele, forward, TGGGCCAGCTAAACATGCTT; reverse, AACAGCATTGCTGTCACTTGGT; Rosa26 locus common primer, forward, TCGTGATCTGCAACTCCAGTCTT; Ai32 allele specific primer, reverse, TGGGCTATGAACTAATGACCCC; Rosa26 WT locus specific primer, reverse, CGCCCACACACCAGGTTAG; Ai39 allele specific primers, forward, CGTGCTGAACGACCCTTTG; reverse, CCTGCAAGTGCGATGTTGAT; Pou4f1 flox gene, forward, TACGAAGTTATTAGGTCCCTCG; reverse, GCAGAGTGGAGCTGAGGGAA; Pou4f1 WT gene, forward, CGCAGCGTGAGAAAATGA AC; reverse, CGATGGAAGTCCGCTTGC; tauLacZ transgene, forward, GAGGAACCAGGCTCTGAAACC; reverse, AAGGGTGCTGTCGCATCTTC.
Immunohistochemistry and in situ hybridization (ISH).
Adult mice were injected intraperitoneally with pentobarbital (120 mg/kg) and then perfused transcardially with PBS followed by 4% PFA in PBS. Brains were fixed in 4% PFA in PBS at room temperature for 4 h while shaking, equilibrated in 15% then 30% sucrose, embedded in OCT medium, and cryo-sectioned at 30 μm for both immunohistochemistry and ISH. Primary antibodies used included the following: guinea pig anti-Brn3a (Quina et al., 2005), goat anti-ChAT (Millipore), and rabbit anti-β-galactosidase (β-Gal; 5Prime-3Prime). Secondary antibodies conjugated to Alexa series fluorophores were obtained from Invitrogen. ISH methods and sources of ISH probes have been previously described (Eng et al., 2004; Quina et al., 2009).
The full extent of the dMHb lesion in dMHbCKO mice compared with controls was examined by immunofluorescent labeling of alternate sections for Brn3a plus ChAT and β-gal (Pou4f1tlacZ) plus ChAT in coronal sections at 120 μm intervals from bregma −0.98 mm to −2.06 mm, which encompass the extent of the habenula as identified in a standard atlas (Paxinos and Franklin, 2001). The Pou4f1tlacZ allele, in which a lacZ transgene replaces the Brn3a coding sequence, allowed the use of β-gal expression to identify any remaining neurons in the dMHb of dMHbCKO mice that may not express Brn3a protein due to Pou4f1 excision. The area of expression of Brn3a+ (control) and β-gal+ (dMHbCKO) in the right and left dMHb was measured using ImageJ. ChAT immunostaining was used to demarcate the boundary between dMHb and vMHb. The extent of the genetic dMHb lesion was also examined in the large behavioral cohort of mice used in the WRA, open field, PhenoTyper cages, CatWalk, rotarod, and sucrose preference tests. Two habenular sections (30 μm sections at bregma −1.6 and −1.9 mm) from the dMHbCKO, and control mice were immunostained for Brn3a, β-gal, and ChAT and analyzed for the area of the dMHb. The dMHbCKO mice showed significant reductions in both the rostral and caudal sections compared with the control mice. Little variation was observed in the small remaining area of the dMHb, and there was no correlation between the extent of the lesion and the behavioral outcomes.
Behavioral analysis of dMHbCKO mice.
Behavioral analysis was performed on dMHbCKO mice (Pou4f1tlacZ/flox, Syt6Cre/+) and their control littermates (Pou4f1flox/+, Syt6Cre/+). All tests were conducted using male mice between 2 and 8 months of age. Three large cohorts of control and dMHbCKO mice were used for preliminary and reported experiments. In Figures 3, 4, and 5, the data reported were obtained from a single cohort of mice, except for the forced swim test (FST), which was performed on a different cohort. For the main cohort, the order of the tests performed was as follows: Sucrose preference test, accelerating rotarod, home cage locomotion, CatWalk gait analysis, open field locomotion, and WRA. Wheel running was performed last such that the differential running of the control and dMHbCKO mice would not impact the other tests. Mice were housed on a 12:12 light and dark (LD) cycle with lights off at 7 PM.
WRA.
Mice were individually housed in cages equipped with a running wheel and maintained in a ventilated, light-tight chamber under a 12:12 LD cycle with 200-lux intensity. WRAs were monitored using Clocklab software (Actimetrics) and analyzed using the El Temps software (Han et al., 2012). Data were analyzed in 10 min or 60 min bins for double-plotted WRA actograms or WRA waveform, respectively.
Accelerating rotarod.
The Rotamex-5 (Columbus Instruments) was used to test rotarod performance. The apparatus is fitted with a 7.0-cm-diameter dowel with a fall height of 44.5 cm from the dowel center. Before testing, all mice were trained on a fixed-speed protocol at 4 rpm until they could stay on the rod for 30 s (4–7 trials). On the same day as the training sessions, mice underwent four 5 min accelerating rotarod trials. In each trial, the rotarod accelerated from 4 to 40 rpm at the rate of 1 rpm every 8 s, then remained at 40 rpm until the end of the trial. The principal outcome was the time (latency) until the mouse fell from the rod. Mice were given at least 15 min of rest in between each trial. A second session of four trials was performed the following day to determine the reproducibility of the test. No training period was conducted on the second day of testing. To calculate the average latency to fall for each mouse, the lowest of four values was discarded.
Open field locomotion.
The open field enclosure consists of a 27.3 cm × 27.3 cm arena (ENV-510, Med Associates) housed in a sound attenuating chamber (ENV-018MD, Med Associates). Mice were recorded by an Ikegami ICD-49 B/W video camera with heliopan E35.5 RG 850 IR lens filter (Ikegami Tsushinki) lit by an IR Illuminator (YY-IR30, YY Trade) for 30 min. 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 8.5 (Noldus Information Technology).
PhenoTyper home cages.
The PhenoTyper home cage (model 3000, floor dimensions: 30 × 30 cm, Noldus Information Technology) is equipped with food, water, and a flat-topped, triangular floor shelter with two entry holes that fits in one corner of the cage. The top unit is equipped with LED units and an IR camera so that the mice can be monitored undisturbed during both the light and dark cycles for several days. The mice were individually housed in the cages, and activities were recorded and tracked by EthoVision XT 8.5 for 3 d. Activities on the third day were analyzed after 2 d of habituation in the cages.
CatWalk automated gait analysis.
The CatWalk (Noldus Information Technology) is a proprietary apparatus for assessing small animal gait (Gabriel et al., 2007). The walking platform consists of a transilluminated glass floor (100 cm × 12 cm × 0.6 cm) fitted with a portcullis at both ends that prevents the mouse from leaving the walkway. Internal reflection inside the glass walkway contains the light, except at points where an object touches the glass, where the light exits, scatters, and illuminates the contact area with intensity depending on the pressure exerted. These areas of contact are video recorded while mice traverse the walkway, and the video is analyzed using proprietary software (CatWalk XT 9.1, Noldus Information Technology), which assigns and categorizes each point of contact as right/left and fore/hindpaw, allowing the footfalls and gait to be assessed quantitatively. The points of contact are also examined by the investigator and corrected manually if necessary. The experiment was performed using the manufacturer's recommended CatWalk parameters, including a threshold value for skin-floor contact of 0.10 and video camera gain of 14.72 db. Each individual mouse had to remain on the walkway until three successful trials with a maximum variance in speed of 60% and maximum run duration of 6.0 s were achieved.
Balance beam test.
The balance beam apparatus consisted of two platforms, 4 inch square, placed 18 inches apart, with the starting platform elevated 4 inches above the base of the apparatus and the ending platform elevated 8 inches. A dark plastic enclosure (“safe box”) was placed on the 8 inch ending platform. The platforms were connected by a 1 inch square beam during training, and a series of graded round beams (3/4, 5/8, 1/2, and 3/8 inch diameter) during actual testing. Testing was conducted in an enclosed hood under ambient room light. A video camera was placed behind the 4 inch starting platform but elevated slightly above the animal such that hindpaw faults could be easily recorded and the head of the animal was always visible (Movies 1 and 2). Preliminary examination from the side of the apparatus revealed that forepaw faults were rare, and these were not recorded.
Each mouse received five training sessions on the 1 inch square beam. During training, the mouse was given up to 1 min to explore the starting platform before being gently coaxed up and across the beam. After the mouse entered the safe box on the ending platform, it was given 30 s to rest inside before returning to the home cage. Cage mates were trained in rotation, through five iterations of training. On the testing day, each mouse was tested five times on each of the four graded beams, in descending order of size, using the same method as the training sessions. Balance beam performance was scored for beam transit time and the number of hindfoot faults observed (Brooks and Dunnett, 2009). Runs in which the mouse left or reversed direction on the bar were excluded as incomplete. The last three completed trials for each animal on a given diameter beam were used to compute the average transit time and the average number of faults for that subject.
Sucrose preference test.
Mice were individually housed in their home cages equipped with two drinking bottles, one containing sucrose solution (1%) and the other containing filtered tap water. Mice were first habituated to their new environment for 3 d. Fluid consumption was measured over 4 d after the habituation period. Mice were allowed to drink freely between the two bottles during both the habituation and testing periods. Both bottles were weighed daily to measure fluid consumption, and the positions of the bottles were switched daily to prevent a side bias. Sucrose preference is calculated as the percentage of the amount of sucrose solution consumed over total fluid consumption [(sucrose solution intake/total fluid intake) × 100]. A second test (including both the habituation and testing periods) was performed to determine reproducibility of the test. All mice were weighed at the end of the second test.
FST.
Mice were individually placed in a 2 L Pyrex beaker filled with ∼12 cm of water (30 ± 1°C) so that they were not able to touch the bottom at any time and left to swim for 15 min. The trials were recorded using a video camera from the side, and recordings were manually scored in 5 s bins for the predominant behavior. An animal was judged immobile when it was floating passively in the water and only making the necessary movements to keep the head above the water.
In vitro electrophysiology recordings.
Electrophysiological studies of the dMHb were performed with dMHbChR2 mice of either sex between 4 and 8 weeks of age. For the preparation of acute brain slices used in loose-seal, cell-attached recording, adult mice were anesthetized with isoflurane, quickly decapitated, and the whole brains were dissected into cold aCSF containing the following (in mm): 118 NaCl, 3 KCl, 25 NaHCO3, 1 NaH2PO4, 1 MgCl2, 1.5 CaCl2, 10 glucose, pH 7.4. Coronal brain slices (300 or 350 μm) were cut with a vibratome (Leica VT1000S) in cold oxygenated aCSF, and slices were allowed to recover for at least 30 min at room temperature in oxygenated aCSF before transferring to the recording chamber. During recording slices were superfused with aCSF at 30°C.
Loose-seal cell-attached recordings from dMHb neurons were obtained using blind-patch technique in slice preparations using bright-field optics to visualize the habenula. Recording electrodes (3–7 mΩ) were filled with 10 mm glucose aCSF, and recordings were obtained with a MultiClamp 700B Amplifier (Molecular Devices), connected to a digitizing interface (Digidata 1440), and analyzed with pClamp 10.2. Neurons expressing ChR2 in dMHbChR2 mice were identified based on the presence of light-induced action potentials. Optical stimulation of brain slice preparations to activate ChR2 was supplied by a 447 nm blue laser (OptoEngine) coupled to a 200 μm, 0.22 NA optical fiber (Thor Laboratories). The optical fiber was placed 5 mm from the habenula (2.0 mW/mm2). Neurons expressing eNpHR3.0 in of dMHbNpHR mice were identified based on light-induced inhibition of spontaneous action potentials. Optical stimulation to inhibit eNpHR3.0 was supplied by a 640 nm red laser (OptoEngine), and the optical fiber was placed 4 mm from the habenula (46 mW/mm2).
dMHb optogenetic stimulation studies.
Optogenetic behavioral studies were performed with dMHbChR2 optogenetic mice and Syt6Cre littermate controls. All tests were conducted using male mice between 4 and 7 months of age. Mice were implanted with a bilateral fiber optic cannula consisting of two 100 μm, 0.22 NA optical fibers (Thor Laboratories) that were spaced 0.7 mm apart and protruded 2.6–2.8 mm from a Delrin guide. The cannulas were placed at (bregma −1.6 mm and ±0.35 mm from the midline). The implanted optical fibers were connectorized with a single zirconia ferrule, linked in turn to a 200 μm fiber optic patch cord (Thor Laboratories), an optical rotary joint (Doric Lenses), and a 100 mW, 473 nm DPSS laser (OptoEngine). Laser pulses were controlled by an optomechanical switch (Avantes) linked to a programmable microprocessor (Arduino, SparkFun Electronics). Duration and timing of the laser pulse trains were controlled and mouse body movements were tracked by EthoVision XT (Noldus Information Technology).
Optogenetic open field locomotion and ICSS of the dMHb.
Open field locomotion was performed in a 2.59 cm L × 17.78 cm W × 12.7 cm H behavioral chamber with a 100 mA house light (Med Associates). Locomotor behavior was recorded, and the distance traveled was calculated using EthoVision XT. To avoid increased locomotion that may result from exposure to a novel environment, mice were habituated to the chamber before the test. The recorded session consisted of 20 min of baseline recording, followed by 10 min of intermittent laser stimulation, consisting of 30 s stimulation periods and 30 s laser off periods. Stimulation periods consisted of 25 ms, 20 Hz pulse trains, with a 473 nm laser delivering 8 mW total power.
To test ICSS reinforcement by the dMHb, mice were placed into an operant chamber (ENV-307W, Med Associates) with two response wheels that record one event for every 90 degrees of rotation (ENV-113AMW, Med Associates). The wheel turning responses were recorded through baseline, stimulated, and reversal sessions as described below.
Baseline sessions.
Spontaneous wheel turning behavior was recorded for three 45 min baseline sessions during which no laser stimulation was delivered. To counterbalance any initial side bias in wheel turning behavior, the wheel recording the larger number of turns in the baseline sessions was designated the “initially preferred side” (P) and the wheel with fewer turns was designated “the initially nonpreferred side” (NP). For the initial sessions with delivery of a laser stimulus, the NP wheel was chosen as the “active” wheel associated with the laser simulation. Thus subjects must overcome any initial bias to develop a preference for the active wheel in stimulated sessions.
Stimulated sessions.
After the baseline sessions, mice underwent three 45 min ICSS sessions in which the active (NP) wheel was associated with the delivery of laser stimulation via an implanted fiber optic cannula (above), and the P wheel was inactive. Each trial consisted of a variable response period during which wheel turning events on the active wheel were recorded until the reward criterion was met, followed by a 20 Hz train of 25 ms, 8 mW laser pulses for 4 s, and a 2 s time out period. Wheel turning events during the period of laser stimulation and time out were recorded but did not count toward a reward. Wheel turning events on the inactive wheel were recorded; during the response period, responses on the inactive wheel reset the event counter to 0. Within each session, the reward criteria were escalated such that the first five trials required one response, the next five trials required five responses (the “training trials”), and subsequent trials required 25 responses to trigger delivery of the laser pulse train. The low fixed ratio of response:reward in the training trials was designed such that all mice received some rewards as a result of spontaneous wheel turning.
Reversal sessions.
The reinforcing effect of ICSS was then confirmed in reversal sessions during which the active wheel was switched to the P side. The trial and session structure in the reversal sessions was the same as the initial stimulation sessions. Mice received three or six 45 min crossover sessions; only the last three sessions were used for analysis. The same cohort of mice was used for all the ICSS experiments (n = 12 dMHbChR2 and 10 control mice).
IP optogenetic inhibition studies.
Optogenetic inhibition studies were performed with dMHbNpHR optogenetic mice and controls that lacked either the Ai39 or Cre driver allele. All tests were conducted using male mice between 6 and 8 months of age. Mice were implanted with a cannula consisting of a single 100 μm, 0.22 NA optical fiber (Thor Laboratories) that protruded 5 mm from the connecting ferrule and Delrin guide. The point of insertion at the skull was at bregma −3.52 mm and lateral −1.57 mm, at a −20 degree angle relative to the midline. The targeted coordinates for the cannula tip were as follows: rostrocaudal, bregma −3.52 mm; dorsoventral, −4.29 mm ventral to the lambda-bregma line; lateral, at midline. The implanted optical fibers were connected as described for the optogenetic activation experiments. The tracking of mouse body movements and event-triggered laser activation were controlled by EthoVision XT (Noldus Information Technology).
Optogenetic open field locomotion and acute place preference/aversion in dMHbNpHR mice.
The open field sessions for dMHb output inhibition in control and dMHbNpHR mice consisted of 5 min of habituation, 10 min of baseline activity, and 10 min of intermittent laser illumination, consisting of 30 s laser-on periods and 30 s laser-off periods (640 nm, 8 mW from the end of the patch cord).
Optogenetic place preference/aversion studies were conducted in a three-chamber place-preference box (ENV-3013, Med Associates). The laser-active lateral chamber was distinguished by vertically striped walls and stainless steel grid rods for flooring, whereas the inactive lateral chamber had horizontally striped walls and metal mesh flooring. The small central entry chamber had a smooth floor. Recording and event triggering were controlled with EthoVision XT using center point tracking. Sessions were initiated by placing the mouse in the center compartment of the apparatus. The mice were then given free access to all three chambers for 15 min. When the mouse entered the active chamber, the 640 nm laser was activated, delivering continuous illumination at 8 mW of total power from the end of the patch cord, which remained on until the mouse entered the inactive chamber. Conversely, after occupancy of the inactive chamber, the laser remained off until the mouse entered the active chamber. Thus, time spent in the central chamber was scored as active chamber time or inactive chamber time depending on the state of the laser linked to the immediately prior occupancy of a lateral chamber. Mean time spent in the center chamber was (give percentage), and inclusion or exclusion of the center time did not significantly affect the analysis. Active and inactive chamber occupancy and total distance traveled were then calculated for each 5 min interval of the 15 min session.
General statistical methods.
Statistical analyses were conducted using unpaired or paired two-tailed t tests, two-way ANOVA, or two-way repeated-measures ANOVA with Bonferroni's post hoc analyses in GraphPad Prism 6 (GraphdPad Software). Results are presented as mean ± SEM.
Results
Conditional deletion of the transcription factor gene Pou4f1 results in nearly complete loss of dMHb neurons
Mice that are homozygous for a constitutive Pou4f1 (Brn3a, Brn3.0) null allele die at birth, probably from the loss of brainstem cholinergic neurons in the nucleus ambiguus (McEvilly et al., 1996). To circumvent this neonatal lethality, and allow the habenula-specific excision of Pou4f1, we generated a line of mice in which exon 2 of the Pou4f1 gene, containing nearly all of the coding sequence, is flanked by loxP sites (Fig. 1A). We then excised the neomycin selection cassette from transgenic founders to yield the allele Pou4f1flox (see Materials and Methods). Pou4f1flox/flox mice are fertile, exhibit no obvious behavioral abnormalities, and express levels of Brn3a protein that are not distinguishable from wild-type by immunofluorescence. Mice lacking tissue-specific expression of Brn3a were generated by interbreeding mice bearing an appropriate Cre driver and a constitutive Pou4f1 null allele with mice bearing the genotype Pou4f1flox/flox. The null allele used, Pou4f1tlacZ, contains a lacZ expression cassette replacing the Pou4f1 coding sequence that allows the identification of neurons normally fated to express Brn3a in the absence of Brn3a protein (Quina et al., 2005). A database screen for potential Cre drivers for the dMHb revealed a synaptotagmin-6 BAC Cre line (Syt6Cre), generated by the Gensat project, as a candidate driver. Mice bearing this allele were recovered from a repository (see Materials and Methods), and its specificity was tested using a reporter line, Ai6, which conditionally expresses the fluorescent protein ZsGreen (Madisen et al., 2010). An examination of ZsGreen expression in the habenula of phenotypically normal mice showed that Cre activity in the MHb was confined almost entirely to the dorsal part, and activation of reporter expression in scattered neurons in the LHb was also noted (Fig. 1B,C). Although Syt6Cre expression is not specific for the habenula, an examination of the sites of induced reporter expression throughout the CNS revealed very little coexpression of ZsGreen with Brn3a outside the habenula (Fig. 1D–I). Conditional knock-out (dMHbCKO) mice used in the behavioral studies described below have the genotype Pou4f1flox/tlacZ, Syt6Cre, and control mice have the genotype Pou4f1flox/+, Syt6Cre. Because Syt6Cre and Brn3a coexpression is largely restricted to the habenula, any loss of function induced by this combinatorial strategy should be habenula-specific.
When we examined adult dMHbCKO mice for the expected loss of Brn3a protein using immunofluorescence, we immediately observed that the dMHb in these mice exhibited profound cell loss. Immunostaining for Brn3a and β-gal protein expressed from the Pou4f1tlacZ allele showed that the dMHb was markedly reduced in size, and also that the small number of surviving β-gal+ dMHb neurons often lacked Brn3a protein (Fig. 2A,B). Because the extent of the dMHb is not markedly reduced in late embryonic mice (Quina et al., 2009), we conclude that Brn3a is required for the postnatal survival of dMHb neurons. Immunostaining for the vMHb-specific marker ChAT showed that the vMHb was of normal size in dMHbCKO mice (Fig. 2C,D). Quantitative analysis of the left/right cross-sectional area of the dMHb along its full rostrocaudal extent (bregma −1.0 to −2.1 mm), defined by cells expressing the β-gal marker but not ChAT, showed that the dMHb was reduced to 25.8% of its normal volume (F(3,114) = 179.4, p < 0.0001; Fig. 2E). Two-way ANOVA instead of repeated measures two-way ANOVA was used because of missing measurements in the dataset (occasional lost sections). Two-way repeated-measures ANOVA was conducted to examine the extent of the dMHb lesion in a large behavioral cohort of mice, which shows a consistent effect on cell loss (F(1,21) = 552.7, p < 0.0001; Fig. 2F).
We also used ISH to characterize the expression of mRNA for Tac1 (encoding SP) and the K+-channel interacting protein Kcnip1, which have been previously described as specifically expressed in the dMHb and vMHb, respectively (Quina et al., 2009). The area of Tac1 expression was markedly reduced, whereas the area of Kcnip1 expression appeared unchanged (Fig. 2G–J).
Ablation of dMHb neurons reduces voluntary wheel running and rotarod performance with relatively minor changes in other motor tests and normal locomotion
To assess the role of dMHb in behavior, we began by measuring voluntary WRA in control and dMHbCKO mice. Most rodents will voluntarily engage in WRA, and in rats it has been shown to specifically act as a reinforcer for other behaviors, such as a lever press, leading to studies of the interaction of WRA with other reward systems (Sherwin, 1998; Novak et al., 2012). Nocturnal animals, such as mice, housed under conditions with an alternating LD condition, exhibit most of their WRA in the dark phase. Both control and dMHbCKO mice entrained to the 12:12 LD cycle and, as expected, exhibited more WRA during the dark period between zeitgeber time (ZT)12 (time of lights off) and ZT0 (time of lights on; main effect of time, F(23,529) = 86.68, p < 0.0001; Fig. 3A). However, the overall WRA of dMHbCKO mice was greatly reduced compared with control mice (F(1,23) = 17.19, p = 0.0004). Post hoc analyses indicated that this effect was driven by reduced activity of the dMHbCKO mice during the first half of the night between ZT12 and ZT17 (Fig. 3B).
The observed reduction in WRA of the dMHbCKO mice could result from a motor deficit leading to decreased locomotor ability or from decreased motivation to run, possibly reflecting an interaction with reward mechanisms. For this reason, we assessed the dMHbCKO mice in a series of motor tests, some of which also involve a sustained motivational component. We first administered the rotarod test, in which mice are required to stay on top of a rotating drum to avoid a harmless but aversive fall. The rod accelerates throughout the test, and the time elapsed until a fall is the principal outcome. A short falling latency indicates a deficit in the animal's ability or drive to maintain an upright position on the device. Although frequently used as a test of motor function, the rotarod includes a motivational component (incentive to avoid falling from the device), and chronic stress, a model of depression, impairs rotarod performance (Mizoguchi et al., 2002). After one training day at fixed speed, mice were tested on the accelerating rotarod for 2 d of four trials each (see Materials and Methods). The dMHbCKO mice exhibited shorter latency to fall compared with control mice across the four trials in the first testing session (F(1,23) = 24.95, p < 0.0001; Fig. 4A) and over the four trials combined on day 1 (t(23) = 5.06, p < 0.0001; Fig. 4B). In repeated testing, the performance on day 1 is predicative of that on day 2 (r = 0.66, p = 0.0003; Fig. 4C). The latency increased across trials on day 1 (F(3,69) = 5.34, p = 0.0023), and post hoc analyses show a significant difference between the control and dMHbCKO mice from trials 2–4 (p < 0.05), which suggests that the control mice also improved their performance as the trials progressed but the dMHbCKO mice did not. Although exposure to WRA has been shown to improve rotarod performance (Clark et al., 2008), the rotarod assessment preceded the WRA trials in this cohort of mice, so the deficit in rotarod performance is not a secondary effect of the differential WRA of control and dMHbCKO mice.
We next assessed the dMHbCKO mice for spontaneous motor activity in an open field enclosure and in their home cages, and we performed a detailed assessment of their gait. In the open field test, dMHbCKO mice traveled a similar distance compared with controls (t(23) = 0.86, p = 0.40; Fig. 4D). In this test, the fractional time spent in the periphery versus the center of the enclosure is frequently used as a model of anxiety, and no significant difference was observed between control and dMHbCKO mice in time spent in these areas (p > 0.05). To assess spontaneous locomotion in a more naturalistic environment, we usedPhenoTyper cages (Noldus Information Technology), which allow extended, undisturbed tracking of motor behavior in a home cage setting. As in the WRA experiments, both genotypes entrained to the 12:12 LD cycle used in the housing area and were more active during the dark period (main effect of time, F(23,529) = 13.03, p < 0.0001). No difference was observed in distance traveled between control and dMHbCKO mice during the 24 h testing phase (F(1,23) = 0.33, p = 0.57; Fig. 4E).
To determine whether the WRA and rotarod performance results could be attributed to ataxia or other impairment of gait, we performed a detailed assessment using the CatWalk (Noldus Information Technology), which allows optical measurement of several gait parameters while the animal traverses a glass track (Table 1) (Gabriel et al., 2007). The dMHbCKO mice showed a nonsignificant trend toward shorter run duration compared with the control mice (t(23) = 0.99, p = 0.33). As expected, the shorter run time may be a result of faster limb swing speed and longer stride length for the dMHbCKO mice compared with controls. Limb swing speed normalized for run duration showed no difference between control and dMHbCKO mice, nor did other measures of the form or regularity of gait (stands, duty cycle, and regularity index; Table 1). Small but significant differences were observed between control and dMHbCKO mice in two other parameters, the forelimb “base of support” (BOS; Table 1), the measured distance between the front paws while walking, and the “print position,” the distance between the imprint of a given hindpaw and the previous position of the ipsilateral front paw during the gait cycle (Table 1). The forelimb BOS was reduced for dMHbCKO mice compared with the control mice (t(23) = 2.51,p = 0.020), and the print position was closer (right paws, t(23) = 2.69, p = 0.013; left paws, t(23) = 2.58, p = 0.017). We note that, like the WRA and rotarod tests, the cohort of dMHbCKO mice showed a fairly broad distribution of individual performance values on these tests. However, no significant correlation was observed between forelimbs BOS to WRA (control, r = −0.057; dMHbCKO, r = −0.11) or rotarod performance (control, r = −0.58; dMHbCKO, r = −0.27). Similarly, no correlation was observed between print position (right paws) and WRA (control, r = −0.15; dMHbCKO, r = −0.12) or latency to fall in the rotarod test (control, r = 0.22; dMHbCKO, r = and −0.15). Thus, the observed small differences in these gait parameters are unlikely to explain, or be related to, the reduced WRA or the impaired rotarod performance seen in the dMHbCKO mice.
Subtle defects in motor function that may be missed in simple gait analysis may be detected by requiring mice to traverse a narrow balance beam (Brooks and Dunnett, 2009). The incentive to complete the task can be increased by requiring that mice cross a beam from an exposed platform to a platform enclosed by a “safe” box. To implement this test, we trained mice of both genotypes to perform the task on a 1 inch square beam, which did not present a physical challenge, then tested them with increasingly difficult round beams of 3/4, 5/8, 1/2, and 3/8 inch diameter. Trials were scored for the time of transit from the exposed platform to the safe box and for “faults” in which a hindpaw lost contact with the beam, resulting in a sudden shift in posture. Although this test is in a sense voluntary, mice are highly motivated to escape to the safe location, and 100% of animals of both genotypes rapidly acquired the behavior. The abilities of control and dMHbCKO mice to traverse the beam were grossly similar, as seen in side-by-side videos of mice of both genotypes that exhibited either zero or one fault in a representative trial (Movies 1 and 2). Despite this, the dMHbCKO mice traversed the bars somewhat more slowly than controls at all diameters (F(1,68) = 11.70, p = 0.0011; Fig. 4F). Both control and dMHbCKO mice exhibited an average of less than one fault per trial, but dMHbCKO mice committed significantly more faults than controls (F(1,68) = 12.62, p = 0.0007; Fig. 4G). Post hoc analyses did not show a difference between beam sizes for either behavioral analysis.
Mice with dMHb lesions exhibit alterations in sucrose preference, a measure of hedonic state, but not in the FST
As a model of hedonic state, we assessed the behavior of dMHbCKO mice in the sucrose preference test. In this two-bottle choice paradigm, mice normally prefer a sweetened source of drinking water, and loss of this preference is observed in models of stress and depression (Willner et al., 1996; Strekalova et al., 2004; Couch et al., 2013). dMHbCKO mice showed a significantly reduced preference toward the sucrose solution compared with control mice averaged over the 4 d testing period (t(23) = 2.14, p = 0.043; Fig. 5A). Individual dMHbCKO mice exhibited sucrose preference ranging from normal (∼80% sucrose consumed) to sucrose aversion (<50% sucrose consumed). To better understand this range of outcomes, we performed a complete replication of the test, which showed retention of individual preferences and a high test–retest correlation (r = 0.82, p < 0.0001; Fig. 5B). Average total fluid consumption across the 4 d testing period (F(1,23) = 0.04, p = 0.85) and average body weight (t(23) = 1.02, p = 0.32) at the conclusion of the test were not significantly different between genotypes. The sucrose preference test was performed before WRA assessment in this cohort of mice, excluding a secondary effect of WRA experience on sucrose preference.
Given the reduced WRA and evidence for an anhedonic state seen in the dMHbCKO mice, we also evaluated these mice using the FST, a model of stress response and depression (Duman, 2010). Initially, mice will swim vigorously if placed in a water-filled cylinder, but if escape is not possible, they will eventually stop swimming and exhibit only the minimum movements necessary to keep their heads above the water. The time that the animals spend immobile is recorded as the principal outcome of the test. Both the control and dMHbCKO mice swam vigorously in the beginning of the test, but immobility increased as the test progressed (F(14,196) = 6.72, p < 0.0001; Fig. 5C). Both genotypes exhibited a similar degree of immobility time during the 15 min of testing (F(1,14) = 1.11, p = 0.31). Collectively, our data indicate that the dMHbCKO lesion has a marked effect on motor behaviors with a significant sustained motivational component, such as WRA and the rotarod test, and also affects hedonic state as measured by sucrose preference. In contrast, this lesion has much less impact on gait and balance beam performance, and no detectable effect on basal locomotion in an open field or home cage. The most straightforward explanation for the marked reduction in WRA in dMHbCKO mice would be a role of the dMHb in mediating the motivational aspects of WRA, perhaps through a link to reward pathways. However, it is also possible that the subtle gait and balance beam deficits detected in these mice could represent a slight motor deficit that could alter the motivation to run indirectly, perhaps by increasing the difficulty (i.e., cost) of the task. To help resolve this question we turned to optogenetic models of dMHb activation and silencing to assess the role of the dMHb in primary reinforcement and aversion. These acute, reversible optogenetic models also help to overcome the intrinsic limitation that the dMHbCKO mice have a permanent developmental lesion, which may be partially subject to compensatory mechanisms.
Activation of the dMHb is intrinsically reinforcing
Because dMHb ablated mice have changes consistent with anhedonia, we wished to test whether this nucleus can mediate primary reinforcement. The habenula complex has been shown to support conventional ICSS (Sutherland and Nakajima, 1981), but this method cannot distinguish the habenula itself from fibers of passage, nor the specific role of the habenula subnuclei in mediating the response. Thus, we used an optogenetic strategy to test whether dMHb neurons can mediate reinforcement. Experimental mice expressing a channelrhodopsin 2-EYFP fusion protein (ChR2-EYFP) specifically in the dMHb were generated by interbreeding the Syt6Cre driver line, also used to generate the dMHbCKO mice, with the transgenic strain Ai32 (Madisen et al., 2012), which conditionally expresses ChR2-EYFP from the Gt(Rosa)26Sor locus (Fig. 6A). As expected, expression of the induced ChR2-EYFP protein was restricted largely to the dorsal part of the MHb, and terminal fibers were present in the IPL (Fig. 6B). In the behavioral studies described below, experimental mice expressing ChR2 (dMHbChR2) have the genotype Ai32, Syt6Cre, whereas control mice have the genotype Ai32.
To test the effectiveness of photostimulation of dMHb neurons in dMHbChR2 mice, we performed loose-seal cell-attached recording from dMHb cell bodies in acute slice preparations. Neurons in the dMHb exhibited slow, irregular spontaneous firing at baseline, consistent with prior reports of spontaneous habenula activity (Kim and Chang, 2005). Trains of blue light pulses elicited 1:1 spiking of action potentials from dMHb neurons at frequencies of up to 20 Hz (Fig. 6C). To perform bilateral light stimulation of the dMHb, we constructed bilateral fiber optic cannulas for permanent implantation (100 μm diameter, 2.4 mm long, spaced at 0.7 mm; Fig. 6D; Materials and Methods). Using stereotaxic surgery (see Materials and Methods), the cannula tips were placed above the dorsal margin of the dMHb, at a point approximately two-thirds of the way from the rostral to the caudal pole of the nucleus (near bregma −1.6 mm). Groups of control and dMHbChR2 mice were cannulated in the same position, which was verified postmortem after the completion of behavioral testing (Fig. 6E).
To determine the effects of dMHb stimulation in vivo, we first used open field locomotion to assess the effects of dMHb activation on motor activity. In the open field test with intermittent laser stimulation, there was a significant interaction between genotype and laser on/off state (F(1,16) = 8.13, p = 0.012). The control mice showed no difference in distance traveled during the laser on and laser off periods (Fig. 6F,H), whereas the dMHbChR2 mice showed greater distance traveled during the laser on periods, especially during the initial intervals of light stimulation (Fig. 6G,H). Post hoc analyses indicated a significant overall effect of light stimulation on locomotion in dMHbChR2 mice (p = 0.0051), whereas control mice did not exhibit a significant change in locomotion with light (Fig. 6H).
To investigate the specific role of the dMHb in primary reinforcement, we implemented an ICSS paradigm using an operant chamber with two response wheels (Fig. 7A,B). ICSS sessions were conducted daily for 9 d. Three baseline sessions were obtained without laser stimulation, during which the mice exhibited low levels of spontaneous wheel turning. The wheel eliciting the larger number of turns in the baseline sessions was designated the initially preferred side (P) and the other wheel the initially nonpreferred side (NP). To counterbalance any preexisting bias, the laser stimulation was initially linked to turns on the NP wheel. Three testing sessions were then performed with light pulses delivered for responses on the NP wheel (sessions 4–6), followed by three crossover sessions with stimulation delivered for responses on the P wheel (sessions 7–9).
Examples of individual responses for one dMHbChR2 mouse and one control mouse are shown in Figure 7C and Figure 7D, respectively. Occasional wheel turning was observed during the baseline sessions in both genotypes (Fig. 7C,D, sessions 1–3). In sessions with stimulation linked to the NP wheel (sessions 4–6), the dMHbChR2 mouse showed a marked increase in wheel turning and a strong preference for the stimulated wheel by session 5. When the stimulation was switched to the P wheel (sessions 7–9), a high response rate on the NP wheel persisted for one session (session 7) but soon extinguished, and the turning behavior switched to the P wheel thereafter. In contrast, the control mouse did not exhibit increased wheel turning or preference (Fig. 7D, sessions 4–9).
We also examined total wheel turns, total stimulations delivered, and the establishment of side preference for the entire control and dMHbChR2 groups. The effect of the laser stimulation was different between control and dMHbChR2 mice, indicated by the significant interaction between genotype and day (F(8,160) = 4.52, p < 0.0001). The dMHbChR2 mice exhibited more wheel turnings compared with the control mice (F(1,20) = 5.31, p = 0.032), and wheel turnings increased over sessions (F(8,160) = 6.34, p < 0.0001). Post hoc analyses indicate a trend for increased wheel turns exhibited by dMHbChR2 mice in the first two stimulated sessions, which reached significance by the third stimulated session (Fig. 7E, session 6; p = 0.013). As expected, the frequency of laser activation (“rewards”) received followed a similar pattern of differential responding by genotype (F(8,160) = 4.54, p < 0.0001; Fig. 7F). The dMHbChR2 mice received significantly more rewards compared with the control mice (F(1,20) = 9.35, p = 0.0062), and the rewards received increased as the sessions progressed (F(8,160) = 29.29, p < 0.0001). Post hoc analyses indicate a trend for increased rewards received by dMHbChR2 mice in the first two stimulated sessions, which reached significance by the third stimulated session (Fig. 7F, session 6; p = 0.0016). Within-subjects analysis for wheel preference for the dMHbChR2 mice was analyzed by combining total wheel turns for the stimulated wheels (NP during 4–6 and P during 7–9) or the nonstimulated wheels (P during 4–6 and NP during 7–9) (Fig. 7G). There was a significant interaction between wheels and sessions (F(5,110) = 3.94, p = 0.0025). The dMHbChR2 mice showed increased wheel turns on the stimulated wheel (F(1,22) = 8.37, p = 0.0085) that increased over sessions (F(5,110) = 2.98, p = 0.015). Post hoc analyses indicate a trend toward increased wheel turning in the first two sessions (Fig. 7G, sessions 4 and 5), which reached significance by session 6 (p = 0.0094). On the initial crossover session (7), wheel preference was mixed as responses to the previously stimulated wheel persisted, but strong preference for the new stimulated wheel was established by session 8 (p = 0.013). Within-subjects analysis for the control mice showed no difference between stimulated versus nonstimulated wheels (p > 0.05, data not shown). Fractional wheel preference expressed as the ratio of NP turns/total turns showed a similar pattern. The dMHbChR2 mice showed a differential relative preference for the NP wheel depending on whether or not it was the stimulated wheel (F(8,160) = 5.03, p < 0.0001). This relative preference for either the stimulated NP wheel (>50% on y-axis) or the stimulated P wheel (<50% on y-axis) persisted across the sessions (F(8,160) = 12.69, p < 0.0001; Fig. 7H). Overall, there was no difference of relative preference between the control and dMHbChR2 mice (F(1,20) = 0.72, p = 0.41), but post hoc analyses revealed a preference for the stimulated wheel by the dMHbChR2 mice that reached significance by the last NP-stimulated session (session 6, p = 0.022) and the last crossover session (session 9, p = 0.0071).
Unlike the mean wheel turn counts for each group, the fractional wheel preference is not skewed by the performances of mice that gave especially large numbers of responses. The establishment and reversal of wheel preference also cannot be explained by a nonspecific increase in motor activity in the dMHbChR2 mice. Together, these results demonstrate that the dMHb mediates primary reinforcement.
Inhibition of tonic dMHb output to the IPL is aversive
MHb neurons exhibit irregular spontaneous firing (Kim and Chang, 2005) and thus provide continuous excitatory input to the IP. If the dMHb mediates primary reinforcement, inhibition of this tonic input would be predicted to be aversive. The two-wheel choice paradigm used to assess dMHb reinforcement is not suitable for assessing aversive responses because the wheel turning response cannot be inhibited below the low baseline rate of this behavior. Instead, we adopted an acute place preference/aversion paradigm, which has been previously used to demonstrate aversive responses to the ChR2-mediated stimulation of the lateral habenula outputs to the tegmentum (Stamatakis and Stuber, 2012a).
Experimental mice expressing the light-activated inhibitory chloride pump/fluorescent reporter eNpHR3.0-EYFP specifically in the dMHb (dMHbNpHR mice) were generated by interbreeding the Syt6Cre driver line (also used to generate the dMHbCKO and dMHbChR2 mice) with the transgenic strain Ai39 (Madisen et al., 2012), which conditionally expresses eNpHR3.0-EYFP from the Gt(Rosa)26Sor locus. Control mice were derived from the same crosses but lacked either the opsin or Cre driver allele. Optogenetic inhibition of the dMHb was first tested in acute brain slice preparations. In cell-attached recordings, tonically firing dMHb neurons were effectively inhibited for short periods by continuous red light (Fig. 8A). However, frequent “escape” of the firing from inhibition was observed for light intervals exceeding a few seconds (arrow), presumably because of the strong pacemaker properties of these neurons. For this reason, we opted to target the fiber optic cannula to the terminals of the dMHb fibers in the IPL rather than the dMHb cell bodies (Fig. 8B). Inhibition of action potential propagation by hyperpolarization of these terminals should be sufficient to prevent release of neurotransmitter signaling from the dMHb to the IP and does not require blocking the spontaneous firing of the dMHb neurons.
We first used open field locomotion to assess the effects of dMHb output inhibition in vivo in control and dMHbNpHR mice. After a 10 min assessment of baseline locomotion, dMHb fibers in the IP were exposed to 10 min of intermittent periods of illumination with 640 nm light (30 s on, 30 s off). No difference in distance traveled was observed between the acute light on/off periods for either genotype (F(1,8) = 0.097, p = 0.76). However, there was an interaction between genotype and laser stimulation (F(1,8) = 7.44, p = 0.026), suggesting that the effect of the laser stimulation was different between control and dMHbNpHR mice. Post hoc analyses revealed that the dMHbNpHR mice showed a decrease in the cumulative distance traveled during the entire 10 min illumination trial compared with baseline (p = 0.0034, Fig. 8C).
To study the effects of dMHb inhibition on acute place preference/aversion, we used a three-chamber apparatus with a small central entry chamber, an active lateral chamber associated with continual laser illumination, and an inactive chamber without laser illumination. Because the laser was turned on or off by entry into the active or inactive lateral chambers, respectively, the small fraction of time spent in the middle entry compartment was assigned to whichever lateral chamber the subject most recently occupied. Examples of individual video activity records for one control and one dMHbNpHR mouse in a 15 min place preference/aversion trial are shown in Figure 8D. During the initial 5 min of exploratory activity, the control and dMHbNpHR mice exhibited approximately equal occupancy of the two compartments. However, over the subsequent 10 min, the dMHbNpHR mouse developed aversion to the active chamber, whereas the control was unaffected.
Analysis of chamber occupancy over the course of the 15 min trial for the entire experimental group indicated a development of aversion to the stimulated chamber in dMHbNpHR mice, whereas control mice consistently divided their time approximately equally between compartments (F(1,9) = 5.50, p = 0.044; Fig. 8E). The difference in chamber occupancy achieved significance between genotypes after 5 min (p = 0.046). The dMHbNpHR mice also showed a trend toward less total distance traveled (t(9) = 1.91, p = 0.089; Fig. 8F), which may result from the fact that they confined their movement to a smaller area for most of the trial.
Some caution is warranted in the interpretation of these optogenetic experiments because Syt6Cre is expressed in brain regions other than the habenula, and thus the induced ChR2 or eNpHR3.0 expression is not confined to the dMHb. Like most optogenetic experiments, the specificity of the effects in this model depends both on the pattern of opsin expression and on locally targeted illumination. In the dMHbChR2 mice, cannulated dorsal to the dMHb, the only ChR2-expressing cell bodies plausibly within the light path are in the habenula. However, Syt6Cre-expressing cortical neurons project to thalamic nuclei that are ventral to the habenula, which lie >1 mm from the cannula tip. Based on in vitro models for the transmission of 473 nm light in brain tissue (Aravanis et al., 2007), <1% of the light intensity supplied to the dMHb should be received by these corticothalamic fibers. Although this calculation does not completely exclude in vivo stimulation of these fibers by a cannula placed above the dMHb, this potential off-target effect is entirely circumvented by the IP inhibition experiments in dMHbNpHR mice, in which the optical fiber bypasses the thalamus. Some ChR2-expressing neurons are also present in the LHb of dMHbChR2 mice and may receive some illumination from a cannula placed dorsal to the dMHb. However, manipulation of LHb input to this area has been shown to have the opposite effects to those described here, with the activation of LHb fibers leading to aversion (Stamatakis and Stuber, 2012a). Thus, it is unlikely that either the reinforcement mediated by dMHb activation or the aversion mediated by IP inhibition can be attributed to the off-target stimulation of LHb neurons or the inhibition of LHb terminal fibers in the tegmentum.
Discussion
The habenula has been implicated in the regulation of diverse behavioral functions, including learning and memory, reward, maternal behaviors, sleep, stress and affective states, and the behavioral effects of nicotine (Klemm, 2004; Lecourtier and Kelly, 2007). However, most prior studies have disrupted habenula function with lesions of the entire nucleus or the fasciculus retroflexus (FR), the output tract for the entire nucleus, and thus have not allowed the specific functions of MHb and LHb, or their subnuclei, to be discerned. The distinct neurotransmitter properties of the dMHb, characterized by the expression of the neuropeptide SP, and the specific projection of the dMHb to the lateral subnucleus of the IP, have been known for >30 years (Swanson and Cowan, 1979; Lecourtier and Kelly, 2007; Quina et al., 2009), yet no unique function has been assigned to this structure. To determine the role of the dMHb in behavior, we have created two mouse models: a developmental/genetic ablation model of the dMHb using Cre-mediated disruption of the Pou4f1 gene and optogenetic models in which we conditionally expressed channelrhodopsin or halorhodopsin in the dMHb. Together, data derived from these models show that the dMHb mediates exercise motivation, plays a role in the regulation of hedonic state, and supports primary reinforcement.
One of the most profound deficits observed in dMHbCKO mice is a reduction in voluntary WRA. A recent study in hamsters has shown that complete lesion of the FR also results in markedly decreased nocturnal WRA, as well as increased home cage locomotion (Paul et al., 2011). Here, using specific ablation of the dMHb, we observed no difference between dMHbCKO and control mice in baseline locomotion in the open field enclosure (30 min, short-term) or in home cages (24 h, long-term). Direct comparison of these results must be made cautiously because of the different species used. However, it appears likely that specific dMHb lesions in mice affect the motivation to engage in WRA without the general effect on locomotion observed in FR-lesioned hamsters, which may be attributable to the function of other habenula subnuclei. Thus, at least part of the motivation to engage in WRA can be assigned to the dMHb pathway.
WRA in rodents affects multiple physiological processes, ranging from physical activity and energy balance to motivation and regulation of hedonic state (Novak et al., 2012). Rats will lever-press for access to a running wheel or develop conditioned place preference in response to its presence (Iversen, 1993; Greenwood et al., 2011), and self-administration of cocaine is reduced in rats with access to running wheels, suggesting that wheel running is reinforcing and interacts with known reward mechanisms (Cosgrove et al., 2002). Access to a running wheel has also been shown to mitigate depression-like phenotypes in several rodent models, including learned helplessness, and the forced swim and tail suspension tests (Greenwood et al., 2003; Duman et al., 2008). The relevance of these findings to depression in humans is supported by clinical studies of major depression, in which structured exercise is one of the best-documented nonpharmacological therapeutic interventions (Blumenthal et al., 2007; Trivedi et al., 2011).
Despite normal baseline locomotion, dMHbCKO mice also exhibited marked performance deficits in the accelerating rotarod compared with controls. Although the rotarod is frequently used to assess the function of motor circuits (Rustay et al., 2003), performance in this test also requires motivation, which is linked to affective state. Rotarod performance is impaired in rats exposed to chronic stress (Mizoguchi et al., 2002), a model of depression in rodents (Willner, 2005; Wiborg, 2013). Thus, together, the WRA and rotarod deficits in dMHb lesioned mice are likely to result from loss of motivation rather than a direct impairment of motor function. However, we cannot rule out some interaction between the subtle locomotor defects detected by the more sensitive gait analysis and balance beam tests and the motivation to engage in WRA.
dMHbCKO mice also show a reduced preference for the sucrose-sweetened drinking water normally preferred by control mice. A reduction in sucrose preference develops in chronically stressed rodents and is a model of anhedonia, the general loss of pleasure often exhibited by patients suffering from major depression (Willner et al., 1996; Willner, 2005). However, the dMHbCKO mice did not exhibit increased immobility time compared with control mice in the FST, another measure of depression-like states in rodents, suggesting that sucrose preference and the FST are mediated by distinct neural pathways. Both serotonin and opioids may be involved in the regulation of affective state by WRA (Greenwood et al., 2003; Duman et al., 2008). Future research will be needed to determine whether hedonic regulation by the dMHb intersects with serotonergic or opioid signaling pathways.
To determine whether dMHb activity could mediate reinforcement, we developed an optogenetic version of ICSS, which in its classic form is widely used to assess the ability of brain regions to support reinforcement (Carlezon and Chartoff, 2007). Conventional ICSS using stimulation electrodes has suggested that habenula stimulation may be reinforcing (Sutherland and Nakajima, 1981; Nakajima, 1984; Morissette and Boye, 2008), but it is difficult to stimulate specific neural populations, such as the dMHb with ICSS, or to distinguish effects on neurons from fibers of passage. In optogenetic ICSS, using a two response-wheel choice paradigm, dMHbChR2 mice show a strong preference for the wheel associated with laser stimulation and markedly increase wheel turning in response to this reinforcement. In contrast, inhibition of dMHb fiber terminals in the IP of dMHbNpHR mice leads to the rapid development of place aversion. These results, together with the loss of WRA motivation and sucrose preference, support a model in which the dMHb mediates signals with a positive valence. However, unlike the ICSS effects seen in classic studies of the medial forebrain bundle, the dMHb pathway is not obviously linked to dopaminergic (DA) systems. The output circuit of the dMHb is not fully known, but there is no evidence that the dMHb or its primary target, the lateral part of the IP, project to DA centers (Groenewegen et al., 1986). Thus, our sucrose preference data on dMHbCKO mice are consistent with the observation that dopamine-deficient mice still develop sucrose preference (Cannon and Palmiter, 2003).
Recent work has shown that lesions of the specific septal inputs to the dMHb and vMHb (rather than the MHb nuclei themselves) show dissociable effects on fear conditioning and anxiety, respectively (Yamaguchi et al., 2013). We did not observe a significant difference between dMHbCKO and control mice in open field center/periphery occupancy, one measure of anxiety, a result consistent with this study. In another recent study, Cre-recombinase expressed from the gene locus Gpr151, an orphan receptor expressed in the vMHb and a subset of LHb neurons, was used to ablate neurons in these nuclei (Kobayashi et al., 2013), sparing the dMHb. These mice exhibited a decrease in open field center occupancy, interpreted as increased anxiety. It will be interesting to determine the exact role of the specific dMHb and vMHb inputs, nuclei, and output pathways in fear and anxiety behaviors.
The finding here that the dMHb signaling has a positive valence contrasts markedly with recent findings on the role of the LHb in reward mechanisms. A subset of LHb neurons has recently been shown to mediate negative motivational values, including the absence of an expected reward and punishment, through inhibition of DA neurons (Jhou et al., 2009; Hikosaka, 2010; Bromberg-Martin and Hikosaka, 2011). Synaptic inputs to the LHb also are increased in learned helplessness, a model of depression (Li et al., 2011). The glutamatergic LHb neurons project to the GABAergic RMTg, which in turn inhibits DA reward-encoding neurons in the VTA (Barrot et al., 2012), and optogenetic activation of LHb inputs to this area elicits aversion/behavioral avoidance (Lammel et al., 2012; Stamatakis and Stuber, 2012b).
This LHb-negative motivational system, however, cannot account for the entire role of the habenula in motivated behavior or the regulation of affective states. The existence of a parallel positive reinforcement system somewhere in the habenula can be inferred from prior reports, which have shown that the habenula supports ICSS reinforcement and that habenula lesions block ICSS of several brain regions (Morissette and Boye, 2008). Because the input and output circuitry of the LHb and MHb are quite distinct, as are the dorsal and ventral MHb (Quina et al., 2009; Hsu et al., 2013; Kobayashi et al., 2013), the function of each of these pathways should be considered separately. Our data demonstrate that some of the positive motivational properties of the habenula reside in the dMHb, including exercise motivation and the maintenance of some aspects of hedonic state. In a pilot case, deep brain stimulation of the habenula has been shown to be therapeutic in depression (Sartorius et al., 2010). Although to date this has been interpreted in the context of the LHb-negative motivational system, the findings presented here indicate that the potential role of dMHb activity in clinical depression and the possible therapeutic manipulation of dMHb function should also be considered.
Footnotes
This work was supported by National Institute of Mental Health Awards R21 MH090478 and R01 MH093667 to E.E.T., National Science Foundation IOS-1145796 to H.O.d.l.I., and National Institute of Child Health and Human Development Award P30HD02274 to the Center on Human Development and Disability. Y.-W.A.H. was supported by the National Institute of Mental Health Training Award F32MH098498. We thank Dr. Toby Cole at the Center on Human Development and Disability Mouse Behavior Core for providing behavioral testing equipment and guidance, Dr. Susan Ferguson for extensive advice on behavioral models, Drs. Ben Greenwood and Elizabeth Glenn Guy for helpful comments on the experimental results and manuscript, and Lynne Tempest and Lely Quina for technical assistance. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Eric E. Turner, Seattle Children's Research Institute, 1900 Ninth Avenue, Mail Stop C9S-10, Seattle, WA 98101. eric.turner{at}seattlechildrens.org