Whole-Brain Neural Connectivity to Lateral Pontine Tegmentum GABAergic Neurons in Mice

The GABAergic neurons in the lateral pontine tegmentum (LPT) play key roles in the regulation of sleep and locomotion. The dysfunction of the LPT is related to neurological disorders such as rapid eye movement sleep behavior disorder and ocular flutter. However, the whole-brain neural connectivity to LPT GABAergic neurons remains poorly understood. Using virus-based, cell-type-specific, retrograde and anterograde tracing systems, we mapped the monosynaptic inputs and axonal projections of LPT GABAergic neurons in mice. We found that LPT GABAergic neurons received inputs mainly from the superior colliculus, substantia nigra pars reticulata, dorsal raphe nucleus (DR), lateral hypothalamic area (LHA), parasubthalamic nucleus, and periaqueductal gray (PAG), as well as the limbic system (e.g., central nucleus of the amygdala). Further immunofluorescence assays revealed that the inputs to LPT GABAergic neurons were colocalized with several markers associated with important neural functions, especially the sleep-wake cycle. Moreover, numerous LPT GABAergic neuronal varicosities were observed in the medial and midline part of the thalamus, the LHA, PAG, DR, and parabrachial nuclei. Interestingly, LPT GABAergic neurons formed reciprocal connections with areas related to sleep-wake and motor control, including the LHA, PAG, DR, parabrachial nuclei, and superior colliculus, only the LPT-DR connections were in an equally bidirectional manner. These results provide a structural framework to understand the underlying neural mechanisms of rapid eye movement sleep behavior disorder and disorders of saccades.


INTRODUCTION
The pontine tegmentum, a brainstem structure, are all the regions from the basilar pons to the fourth ventricle. It plays important roles in perception, movement, vigilance, respiration, and the sleep-wake regulation (Plazzi et al., 1996;Alheid et al., 2004;Tziridis et al., 2012). The lateral pontine tegmentum (LPT) is the dorsolateral part of the pontine tegmentum, which is included in the deep mesencephalic nucleus (DpMe). The LPT is identified as the site between the ventral lateral periaqueductal gray (vlPAG) and the pedunculopontine tegmental nucleus (PPT) on the ventro-dorsal localization, and is a trip lying the region medial to the PPT on the rostro-caudal extent in rodents (approximately −4.04 to −5.02 mm corresponding to Bregma in the mouse brain atlas) (Franklin and Paxinos, 2001;Lu et al., 2006;Sherman et al., 2015). The LPT mainly contains glutamatergic and GABAergic neurons .
The LPT is involved in voluntary locomotion, general anesthesia, and even incentive learning (Sukhotinsky et al., 2005;Fogel et al., 2010;Le Ray et al., 2011). Dysfunction of the LPT is associated with neurological disorders such as rapid eye movement (REM) sleep behavior disorder (RBD), saccade disorders and temporal lobe epilepsy (TLE). Using magnetic resonance imaging, diffusion changes of visualized structural abnormalities were found in the areas including the LPT in patients with idiopathic RBD (Scherfler et al., 2011). Ocular flutter was found to be related to the right upper pontine tegmentum damage, which contained LPT neurons (Tsutsumi et al., 2009). LPT dysfunction in patients with TLE was found to contribute to neurocognitive disturbances (Englot et al., 2018). In addition to the above, LPT neurons are especially known to have a vital function in REM sleep regulation in rodents (Sapin et al., 2009) and demyelination of LPT neurons led to dream-like hallucinations during wake, relating to the impairment of REM sleep inhibitory mechanisms (Vita et al., 2008). Furthermore, c-fos, a marker of activated neurons, was strongly expressed in the LPT after REM sleep deprivation and colocalized with glutamate decarboxylase, a key enzyme for the production of GABA, suggesting that LPT GABAergic neurons are involved in REM sleep regulation (Lu et al., 2006;Sapin et al., 2009).
Recent genetic tracing methods using the retrograde trans-synaptic rabies virus (RV) allow us to identify the presynaptic connections of a genetically well-defined neuronal population (Wickersham et al., 2007;Yuan et al., 2018). Moreover, the use of adeno-associated virus (AAV) vectors allows cell-type-specific, non-toxic, long-term expression of transgenes for anterograde tracing of neuronal pathways (Kuhlman and Huang, 2008;Holloway et al., 2013;Zhang et al., 2013). Here, we employed a modified RV and AAV in combination with the Cre/LoxP system to map presynaptic inputs and output projections of LPT GABAergic neurons. We identified several monosynaptic inputs of LPT GABAergic neurons by double labeling with antibodies for molecular markers related to sleep-wake regulation. Our results revealed the wide-ranging connectivity of LPT GABAergic neurons and provided a structural framework to understand the underlying neural circuits of important physiological functions.

