GABAergic Control of Nigrostriatal and Mesolimbic Dopamine in the Rat Brain

Purpose: The present study assessed the effects of the GABAA receptor (R) agonist muscimol (MUS), and the GABAAR antagonist bicuculline (BIC) on neocortical and subcortical radioligand binding to dopamine D2/3Rs in relation to motor and exploratory behaviors in the rat. Methods: D2/3R binding was measured with small animal SPECT in baseline and after challenge with either 1 mg/kg MUS or 1 mg/kg BIC, using [123I]IBZM as radioligand. Motor/exploratory behaviors were assessed for 30 min in an open field prior to radioligand administration. Anatomical information was gained with a dedicated small animal MRI tomograph. Based on the Paxinos rat brain atlas, regions of interest were defined on SPECT-MRI overlays. Estimations of the binding potentials in baseline and after challenges were obtained by computing ratios of the specifically bound compartments to the cerebellar reference region. Results: After MUS, D2/3R binding was significantly reduced in caudateputamen, nucleus accumbens, thalamus, substania nigra/ventral tegmental area, and posterior hippocampus relative to baseline (0.005 ≤ p ≤ 0.012). In all these areas, except for the thalamus, D2/3R binding was negatively correlated with grooming in the first half and positively correlated with various motor/exploratory behaviors in the second half of the testing session. After BIC, D2/3R binding was significantly elevated in caudateputamen (p = 0.022) and thalamus (p = 0.047) relative to baseline. D2/3R binding in caudateputamen and thalamus was correlated negatively with sitting duration and sitting frequency and positively with motor/exploratory behaviors in the first half of the testing time. Conclusions: Findings indicate direct GABAergic control over nigrostriatal and mesolimbic dopamine levels in relation to behavioral action. This may be of relevance for neuropsychiatric conditions such as anxiety disorder and schizophrenia, which are characterized by both dopaminergic and GABAergic dysfunction.

The DAergic system is under inhibitory γ-amino butyric acidergic (GABAergic) control (Precht and Yoshida, 1971). From the CP, GABAergic projections run either directly or via the external part of the globus pallidus (GPe) and the subthalamic nucleus (STN) to the internal part of the globus pallidus (GPi) and to the pars reticulata of the SN (SNr). From GPi and SNr, further GABAergic efferents project to THAL (ventral anterior, ventrolateral, dorsomedial, and centromedian nucleus), pedunculopontine nucleus, inferior and superior colliculus, and periaqueducatal gray (Kuo and Carpenter, 1973;Graybiel and Ragsdale, 1979;Herkenham, 1979;Coimbra and Brandao, 1993). In addition, the CP is inhibited by GABAergic mircrocircuits, which are formed by fast-spiking interneurons and collaterals of descending projection neurons (Groves, 1983). The NAC receives GABAergic input from the PFC (Lee et al., 2014), and, in turn, sends GABAergic efferents back to the VTA as well as to the ventral GP and from the GP to the dorsomedial thalamic nucleus (Ueki et al., 1977;Yamamoto et al., 1983). DA neurons serve to modulate the GABAergic system: in the direct pathway, DA exerts an inhibitory effect on those GABAergic neurons, which project to the GPe, leading to a net disinhibition of the THAL, whereas in the indirect pathway, DA stimulates GABAergic projections to GPi and SNr with a net inhibitory effect on the THAL. The THAL, in turn, sends GLUergic efferents to MC, FC, PFC, and CING.
We have previously shown that MUS but not BIC reduced radioligand binding to the neostriatal D 2 R-like subtype relative to baseline (BAS) and to vehicle. Moreover, motor/exploratory behaviors were diminished after MUS, but elevated after BIC relative to vehicle (Nikolaus et al., 2017). These effects of MUS and BIC prompted us to assess binding to the D 2 R-like subtype separately in CP and NAC as well as in other cortical and subcortical regions including THAL, SN/VTA, FC, MC, parietal cortex (PC), anterodorsal hippocampus (aHIPP), and posterior hippocampus (pHIPP). In the previous study, for each rat, small volumes of interest (VOIs) had been defined around the neostriatal hot spots of maximum radioactivity accumulations (Nikolaus et al., 2017). This mode of analysis has proven sufficient for a comparatively large region with a high density of D 2 Rlike binding sites such as the striatum. However, it may be not valid to determine radioligand binding in smaller regions such as SN/VTA, or in cortical or limbic regions with lower amounts of D 2 R-like binding sites. Therefore, in the present investigation, images of D 2 R-like binding were coregistered with morphological images obtained with a dedicated small animal MRI tomograph, which allowed the transference of the predefined cortical and subcortical VOIs of the Paxinos standard rat brain MRI (Schiffer et al., 2006) to the functional SPECT images and, thus, permitted the first-time in vivo imaging analysis of radioligand binding to D 2 R-like binding sites in relevant regions of the rat nigrostriatal and mesolimbic system such as NAC, THAL, SN/VTA, FC, MC, PC, aHIPP, and pHIPP. Additionally, the relation between D 2 R-like binding in these regions and motor/exploratory behaviors was assessed with both correlation analysis and cluster analysis.

