Introduction

The lateral habenula (LHb) is an epithalamic brain structure, which acts as an inhibitory modulator of the midbrain reward circuitry, including the ventral tegmental area (VTA)1,2,3 and the raphe nuclei4. Recently, the LHb has drawn renewed attention5 due to the revelation that it becomes pathologically hyperactive in major depressive disorder (MDD)6,7,8,9,10,11,12, thus providing excessive inhibition to these reward centres and silencing the associated positive emotions. Concurrently, this makes the LHb an intriguing target for novel antidepressant therapies6,7,9,11,12,13,14. Notably, multiple studies have shown that potentiation of inhibition within the LHb yields an antidepressant effect9,10,13,15, likely by disinhibition of the reward circuitry. Despite this, however, the neural circuitry which mediates inhibitory signalling within the LHb is not well understood. To date, very few populations of local inhibitory neurons have been identified within the LHb16,17,18, and external inhibitory input appears to arise primarily from GABA/glutamate co-releasing fibres8,19,20,21. Yet the LHb does receive a large inhibitory input9,17, and recently some exclusively GABAergic fibres have been reported to project to it15,22,23. It is hence likely that other as-of-yet undefined circuitry components mediate inhibitory signalling within the LHb, and studying these may yield further insights into the underlying mechanisms by which it is implicated in the pathogenesis of MDD.

We therefore sought to address this issue of inhibitory control within the LHb. Both parvalbumin (PV) and somatostatin (SOM) have long been accepted as neuronal markers associated with particular sub-populations of inhibitory interneurons within many structures throughout the brain including the neocortex24, hippocampus25 and striatum26. Referring to the question of inhibitory control within the LHb, we therefore asked if these markers also represented distinct sub-populations of inhibitory neurons within this structure, and aimed to characterise the circuitry formed by neurons expressing them. We report three sources of inhibitory input to the LHb arising from locally-targeting PV-positive neurons within the LHb, from PV-positive neurons within the medial dorsal thalamic nucleus, and from SOM-positive neurons in the ventral pallidum.

Results

Excitatory and inhibitory transmission within the lateral habenula

GABAergic signalling has previously been reported in the LHb9,20. To assess the balance of excitation and inhibition within the LHb, we simultaneously recorded spontaneous excitatory and inhibitory currents in LHb neurons using a low chloride intracellular recording solution (Fig. 1A; n = 10 neurons; N = 5 mice). For comparison, we recorded spontaneous events in L2/3 pyramidal neurons in the somatosensory cortex (n = 9; N = 4 mice). While the ratio of excitatory to inhibitory events was far larger in LHb neurons than in cortical neurons (Fig. 1B; 22.9 vs 2.3, respectively), we could clearly observe inhibitory currents with comparable frequency (Fig. 1C; 0.1 ± 0.0 Hz vs 0.5 ± 0.3 Hz, respectively; p = 0.23; two-tailed unpaired t-test), amplitude (Fig. 1C; 13.8 ± 1.3 pA vs 12.6 ± 1.1 pA, respectively; p = 0.49; two-tailed unpaired t-test) and kinetics (Fig. 1D; rise time 1.3 ± 0.2 ms vs 1.0 ± 0.1 ms, respectively; p = 0.11; half-width 2.1 ± 0.3 ms vs 1.7 ± 0.4 ms, respectively; p = 0.45; decay 16.1 ± 3.4 ms vs 13.4 ± 4.1 ms, respectively; p = 0.64; two-tailed unpaired t-test) to those in cortical neurons. Hence we asked whether these events could be mediated by similar inhibitory neurons as those in the neocortex.

Figure 1
figure 1

Spontaneous inhibitory currents within the lateral habenula are similar to those within the somatosensory cortex. (A) Example traces of spontaneous currents in an LHb neuron recorded in voltage clamp configuration. AMPA and NMDA-mediated excitatory currents and GABAA-mediated inhibitory currents could be observed simultaneously using a low-chloride intracellular solution. (B) Relative frequency of excitatory vs inhibitory currents in both LHb neurons (n = 10 neurons from 5 mice) and somatosensory cortex L2/3 pyramidal neurons (n = 9 neurons from 4 mice). (C) Comparison of frequency (left) and amplitude (right) between inhibitory currents in LHb neurons and L2/3 pyramidal neurons. Data are mean ± SEM. (D) Comparison of kinetics (rise time, half-width and decay) between inhibitory currents in LHb neurons and L2/3 pyramidal neurons. Note for 4 LHb neurons, and 1 cortical neuron, no inhibitory currents were observed. For these neurons, frequency value was 0 Hz. However, as there were no measurable currents, these neurons have been excluded from analyses of current amplitude and kinetics.

Characterisation of inhibitory sub-populations of LHb neurons

To address this question, we implemented the use of two well-characterized markers of inhibitory identity within the neocortex: PV and SOM24. To test if these markers were also representative of inhibitory identity within the LHb, we used either transgenic mouse lines where fluorescent reporter proteins are expressed within neurons expressing the respective marker, or an antibody against the respective marker, and co-stained sections from these animals with an antibody against GABA.

We stained slices from C57BL/6 mice (N = 3) with an antibody against PV. While PV-positive neurons were also clearly visible within the LHb (Fig. 2A,B), only a small fraction of these (8.8%; N = 3 mice) were co-labelled with GABA (Fig. 2C). Interestingly, these formed two distinct clusters; one within the medial LHb (Fig. 2Aa) and one within the lateral LHb (Fig. 2Ab) along the rostral-caudal axis. Moreover, the GABAergic PV-positive neurons appeared to be exclusively confined to the lateral LHb (Fig. 2Ac,B), and were very brightly labelled with GABA (Fig. 2Ac), thus indicating a sub-population of PV-positive GABAergic neurons.

Figure 2
figure 2

A sub-class of PV-positive LHb neurons are GABAergic. (Aa) 20X representative confocal micrographs displaying PV-immunoreactivity (left), GABA-immunoreactivity (middle) and merge of both (right) within the LHb. (Ab) Zoom of boxed region in (Aa) depicting PV-immunoreactive neurons which did not co-localise with GABA. (Ac) Zoom of boxed region in (Aa) depicting a PV-immunoreactive neuron which did co-localise with GABA. (B) Schematic illustrating location of PV-immunoreactive only, or PV/GABA co-labelled neurons throughout the LHb in the rostral-caudal plane from one mouse, in which every second 60 µm section was analysed. Sub-nuclear boundaries as defined by Andres et al., (1999) are indicated by dashed lines. Approximate rostral-caudal distances from Bregma (in mm) are indicated. (C) Bar graphs showing total number of PV-immunoreactive and GABA/PV co-localising (Co-L) LHb neurons (left) and fractions of co-localising neurons as a percentage of total PV-immunoreactive and of GABA-immunoreactive neurons (right; N = 3 mice). (D) Breeding scheme for generating SOM-IRES-Cre::Ai9 mice. (Ea) Confocal micrograph from SOM-IRES-Cre::Ai9 slices depicting SOM-positive neurons in the LHb. (Eb) Zoom of boxed region in (Ea) depicting SOM-positive LHb neurons localized within the superior sub-nucleus (LHbS) of the LHb. (F) Bar graphs showing total number of SOM-positive and GABA/SOM co-localising (Co-L) LHb neurons (left) and fractions of co-localising neurons as a percentage of total SOM-positive and of GABA-immunoreactive neurons (right; N = 2 mice).