Animals
The viral-based tracing experiments were all performed in male GAD2-IRES-Cre mice obtained from Miao He and Z. Josh Huang's laboratory, which generated the strain on a mixed genetic background (129SVj/B6; with black coat color) to genetically target almost whole-brain GABAergic neurons with a high level of specificity (Taniguchi et al., 2011). The mice were bred and housed under an automatically controlled 12-h light/12-h dark cycle (lights on at 7 a.m.; 100 lux intensity) . After surgery, the mice were allowed to recover for at least 2 weeks before further experiments.

Surgery and Viral Injections
Totally, 28 mice were injected with viruses. Only eight animals were correctly injected after histological confirmation and were included in the results. For virus injection, adult GAD2-IRES-Cre mice were anesthetized with ∼1.5% isoflurane in oxygen (flow rate of 1 L/min). Using the Nanoject II (Drummond Scientific, Broomall, PA), 23 nL AAV9-EF1α-DIO-TVA-GFP and AAV9-EF1α-DIO-RvG was stereotaxically injected (4.6 nL/s) via a micropipette into the unilateral LPT (anterior-posterior: −4.2 mm, medial-lateral: +1.0 mm, dorsal-ventral: −3.5 mm) (Krenzer et al., 2011). Following the injection, the pipette was held in place for an additional 10 min to allow diffusion of viral particles away from the injection site before being slowly withdrawn. Two weeks later, 46 nL volume of RV-EnvA-RG-DsRed was prepared to inject into the previous site of the LPT (n = 4 mice). For anterograde tracing, AAV-EF1α-FLEX-ChR2-mCherry was injected into the LPT following the same procedures described above. After 3 weeks, these mice were perfused (n = 4 mice).

Histology and Immunostaining
After being deeply anesthetized, adult mice were transcardially perfused with cold normal saline followed by 4% paraformaldehyde in 0.1 M phosphate buffer (PB). The brains were then extracted, post-fixed for 6 h, and then incubated in 20% sucrose in PB at 4 • C until they sank. Coronal sections (30 µm) were cut from a fixed brain on a freezing microtome (CM1950, Leica, Germany) into four series to obtain groups of tissue for multiple manipulations, and the distance between sections in each group was 90 µm. Tissue groups were restored in a cryopreservative solution.

Imaging and Analysis
Whole-brain sections were imaged by a 10 or 20× objective on the VS120 virtual microscopy slide scanning system (Olympus) and magnified images of brain sections were captured using a 20 or 40× objective by a confocal microscope, to obtain more detail (Olympus Fluoview 1000, Tokyo, Japan).
For cell mapping of neurons, neural bodies were quantified semi-automatically using ImageJ software. Based on mouse brain map, ImageJ was used to manually depict the boundaries of specific brain regions (Franklin and Paxinos, 2001). For starter cell mapping, we first applied ImageJ to distinguish the cells coexpressing DsRed and GFP as starter cells, then generated cell representation by applying the automatic wand (tracing) tool and bicubic interpolation to maximize neuronal fidelity. Next, we inverted the colorless areas to white and matched the cellular outlines to the appropriate brain regions based on the mouse brain atlas (Figure 2A). Starter cells were binned at 0.12 mm along the anterior-posterior axis, centered at the brain slice coordinate for each coronal section image.
For axonal varicosity counting, the images were captured by a 20× objective on the Olympus VS120 system. The axonal varicosity values of overall brain were calculated semi-automatically by particle analyzing plugin in ImageJ. If the axons had a transverse diameter greater than 0.5 µm, then varicosities were defined (Mechawar et al., 2000;Li et al., 2002). The outlines of brain regions were also depicted by ImageJ based on the reference brain atlas.
All data values were presented as the mean ± standard error of the mean (SEM). The strength and direction of the linear relationship between subregions and cells or varicosity proportion were measured by the Pearson product-moment correlation coefficient (CC). Two-tailed Student's t-tests were used to compare the monosynaptic inputs and axonal projections of LPT GABAergic neurons.