Animals
Imaging studies of D 2 R-like binding sites were conducted on 32 adult male Wistar rats (ZETT, Heinrich-Heine University, Düsseldorf, Germany), weighing 429 ± 42 g [mean ± standard deviation (SD)]. Thereby, animals underwent one SPECT measurement in BAS and one SPECT measurement plus behavioral testing after challenge with MUS or BIC (n = 16, respectively). Morphological images were obtained on eight further male Wistar rats (ZETT, Heinrich-Heine University, Düsseldorf, Germany; weight: 434 ± 37 g). Findings on striatal D 2 R-like binding and behavioral data were previously published, with binding data, however, having been obtained with a different mode of analysis (see above; Nikolaus et al., 2017). Rats were maintained in standard macrolon cages (590 × 380 × 200 mm; three animals per cage) in a climate cabinet (Scantainer, Scanbur BK, Karslunde, Denmark; temperature: 20 • C; air humidity; 70%) with an artificial ligh-dark cycle (lights on at 6:00 a.m., lights off at 6:00 p.m.) and food and water freely available. The study was carried out in accordance with the recommendations of the "Principles of laboratory animal care" (NIH publication No. 86-23, revised 1985) and the German Law on the Protection of Animals. The protocol was approved by the regional authority (Landesamt für Natur, Umwelt und Verbraucherschutz, Nordland-Westfalen, Recklinghausen, Germany).