We next crossed SOM-IRES-Cre mice to Ai927 reporter mice, so as to generate SOM-IRES-Cre::Ai9 (N = 2), which expressed TdTomato in SOM-positive neurons (Fig. 2D). Strikingly, the majority of SOM-positive LHb neurons were clustered in the superior sub-nucleus of the LHb (Fig. 2Eb)28. However, SOM-positive LHb neurons were a primarily non-GABAergic population, as only a small sub-population of these (3.65% respectively) co-localised with GABA (Fig. 2F). Thus we conclude that while a clear sub-population of GABAergic PV-positive neurons is located within the LHb, SOM-positive LHb neurons are largely non-GABAergic.

Despite the clear GABA co-labelling we observed with some PV-positive neurons, however, the antibody background staining was relatively high (Fig. 2A). We therefore also performed more sensitive in situ hybridizations with probes for PV and vesicular GABA transporter, VGAT (Fig. 3; N = 3 mice); and PV and vesicular glutamate transporter 2, VGLUT2 (Supplementary Fig. S1; N = 3 mice). Consistently, we observed that all PV-positive neurons within the lateral LHb (Fig. 3B; n = 10) were VGAT-positive. PV-positive neurons within the medial LHb were VGAT-negative (Fig. 3C; n = 30), but VGLUT2-positive (Supplementary Fig. S1). Together with our immunohistochemistry data (Fig. 2), these results indicate that the majority of PV-positive LHb neurons are non-GABAergic, but they also suggest the existence of a unique sub-class of inhibitory PV-positive neurons located within the lateral LHb.

Figure 3
figure 3

A population of GABAergic PV-positive neurons is located within the lateral LHb. (A) 20X in situ hybridization overview image displaying the LHb co-stained with probes for PV and VGAT. (Ba) Zoom of the left boxed region in (A) displaying VGAT-positive PV-positive neurons in the lateral LHb. (Bb) Bar chart quantifying the number of PV-positive neurons in the lateral LHb which were VGAT-positive. (Ca) Zoom of the right boxed region in (A) displaying VGAT-negative PV-positive neurons in the medial LHb. (Cb) Bar chart quantifying the number of PV-positive neurons in the medial LHb which were VGAT-positive.

PV-positive and SOM-positive LHb neurons form physiologically distinct sub-classes

We next sought to further characterise PV and SOM-positive LHb neurons by assessing their physiological properties. We crossed each Cre-driver line with the Ai927 reporter line to generate PV-IRES-Cre::Ai9 (N = 8) and SOM-IRES-Cre::Ai9 (N = 3) transgenic offspring (Fig. 4A), and used fluorescence-assisted patch-clamp recordings to record from TdTomato-expressing neurons in acute slices from each line (n = 19 PV neurons; n = 24 SOM neurons). We also recorded from a control sample of neurons from the general population of LHb neurons (n = 16 from 5 C57BL/6 mice), and compared passive physiological properties between all groups (Fig. 4B). Resting membrane potential was comparable between all groups (p = 0.72; one-way ANOVA); however clear differences could be observed between input resistances (p < 0.0001; one-way ANOVA). SOM-positive neurons had a far larger mean input resistance than any other group (1364.0 ± 111.3 vs 512.1 ± 67.4 and 366.2 ± 39.1 MΩ for general population and PV neurons respectively; p < 0.0001; Bonferroni’s multiple comparison test). Interestingly, SOM-positive neurons are generally located near the border with the medial habenula (Figs. 2E, and 4J), where neurons are known to have very large input resistances29, and as such these findings lend support to the recently proposed idea of an area of overlap between the lateral and medial habenulae30.

Figure 4
figure 4

PV-positive and SOM-positive LHb neurons form physiologically distinct sub-classes. (A) Breeding scheme for generating PV and SOM-IRES-Cre::Ai9 mice. (B) Comparison of passive physiological properties between all groups. (C) Confocal micrograph depicting localization of PV-positive neurons within the LHb. (D) Initial action potential firing frequency versus current injection for both PV-positive (n = 18; N = 8 mice) and general population (n = 16; N = 5 mice) LHb neurons recorded in C57 mice. Data are mean ± SEM. (E) Schematic illustrating location of patched neurons within the habenular complex. (F) Example traces from PV-positive neurons upon current injection. (Fa): traces from a neuron in the medial LHb. Asterisk indicates a section from the recording displaying characteristic sub-threshold voltage oscillations. (Fb): traces from a neuron in the lateral LHb. Current steps: −50 pA and 50 pA; 1 s duration. (G) Comparison of passive properties between PV positive neurons in the lateral LHb (n = 11) and medial LHb (n = 8). Data are mean ± SEM. (H) Example reconstructions from two PV-positive neurons. (I) As for (D), for SOM-positive LHb neurons (n = 22; N = 3 mice). (J) As for (E), in for SOM-positive neurons. (K) Example traces from one neuron when held at −50 mV and −70 mV upon 1 s injection of 50 pA current. Zoom: comparison of the first induced action potential at each holding potential. (La) Inter-event interval (IEI) between the first and second action potential induced upon depolarising current injection (20–150 pA; 1 s) for neurons at different holding potentials (n = 11). (Lb) Comparison of amplitude of the first (‘early spike’) action potential induced upon depolarising current injection (20–150 pA; 1 s) for neurons at different holding potentials (n = 11). (Lc) Comparison of kinetics (rise time and half-width) of the first action potential induced upon depolarising current injection (20–150 pA; 1 s) for neurons at different holding potentials (n = 11). (M) Example reconstructions from a SOM-positive neuron within the superior sub-nucleus (top; n = 6), and from another SOM-positive neuron out with the superior sub-nucleus (bottom; n = 4).

We next assessed the active physiological properties of each class of neuron by injection of a series of current steps. Consistent with our histological data, most PV-positive neurons were clustered in either the medial or lateral LHb (Fig. 4C,E). There was also a clear reduction in the firing frequency of the first induced action potential in response to depolarising current injection in comparison to neurons in the general population (Fig. 4D; p < 0.0001; two-way ANOVA). This was as a result of the fact that only a minority (4 of 19) of PV-positive neurons exhibited any kind of high-frequency bursting behaviour (Fig. 4Fb), and this was only observed upon larger current injections (Fig. 4D). This was a striking observation as high-frequency bursting has long been considered a hallmark physiological phenotype of most LHb neurons12,16,31. Furthermore, the clusters of PV-positive neurons in the medial and lateral LHb could clearly be differentiated based on their physiological profile (Fig. 4F,G). PV-positive neurons in the medial LHb frequently exhibited sub-threshold voltage oscillations (6 of 8 medial LHb neurons) while their counterparts in the lateral LHb did so very rarely (1 of 11 lateral LHb neurons, p = 0.006; Fisher’s exact test; Fig. 4G). Moreover, a second distinctive population of PV-positive neurons appeared in the lateral LHb, identifiable by their hyperpolarized resting membrane potential (5 of 11 lateral LHb neurons; −77.0 ± 1.3 vs −59.4 ± 2.6 mV for all 19 PV-positive neurons; p = 0.001; two-tailed unpaired t-test) and lack of rebound action potential discharge (Fig. 4Fb), a hallmark phenotype of LHb neurons16,32,33. We also performed morphological reconstruction of these neurons (Fig. 4H). Ten neurons were sufficiently reconstructed to visualize a prolonged section of the axon. In nine of these, this was an unbranching axon, possibly indicative of a projection neuron (Fig. 4H). In the one remaining neuron, we did observe the axon to branch locally and extensively, suggesting this to be a locally-targeting neuron. Taken together, these results indicate that PV-positive LHb neurons form multiple distinct sub-populations based on physiological profile and location within the LHb.