Applying the Rabies Virus System to Identify Monosynaptic Inputs to LPT GABAergic Neurons
In order to identify monosynaptic inputs to LPT GABAergic neurons, we used a cell-type-specific, RG-deleted RV strategy in GAD2-IRES-Cre mice. The approach has been reported to mark monosynaptic inputs to specifically selected starter cells and has been successfully used in previous studies by our group . GAD2-IRES-Cre mice were injected in the LPT with AAV helper viruses fused to green fluorescent protein (GFP). Two weeks later, the modified RV was injected in the same area. One week after this second injection, the mice were perfused, and the brains were processed (Figures 1A,B).
Rabies virus infection was clearly detectable at the injection site in GAD2-IRES-Cre mice compared with wild-type littermates, indicating no leakage of viral infection ( Figure 1C). The location of the injection site was verified by identifying virus-infected neurons in the dorsomedial corner of the PPT and the ventrolateral corner of the PAG. The yellow staining in neurons indicated that the starter cells were co-infected by AAV helper viruses and RV, and the purple staining in neurons presented ChAT immunoreactive signals for the PPT, which was the border of the LPT (Figures 1D,E). There were also some DsRed-positive neurons not expressing GFP in the LPT, revealing the existence of local input directly to LPT GABAergic neurons ( Figures 1D,E, 2A). In addition, we mapped the GABAergic starter cells as described in Section "Materials and Methods" (Figure 2A) and revealed that they were mostly restricted to the LPT. Most of the starter neurons were found to be located in the coronal planes between −4.16 and −4.84 mm from bregma, especially between −4.36 and −4.60 mm from bregma ( Figure 2B).

Input Patterns of LPT GABAergic Neurons
To determine the presynaptic connections of LPT GABAergic neurons, we analyzed serial coronal sections from mice injected with TVA-based trans-synaptic tracer. Sections from a representative GAD2-IRES-Cre brain revealed RV-DsRed-labeled presynaptic input neurons located only in specific brain nuclei in a unilateral manner (Figure 3). Figure 4 depicts input tracing to LPT GABA neurons in greater detail, with enlarged images of inputs from typical subregions such as the bed nucleus of the stria terminalis (BNST), central nucleus of the amygdala (CeA), parasubthalamic nucleus (PSTN), substantia nigra pars reticulata (SNr), SC, dorsal raphe nucleus (DR), PAG, and tuberomammillary nucleus (TMN). Further immunofluorescence assays showed that the inputs to LPT GABAergic neurons were colocalized with several markers associated with important neural functions, especially the sleep-wake cycle. In the hypothalamus, the inputs to LPT GABAergic neurons were found to be partly colocalized with orexin (0.14 ± 0.04 in Figure 5A) or MCH (0.11 ± 0.02 in Figure 5B) in the LHA, which was shown to inhibit or promote REM sleep (Adamantidis et al., 2007;Sasaki et al., 2011;Vetrivelan et al., 2016). The input neurons from the PSTN to LPT GABAergic neurons were mostly colocalized with calretinin (CR, 0.44 ± 0.04 in Figure 5C), where the inputs might to be related to the sleep-wake cycle (Steininger et al., 1999). In addition, the inputs from the limbic system to the LPT were partly colocalized with GABA neurons, such as the BNST (0.32 ± 0.03 in Figure 5D) and CeA (0.41 ± 0.03 in Figure 5E), which may be associated with sleep-wake regulation (Kodani et al., 2017;Mahoney et al., 2017) The inputs from the DR to the LPT were found partly colocalized with serotonin (0.20 ± 0.04 in Figure 5F), which was proved to promote wakefulness (Ito et al., 2013). Beyond our expectation, numerous inputs from the LPT were colocalized with parvalbumin in the SNr (0.65 ± 0.03 in Figure 5G), a region that was reported to mostly contain GABAergic neurons and to be possibly involved in REM sleep (Pal and Mallick, 2009).  After identifying the brain regions with monosynaptic input to LPT GABAergic neurons, the distribution of these DsRed-labeled afferent neurons was analyzed based on detailed statistics. We divided each brain into six general structures, namely the midbrain, diencephalon, telencephalon, cerebellum, medulla, and pons, which together encompassed 48 specific brain regions containing the DsRed-labeled neurons in the whole brain (n = 4 mice). Then we calculated the proportion of input from each region against the total number of input neurons (Figure 6, left). We found that the midbrain provided the highest numbers FIGURE 4 | Schematic presenting the typical regions with monosynaptic inputs to LPT GABAergic neurons. Primary inputs originated from regions involved in sleep-wake regulation (e.g., LHA, TMN, and SLD), locomotion (e.g., SC and SNr), emotion (e.g., CeA, and DR), and reward control (e.g., VTA and NAc). of inputs to the LPT, while few labeled neurons were found in the medulla or cerebellum. The largest number of inputs to LPT GABAergic neurons was found to arrive from the SC (15.10 ± 3.44%). The brain regions sending the second and third largest number of inputs to LPT GABAergic neurons were the SNr (10.00 ± 1.62%) and the DR (6.13 ± 1.07%) in the midbrain. In addition, the DpMe (5.17 ± 0.67%), PAG (3.64 ± 0.40%), and the substantia nigra pars compacta (SNc, 2.83 ± 0.68%) in the midbrain, the LHA (4.41 ± 1.05%) and PSTN (3.12 ± 1.08%) in the diencephalon, and the CeA (4.62 ± 1.71%) and the BNST (3.29 ± 0.48%) in the telencephalon also had strong projections to LPT GABAergic neurons.