MRI Studies
Upon anesthization with ketaminehydrochloride (dose: 0.45 ml/kg i.p., concentration: 100 mg/ml) and xylazinehydrochloride (dose: 0.2 ml/kg i.p., concentration: 0.02 mg/ml), rat heads were scanned with a dedicated small animal MRI tomograph (MRS3000 Pre-clinical MRI, 3.0 T, MR Solutions, Guildford, UK). A rat body volume coil with an inner diameter of 54 mm was used to transmit the radio frequency pulses and to receive the MR signals. A standard gradient echo pilot scan in three orthogonal directions was used for the positioning of a 3D fast low angle shot (FLASH) MR sequence to acquire high-resolution anatomical images (Haase et al., 1986). The measurement parameters for the 3D FLASH sequence were as follows: image matrix: 192 × 192 × 96, interpolated by zero-filling before reconstrution to a 256 × 256 × 128 matrix in coronal slice orientation; field of view: 64 × 64 × 44 mm; spatial resolution: 0.25 × 0.25 × 0.69 mm; repetition time: 30 ms; echo time: 4.87 ms; excitation flip angle: 30 • ; total acquisition time: 553 s.
Evaluation of SPECT Imaging Studies D 2/3 R imaging data were evaluated with PMOD (version 3.5, PMOD Technologies Ltd., Zürich, Switzerland). Firstly, the SPECT image of each rat was coregistered with a MR image of an animal of the same weight. Then, the respective MR image was coregistered with the Paxinos standard rat brain MRI (Schiffer et al., 2006) provided by PMOD. The necessary mathematical transformations were saved. The SPECT image as coregistered with the fitting MRI was imported using these transformations, which allowed creation of an overlay with the Paxinos standard rat brain MRI. On these overlays, the following volumes of interest (VOIs) were defined: CP, NAC, THAL, SN/VTA, FC, MC, PC, aHIPP, and pHIPP. The tomographic resolution of the employed SPECT camera amounts to 2.8 and 3.4 mm for 99m Tc and 123 I, respectively (Schramm et al., 2000). According to the rat brain atlas, all these regions have maximum craniocaudal (CC) and one-sided mediolateral (ML) and dorsoventral (DV; vertical or oblique) dimensions in the range of or beyond the spatial resolution of  (Paxinos and Watson, 2014). Moreover, in autoradiographic studies performed with a variety of radioligands including [ 123 I]IBZM, these regions have been shown to express D 2 R-like binding sites (Bouthenet et al., 1987;Verhoeff et al., 1991). Since [ 123 I]IBZM accumulation in the cerebellum (CER) is non-specific, the CER was used as reference region (REF; see Supplementary Figure). Estimations of regional binding potentials (BPs) for BAS and challenges were obtained according to the simplified reference tissue model by computing ratios of radioactvity counts obtained in the specifically-bound compartments (CP, NAC, THAL, SN/VTA, FC, MC, PC, aHIPP, and pHIPP) to radioactivity counts in the CER (Ichise et al., 2001).

Statistical Analysis
Distributions of both regional BPs and behavioral data were tested for normality with the non-parametric Kolmogorov-Smirnov test (α = 0.05). Neither in BAS, nor after MUS or BIC, regional BPs were uniformly normally distributed (0.021 ≤ p ≤ 0.20). This also held for behavioral parameters after both MUS and BIC (0.0001 ≤ p ≤ 0.20).
Medians with 25-/75-percentiles were computed for regional BPs. Regional BPs were compared between pre-treatment conditions (BAS MUS vs. MUS, BAS BIC vs. BIC) with the Wilcoxon signed rank test for paired samples (two-tailed, α = 0.05). Moreover, percentual differences of BPs relative to BAS were computed for both MUS and BIC. No corrections of the alpha value were made for multiple comparisons. Calculations were performed with IBM SPSS Statistics 23 (IBM SPSS Software Germany, Ehningen, Germany).
In addition, cluster analyses were performed for regional BPs and behavioral data (traveled distance and durations and frequencies of ambulation, sitting, rearing, head-shoulder motility, and grooming) in the individual time frames (1-5, 6-10, 11-15, 16-20, 21-25, and 26-30 min) after MUS and after BIC. The individual data were standardized using the Ztransformation [(individual value -variable mean)/SD from the variable mean]. The number of clusters and the centroid mean of each variable in them were determined with an Xmeans algorithm (Pelleg and Moore, 2000). Cluster analyses were calculated with Rapid Miner (version 5.3., Rapid-I GmbH, Dortmund, Germany). The individual values of the behavioral parameters, which entered correlation analysis and cluster analysis, are given for each animal in the Supplementary Table.