When recording active properties of SOM-positive neurons, a striking pattern quickly emerged. These neurons had a far lower firing frequency than the general population, regardless of input current (Fig. 4I; two-way ANOVA; p < 0.0001; and Bonferroni’s multiple comparison test; p < 0.05), and many of these neurons (15 of 24) displayed a prominent ‘early spike’ upon depolarising current injection (Fig. 4K; defined as spiking latency <20 ms after current injection; current injection 20–150 pA, 1 s). Furthermore, we observed that this spike became more pronounced upon hyperpolarization (Fig. 4K,L; tested in 11 neurons). Upon hyperpolarization as far as −70 mV, we could observe a greater interval between this spike and the next spike in the train (Fig. 4La; p = 0.0014; one-way ANOVA), and that this early spike had faster kinetics than when depolarized (Fig. 4Lc; rise time p = 0.0011; half-width p < 0.0001; two-way ANOVA). Additionally, we reconstructed a subset of these neurons (Fig. 4M; n = 10), and observed that of those in the superior sub-nucleus (n = 6), 5 of these displayed the short stubby dendrites associated with medial habenular neurons33, while those outside of the superior sub-nucleus (n = 4) had the elongated dendrites known to be far more conventional of LHb neurons16,33. Membrane potential-mediated change in firing modality is a hallmark physiological characteristic of LHb neurons12,16,31. Yet these SOM-positive neurons discharge only one action potential as opposed to the bursting discharge commonly displayed by LHb neurons12,16,31, and then continue to discharge action potentials in a tonic train akin to that described in the medial habenula33. These neurons also have huge input resistances (Fig. 4B) and similar morphological properties (Fig. 4M) more comparable to that of medial habenular neurons. Thus, while we conclude that SOM-positive neurons are not GABAergic (Fig. 2), we provide physiological evidence for a sub-class of habenular neurons in the superior sub-nucleus (Figs. 2E and 4J) which possess intermediate characteristics of neurons from both the lateral and medial habenula, adding support to transcriptomic data proposing the existence of an area of overlap between the two regions30.

PV-positive and SOM-positive neurons provide inhibitory input to the LHb

Our data indicate that while some PV-positive LHb neurons are inhibitory (Fig. 2), PV and SOM-positive LHb neurons are likely not primarily inhibitory populations as they are known to be in the neocortex24,34. Additionally, previous studies have reported excitatory PV-positive neurons which project to the LHb21,35 and GABA/glutamate co-releasing SOM-positive neurons which also target the LHb21,36. Considering this information, we next asked whether each of these classes of neurons mediate primarily excitatory or inhibitory transmission within the LHb.

We therefore crossed mice from each Cre-driver line to Ai3237 reporter mice, to generate offspring that express Channelrhodopsin-2 (ChR2) and eYFP in PV or SOM-positive neurons respectively (Fig. 5Aa). We cut acute slices from these mice and recorded postsynaptic potentials in LHb neurons (Fig. 5Ab) in response to photostimulation. In slices from PV-IRES-Cre::Ai32 transgenic mice imaged in the rostral-caudal axis (N = 2), we could visualize both PV-positive somata and fibres within the LHb (Fig. 5C and Supplementary Fig. S2). As with the staining data (Fig. 2), and data from the PV-IRES-Cre::Ai9 line (Fig. 4C), PV-positive somata were mostly confined to two clusters in either the medial or lateral LHb (Fig. 5C and Supplementary Fig. S2). These clusters also appeared densely enriched with fibres, and appeared to roughly correlate to the previously described central sub-nucleus of the medial LHb, or oval sub-nucleus of the lateral LHb28. While this dense enrichment of fibres could be processes from local PV-positive neurons, we speculated that they may also be from upstream PV-positive projection neurons, as these are known to specifically target the oval sub-nucleus of the LHb21.

Figure 5
figure 5

Optogenetic dissection of signalling pathways mediated by PV and SOM-positive neurons within the LHb (Aa) Schematic illustrating breeding scheme for generating PV or SOM-IRES-Cre::Ai32 mice. (Ab) Schematic of recording scheme. LHb neurons were patched and blue light was used to stimulate both local and long-distance PV or SOM-positive neurons. (B) Pie chart quantifying fraction of neurons responsive to photostimulation in PV-IRES-Cre::Ai32 slices. (C) Confocal micrograph depicting location of both PV-positive somata and fibres within the LHb. Arrowheads indicate PV-positive somata. (Da) Example traces from one neuron in which photostimulation elicited an EPSP, and (Db) one in which it elicited an IPSP. Traces are multiple trace averages superimposed over the traces from which the average was derived. Blue square denotes 2 ms LED pulse. (E) Before-after plot showing PSP size for neurons responsive to photostimulation which exhibited EPSPs (n = 16) and IPSPs (n = 6) at −50 mV and with application of either GABAzine (2 µM) or NBQX (10 µM). Two neurons which displayed both excitatory and inhibitory components are highlighted in red. (F) Schematic illustrating location of patched neurons within the habenular complex in PV-IRES-Cre::Ai32 slices. Responses are indicated as follows: non-response highlighted in grey; + EPSP in green; + IPSP in red; + dual response in light blue (GABA/glutamate co-release); + ChR2-expressing neuron in dark blue. (Ga) Example traces overlaid from a neuron in which photostimulation was sufficient to inhibit spontaneous action potential firing (top), and which was completely blocked by GABAzine application (bottom). (Gb) Raster plot of (Ga). Grey area denotes presence of GABAzine. Blue bar denotes light 2 ms light stimulation. (H) As for (B), in SOM-IRES-Cre::Ai32 slices. (Ia) Example traces from a neuron at two different holding potentials in which photostimulation elicited an IPSP. (Ib) Before-after plot from neurons receiving inhibitory input in which the IPSP was pharmacologically blocked (n = 4 neurons). (Ja) Example traces from a neuron in which photostimulation elicited GABA/glutamate co-release. (Jb) Before-after plot from neurons receiving GABA/glutamate co-releasing input in which both excitatory and inhibitory components of the postsynaptic potential were pharmacologically blocked (n = 3 neurons). (K) As for (F), in SOM-IRES-Cre::Ai32 slices.

In acute slices from these mice (N = 19), photostimulation-induced postsynaptic potentials were observed in 29 of 76 (38.2%) neurons recorded (Fig. 5B and D–F). The majority of these (18 of 29; 62.1%) displayed a solely excitatory response (Fig. 5Da). However, solely inhibitory responses were also observed with relative frequency (6 of 29 responsive cells; 20.7%; Fig. 5Db) while on two occasions, we recorded postsynaptic potentials consisting of both an excitatory and inhibitory component (Fig. 5B). Furthermore, inhibitory postsynaptic potentials (IPSPs) were comparably larger than excitatory postsynaptic potentials (EPSPs) (mean 5.6 ± 1.3 vs 2.3 ± 0.6 mV, respectively; p = 0.013; unpaired t-test; Fig. 5E). In three of these neurons, the induced IPSP was sufficiently large to momentarily silence spontaneous action potential discharge of the recorded neuron (Fig. 5G) in a manner which could be blocked by application of GABAzine. In contrast, we report that PV-positive neurons mediate both excitatory and inhibitory input to the LHb.