Output Patterns of LPT GABAergic Neurons
We next mapped the output of LPT GABAergic neurons using an AAV-mediated virus, similar to our previous work with viral anterograde tracing (Zhang et al., 2013;Chen et al., 2016;Oishi et al., 2017;Wang et al., 2017). To label the axonal projections, we injected AAV expressing Cre-dependent ChR2-mCherry ( Figure 7A) into the LPT of GAD2-IRES-Cre mice. After 3 weeks, mCherry expression was observed to be restricted to the LPT, between the PPT and vlPAG (Figure 7B).

Reciprocal Connections Between the LPT and Other Nuclei
Comparing the broad distribution of the inputs and outputs of the LPT, we found that there was a high correlation between LPT afferents and efferents on a large scale of the whole brain. We found the CC was high between the afferent and efferent spatial distributions in all six major brain subdivisions (R = 0.867). Among the six major brain subdivisions, the diencephalon, midbrain, and pons revealed strong reciprocal connections with the LPT; they heavily received projections from and sent axonal fibers to LPT GABAergic neurons. Specifically, the midbrain provided a significantly higher proportion of inputs to the LPT, compared with the proportion of projections from the LPT to the midbrain (53.31 ± 6.33 versus 34.87 ± 2.84%, P = 0.038; Individual P-values indicate two-tailed Student's t-test comparisons of inputs and outputs for each brain region. In contrast, the pons and diencephalon had higher proportions of fibers arriving from the LPT than inputs going into the LPT (Figure 10A, pons: 19.94 ± 1.36 versus 12.62 ± 1.17%, P = 0.006; diencephalon: 24.03 ± 2.59 versus 14.22 ± 2.13%, P = 0.027 by two-tailed Student's t-test).
Our analysis of the proportion of inputs versus outputs in each region revealed that LPT GABAergic neurons had equally interactive connections to structures within the DR in the midbrain ( Figure 10B; the proportion of the inputs to the LPT versus the outputs from the LPT, 6.04 ± 1.07 versus 6.39 ± 0.50%). However, there were many brain regions where connections with the LPT dominated in a unidirectional manner, such as the PAG and SC (Figure 10B). Figures 10C,D summarizes the afferents and efferents of LPT GABAergic neurons in schematic form, reflecting the anatomical distribution of these innervations.

DISCUSSION
By applying advanced viral tracing methods, we identified for the first time the whole-brain distribution of input and output projections of LPT GABAergic neurons. The afferents and efferents of LPT GABAergic neurons were mostly asymmetric in major nuclei such as the LHA, SC, PAG, PBN, and SNr. Only the DR had equally reciprocal connections with the LPT. In addition, LPT GABAergic neurons received projections from several brain areas related to sleep-wake regulation and motor activity. Finally, we found a strong connection between LPT GABAergic neurons and the limbic system, suggesting that emotional processing may be associated with LPT-mediated REM sleep and motor regulation.