D 2/3 R Binding
In Figure 1, characteristic coronal images of regional [ 123 I]IBZM accumulations in the BAS condition and after challenge with MUS and BIC, respectively, are presented at different positions from Bregma (Paxinos and Watson, 2014). SPECT images in the conditions BAS MUS and MUS were obtained in the same rat. This also holds for BAS BIC and BIC.
After BIC (Figure 3), the thalamic BP was significantly elevated (+17%, p = 0.047) compared to BAS BIC . Likewise, the post-challenge BP in the CP was significantly increased relative to the BP in BAS (+8%, p = 0.022; Figure 3). There were no differences between BAS BIC and BIC in NAC, SN/VTA, FC, MC, PC, aHIPP, pHIPP, and CER (0.130 ≤ p ≤ 0.678). Correlations Between Regional D 2/3 R Binding and Behaviors Table 1 shows the significant (0.0001 ≤ p ≤ 0.05) positive and negative correlations between regional BPs and behavioral parameters in the individual time frames after treatment with MUS or BIC.
After MUS, lower D 2/3 R binding in CP, NAC, SN/VTA, FC, MC, aHIPP, and pHIPP was associated with a reduction of motor/exploratory parameters in the second half of the testing session and an increase of sitting and grooming primarily between 11 and 15 min.
For BIC, the correlation analysis showed an association between lower D 2/3 R binding in CP, NAC, THAL, MC, and aHIPP and more sitting and head-shoulder motility as well as less traveled distance and rearing primarily during the first 5 min of the testing session.

MUS
The cluster analysis (Table 2) of regional D 2/3 R binding data and behavioral variables in the individual time frames yielded two clusters (cluster 1: n = 7, cluster 2: n = 9). Relative to the second cluster, the first one was characterized by lower centroid means of D 2/3 R binding in CP, NAC, THAL, and FC and higher centroid means of D 2/3 R binding in SN/VTA, MC, PC, aHIPP, and pHIPP. Moreover, relative to the second cluster, centroid means in the first one were lower for (1) sitting duration and frequency in the fifth time frame, (2) traveled distance, ambulation duration, ambulation frequency, and frequency of head-shoulder motility in the first, second, third, fourth, and sixth time frame, and (3) rearing duration, rearing frequency, and duration of head-shoulder motility in all time frames. In contrast, centroid means in the first cluster were higher for (1) traveled distance, ambulation duration, ambulation frequency, and frequency of head-shoulder motility in the fifth time frame, (2) sitting duration and sitting freuency in the first, second, third, fourth, and sixth time frame, and (3) both grooming duration and grooming frequency in all time frames.