In SOM-IRES-Cre::Ai32 slices (N = 5 mice), we observed solely inhibitory responses in over half of the neurons tested (Fig. 5H and I; n = 13 of 21 neurons), confirmed by complete blockade upon GABAzine application (Fig. 5Ib; n = 4 neurons) We also frequently observed postsynaptic potentials consisting of both an excitatory and inhibitory component (Fig. 5H,J; n = 5 of 21 neurons); most likely arising from those previously described GABA/glutamate co-releasing neurons in the entopeduncular nucleus21,36, which appeared to be specifically confined to the caudal portion of the LHb (Fig. 5K). Therefore, we conclude that input to the LHb from SOM-positive neurons is primarily inhibitory.

PV-positive LHb neurons mediate local inhibitory transmission

The above data (Fig. 5) shows that PV-positive neurons provide both excitatory and inhibitory input to the LHb, while SOM-positive neurons provide inhibitory input, and also co-release GABA and glutamate. Previous work can explain the origin of these excitatory21,35, and co-releasing21,36 inputs. So where then do these inhibitory inputs arise from?

Our histological data (Figs. 2 and 3) suggests that PV-positive, but not SOM-positive inhibitory neurons, exist within the LHb. To verify this, and that these neurons make local synaptic connections, we used a viral approach. We injected PV-IRES-Cre mice and SOM-IRES-Cre mice with Cre-dependent AAVs encoding ChR2 and eYFP directly into the LHb to drive expression of ChR2 and eYFP in local PV and SOM-positive neurons (Fig. 6Aa and Supplementary Fig. S3). For PV-IRES-Cre mice, we initially injected an AAV containing capsid protein 9 (N = 6), however this occasionally resulted in some infected neurons outside of the LHb, likely due to the high neuronal tropism of AAV938 (data not shown). We therefore also injected mice with an AAV containing capsid proteins 1 and 2 which allowed us to confine eYFP expression to the LHb (Fig. 6B), with the exception of some sparse fibres which we assume to be processes of glutamatergic PV-positive LHb neurons (Supplementary Fig. S1). For SOM-IRES-Cre mice (N = 4) we used the serotype 9 virion. We then created acute slices from these mice and recorded postsynaptic events while photostimulating (Fig. 6Ab and Supplementary Fig. S3).

Figure 6
figure 6

PV-positive neurons mediate local inhibitory signalling within the LHb. (Aa) Schematic illustrating stereotaxic injection protocol of AAV9 or AAV1/2 into the LHb of PV-IRES-Cre mice (N = 8). (Ab) Schematic illustrating electrophysiology recording protocol for LHb neurons following stereotaxic viral injection. Transduced neurons are photostimulated while recording from nearby LHb neurons. (B) Confocal micrograph depicting PV-positive LHb neurons transduced to express eYFP following viral injection. (C) Pie chart quantifying fraction of neurons responsive to photostimulation of PV-positive neurons. (Da) Example traces from one neuron in which a GABAzine-sensitive IPSP could be elicited following photostimulation. Blue bar denotes 2 ms photostimulation. (Db) Before-after plot of the IPSP amplitude from 5 different neurons at −50 mV before and after application of GABAzine (4 neurons). (E) Schematic illustrating location of patched neurons within the habenular complex. Responses are indicated as follows: non-response highlighted in grey; + EPSP in green; + IPSP in red; + dual response in light blue (GABA/glutamate co-release); + ChR2-expressing neuron shown in dark blue.

Upon photostimulation of PV-positive neurons, we could occasionally observe GABAzine-sensitive IPSPs (Fig. 6C,D; n = 5 of 52 neurons tested), which were confined to the lateral portion of the LHb (Fig. 6E). These were observed in slices from mice injected with either AAV variant. Interestingly, we also observed one NBQX-sensitive excitatory response, and one response featuring both a GABAergic and glutamatergic component (Fig. 6C), possibly indicative that glutamatergic PV-positive neurons (Supplementary Fig. S1) make local contacts. In the case of SOM-positive neurons, as expected we observed no synaptic responses (Supplementary Fig. S3), hence confirming that SOM-positive neurons do not make local inhibitory or excitatory contacts. Altogether, these data hence show that some PV-positive, but not SOM-positive, LHb neurons are locally-targeting inhibitory neurons.

Distinct extrinsic inhibitory inputs to the LHb arise from PV-positive neurons in the medial dorsal thalamus and SOM-positive neurons in the ventral pallidum

If not local, then the large inhibitory responses we observed mediated by SOM-positive neurons (Fig. 5H,I) remained unaccounted for. We hence referred to the Allen Brain Atlas, and by searching for both SOM and GAD65/67 expression in the main afferent input regions to the LHb39, we observed that the ventral pallidum (VP) stood out as a region enriched with both SOM-positive and GAD-positive neurons, and speculated that this area may be the origin of the observed input.

Interestingly, recent work has shown that the lateral geniculate nucleus of the thalamus projects to and inhibits the LHb15, and as such we speculated on the possibility that other thalamic neurons could also provide inhibitory input. Indeed, double in-situ hybridizations for PV and VGAT showed neurons positive for both these markers in the MDT (Supplementary Fig. S4), and we therefore sought to address if both the MDT and VP were providing inhibitory input to the LHb.

We performed stereotaxic injection of Cre-dependent AAV9 encoding ChR2 and eYFP into the MDT of PV-IRES-Cre mice (Fig. 7A; N = 6), and into the VP of SOM-IRES-Cre mice (Fig. 7G; N = 3). By targeting injections to the ventral MDT, we could confine injections to this region without infecting the LHb (Fig. 7Aa and D). We recorded from LHb neurons (n = 47) while photostimulating MDT PV-positive neurons and observed inhibitory events in seven neurons (Fig. 7B,C; note in three these were only visible when the neuron was strongly depolarized). Consistently, upon post-hoc confocal imaging, we could visualize fibres which appeared to be projecting dorsally from neuronal somata located in the MDT to the LHb (Fig. 7D). Strikingly, and consistent with our observation of fibre enrichment (Supplementary Fig. S2), these fibres appeared to be exclusively targeting the lateral LHb; particularly the oval sub-nucleus28 where all responsive neurons were recorded (Fig. 7E). We also filled PV-positive MDT neurons (n = 8) with biocytin in slices from PV-IRES-Cre::Ai9 mice (N = 2), and upon reconstruction could observe fibres penetrating the LHb in 5 of 8 neurons (Fig. 7F).

Figure 7
figure 7

Distinct extrinsic inhibitory inputs to the LHb from the MDT and VP. (Aa) Schematic illustrating stereotaxic injection protocol of AAV9 into the MDT of PV-IRES-Cre mice (N = 6). (Ab) Schematic illustrating electrophysiology recording protocol for LHb neurons following stereotaxic viral injection. Transduced PV-positive neurons are photostimulated while recording from nearby LHb neurons. (B) Pie chart quantifying fraction of neurons responsive to photostimulation. (C) Example traces from one neuron in which an IPSP could be elicited following photostimulation. Blue bar denotes 2 ms photostimulation. (Da) Confocal micrograph depicting eYFP expression within the MDT following stereotaxic viral injection. (Db) Zoom of the outlined region in (Da). Note the presence of fibres originating from neurons in the MDT which penetrate the LHb. (E) Schematic illustrating location of patched neurons throughout the habenular complex. All neurons responsive to photostimulation were located in the oval sub-nucleus of the lateral habenula. Sub-nuclear boundaries as defined by Andres et al., (1999) are indicated by dashed lines. Approximate rostral-caudal distances from Bregma (in mm) are indicated. (F) Example reconstructions from 2 PV-positive MDT neurons which had axonal fibres projecting to the LHb. (Ga) Schematic illustrating stereotaxic injection protocol of AAV9 into the VP of SOM-IRES-Cre mice (N = 3). (Gb) Schematic illustrating electrophysiology recording protocol for LHb neurons following stereotaxic viral injection. Transduced SOM-positive neurons are photostimulated while recording from postsynaptic LHb neurons. (Ha) Confocal micrograph depicting transduced neurons following injection into the VP. (Hb) Confocal micrograph depicting terminals from transduced VP neurons in the LHb. (I) As for (B), in SOM-IRES-Cre slices. (Ja) Example traces from one neuron in which a GABAzine-sensitive IPSP could be elicited following photostimulation. Blue bar denotes 2 ms photostimulation. Traces are averages taken from multiple sweeps. (Jb) Before-after plot showing amplitude of photostimulation-induced IPSPs before and after application of GABAzine (n = 6 neurons). (K) As for (E), in SOM-IRES-Cre slices.