Technical Advances and Considerations
Genetically engineered RV has been widely used in neuroscience for its cell type-specific infection of neurons that does not affect passing neuronal tracts and its effective passage across known synaptic connections (Wickersham et al., 2007;Weissbourd et al., 2014). Accordingly, our data obtained using RV-mediated retrograde tracing in GAD2-IRES-Cre mice demonstrated that there are, in fact, direct projections from the SLD and LHA to LPT GABAergic neurons. This result not only corroborated a similar finding using CTB retrograde tracing (Lu et al., 2006;Clement et al., 2012), but also firstly identified GABAergic neurons in the LPT that received the inputs from the SLD and LHA. In addition, there are some differences between the distributions of the inputs to LPT GABAergic neurons compared with the inputs to LPT non-specific neurons. Clement and colleagues report that the ZI sends numerus projections to the LPT by CTB tracing (Clement et al., 2012). However, we found that the specific connection from ZI to LPT GABAergic neurons by RV-based tracing was weak, suggesting that the ZI may send most projections to glutamatergic neurons in the LPT.
On the other side, we built a detailed and quantitative map of LPT GABAergic neuron projections to the entire brain by AAV-mediated anterograde tracing, which was more specific and more restricted than the traditional methods. For example, although the non-specific projections from the LPT to the subthalamic nucleus (STN) were found by biotinylated dextran amine (BDA) (Bevan and Bolam, 1995), they didn't give the statistics of GABAergic projections from the LPT or full view of any other projecting nuclei. In our results, we quantitatively showed that LPT GABAergic neurons sent GABAergic varicosities in the STN (<0.1% of whole-brain varicosities). By visualizing the AAV reporter mCherry with DAB immunostaining, we avoided false positive FIGURE 8 | Schematic presenting the typical regions with LPT GABAergic neuronal outputs. These areas received LPT GABAergic innervation involved in sleep-wake regulation (e.g., LHA, LC, and SLD), locomotion (e.g., SC and SNc), emotion (e.g., CeA, and BNST), and reward control (e.g., VTA and NAc). Scale bar: 20 µm. STRv, stratum ventral region; NAc, nucleus accumbens; BNST, bed nucleus of the stria terminalis; CeA, central amygdaloid nucleus; PALv, pallidum, ventral region; VP, ventral pallidum; ZI, zona incerta; STN, subthalamic nucleus; PSTN, parasubthalamic nucleus; LHA, lateral hypothalamic area; MHA, medial hypothalamic area; PF, parafascicular thalamic nucleus; MD, mediodorsal thalamic nucleus; PVT, paraventricular thalamic nucleus; SNc, substantia nigra pars compacta; DpMe, deep mesencephalic nucleus; vlPAG, ventrolateral periaqueductal gray; SC, Superior colliculus; DR, dorsal raphe nucleus; PBN, parabrachial nucleus; SLD, sublateral dorsal nucleus; PCRt, parvicellular reticular nucleus, alpha part; IRT, intermediate reticular nucleus; LC, locus coeruleus; fr, fasciculus retroflexus; acp, anterior commissure, posterior; scp, superior cerebellar peduncle. signals, which are commonly observed after fluorescence immunohistochemistry (Oh et al., 2014). We chose to quantify projections by counting axonal varicosities instead of measuring fluorescence intensity, in order to avoid differences between thickness of brain sections or light exposure by the microscope.