BIC
For treatment with BIC, the cluster analysis (Table 2) of regional D 2/3 R binding data and behavioral variables in the individual time frames also yielded two clusters (cluster 1: n = 6, cluster 2: n = 10). Relative to the second cluster, the first one was characterized by lower centroid means of D 2/3 R binding in CP, NAC, THAL, and MC and higher centroid means of D 2/3 R binding in SN/VTA, FC, aHIPP, and pHIPP. Moreover, relative to the second cluster, centroid means of the first one were lower for (1) sitting frequency in the first and fourth time frame, (2) duration of head shoulder motility in the first and second time frame, and (3) sitting duration in the first, second, fourth, fifth, and sixth time frame, but higher for (1) sitting duration in the third time frame, (2) sitting frequency in the second, third, fifth, and sixth time frame, (3) grooming duration and grooming frequency in the fourth, fifth, and sixth time frame, (4) traveled distance, ambulation duration, and ambulation frequency, duration of head-shoulder motility in the second, third, fourth, fifth, and sixth time frame, and (5) rearing duration, rearing frequency, and frequency of head-shoulder motility in all time frames.
Former investigations have shown that MUS injections into VTA (Klitenick et al., 1992) and NAC (Yoshida et al., 1997;Aono et al., 2008) elevated DA concentrations in these regions. Moreover, MUS increased neostriatal and intranigral DA levels when applied into the SN  and into the GP Abercrombie, 2002, 2003). In a precedent study, we had presented first evidence that MUS also after systemic application elevated neostriatal DA levels, leading to a competition between endogenous DA and the exogenous radioligand and a subsequent reduction of radioligand binding to the neostriatal D 2 R-like subtype (Nikolaus et al., 2017). The present analysis of these imaging data, using MRI overlays additionally yielded evidence of diminished D 2/3 R binding in NAC, THAL, SN/VTA, and pHIPP, reflecting increased availability of DA also in these regions.
Former in vivo microdialysis studies have shown elevated DA levels in CP (Smolders et al., 1995 andNAC (Yan, 1999;Aono et al., 2008) upon administration of BIC into these regions. Furthermore, BIC augmented DA concentrations in the CP, when injected into PFC (Karreman and Moghaddam, 1996;Matsumoto et al., 2003) and SN Westerink et al., 1992) and in the NAC, when injected into the VTA (Ikemoto et al., 1997). This is contrasted by the results of our in vivo imaging study, which revealed increased D 2/3 R binding in CP and THAL upon systemic challenge with BIC, indicating decrements of synaptic DA in these regions. This disagreement, firstly, may be due to the applied dose, which, in the present study (2.7 mM) was considerably higher than in the precedent investigations (10-100 µM; Santiago and Westerink, 1992;Westerink et al., 1992;Smolders et al., 1995;Karreman and Moghaddam, 1996;Ikemoto et al., 1997;Yan, 1999;Matsumoto et al., 2003;Aono et al., 2008). A second reason may be the fact, that in the present study the pharmacological challenges were administered systemically in contrast to the localized injection in the other investigations. It may be inferred that the simultaneous targeting of all regions expressing GABA A R binding sites may have influenced the GABAergic (and DAergic) afferents and efferents, regulating nigral, tegmental, striatal, thalamic, and cortical function differently from the localized approach in the in vivo microdialysis studies.
In the present study, systemic application of the GABA A R agonist MUS reduced D 2/3 R binding in striatum (CP, NAC), THAL, SN/VTA, and pHIPP, reflecting increased availability of DA in these regions. The D 2 R has two interconvertible binding states for DA, which are referred to as high-affinity (Gprotein-coupled) and low-affinity (G-protein-uncoupled) state (De Lean et al., 1982). To our knowledge, the individual affinities of [ 123 I]IBZM for the high-and the low-affinity D 2 R configuration have not yet been determined. However, studies on [ 11 C]raclopride-a further D 2 R antagonistic benzamide analog widely used for the assessment of D 2 R binding and competion with endogenous DA (Laruelle, 2000)-have shown similar affinities for both types of D 2 Rs (Seneca et al., 2006). From this may be inferred that also [ 123 I]IBZM binds to D 2 R-like binding sites in both configurations and that the regional BPs obtained in the present investigation may be considered to reflect the regional densities of D 2/3 Rs as such, irrespective of the individual contributions of either affinity state. This not only holds for the BPs in BAS, but also for BPs obtained after challenge with either MUS or BIC, and exempts us from the necessity to differentiate between the individual effects exerted by D 2/3 Rs in the highand in the low-affinity configuration. With this simplification in mind, the following actions may be hypothesized to occur in the individual brain regions: firstly, within the neostriatal microcircuits (Groves, 1983), the GABA A R agonist action of MUS can be assumed to lead to a decline of DA efflux. In the DAergic system, DA concentrations undergo regulation by autoreceptors of the D 2 R subtype, which are situated at the presynaptic terminal (Langer, 1974). Consequently, the decrement of striatal DA efflux elicited by MUS is likely to diminish DA binding to presynaptic D 2 autoreceptors, leading to a reduction of feedback inhibition, subsequent elevation of DA efflux and the observed reduction of radioligand binding to the D 2/3 R in the CP.
The increased GABA A R agonistic action in the CP can be inferred to augment the inhibition of the SN via the direct pathway. The same effect is exerted by GABAergic thalamonigral and DAergic striatonigral efferents (the latter via inhibitory D 2/3 R binding sites). The consequence of both GABAergic and DAergic inhibition is the decrease of nigral DA levels and the subsequent reduction of inhibitory D 2 autoreceptor action. Via the indirect pathway, the SN is disinhibited, which-together with the increased input of excitatory striatonigral efferents (via D 1 R binding sites) and the mentioned decline of D 2 autoreceptor action-can be hypothesized to cause the net increase of nigral DA reflected by the reduction of D 2/3 R under the present experimental conditions. The increased GABA A R agonistic action in the indirect pathway disinhibits the THAL. Furthermore, elevated DAergic input from the CP both disinhibits (via D 2/3 heteroreceptors) FIGURE 3 | Binding potentials in baseline (white) and after challenge with 1 mg/kg bicuculline (gray). Rendered are medians and 25-/75-(boxes) and 9-/95-quartiles (whiskers). The circles represent the individual animals. For significant between-group differences, the respective p-values are given (Wilcoxon signed rank test for paired samples, two-tailed, α = 0.05). and inhibits (via D 1 heteroreceptors) the THAL in the direct and indirect pathway, respectively (Bolam et al., 2000). As indicated by the observed reduction of thalamic D 2/3 R binding, these effects altogether appear to result in a net increase of thalamic DA.
From the THAL, GLUergic efferents run to the neocortex, which in turn sends GLUergic projections back to the CP (Rouse et al., 2000). Thalamic excitation of the neocortex can be expected to enhance the excitatory input exerted by the descending corticostriatal fibers. The stimulation of striatal N-methyl-Daspartate (NMDA) receptors (Clow and Jhamandas, 1989) likely joins the action of both nigral and striatal D 2 autoreceptors in augmenting striatal DA levels in response to the GABA A R agonist action of MUS.
In line with the action of MUS on the CP, leading to an elevated inhibition of the SN, the action of MUS on the NAC may be conceived to enhance the inhibition of the VTA, resulting in a reduction of DA release in the NAC. Again, the diminished DA efflux can be hypothesized to decrease DA binding to presynaptic D 2 autoreceptors, leading to a decrease of feedback inhibition, subsequent elevation of DA efflux and the observed reduction of D 2/3 R binding in the NAC.
Moreover, elevated GABA A R agonist action in the NAC can be assumed to increase the inhibition of the GP, leading to disinhibition of the THAL. This likely adds to the increased GABAergic and DAergic action in the indirect and direct pathways, respectively, causing a net increase of thalamic DA and the observed reduction of thalamic D 2/3 R binding.
The NAC receives inhibitory GABAergic and DAergic efferents from the PFC (Lee et al., 2014) and the THAL (Hara et al., 1989), respectively. The increase of DA in the THAL probably acts jointly with the inhibition of VTA and NAC to reduce DA levels in the latter region, leading to the decrease of feedback inhibition and subsequent elevation of DA efflux hypothesized to underlie the observed decline of D 2/3 R binding.
The HIPP receives DAergic neurons originating in the VTA (Nazari-Serenjeh et al., 2011). Hence, the MUS-induced inhibition of the VTA is likely to diminish the DA release in the HIPP, again resulting in reduced feedback inhibition, subsequent enhancement of DA release, and the observed decline of radioligand binding to the D 2/3 R in the pHIPP. The HIPP sends GLUergic efferents to the NAC (Nazari-Serenjeh et al., 2011), with the stimulation of NMDA receptors likely joining the action of D 2 autoreceptors in augmenting DA levels in the NAC and subsequently also in the THAL.
In the present study, systemic application of the GABA A R agonist BIC increased D 2/3 R binding in CP and THAL, reflecting decreased availability of DA in these regions. Firstly, it may be hypothesized that, within the microcircuits in the CP, the GABA A R antagonistic action elicited an increase of DA efflux,   which is compensated by an enhancement of feedback inhibition and subsequent reduction of DA levels, as indicated in the present study by the elevation of D 2/3 R binding in the CP.
In vivo imaging studies did not show alterations of D 2/3 R binding (and DA) in SN/VTA and NAC. Therefore, actions of BIC on thalamic DA via the NAC can be dismissed. It may be inferred, however, that the increased GABA A R antagonistic action in the indirect pathway leads to an inhibition of the THAL. Furthermore, also the BIC-induced reduction of DAergic input from the CP decreases both thalamic inhibition (via D 2/3 heteroreceptors) and excitation (via D 1 heteroreceptors) in the direct and indirect pathway, respectively. Altogether, these effects may be assumed to result in a net decrease of thalamic DA and the observed elevation of thalamic D 2/3 R binding.