More than three weeks post-injection in SOM-IRES-Cre mice (N = 3), we could clearly see virally-transduced eYFP-expressing neurons in the ventral pallidum (Fig. 7Ha), which had dense terminals in the anterior portion of the LHb (Fig. 7Hb and K). Consistently, recording from LHb neurons in acute slices while photostimulating, we observed large (−8.8 ± 2.3 mV) GABAzine-senstive IPSPs in 6 of 13 recorded neurons (Fig. 7I,J), consistent with the solely inhibitory responses observed in data from our SOM-IRES-Cre::Ai32 transgenic line (Fig. 5H,I). Hence, taking this data altogether, we report two distinct source of extrinsic inhibitory input to the LHb, arising from PV-positive neurons in the MDT, and SOM-positive neurons in the VP.

Discussion

In this study we sought to address the issue of inhibitory control within the LHb. We implemented the use of two markers known to represent distinct sub-populations of inhibitory neurons within the neocortex24, hippocampus25 and striatum26. We provide evidence for three sources of inhibitory input to the LHb mediated by separate populations of PV-positive neurons within the LHb itself, and in the medial dorsal thalamus; and SOM-positive neurons in the ventral pallidum.

PV-positive neurons provide input to the LHb from two distinct sources

The subject of locally-targeting inhibitory neurons within the LHb has classically been a topic of debate16,40. In this work we show that at least two distinct populations of PV-positive neurons provide inhibitory input to the LHb; one locally-targeting and one located ventrally in the MDT. Notably, LHb innervation from each of these sources of inhibition does appear relatively sparse (Figs. 6 and 7). Thus, one could question whether each source individually can provide sufficient inhibition to be functionally relevant. However, prefrontal cortical innervation of the LHb could also be considered sparse41, and yet has been shown to influence working memory42. Perhaps then, each of these sources could indeed be individually capable of influencing LHb function, or perhaps their inhibitory input could converge upon the same postsynaptic neurons to result in an overall greater inhibition of LHb activity. Strikingly, both of these neuronal populations appear to selectively target neurons located within the lateral LHb; particularly the oval sub-nucleus (Figs. 6E and 7E). It is a relatively novel idea that topographical organisation of LHb outputs exists43, and interestingly it is believed that neurons in the lateral LHb specifically project to the rostromedial tegmental nucleus (RMTg)43, the primary inhibitory modulator of the VTA44,45. Thus it is interesting to speculate whether these two distinct sources of inhibition converge specifically on RMTg-projecting neurons in the lateral LHb. If this were to be the case, theoretically these neurons would be well poised to reduce excitatory output to the RMTg and reduce the ‘reward aversion’ signalling from the LHb3,7. Further work employing in vivo optogenetic and chemogenetic manipulations in combination with behavioural testing could serve to answer these questions.

In addition to these two populations of PV-positive inhibitory neurons, our results also imply the existence of several physiologically diverse sub-populations of PV-positive neurons within the LHb (Fig. 4), which we assume to be projection neurons. These neurons possess clearly distinct physiological properties from the generic LHb neuron physiology16,32,33, and also appear to be roughly obey sub-nuclear boundaries28. Thus we also speculate that these neurons may have differing projection targets and consequently differing functions. Further studies employing the use of Cre-dependent viral tracers are necessary to delineate the circuitry that these projection neurons comprise.

Possible implications for inhibitory SOM-positive ventral pallidal neurons

We also report SOM-positive inhibitory projection neurons in the ventral pallidum which target the LHb (Fig. 7). Although inhibitory LHb-targeting ventral pallidal neurons have very recently been described23,46, their function has remained largely elusive. Recent work suggests that LHb-projecting inhibitory and excitatory ventral pallidal neurons act oppositely to encode positive and negative motivational valence46. By identifying that SOM acts as a specific molecular marker for these inhibitory neurons, we can thus facilitate further study of this pathway and more thoroughly investigate the role of the ventral pallidum to LHb pathway in controlling motivation.

A need for sub-classifying LHb neurons based on physiological properties and molecular marker expression

Much invaluable information regarding the LHb has come from histological28,40,47 and transcriptomic studies17,30. These works have provided great insight into the organization of the LHb neuronal circuitry with reference to particular protein and gene expression patterns. However, these approaches do not permit these findings to be correlated with neuronal physiology at the single cell level. Indeed, recent reports have highlighted that individual LHb neurons which exhibit ‘bursting’ physiological phenotypes are associated with depression11,12, and thus a greater understanding of LHb neuronal physiology, and how this physiology links to molecular protein expression is of great value.

Within this work, we have approached this problem by investigating the physiological properties of neurons expressing defined molecular markers; indeed a highly relevant approach as recent work indicates that LHb neurons can be divided into 4 major sub-populations based on genetic profile48. This study failed to identify GABAergic neurons; however, a similar transcriptomics study in zebrafish revealed a small population of GABAergic neurons49.

In the case of PV-positive LHb neurons alone, we show that these can be sub-classified into inhibitory neurons, and at least two distinct classes of non-inhibitory neurons (Fig. 4). We also show that SOM-positive LHb neurons represent a physiological sub-class with hybrid properties of LHb and MHb neurons (Fig. 4). What could be the functional importance of these neurons then? We now know that calcium-channel mediated bursting acts as a synchronization mechanism31,50, which in the case of the LHb acts to increase overall output to downstream structures and is thus critical in the pathogenesis of depression11,12. Assuming that the hyperpolarization-induced change in firing modality in these neurons is at least in part mediated by similar low-threshold channels, one could then speculate that this may also be some form of synchronization mechanism specific for these SOM-positive LHb neurons. It has also been proposed that this superior sub-nucleus of the LHb where these SOM-positive neurons are mostly confined (Fig. 2) may act as an area of interaction between the LHb and MHb30. Assuming this to be true, perhaps then these neurons may act to synchronize activity between these two regions and increase overall habenular output, or perhaps they may have different downstream targets altogether from both the LHb and MHb.

As of yet, the specific functions of these PV-positive and SOM-positive neurons remain unclear. However, implementing the use of molecular markers as tools to study particular sub-populations of LHb neurons will most likely facilitate a greater understanding of habenular circuitry at the cellular level.

The functional importance of inhibitory signalling within the LHb

It is well established that the LHb is hyperactive in MDD6,7,8,9,10,11,12, and that reducing LHb hyperexcitability by altering the balance of excitatory and inhibitory signalling has an antidepressant effect7,8,9,20,21. A striking feature of many neurons innervating the LHb is that they co-release GABA and glutamate and hence much work has focussed on the role that these neurons play in modulating LHb excitability8,19,20,21,36. Yet, both local18 and projection15,23 GABAergic neurons are also known to innervate the LHb.