Morphological Findings
Due to the technical shortcomings of classical tracing methods, the complete neuronal connectivity of GABAergic neurons in the LPT was previously unknown. We used cell-type-specific anterograde and retrograde viral tracing to map the connectivity of LPT GABAergic neurons and found that their main inputs (≥20%) originated from the midbrain, whereas the brain regions that received projections from LPT GABAergic neurons were mainly distributed (≥20%) between the telencephalon, midbrain, and pons.
From comparison of distribution among the major nuclei more unidirectionally connecting to LPT GABAergic neurons, we found that afferents to LPT GABAergic neurons concentrated on the SC and SNr, whereas GABAergic efferents from the LPT were mostly distributed in the PAG, PBN, and LHA. Most afferents to LPT GABAergic neurons originated from the SC, particularly the intermediate and deep layers of the SC (Figures 3, 6). Although the projections from the SC to the LPT have been demonstrated in monkeys by anterograde tracing with BDA (Wang et al., 2013), the fact that these projections are predominantly GABAergic was unknown. Although previous studies indicated that the LPT received inhibitory projections from the SNr (Sherman et al., 2015), we identified the SNr as the second largest input to LPT GABAergic neurons. In addition, we showed that the main outputs of LPT GABAergic neurons project to the PAG, especially the vlPAG (Figures 8, 9), while LPT-projecting PAG neurons were already known (Weissbourd et al., 2014). It was previously suggested that the majority of LPT cells provide GABAergic inhibition of PPT cholinergic neurons, including those cholinergic neurons that provide a major ascending pathway into the posterior LHA (Ford et al., 1995). However, our data revealed that the LHA directly received strong projections from LPT GABAergic neurons.
The proportion of inputs to the LPT among whole-brain inputs is different from that of the LPT. A pathway from the PSTN to the LPT was previously reported using phaseolus vulgaris leucoagglutinin anterograde tracing (Goto and Swanson, 2004), and by this method it is not possible to determine whether a majority or a minority of PSTN neurons project to LPT. Here, we reported the proportion of the inputs from the PSTN projecting to LPT GABAergic neurons among whole-brain inputs was 3.02 ± 1.09%, which was about 44% colocalized with CR neurons (Figure 5C). Combined with dense inputs from the CeA in the forebrain, our results suggested a strong ability of the LPT to integrate behavioral information from the PSTN and CeA such as feeding and fear behaviors (Zimmerman et al., 2007;Chometton et al., 2016). Although the LPT has been shown to provide strong GABAergic projections to the dorsocaudal region of the pontine tegmentum, including the PBN, LC, and SLD (Hayashi et al., 2015), our finding of dense, LPT-originating, GABAergic axonal varicosities in the SNc of the midbrain suggested that the LPT also transmitted strong information to the SNc.

Functional Implications
LPT GABAergic neurons have strong connections with brain regions related to sleep-wake regulation. Firstly, the output to the vlPAG is known to gate REM sleep and the ultradian rhythm of REM/NREM alternation (Weber et al., 2018). In addition, vlPAG/DpMe GABAergic neurons were shown to be important in sleep regulation (Hayashi et al., 2015;Weber and Dan, 2016) and thus, bidirectional projections between the LPT and FIGURE 10 | Comparison between inputs and outputs of LPT GABAergic neurons. (A) LPT GABAergic inputs and outputs were compared in six major structural subdivisions of the whole brain. LPT GABAergic inputs are colored red and the outputs are shaded in orange. The data represent the mean ± SEM, with n = 4 for each group. Two-tailed Student's t-tests indicated significant ( * P < 0.05, * * P < 0.01) differences between the inputs and outputs of LPT GABAergic neurons in the same subdivision.  vlPAG may be important for sleep regulation. Strong inhibitory projections from the LPT to the arousal-promoting PBN may also indicate a role of the LPT in sleep/wake control and thermoregulation (Fuller et al., 2011;Kaur et al., 2013;Nakamura, 2018;Wang et al., 2018). Furthermore, we and others identified SLD-projecting LPT GABAergic neurons (Boissard et al., 2003), which may play important roles in REM sleep regulation and related sleep disorders, such as narcolepsy or REM sleep behavior disorder. The direct connections we observed between the LPT and the hypothalamus suggested that LPT GABAergic neurons might be involved in regulation of hypothalamus functions, such as REM sleep, arousal, feeding behaviors, body temperature, and sense of pain (Clement et al., 2012;Li et al., 2014). Finally, another study showed that the neuronal activity of the posterior lateral hypothalamus containing the PSTN was related to the wake and REM sleep (Steininger et al., 1999). In the present study, the PSTN sent abundant projections to the LPT GABAergic neurons, suggesting that the pathway of the PSTN to LPT GABA neurons might be involved in sleep-wake regulation. Besides that, chemical activation of PSTN neurons elicited depressor as well as bradycardia (Ciriello et al., 2008), and hedonic tastes increased c-Fos expression of the PSTN in rats (Chometton et al., 2016), suggesting that cardiovascular regulation and food intake were probably mediated by the PSTN innervated LPT GABAergic neurons.
Anatomical connectivity of LPT GABAergic neurons with the SC, SN, and PAG in the midbrain may implicate a role of LPT in the regulation of motor behaviors. Our data revealed the existence of asymmetrically reciprocal circuits between the SC and LPT GABAergic neurons for the transmission of saccadic information. The LPT was suggested to be involved in the control of eye movements, due to a causal link that was reported between ocular flutter and small lesions in the right upper pontine tegmentum containing the LPT (Tsutsumi et al., 2009). Moreover, strong connections between the LPT and SN suggested a role in regulating locomotor behavior (Deng et al., 2016;Lintz and Felsen, 2016;Barchini et al., 2018;Weber et al., 2018). A study previously proposed that projections of SNr neurons to PPT glutamatergic neurons primarily regulated motor activity (Roseberry et al., 2016); however, the strong projections from the SNr to LPT GABAergic neurons may also implicate inhibitory LPT neurons in some aspects of motor control such as saccadic eye movement (Sato and Hikosaka, 2002). In contrast, SNcprojecting LPT GABAergic neurons were shown to be involved in motor and reward control (Pioli et al., 2008). Our finding that SNc neurons also modestly projected to LPT GABAergic neurons could indicate a sophisticated feedback mechanism between SNc and LPT GABAergic neurons.
Finally, the inputs from the limbic system, including the CeA and BNST, suggested that LPT GABAergic neurons may be regulated by neural circuits involved in stress and fear processing (Zhao and Davis, 2004). Cataplexy, a common symptom of narcolepsy, is known to be triggered by strong emotions associated with GABAergic neurons in the limbic system and may also result from the disturbance of REM sleep atonia (related to the LPT/vlPAG/LC) into wakefulness (Snow et al., 2017). Our data provided evidence that LPT GABAergic neurons directly received projections from BNST GABAergic neurons (Figures 4, 5D) and CeA GABAergic neurons (Figures 4, 5E). Therefore, we confirmed the anatomical connectivity of the limbic system with LPT GABAergic neurons that control REM sleep atonia, as hypothesized by Fuller (Peever and Fuller, 2017).