Rat Behavior
Imaging and behavioral data obtained after MUS can be grouped into two clusters with the first one characterized by lower D 2/3 R binding in CP, NAC, THAL, and FC and higher D 2/3 R binding in SN/VTA, MC, PC, aHIPP, and pHIPP and the second one displaying the opposite. Moreover, the animals of the first cluster showed less motor/exploratory behaviors and more sitting and grooming, whereas the animals of the second cluster exhibited more motor/exploratory activity and less sitting and grooming. Hence, the cluster analysis implies an inverse relationship between DA levels in the nigrostriatal/mesolimbic system and motor/exploratory activity after GABA A R agonistic treatment. In addition, correlation analysis revealed that, after MUS, lower D 2/3 R binding (and higher DA) throughout the entire DAergic system (CP, NAC, SN/VTA, FC, MC, aHIPP, and pHIPP) was associated with a reduction of motor/exploratory parameters in the second half of the testing time and an increase of sitting and grooming primarily from 11 to 15 min. From this, we infer that, after systemic MUS, DA concentrations in the nigrostriatal and mesolimbic system started to rise around 11 min post-injection. Furthermore, the reductions of D 2/3 R binding in the imaging studies indicate that the synaptic DA levels were still elevated in 75-135 min after MUS challenge. 2 | Centroid means of D 2/3 R binding in CP, caudateputamen; NAC, nucleus accumbens; THAL, thalamus; SN/VTA, substantia nigra/ventral tegmental area; FC, frontal cortex; MC, motor cortex; PC, parietal cortex; aHIPP, anterodorsal hippocampus; pHIPP, posterior hippocampus; as well as traveled distance and ambulation and frequency of ambulation, sitting, rearing, head-shoulder motility, and grooming in the individual time frames (min 1-5, 6-10, 11-15, 16-20, 21-25, and 26-30) in rats pre-treated with 1 mg/kg muscimol (MUS) or 1 mg/kg biculline (BIC) after Z-transformation of the individual data.  For imaging and behavioral data obtained after BIC, the cluster analysis yielded two clusters with the first one characterized by lower D 2/3 R binding in CP, NAC, THAL, and MC and higher D 2/3 R binding in SN/VTA, FC, aHIPP, and pHIPP and the second one displaying the opposite. Moreover, the animals of the first cluster showed less sitting, but more motor/exploratory behaviors as well as grooming throughout the testing time, whereas those of the second cluster exhibited more sitting and less motor/exploratory activity and grooming. Thus, interestingly, after GABA A R antagonistic treatment, leading to decreased DA in CP and THAL but normal DA in the other brain regions, a direct relationship emerges between DA concentration in the nigrostriatal/mesolimbic system and motor/exploratory activity.