How these GABAergic neurons modulate the LHb in MDD remains to be fully established. Indeed, one could speculate that if excitatory inputs to the LHb are potentiated in MDD7, then GABAergic tone may be reduced. If this were to be the case, then one may also speculate that increasing GABAergic tone should also have an antidepressant effect. For this reason, we believe GABAergic innervation of the LHb to be an important topic of study. Interestingly, our results seem to point to a target selectivity of inhibitory PV-positive neurons at the sub-regional, possibly even single cellular level (Figs. 5E and 7E). We also show that these inhibitory neurons are capable of silencing the spontaneous activity of LHb neurons (Fig. 5G). Therefore, it is reasonable to assume that the role played by these inhibitory neurons in modulating LHb neuronal excitability is not negligible, and hence a greater understanding of this inhibitory control is an exciting concept.

Methodological limitations

The LHb is relatively small structure, and as such it is difficult to confine transgene expression exclusively to this region following stereotaxic AAV injection. This problem was more prominent with our serotype 9 vector, as vectors of serotype 9 are known to have a large spread38. We approached this problem by using AAV 1/2 vectors, which transduced smaller volumes (Fig. 6B) and thus allowed us to confine transfected neuronal somata to the LHb. As we observed inhibitory responses under both conditions, we therefore conclude that there are indeed locally-targeting inhibitory PV-positive neurons within the LHb. However, future studies using paired recordings will have to investigate the properties of synaptically coupled neurons in detail.

Conclusions

We have investigated the mechanisms of inhibitory control within the LHb and have defined three sources of inhibitory input from both locally-targeting PV-positive LHb neurons, and from those in the MDT; and from SOM-positive neurons in the VP. We also report multiple physiologically distinct populations of PV-positive and SOM-positive LHb neurons (Fig. 4), and provide further evidence21,35 for excitatory input to the LHb from PV-positive projection neurons (Fig. 5). These results also indicate that these markers represent broadly diverse populations of neurons on a region-to-region basis and therefore these populations must be validated in each case.

Materials and Methods

Animals

All procedures were approved by the Ethics committee of the University of Strathclyde, Glasgow, and carried out in accordance with the relevant UK legislation (the Animals (Scientific Procedures) Act, 1986). Male and female mice from each strain were used in this work, and all data from both genders were pooled for each experiment. All animals were maintained on a C57BL/6 background, and kept on a 12:12 light/dark cycle under standard group housing conditions with unlimited access to water and normal mouse chow. New-born pups were housed with parents until weaning at P21. To generate transgenic reporter-bearing offspring, transgenic mice of the PV-IRES-Cre (Jax. ID 017320)51 or SOM-IRES-Cre (Jax. ID 018973)52 driver lines were crossed with either Ai32 (Jax. ID 025109)37 or Ai9 (Jax. ID 007909)27 reporter mice driving expression of Channelrhodospin-2 (ChR2) and enhanced yellow fluorescent protein (eYFP), or the enhanced red fluorescent protein variant tdTomato in a Cre-dependent manner, respectively. The resulting offspring strains are hence referred to as: PV-IRES-Cre::Ai32, SOM-IRES-Cre::Ai32, PV-IRES-Cre::Ai9 and SOM-IRES-Cre::Ai9.

Production of recombinant AAV vectors

pAAV-Ef1a-DIO hChR2(E123T/T159C)-EYFP plasmid (a gift from Karl Deisseroth; RRID:Addgene_35509) was packaged into AAVs as described previously53. Briefly, virions containing a 1:1 ratio type 1 and type 2 capsid proteins were produced by transfecting human embryonic kidney (HEK) 293 cells with the rAAV backbone plasmid pAAV-Ef1a-DIO hChR2(E123T/T159C)-EYFP along with AAV1 (pH21), AAV2 (pRV1) and adenovirus helper plasmid pFdelta6 using the calcium phosphate method. 48 hours post transfection, cells were harvested and rAAVs were purified using 1 mL HiTrap heparin columns (GE Healthcare Bio-Sciences, Uppsala, SwedenSigma) and concentrated using Amicon Ultra centrifugal filter devices (Merck Millipore, Tullagreen, Ireland). Infectious rAAV particles (viral titre) were calculated by serially infecting HEK293 cells stably expressing Cre-recombinase and counting GFP-positive cells.

Stereotaxic viral injections

PV-IRES-Cre mice (P31–48) or SOM-IRES-Cre mice (P34–44) were deeply anaesthetized via inhaled isoflurane (5% for induction; 1–2% for maintenance), transferred to a stereotaxic frame (Narishige, Tokyo, Japan) and were subcutaneously injected with the analgesics carprofen (5 mg/kg) in the nape and lidocaine (4 mg/kg) under the scalp. Intracranial injections were made using a glass micropipette pulled using a PC-100 vertical puller (Narishige, Tokyo, Japan). Under aseptic conditions, the skull was exposed and a small burr hole was drilled above the habenula within the left hemisphere. The injection capillary was then advanced, and a Cre-dependent viral vector was injected into the LHb, the medial dorsal thalamic nucleus (MDT) or the ventral pallidum (VP) at a rate of 25 nL/min using a pressure microinjector (Narishige, Tokyo, Japan). Two different vector solutions were used in this study: one containing a pAAV-EF1α-Switch:NLSmRuby2/ChR2(H134R)-EYFP41 HGHpA54 (titre: 2 × 1013 of viral particles containing viral genomes/mL), an adeno-associated virus containing capsid protein 9 which drives expression of ChR2 and eYFP in Cre-expressing cells, or the red fluorescent protein variant mRuby2 in the absence of Cre; or AAV Ef1α-DIO-hChR2-eYFP (in-house titre: 1.6 × 108 infectious particles/mL (as mentioned above) corresponding to about 1.6 × 1012 genomic particles), containing a 1:1 ratio of capsid proteins type 1 and 2. Parameters for LHb-targeted injections were: coordinates (from Bregma, in mm) AP −1.5, ML −0.4, depth 3; volume 50–100 nL; N = 8 PV-IRES-Cre mice, N = 4 SOM-IRES-Cre mice. Parameters for MDT-targeted injections were: AP −1.355, ML −0.75 depth 3.5; volume 200 nL; N = 2 mice or AP −1.155, ML −0.45 depth 3.7; volume 200 nL; N = 4 mice. Note that for the MDT-targeted injection coordinates were adjusted following data acquisition from the first round of experiments. Parameters for VP-targeted injections were: AP 0.145, ML −1.55, depth 5.65; volume 100–250 nL; N = 3 mice. Following injection, the needle was left for 10 minutes to allow the virus to diffuse before being slowly withdrawn. Animals were allowed to recover from anaesthesia on a heat pad. Following completion of surgery, animals were given at least two weeks to allow expression of the virus before acute slice preparation for electrophysiology. Upon completion of electrophysiology, viral spread was assessed by imaging on a Leica SP5 or SP8 confocal microscope.

Acute brain slice preparation