Comparison With Areas Near the LPT
Although the vlPAG and PPT are in close proximity to the LPT, the connectivity of these GABAergic neurons is different from that of the LPT. Whereas vlPAG GABAergic neurons primarily receive inputs from the ventral medulla, the SLD, and the LHA (Boissard et al., 2003;Clement et al., 2012;Weber et al., 2015), the LPT is mainly innervated by the SC and SNr. In contrast to LPT GABAergic neurons, the vlPAG has major projections to the RMg and nearby areas, commonly known as the rostral ventral medulla (Bowman et al., 2013). Moreover, PPT GABAergic neurons are rarely innervated by the SNr (Roseberry et al., 2016), while LPT GABAergic neurons are strongly innervated by this area.
There are also clear differences in the physiological and pathological functions of the LPT, vlPAG, and PPT. GABAergic neurons in the vlPAG are important for the sleep-wake cycle, micturition, and nociception. In addition, optogenetic studies revealed that vlPAG GABAergic neurons regulated non-REM and REM sleep (Weber et al., 2018). Other studies based on genetically engineered systems showed that activation of vlPAG GABAergic neurons delayed detrusor contraction and inhibited voiding (Zare et al., 2018). In contrast, PPT GABAergic neurons are mostly involved in gait and balance regulation, a function that is not identified in LPT GABAergic neurons. For example, loss of rostral PPT GABAergic neurons was reported in Parkinson's disease patients with gait and balance disorders (Sebille et al., 2018). In spite of the close proximity of the vlPAG, PPT, and LPT GABAergic neurons, it is likely that these areas have distinguishing physiological and pathological roles.

CONCLUSION
We mapped, for the first time, the afferents and efferents of LPT GABAergic neurons and found them to be extensively interconnected with other brain areas. This suggested a vital role of LPT neurons in a wide range of physiological and pathological functions, especially sleep-wake regulation and locomotor control. Our anatomical data could be useful for future functional studies of the brain and also provide a structural basis to understand neurological disorders.

ETHICS STATEMENT
This study was carried out in accordance with the recommendations of the China Regulations on the Administration of Laboratory Animals, the Decree No.

AUTHOR CONTRIBUTIONS
Z-KC, X-SY, and HD designed and performed the experiments, analyzed the data, and wrote the manuscript. Y-FW performed the experiments and analyzed the data. W-MQ and G-HC conceived the experiments and wrote the manuscript. MH provided GAD2-IRES-Cre mice. Z-LH conceived the experiments, analyzed the data, and wrote the manuscript.

FUNDING
This study was supported in part by grants-in-aid for scientific research from the National Natural Science Foundation of China (Grant Nos. 81420108015, 31530035, 31871072, 31671099, and 31571103) and the National Basic Research Program of China (Grant No. 2015CB856401).