MUS
Yet, the correlation analysis revealed an association between lower D 2/3 R binding (and higher DA) throughout the nigrostriatal/mesolimbic system (CP, NAC, THAL, MC, and aHIPP) and more sitting and head-shoulder motility as well as less traveled distance, and rearing primarily during the first 5 min of testing time. This, for one, indicates that the GABA A R antagonist elicited an almost immediate decline of synaptic DA concentrations, which, in NAC and THAL, was still visible at the time of in vivo imaging studies, and which is reflected by the immediate reduction of motor/exploratory behaviors. Secondly, however, it seems that the association between nigrostriatal/mesolimbic DA and motor/exploratory behaviors is different for generally high (as after L- DOPA Nikolaus et al., 2016 or MUS) and generally low or normal levels of DA (as after BIC) with the former inversely and the latter directly related to motor/exploratory activity. Further in vivo imaging studies of D 2/3 R binding after increasing doses of L-DOPA, MUS, and BIC are required to gain more information on this regulation mechanism, which is likely to involve tight region-specific D 2 autoreceptor action.

Appraisal
Against the present findings the objection can be raised that D 2/3 R SPECT and MR images were not obtained on the same rats, which may have led to misalignments not only between D 2/3 R SPECT and MRI of individual animals, but also between D 2/3 R SPECT and the standard Paxinos rat brain MRI. However, the BPs obtained for the CP in the present study (BAS MUS : 2.845, MUS: 2.37) were in agreement with the previous mean neostriatal equilibrium ratios (V 3 ") of 1.818 (BAS MUS ) and 1.397 (MUS), corresponding to BPs of 2.818 and 2.397, respectively (Nikolaus et al., 2017). Also after BIC, the BPs obtained for the CP (BAS BIC : 2.380, BIC: 2.575) were consistent with the fromer findings (BAS BIC : V 3 " = 1.532, BP = 2.532; BIC: V 3 " = 1.646, BP = 2.646; Nikolaus et al., 2017), which argues in favor of the present method. Although also the previous data had indicated an increase of neostriatal D 2/3 R binding after BIC, this difference had not reached statistical significance. It may be argued, thus, that the MRI-based mode of analysis is superior to the former method of defining a small VOI around the neostriatal hot spot of maximum radioactivity accumulation, since it allows the definition of exact anatomical VOIs within the striatum. Besides, also the definition of VOIs other than striatal ones is rendered possible, resulting in the present first-time evidence of GABAergic effects on thalamic, nigral/ventral tegmental, and hippocampal DA obtained with in vivo imaging methods.
The maximum VOI diameters are either in the range or beyond the spatial resolution of the employed imaging system. It must be borne in mind, however, that the quantification of D 2/3 R binding in those portions of VOIs, which fall short of the full width at half maximum, may be hampered by partial volume effects, leading to underestimations of radioligand accumulations. Another pitfall may be overestimation of radioligand binding due to spill-over from regions with high radioligand accumulation such as the extraorbital Harderian glands to the adjacent VOIs of FC, CP, and NAC, or from the CP to NAC, THAL, and aHIPP. It can be maintained, however, that (semi)quantitative data-also of regions such as NAC, THAL, SN/VTA, and HIPP-are comparable between BAS and challenge in the present study as well as between investigations on other rats perfomed with the same imaging tool (Nikolaus et al., 2011(Nikolaus et al., , 2013(Nikolaus et al., , 2014a(Nikolaus et al., , 2016(Nikolaus et al., , 2017. As soon as further imaging studies on D 2 R-like binding sites in these regions will be available, it will be interesting to see, in which range the present binding ratios lie in comparison to binding ratios obtained with new-generation SPECT systems.
Dysfunctions of GABA A Rs and D 2 R-like binding sites have been implied in neuropsychiatric disorders including anxiety disorders and schizophrenia (Nikolaus et al., 2010(Nikolaus et al., , 2014b. Interestingly, the decline of GABA A R function in anxiety disorders involves the whole nigrostriatal and mesolimbic system, while the function of D 2 R-like binding sites is merely impaired in the CP. Contrarily, in schizophrenia, the reduction of GABA A R binding is confined to the neocortex, while D 2 Rlike binding sites are dysfunctional throughout the nigrostriatal and mesolimbic system. This implies that the emotional and behavioral changes characteristic for the individual diseases are related to regional neurochemic alterations of receptor function. The present study in the rat constitutes an important step toward unraveling the complex interdependencies of behaviors and the neurochemistry of DA and GABA using a two-modality in vivo imaging approach.

CONCLUSION
Findings indicate direct GABAergic control over synaptic DA levels in the migrostriatal and mesolimbic system in relation to behavioral action. The may be of relevance for neuropsychiatric conditions such as anxiety disorder and schizophrenia, which are characterized by both DAergic and GABAergic dysfunction.