C57BL/6 (N = 5; P21–28), PV-IRES-Cre::Ai32 (N = 19; P23–40), PV-IRES-Cre::Ai9 (N = 8; P23–33), SOM-IRES-Cre::Ai32 (N = 5; P19–37), SOM-IRES-Cre::Ai9 (N = 3; P24–28) or surgically injected PV-IRES-Cre (P46–70 for local LHb injections; P62–70 for MDT injections) or SOM-IRES-Cre mice (P49–63 for local LHb injections; P60–68 for VP injections) were humanely euthanized by cervical dislocation and immediately decapitated, and brains were rapidly removed and transferred to ice-cold oxygenated (95% O2; 5% CO2) sucrose-based artificial cerebro-spinal fluid (ACSF) solution containing (in mM): sucrose 50, NaCl 87, NaHCO3 25, KCl 3, NaH2PO4 1.25, CaCl2 0.5, MgCl2 3, sodium pyruvate 3 and glucose 10. Brains sections containing the habenula were then cut in the coronal plane at 250–300 µm on a Leica VT1200S vibratome (Leica Biosystems, Newcastle-upon-Tyne, UK). In order to ensure slices contained the habenula, the hippocampus was used as a visual guidance due to the easily identifiable structure and immediate proximity to the habenula. Following sectioning, slices were incubated in oxygenated sucrose-based ACSF at 35 °C for 30 minutes, and then incubated for a further 30 minutes at room temperature in ACSF containing (in mM) NaCl 115, NaHCO3 25, KCl 3, NaH2PO4 1.25, CaCl2 2, MgCl2 1, sodium pyruvate 3 and glucose 10. Following the incubation period, slices were stored at room temperature in oxygenated ACSF.

Electrophysiological recordings

Individual slices were transferred to a recording chamber and continually perfused with oxygenated ACSF at a flow rate of 2–3 mL/min, and visualized with a Luigs and Neumann LN-Scope System (Luigs and Neumann, Ratingen, Germany). The habenula is easily identifiable under differential interference contrast microscopy even at low magnification and hence a 4X objective was used to locate the lateral habenular nucleus. A 60X objective was then used to identify suitable cells for whole-cell recordings. In the case of transgenic PV-IRES-Cre::Ai9 or SOM-IRES-Cre::Ai9 slices, TdTomato-expressing cells could be selectively visualized with an Olympus XM10 fluorescent camera (Olympus, Southend-on-Sea, UK) upon photostimulation with a blue LED (pE-300ultra, Cool LED, Andover, UK). Recordings were made with a Multiclamp 700B Amplifier (Molecular Devices, California, USA). Glass micropipettes were filled with a solution containing (in mM) potassium gluconate 125, HEPES 10, KCl 6, EGTA 0.2, MgCl2 2, Na-ATP 2, Na-GTP 0.5, sodium phosphocreatine 5, and with 0.2% biocytin. pH was adjusted to 7.2 with KOH. For spontaneous current measurement experiments, a reduced chloride intracellular solution was used consisting of (in mM) potassium gluconate 140, potassium chloride 2, EGTA 0.2, Hepes 10, NaATP 2, NaGTP 0.5 and sodium phosphocreatine 5.

Once in whole-cell patch mode, the intrinsic properties of LHb neurons were assessed in current-clamp configuration using a stepping protocol consisting of 1 s long injections of increasing current (range: −250–250 pA; step size: 5–50 pA for LHb neurons and −500–1000 pA; step size 100 pA for cortical neurons). Action potential firing pattern was assessed in response to depolarizing current injection, while hyperpolarizing current injection allowed the characterisation of rebound action potential firing of neurons. Resting membrane potential (RMP) was assessed by recording the spontaneous activity of each neuron with no current injection for at least 30 seconds, while membrane input resistance was monitored by injecting a small hyperpolarizing pulse (100 ms; −10 to −100 pA) and measuring the voltage change. Spontaneous currents were observed in voltage clamp at a holding potential of −60 mV. Series resistance was monitored throughout. All neuronal voltage and current signals were low pass-filtered between 2 and 10 kHz and acquired between 10 and 25 kHz using an ITC-18 digitizer interface (HEKA, Pfalz, Germany). The data acquisition software used was Axograph X.

Optogenetic experiments and pharmacology

For optogenetic experiments, acute brain slices were prepared from transgenic PV-IRES-Cre::Ai32 and SOM-IRES-Cre::Ai32 offspring as above in darkness. Whole-cell patch configuration was achieved and neuronal recordings were obtained at varying holding potentials as slices were illuminated with a blue LED pulse (a single pulse of 2 to 200 ms, or a train of 2 ms pulses at 10 to 100 Hz, power 11.5 mW) to elicit postsynaptic events. Where required, SR-95531 (2 µM; henceforth referred to as GABAzine), NBQX (10 µM) or CGP-52432 (10 µM) (all from Tocris, Bristol, UK) were washed into the perfusion bath via the perfusion pump.

Immunohistochemistry and neuronal recovery

Following electrophysiological recordings, slices containing neurons which had been patched and filled with biocytin were processed as previously described55. Briefly, slices were fixed overnight in 4% paraformaldehyde (PFA) dissolved in 0.1 M sodium-based phosphate buffered saline (PBS). After fixation, slices were washed 3 ×5 minutes in 0.1 M PBS, and then incubated for 1 hour in a blocking solution consisting of 5% normal goat serum (NGS) and 1% Triton X-100. Slices were then allowed to incubate on a shaker at room temperature overnight in a primary antibody mixture containing 2.5% NGS and 1% Triton in PBS along with the required primary antibodies. Primary antibodies and dilutions used in this study were: mouse anti-PV (1/4000; Swant, Marly, Switzerland) and rabbit anti-GABA (1/200; Sigma-Aldrich, Dorset, UK). Upon completion of the primary incubation step, slices were washed 2 ×5 minutes in 0.1 M PBS and incubated for 2–3 hours in a secondary antibody cocktail containing the relevant secondary antibodies along with streptavidin (conjugated to Alex Fluor 488 or 647; 1/500 dilution; Life Technologies, Paisley, UK), in order to recover neurons which had been patched and filled with biocytin. The secondary antibodies used in this study were: donkey anti-mouse conjugated to Alexa Fluor 488 (1/500 dilution; Life Technologies, Paisley, UK), goat anti-rabbit conjugated to Alexa Fluor 555, 633 or 647 (1/500 dilution; Life Technologies, Paisley, UK) and supplemented with 1% Triton in PBS. Fluorophores excitable at differing wavelengths were implicated depending on whether the slice expressed YFP (Ai32 animals) or TdTomato (Ai9 animals) to minimize crosstalk. Where only neuronal recovery was required, slices were blocked as above and incubated in a solution containing streptavidin supplemented with 1% Triton in PBS. After secondary antibody incubation, slices were washed for 3 ×5 minutes in 0.1 M PBS and mounted on glass slides using Vectashield medium (containing DAPI as required, Vector Labs, Peterborough, UK) and cover-slipped.

Double fluorescentin situhybridization

Fresh, unfixed C57BL/6 mouse brains (N = 3) were embedded in tissue freezing medium (Leica Biosystems, Richmond, UK), frozen on a dry ice ethanol bath, sectioned at 20 μm on a cryostat, and then mounted onto Polysine Adhesion Slides (Thermo Fisher Scientific, Waltham, USA). RNA probe hybridization and subsequent washes were performed as described previously56. Fluorescein-labeled probes were detected using peroxidase-conjugated anti-fluorescein antibodies (Roche Diagnostics, Mannheim, Germany). Peroxidase activity was detected using Cy3-tyramide conjugate. After the detection, peroxidase activity was blocked by 100 mM glycine-HCl (pH 2.0) solution containing 0.1% Tween. Sections were washed with 0.1 M Tris-HCl buffer (pH 7.5) containing 0.15 M NaCl and 0.05% Tween and then blocked by 10% goat serum and 1% blocking reagent powder (Roche Diagnostics, Mannheim, Germany) in 0.1 M Tris-HCl (pH 7.5) with 0.15 M NaCl. Dioxygenin (DIG)-labeled probes were detected using peroxidase-conjugated anti-DIG antibodies (Roche Diagnostics, Mannheim, Germany). Peroxidase activity was detected using fluorescein-tyramide conjugate. Sections were counterstained with 4′,6-diamidino-2-phenylindole (DAPI) (Sigma-Aldrich, Munich, Germany).

Tyramide conjugate synthesis

Fluorescein- and Cy3-tyramide conjugates were synthetized as described previously57. Briefly, the succinimidyl esters of fluorescein (Thermo Fisher Scientific, Schwerte, Germany) and Cy3 (GE Healthcare, Little Chalfont, United Kingdom) were coupled to tyramine (Sigma-Aldrich, Munich, Germany) in dimethylformamide (Carl Roth, Karlsruhe, Germany) adjusted to a pH of 7.0–8.0 with triethylamine (Sigma-Aldrich, Munich, Germany).

Details of probes designed forin situhybridization

DNA fragments used as template for RNA probe transcription were PCR-amplified from C57BL/6 mouse genomic DNA using the following primers:

VGAT:                                                                                               GACCTCGAGCTACCTGGGGTTGTTCCTCA, AATTAACCCTCACTAAAGGGACTAGTCGAAGTGTGGCACGTAGATG

VGLUT2:                                                                                               GACGAATTCTATTCGTTGGACCCATCACC,gacAATTAACCCTCACTAAAGGGCGGCCGCAGAAATTGCAATCCCCAAAC

PV:

Exon3:                                                                                               GACGAATTCCCTCTCCCCTGTCCTTCTTT, AATTAACCCTCACTAAAGGGCGGCCGCaTGGGAACTTTGGGTGCTATC

Exon4:                                                                                               GACGAATTCAGGTTCTGCCTGTGACCTTG, AATTAACCCTCACTAAAGGGCGGCCGCtAAGCTTTGACAGCCGCATAC

Exon5:                                                                                               GACGAATTCCTCCACTCTGGTGGCTGAA, AATTAACCCTCACTAAAGGGCGGCCGCTTTCTCTTTTCAGGTATTTTATCACA

Amplified DNA fragments were cloned into the pBSK backbone and transcribed with T3 RNA polymerase using DIG RNA Labeling Mix and Fluorescein RNA Labeling Mix according to the manufacturer protocol (all from Roche Diagnostics, Mannheim, Germany).

Intracardial perfusion and serial sectioning

To prepare tissue for serial sectioning, PV-IRES-Cre::Ai32 (N = 2; P31–35), SOM-IRES-Cre::Ai9 (N = 2; P28), or C57BL/6 (N = 3; P21–22) mice were terminally anaesthetized by subcutaneous injection with an overdose cocktail of 50% lidocaine and 50% euthatal. Once anaesthetized sufficiently to be non-responsive to noxious tail and toe pinch stimuli, mice were perfused through the left ventricle with 0.1 M PBS followed by perfusion with 4% PFA dissolved in PBS. Brains were then removed and fixed overnight in 4% PFA in PBS, after which they were cryoprotected in a solution containing 30% (w/v) sucrose in PBS for storage until required for serial sectioning.

For these experiments, brains were embedded in OCT compound (VWR, Leicestershire, UK) and sectioned on a Leica SM2010 R microtome (Leica Biosystems, Newcastle-upon-Tyne, UK) at 60–80 µm. Upon completion of sectioning, slices were washed 3 ×5 min in 0.1 M PBS. Where further staining was required, this was carried out as above (see immunohistochemistry and neuronal recovery), however, 0.3% Triton X-100 was used in place of 1% and slices were incubated overnight in primary antibody cocktails to minimize tissue damage. Slices were then mounted using Vectashield medium (Vector Labs, Peterborough, UK) and cover-slipped.

Image acquisition and neuronal reconstructions

For immunohistochemistry-stained sections and biocytin-filled neurons; mounted sections were scanned on either a Leica SP5 or SP8 confocal microscope, imaging z-stacks of each slice at 2–4 µm steps. Confocal laser excitation wavelengths (in nm) were 405, 488, 514, 552 and 503. Objectives used were 10X (dry), 20X (oil immersion), 40X (oil immersion) and 63X (oil immersion) for Leica SP5, or 10X (dry), 20X (dry) and 63X (oil immersion) for SP8. A zoom of up to 2X was applied as required to occasionally visualize soma in enhanced detail. Sections were scanned to ensure that all visible streptavidin-stained cells and their neurites were included in the z-stack. 3D reconstructions of neurons were carried out using NeuTube 3D reconstruction software, with processes which branched extensively, had a thicker diameter and tapered being designated dendrites, while processes with a thinner diameter which remained constant along the length of the process being designated an axon58.

For in situ hybridizations; the two-dimensional overview images were acquired with an Axio Examiner microscope with Axiocam 506 camera and LED Light Source Colibri 7 system using a Plan-Apochromat 20×/0.8 M27 objective (all from Carl Zeiss, Jena, Germany). The illumination wavelengths of 450–488 nm, 540–570 nm and 370–400 nm and emission filter wavelengths of 500–550 nm, 570–640 nm and 420–470 nm were used for Fluorescein, Cy3 and DAPI, respectively. The light source intensity was set to 40%, 50% and 20% for Fluorescein, Cy3 and DAPI, respectively. The three-dimensional images were acquired with a Zeiss LSM880 confocal laser scanning microscope using a Plan-Apochromat 20×/0.8 M27 objective (all from Carl Zeiss, Jena, Germany) with a 2X optical zoom. Laser wavelengths of 561 nm, 488 nm and 405 nm were used for Cy3, Fluorescein and DAPI, respectively. Stacks of 6 to 12 optical slices were captured with a z-step size of 0.785 μm.

Data analysis

Analysis of electrophysiological recordings was carried out using Axograph X. Passive intrinsic properties were calculated as described above, while active intrinsic properties (action potential initial frequency, amplitude, rise-time and half-width) were calculated by subtracting the baseline and then using the event detection feature to analysis the first action potential elicited in response to a 50 pA depolarizing pulse. Neuronal spontaneous activity was classified as either bursting (a clearly distinctive behaviour16), tonic or silent (where spontaneous action potential frequency was < 1 Hz). For optogenetically evoked events, peak size was measured at various holding potentials. For spontaneous current measurements, representative example traces of postsynaptic currents were first generated and currents were detected and measured using the event detection feature.

Image analysis was carried out using ImageJ. Confocal z-stacks were compressed onto a single image and brightness and contrast were occasionally adjusted to enhance cellular visualization. Cell counts were quantified using the cell-counter plugin. For these experiments, serial sections containing the whole habenula were imaged and analysed from one animal for each strain, and for remaining animals every second or third section was imaged and analysed to allow quantification of markers with fair representation of the habenular sub-nuclei. Images were then transferred to PowerPoint (Microsoft), where cells of interest were marked.

Graphs were generated and statistical analysis was performed using GraphPad Prism 5 (California, USA). Statistical tests used were: two-tailed unpaired t-test for single comparisons of passive physiological properties; one-way ANOVA with Tukey’s multiple comparison test for comparison of physiological properties between multiple groups; two-way ANOVA with Bonferroni’s multiple comparison for assessing relationship between input current and action potential discharge (f-I curves), or Fisher’s exact test. Once graphs were generated, they were transferred to PowerPoint 2013 for formatting and assembly into figures. Statistical significance thresholds for all tests were: *p < 0.05; **p < 0.01; and ***p < 0.001.