Endocannabinoid 2-Arachidonoylglycerol Levels in the Anterior Cingulate Cortex, Caudate Putamen, Nucleus Accumbens, and Piriform Cortex Were Upregulated by Chronic Restraint Stress

Abstract

Endocannabinoid 2-arachidonoylglycerol (2-AG) has been implicated in habituation to stress, and its augmentation reduces stress-induced anxiety-like behavior. Chronic restraint stress (CRS) changes the 2-AG levels in some gross brain areas, such as the forebrain. However, the detailed spatial distribution of 2-AG and its changes by CRS in stress processing-related anatomical structures such as the anterior cingulate cortex (ACC), caudate putamen (CP), nucleus accumbens (NAc), and piriform cortex (PIR) are still unclear. In this study, mice were restrained for 30 min in a 50 mL-centrifuge tube for eight consecutive days, followed by imaging of the coronal brain sections of control and stressed mice using desorption electrospray ionization mass spectrometry imaging (DESI-MSI). The results showed that from the forebrain to the cerebellum, 2-AG levels were highest in the hypothalamus and lowest in the hippocampal region. 2-AG levels were significantly (p < 0.05) upregulated and 2-AG precursors levels were significantly (p < 0.05) downregulated in the ACC, CP, NAc, and PIR of stressed mice compared with control mice. This study provided direct evidence of 2-AG expression and changes, suggesting that 2-AG levels are increased in the ACC CP, NAc, and PIR when individuals are under chronic stress.

1. Introduction

Endocannabinoid 2-arachidonoylglycerol (2-AG) is one of the main endogenous cannabinoid ligands in the central nervous system [1]. It is a neutral lipid that has the intrinsic tendency to associate with membranes, synthesize “on demand” [2], release from postsynaptic neurons, and travel backward across synapses, activating the CB1 receptor on presynaptic axons and suppressing neurotransmitter release, functioning as a retrograde synaptic messenger [3]. 2-AG could be biosynthesized from 1-oleoyl-2-arachidonoyl-sn-glycerol (OAG) [4] and 1-stearoyl-2-arachidonoylglycerol (SAG) [5] by diacylglycerol lipases (DAGLs) in the brain tissues. DAGLs are considered the most important enzymes for 2-AG biosynthesis, and monoacylglycerol lipase (MGL) is a key enzyme for 2-AG hydrolysis [2]. 2-AG displays anti-inflammatory and neuroprotective properties [6]. It also participates in responses to stress stimuli and is becoming a potential new therapeutic target in treating major depressive disorder (MDD) [7]. It might act to mitigate depressive symptoms [8], reduce stress-induced anxiety-like behavior [9], and participate in habituation to repeated stress [10]. One study has suggested an association between stress-related disorders and an increased risk of neurodegenerative diseases, such as Alzheimer’s disease [11]. All the information suggests that 2-AG is an essential lipid messenger in the brain. Revealing the spatial distribution of 2-AG in brain tissue will be beneficial for the study of related diseases, such as stress-related disorders and neurodegenerative diseases.

Stress is a risk factor for the development and exacerbation of many diseases, such as neuropsychiatric disorders [12]. Chronic restraint stress (CRS) is a preclinical chronic stress exposure model for investigating stress-relevant disorders [13]. It has been reported that 2-AG levels were significantly increased in the forebrain after 5 days of CRS exposure [14] and in the medial prefrontal cortex (mPFC) after 10 days of CRS exposure [10]. These results were measured by liquid chromatography-mass spectrometry (LC-MS). However, LC-MS could not provide the spatial information of 2-AG, and it is difficult to apply to small brain anatomical structures. The mPFC is a part of the forebrain, containing the anterior cingulate cortex (ACC), prelimbic cortex, and infralimbic cortex. ACC is associated with stress-related illnesses, such as MDD [15], anxiety, and mood disorders [16]. Extensive preclinical evidence has also established the mediating role of the caudate putamen (CP) [17], nucleus accumbens (NAc) [18], and piriform cortex (PIR) [19] in stress processing or stress-related disorders. However, the detailed spatial distribution of and changes in 2-AG after CRS exposure in these regions are still unclear.

Mass spectrometry imaging (MSI) can provide complementary information for molecular analysis and is applied by the scientific community to study proteins, peptides, small molecules, and metabolites [20]. Desorption electrospray ionization MSI (DESI-MSI) is a technique that allows the spatial intensity distribution to be recorded directly from histological sections without prior chemical treatment under ambient conditions. It can simultaneously visualize target molecules, such as metabolites [21] and lipids [22]. Our lab has already measured different molecules applying this technique [23], including 2-AG, and revealed that 3-day water-immersion stress increased 2-AG levels in the hypothalamus, midbrain, and hindbrain of senescence-accelerated mouse prone 8 mice [24].

In this study, we used DESI-MSI to investigate and reveal the 2-AG distribution in wild-type mice brain tissues and the effect of CRS on 2-AG levels in the ACC, CP, NAc, and PIR. Doing so will allow for a clearer understanding of 2-AG expression and changes in the brain and extend the horizon for future etiological and therapeutic studies of related diseases.

2. Materials and Methods

2.1. Animals

Experiments were performed on male C57BL/6JJmsSlc mice (25–30 g, 4–6 months old). They were bred in a climate-controlled standard animal house at 22 ± 2 °C with a light/dark cycle of 12 h (light period: 7 a.m. to 7 p.m.). Mice were housed two to three per cage in standard microisolator polycarbonate cages with ad libitum access to food and water. Subjects were randomly assigned to experimental groups, and the testing sequence was randomized throughout the experiment. All procedures were conducted according to the guidelines of the Institutional Animal Care and Use Committees of Hamamatsu University School of Medicine (HUSM), and the ethical permission number is 2022021. All efforts were made to minimize the number of animals used and their suffering.

2.2. Stress Protocol

Mice were subjected to stress using a published procedure [10]. In brief, mice were marked for identification on the ear as needed and randomly grouped into stress and control groups. Mice in the stress group were separately restrained for 30 min in a 50 mL-centrifuge tube with several small holes (for ventilation) in the morning for 8 consecutive days. They were placed on the bench top during the restraint period. Then we put them back in their home cage after finishing the stressor. Mice in the control group were left undisturbed in their home cages. To ensure the reproducibility of the results while complying with the requirements of the three Rs ethical principle, we used three mice in each group. All mice were weighed before inducing restraint stress on the first and last day.

2.3. Tissue Preparation

Mice in the stress group were immediately sacrificed by cervical dislocation after the last restraint exposure. Control animals were sacrificed immediately after removal from their home cages at the same time of day and in the same room as the stressed mice. All mice brain tissues were quickly removed and rapidly frozen in powdered dry ice, then stored at −80 °C until cryo-sectioning. The brain tissues were mounted on a sample holder with optimal cutting temperature compound (Sakura Finetek Japan, Tokyo, Japan) and coronally sectioned to a thickness of 10 µm using a cryostat system (CM1950; Leica Biosystems, Wetzlar, Germany) at −20 °C. The sliced samples were mounted on the glass slide (Matsunami, Osaka, Japan, thickness: 0.8 to 1 mm, size: 76/26 mm) for DESI-MSI and immunohistochemistry (IHC) analysis. Hematoxylin and eosin (H&E) staining was conducted after DESI-MSI measurement to facilitate anatomical localization.

2.4. DESI-MSI

Data were acquired within 3 days after mice were sacrificed by a Xevo G2-XS quadrupole time-of-flight (Q-TOF) mass spectrometer (Waters, Milford, MA, USA) equipped with a 2D DESI source in positive ionization mode. Prior to acquiring data, in order to obtain better signal intensity, the detector setup was performed by leucine enkephalin solution (500 µM), and the mass spectra calibration was performed by sodium formate solution (500 µM) in 2-propanol: water (90: 10, v/v). Compounds with m/z 443 ions from the ink of a red sharpie pen (Sharpie, USA) were used to adjust the signal intensity to over 10⁶ counts (cts) and m/z 798 from brain tissues over 10⁴ cts, respectively. The final optimized parameters of DESI source were as follows. The capillary voltage was set to 4.0 kV and the sampling cone voltage was set to 50 V. The source temperature was set to 130 °C and nebulizing nitrogen gas was set at a pressure of 0.5 MPa. The spray solvent (methanol: water, 98: 2, v/v) was sprayed continuously at a flow rate of 3 µL/min. The spray impact angle was 75°; the emitter tip to tissue, emitter tip to ion transfer capillary orifice, and ion transfer capillary orifice to tissue distance were approximately 2, 6, and 0.5 mm, respectively. The samples were scanned at a velocity of 200 µm/sec with a pixel size of 100 µm × 100 µm. The data were acquired in the m/z 100 to 1000 range; the speed of acquisition was 1 spectrum sec⁻¹. The analyzer mode was set as “sensitivity” and the data type was set as “continuum”. For data processing using HDImaging software (Waters, Milford, MA, USA; version 1.4), the mass spectrum of the 1000 highest intensity peaks was collected in a mass range of m/z 100 to 1000, with the mass resolution set at 20,000 and a mass window of 0.02 Da.

Ions of m/z 309.2036 were used for lock mass correction, with a tolerance and minimum signal intensity of 0.05 amu and 500 cts, respectively. Standard 2-AG (0.1 μL) (Sigma-Aldrich, St. Louis, MO, USA) dissolved in acetonitrile at different concentrations (0 µg/mL, 1 µg/mL, 5 µg/mL, 10 µg/mL) was used to confirm the mass accuracy.

2.5. IHC Staining

IHC staining was conducted according to a previously described method [25], with some modifications. In brief, we proceeded step by step as follows: (1) The frozen brain sections were removed and allowed to return to room temperature (RT) and dry. (2) Then, they were fixed for 10 min in 4% paraformaldehyde and (3) blocked in 10% horse normal serum (HNS), 0.1% sodium azide, and 0.5% Triton X-100 in Tris-HCl-buffered saline (TBS) (pH 7.4) for 30 min at RT. (4) Primary antibody incubation was conducted with rabbit anti-MGL (1:500, Frontier Institute catalog #MGL-Rb, RRID: AB_2571798) on a shaker for 2 days at 4 °C. (5) Endogenous peroxidase was blocked with 3% (v/v) H₂O₂ for 20 min at RT. (6) The samples were incubated with goat anti-rabbit IgG (1:200, Vector Laboratories catalog #BA-1000, RRID: AB_2313606) for 1 h on a shaker at RT. (7) They were then incubated in the avidin-biotin complex (ABC, 1:50, Vector Laboratories catalog #PK-6100, RRID: AB_2336819) and prepared in a washing solution (1:50) for 30 min at RT. Between all steps, sections were washed in 1% HNS and 0.5% Triton X-100 in TBS. Visualization was performed using 0.05% diaminobenzidine (DAB) in a 0.1 M phosphate buffer solution (pH 7.6) with 0.5% Triton X-100 and 0.01% hydrogen peroxide for 3 min 30 s at RT. Finally, tissues were subsequently counterstained with hematoxylin, dehydrated in graded alcohols (80, 90, 100%), transparentized with xylene, and coverslipped with PathoMount (FUJIFILM). Controls were set up to check the accuracy.

2.6. Data Analysis

The general information of control and CRS mice was analyzed by the Mann-Whitney U test. Quantification of the immunoreactivity of representative areas was determined using the analysis software ImageJ (National Institutes of Health, Bethesda, MD, USA). MassLynx (Waters, Milford, MA, USA; version 4.1) and HDImaging software were used to acquire DESI-MSI data. In HDImaging software, raw data were normalized using total ion current, and ion images were constructed to visualize the spatial distributions of the detected molecules after setting parameters in visualization, including scale: linear, image smoothing: linear interpolation. The normalized spectra were exported to the MassLynx software to further normalize to the largest peak on display. MSI raw data were converted into .imzML by HDImaging software and then further converted to .imdx using an IMDX converter (Shimadzu, version 1.20.0.10960). Finally, the .imdx data were used to analyze the intensity of target m/z by IMAGE REVEAL software (Shimadzu, Kyoto, Japan; version 1.20.0.10960). A thresholding value was determined by the mean + SD of the control mice when comparing the 2-AG and its precursor levels between the control and stressed mice. The lower and higher intensities than the threshold value were counted and set for a 2 × 2 contingency table. Pearson’s Chi-square test was performed to analyze the MS intensity of molecules. Two-tail t-test was performed to evaluate the differences in fold changes in the average intensity of 2-AG and its precursors between the two groups. All statistical analyses were performed using SPSS version 22 (IBM) and GraphPad Prism version.8.0.2.263 (GraphPad Software, GraphPad Software, La Jolla, CA, USA). Differences with p-value < 0.05 were considered significant.

3. Results

3.1. Detection of 2-AG by DESI-MSI

After optimizing DESI-MSI parameters, the standard 2-AG was detected as potassium adduct [M + K]⁺ at m/z 417.2409. Consistently, the endogenous 2-AG in the mice brain sections at m/z 417.2409 was also detected (Figure 1A). The error in the determination of the monoisotopic mass of 2-AG was 0.48 ppm (Table 1). Ion images of different concentrations of the 2-AG standard (0 µg/mL, 1 µg/mL, 5 µg/mL, 10 µg/mL) applied to mice brain tissues were acquired (Figure 1B). We confirmed that the signal intensity of 2-AG varied depending on its concentration.

3.2. 2-AG Levels Are Highest in the Hypothalamus (HY) and Lowest in the Hippocampal Region (HIP)

We visualized 2-AG in the coronal brain sections of non-stressed male mice via DESI-MSI. The approximate localization of the sections from the forebrain to the cerebellum (sections I to VIII) is shown in Figure 2A. To evaluate the distribution of 2-AG in the mouse brains, we annotated the H&E stained coronal sections into detailed regions based on the Allen Brain Atlas, including the anterior commissure of the olfactory limb (aco), isocortex, striatum ventral region (STRv), caudate putamen (CP), pallidum (PAL), HY, amygdala, HIP, motor-related superior colliculus (SCm), midbrain reticular nucleus (MRN), periaqueductal gray (PAG), pons, medulla, cerebellar cortex (CBX), retrohippocampal region (RHP), and forebrain bundle system (Figure 2B). Different 2-AG signal intensities were detected in different functional areas of the brain. Moderate 2-AG signal intensity was detected in the aco, PAL, MRN, SCm, pons, and medulla. The remaining regions except for HY showed weaker signal intensity than the above-mentioned egions. The 2-AG levels were highest in the HY and lowest in the HIP (Figure 2C).

3.3. 2-AG Levels Were Upregulated and 2-AG Precursor Levels Were Downregulated in the ACC, CP, NAc, and PIR by CRS

The physical profiles of the control and CRS mice showed no differences in age and body weight at baseline, but the body weight on the last day of CRS exposure was significantly different (p = 0.046) (Table 2). Compared with the control group, CRS tended to inhibit body weight growth.

OAG and SAG are likely the predominant diacylglycerol species as 2-AG precursors in neuronal cells. They have been identified in the mice brain tissues and used for a 2-AG-related study by LC/MS [26]. The OAG or SAG can be hydrolyzed by DAGLs at the sn1 position to generate 2-AG (Figure 3A). Therefore, we also checked the distribution of these two precursor molecules. [OAG + K]⁺ (m/z 681.48) and [SAG + K]⁺ (m/z 683.50) were detected in the mice brain sections (Figure S1). Errors in the determination of the monoisotopic mass of OAG and SAG were −5.72 ppm and −1.17 ppm, respectively (Table 1).

To elucidate the effect of CRS on the distribution of 2-AG and 2-AG precursors in the different brain regions, we acquired DESI-MSI data from the coronal brain sections of control and stressed mice. Histograms were used to show the signal intensity of each pixel in the defined regions. We found that 2-AG levels were significantly (p < 0.05) upregulated, while OAG and SAG were significantly (p < 0.05) downregulated by CRS in the ACC, CP, NAc, and PIR (Figure 3B,C; Figures S2A,B and S3A,B).

3.4. The Quantification of the Immunohistochemical Expression of MGL in the ACC, CP, NAc, and PIR Was Not Affected by CRS

As CRS caused changes in 2-AG levels in the ACC, CP, NAc, and PIR, the expression of MGL was also investigated by IHC. Antibody validation (Figure 4A) was performed in the HIP, where the expression of MGL was known [25]. MGL immunoreactivity was observed in the ACC, CP, NAc, and PIR of control and stressed mice (Figure 4B). The quantification of the immunohistochemical expression of MGL in ACC (p = 0.130), CP (p = 0.360), NAc (p = 0.432), and PIR (p = 0.5178) between control and stress groups were not different (Figure 4C).

3.5. CRS Upregulated 2-AG Levels in the Non-LI Area but Not in the LI of ACC

The ACC consists of layers I, II–III, and V–VI in rodents and humans [27]. LI has less cell density than non-LI and receives strong dendritic tuft branches from other layers [28]. To explore the detailed changes of 2-AG and 2-AG precursors in different neural cell density areas after CRS exposure, we analyzed the changes in these molecules in LI and non-LI areas of ACC by combining DESI-MSI image results with MGL IHC-stained results. The magnification of DESI-MSI results (Figure 5) showed that 2-AG was significantly upregulated by 12% in the non-LI areas after CRS exposure (p = 0.033). However, the difference between control and stress in the LI was not obvious (p = 0.161). OAG and SAG were significantly downregulated in both non-LI and LI of the ACC after CRS exposure. OAG decreased by 35% in non-LI (p < 0.001) and 33% in LI (p = 0.003) areas. SAG decreased by 52% in non-LI (p = 0.003) and 39% in LI (p = 0.016) areas.

4. Discussion

Our study shows the spatial distribution of 2-AG in control male mouse brains and reveals a moderate 2-AG signal intensity in the aco, PAL, MRN, SCm, pons, and medulla. Except for HY, the remaining regions showed weaker signal intensity than the above-mentioned regions. The HY has the highest 2-AG levels, while the HIP has the lowest. 2-AG signaling would exert different effects in different brain functional areas by targeting corresponding receptors. CB1 receptors expressed in the amygdala, cingulate cortex, prefrontal cortex, ventral PAL, CP, and NAc [29] could also be detected in the anterior commissure [30] and the majority of hypothalamic nuclei [31]. 2-AG within HY regulates the effects of time of day and stress on short-term memory [32] and plays a role in maintaining food intake [33]. In PAG, 2-AG-CB1 signaling promotes stress-induced analgesia [34]. In the hippocampus, the release of 2-AG leads to the transient suppression of gamma-aminobutyric acid-mediated transmission by retrograde signaling [35]. DAGLα [36], MGL, and CB1 receptors [25] were detectable in the hippocampus. Our results show that 2-AG distributes at a low level in the HIP, which belongs to the hippocampus, suggesting that MGL may contribute to low levels of 2-AG.

This study is the first to show that 2-AG levels were upregulated in the ACC, CP, NAc, and PIR after CRS exposure. Others have reported that 2-AG augmentation reduces stress-induced anxiety-like behavior [9] and participates in habituation to repeated stress [10]. Our DESI-MSI results provide direct evidence that 2-AG in the ACC, CP, NAc, and PIR is upregulated under CRS conditions and may be involved in habituation to stress. The CRS not only upregulated 2-AG but also downregulated OAG and SAG. DAGLα immunoreactivity has been reported in the ACC, CP, NAc, and PIR [36]. We speculate that the demand for 2-AG in the ACC, CP, NAc, and PIR increases after CRS exposure; thus, OAG and SAG are hydrolyzed more by DAGLα, leading to a decrease in them and an increase in their degradation product 2-AG. Therefore, the downregulation of OAG and SAG may be a reason for the upregulation of 2-AG after CRS exposure.

The results show that CRS tended to inhibit body weight growth. An analysis of cumulative feeding following 2 h restraint stress revealed that it has a significant effect on body weight; restraint significantly decreased intake during early post-stress time points (1 h) and overall intake at the 22 h time point [37]. It also showed that increased 2-AG following the systemic administration of an MGL inhibitor decreased energy intake. Therefore, we speculate that CRS reduced food intake to inhibit body weight growth, and CRS-induced upregulated 2-AG may play a role in this process.

Moreover, our results show that MGL expression levels in the ACC, CP, NAc, and PIR were not different between the control and stressed groups. Therefore, we suggest that the CRS-induced 2-AG increase was not caused by the change in MGL expression levels. Previous studies have reported that CB1 receptors were detected in the ACC [38,39], CP, NAc, and PIR [40]. Therefore, a functional 2-AG signaling system composed of cannabinoid receptors and the complete machinery for the synthesis and degradation of 2-AG existed in the ACC, CP, NAc, and PIR. 2-AG contributes to stress response termination by inhibiting glutamate release and restraint following anxiety arousal by stimulating presynaptic CB1 receptors [41]. 2-AG can also act on astrocytic CB1 to indirectly regulate glutamate release via the release of gliotransmitters [42]. These reports and our results suggest that CRS-induced upregulated 2-AG in the ACC, CP, NAc, and PIR plays a role in contributing to stress response termination by regulating glutamate release through targeting presynaptic or astrocytic CB1 receptors.

The DESI-MSI results also show that CRS significantly upregulated 2-AG levels in the non-LI area but not in LI of the ACC. The MGL result shows that LI has fewer neural cells than non-LI areas in ACC. LI contains interneurons; non-LI contains both pyramidal neurons and interneurons [27]. 2-AG changes in the ACC caused by CRS may be more concentrated in pyramidal neurons. 2-AG might be involved in antidepressant mechanisms by targeting CB1 receptors in neurons and astrocytes and CB2 receptors in microglial cells [7]. 2-AG could also avoid the detrimental effects of inflammation through microglia [43]. Thus, we supposed that the increased 2-AG by CRS might be involved in antidepressant and anti-neuroinflammation.

In addition, it is important to note that excessive 2-AG inevitably increases its degradation by MGL to generate neuroinflammatory prostaglandins that promote neuroinflammation [44]. Therefore, a prolonged stressor with persistently elevated 2-AG degradation products is not beneficial for neurodegenerative disorders. Interestingly, the inhibition of MGL would promote the augmentation of the endogenous 2-AG that can be therapeutically useful for suppressing neuroinflammation and maintaining the balance between excitatory and inhibitory neurotransmission [45]. Therefore, our results on MGL expression and 2-AG changes also lay the foundation for the future use of MGL inhibitors to achieve antidepressant and anti-neuroinflammation effects.

Furthermore, the bilateral injection of 2-AG into the NAc shell produced a clear-cut, short-term stimulatory action on feeding behavior [46]. The increase of 2-AG in NAc due to stress may explain stress-induced overeating. 2-AG in the ACC likely modulates fear-conditioned analgesia, and 2-AG-CB2 receptor signaling may suppress this form of endogenous analgesia [47]. 2-AG also targets GPR55, and the activation of GPR55 signaling in the ACC facilitates inflammatory pain via a top-down modulation of descending pain control [48]. Thus, the CRS-induced elevated 2-AG in the ACC may also contribute to stress-induced analgesia or hyperalgesia processes through CB2 or GPR55 receptors.

To ensure the accuracy of the 2-AG measurements, we not only calibrated each measurement with standard 2-AG on each slide but also took special care to minimize post-mortem changes and degradation of 2-AG over time throughout the process. The steps and times were consistent from sample to sample, and the exposure time from −80 °C to −20 °C was consistent across samples and as minimal as possible during cryo-section. Moreover, we avoided any interference of drugs or injection stress because it has been indicated that morphine [49] and fentanyl [50] downregulate the contents of 2-AG. Hence, we chose cervical dislocation to sacrifice mice, a method already used in 2-AG-related animal experiments [51]. Neuronal localization of putative 2-AG was achieved by space and time coherent mapping MSI [52]. The three-dimensional distribution of 2-AG in the brain tissue can be achieved in the future by DESI-MSI, as with previous research on the three-dimensional distribution of monoamines in mouse brains [53]. DESI-MSI is a powerful method that has the potential to achieve further progress in 2-AG-related studies.

In summary, we suggest that the CRS-induced elevated 2-AG in ACC CP, NAc, and PIR constitutes self-regulation during habituation to stress. The elevated 2-AG in these regions may play a role in stress response termination, antidepressant and anti-neuroinflammation effects, and stress-induced analgesia or hyperalgesia. An imbalance of 2-AG in the brain seems to trigger or promote the development of stress-related neuropsychiatric diseases and neuroinflammation-related neurodegenerative disorders. The results of this study extend the horizon for the etiological investigation or treatment of these related diseases.

5. Conclusions

The results presented in this study show that 2-AG levels are highest in the HY and lowest in the HIP. Moreover, this study reveals that CRS upregulates 2-AG levels and downregulates its precursors’ levels in the ACC, CP, NAc, and PIR. The expression levels of MGL were not affected by CRS. Our findings extend the horizon for future etiological and therapeutic studies of 2-AG-related diseases, such as stress-related neuropsychiatric diseases.

Author Contributions

Conceptualization, T.K. and M.S.; Methodology, Q.Z., A.I., B.C., T.K. and M.S.; Validation, Q.Z.; Data analysis, Q.Z., A.I., B.C., H.Z. and Y.T.; Investigation, Q.Z. and H.Z.; Resources, D.H.C.; Writing—original draft, Q.Z.; Writing—review & editing, A.I., B.C., D.H.C., M.A.M., Y.T., N.S., H.Y., Y.N. (Yoshiki Nakajima), Y.N. (Yu Nagashima), F.S., T.K. and M.S.; Supervision, T.S., T.K. and M.S.; Project administration, M.S. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Animal Care and Use Committees of Hamamatsu University School of Medicine (protocol code: 2022021, date of approval: 30 May 2022).

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data were reported within the article. Further supporting data will be provided upon a written request addressed to the corresponding author.

Acknowledgments

We thank all the members of Mitsutoshi Setou’s lab for their support and help.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Table 1

m/z accuracy of targeted molecules from standard 2-AG and mice brain tissues.

Assigned Molecule	Theoretical m/z	Observed m/z	Mass Accuracy (ppm)
[2-AG + K]+	417.2407	417.2409	0.48
[OAG + K]+	681.4860	681.4821	−5.72
[SAG + K]+	683.5017	683.5009	−1.17

Table 2

Summary of the general information of control and stressed mice.

Group	Control	Stress	p-Value
Age (day)	165 ± 13	160 ± 13	NS
BW (g) of day 1	27.03 ± 0.28	26.40 ± 0.21	NS
BW (g) of day 8	28.03 ± 0.43	25.73 ± 0.32	0.046

The number of mice in each group was three, NS: p-value > 0.05, BW of day 1: body weight on the first day of stress exposure, BW of day 8: body weight on the last day of stress exposure. Values are mean ± SEM. Mann-Whitney U test was performed.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells12030393/s1, Figure S1: DESI-MSI mass spectra showed that OAG and SAG were detected in mice brain tissues by DESI-MSI in positive ion mode; Figure S2: CRS upregulated the 2-AG levels and downregulated OAG and SAG levels in the ACC, CP, NAc, and PIR in group #2 mice; Figure S3: CRS upregulated 2-AG levels and downregulated OAG and SAG levels in the ACC, CP, NAc, and PIR in group #3 mice.

References

1. Stella, N.; Schweitzer, P.; Piomelli, D. A second endogenous cannabinoid that modulates long-term potentiation. Nature; 1997; 388, pp. 773-778. [DOI: https://dx.doi.org/10.1038/42015]

2. Baggelaar, M.P.; Maccarrone, M.; van der Stelt, M. 2-Arachidonoylglycerol: A signaling lipid with manifold actions in the brain. Prog. Lipid Res.; 2018; 71, pp. 1-17. [DOI: https://dx.doi.org/10.1016/j.plipres.2018.05.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29751000]

3. Wilson, R.I.; Nicoll, R.A. Endocannabinoid Signaling in the Brain. Science; 2002; 296, pp. 678-682. [DOI: https://dx.doi.org/10.1126/science.1063545]

4. Bisogno, T.; Howell, F.; Williams, G.; Minassi, A.; Cascio, M.G.; Ligresti, A.; Matias, I.; Schiano-Moriello, A.; Paul, P.; Williams, E.J. et al. Cloning of the first sn1-DAG lipases points to the spatial and temporal regulation of endocannabinoid signaling in the brain. J. Cell Biol.; 2003; 163, pp. 463-468. [DOI: https://dx.doi.org/10.1083/jcb.200305129] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14610053]

5. Jung, K.M.; Astarita, G.; Zhu, C.; Wallace, M.; Mackie, K.; Piomelli, D. A key role for diacylglycerol lipase-alpha in metabotropic glutamate receptor-dependent endocannabinoid mobilization. Mol. Pharmacol.; 2007; 72, pp. 612-621. [DOI: https://dx.doi.org/10.1124/mol.107.037796]

6. Chen, C. Endocannabinoid metabolism and Alzheimer’s disease. Neural Regen. Res.; 2022; 17, pp. 1987-1988. [DOI: https://dx.doi.org/10.4103/1673-5374.335153]

7. Silveira, K.M.; Wegener, G.; Joca, S.R.L. Targeting 2-arachidonoylglycerol signalling in the neurobiology and treatment of depression. Basic Clin. Pharmacol. Toxicol.; 2021; 129, pp. 3-14. [DOI: https://dx.doi.org/10.1111/bcpt.13595] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33905617]

8. Bersani, G.; Pacitti, F.; Iannitelli, A.; Caroti, E.; Quartini, A.; Xenos, D.; Marconi, M.; Cuoco, V.; Bigio, B.; Bowles, N.P. et al. Inverse correlation between plasma 2-arachidonoylglycerol levels and subjective severity of depression. Hum. Psychopharmacol. Clin. Exp.; 2021; 36, e2779. [DOI: https://dx.doi.org/10.1002/hup.2779]

9. Bluett, R.J.; Báldi, R.; Haymer, A.; Gaulden, A.D.; Hartley, N.D.; Parrish, W.P.; Baechle, J.; Marcus, D.J.; Mardam-Bey, R.; Shonesy, B.C. et al. Endocannabinoid signalling modulates susceptibility to traumatic stress exposure. Nat. Commun.; 2017; 8, 14782. [DOI: https://dx.doi.org/10.1038/ncomms14782]

10. Rademacher, D.J.; Meier, S.E.; Shi, L.; Vanessa Ho, W.S.; Jarrahian, A.; Hillard, C.J. Effects of acute and repeated restraint stress on endocannabinoid content in the amygdala, ventral striatum, and medial prefrontal cortex in mice. Neuropharmacology; 2008; 54, pp. 108-116. [DOI: https://dx.doi.org/10.1016/j.neuropharm.2007.06.012]

11. Song, H.; Sieurin, J.; Wirdefeldt, K.; Pedersen, N.L.; Almqvist, C.; Larsson, H.; Valdimarsdottir, U.A.; Fang, F. Association of Stress-Related Disorders with Subsequent Neurodegenerative Diseases. JAMA Neurol.; 2020; 77, pp. 700-709. [DOI: https://dx.doi.org/10.1001/jamaneurol.2020.0117] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32150226]

12. Sanacora, G.; Yan, Z.; Popoli, M. The stressed synapse 2.0: Pathophysiological mechanisms in stress-related neuropsychiatric disorders. Nat. Rev. Neurosci.; 2022; 23, pp. 86-103. [DOI: https://dx.doi.org/10.1038/s41583-021-00540-x]

13. Fee, C.; Prevot, T.; Misquitta, K.; Banasr, M.; Sibille, E. Chronic Stress-Induced Behaviors Correlate with Exacerbated Acute Stress-Induced Cingulate Cortex and Ventral Hippocampus Activation. Neuroscience; 2020; 440, pp. 113-129. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2020.05.034] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32473277]

14. Patel, S.; Roelke, C.T.; Rademacher, D.J.; Hillard, C.J. Inhibition of restraint stress-induced neural and behavioural activation by endogenous cannabinoid signalling. Eur. J. Neurosci.; 2005; 21, pp. 1057-1069. [DOI: https://dx.doi.org/10.1111/j.1460-9568.2005.03916.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15787710]

15. Peng, X.L.; Wu, X.P.; Gong, R.X.; Yang, R.; Wang, X.; Zhu, W.Z.; Lin, P. Sub-regional anterior cingulate cortex functional connectivity revealed default network subsystem dysfunction in patients with major depressive disorder. Psychol. Med.; 2021; 51, pp. 1687-1695. [DOI: https://dx.doi.org/10.1017/S0033291720000434]

16. Misquitta, K.A.; Miles, A.; Prevot, T.D.; Knoch, J.K.; Fee, C.; Newton, D.F.; Ellegood, J.; Lerch, J.P.; Sibille, E.; Nikolova, Y.S. et al. Reduced anterior cingulate cortex volume induced by chronic stress correlates with increased behavioral emotionality and decreased synaptic puncta density. Neuropharmacology; 2021; 190, 108562. [DOI: https://dx.doi.org/10.1016/j.neuropharm.2021.108562]

17. Palacios, R.D.Y.; Verhoye, M.; Henningsen, K.; Wiborg, O.; Van der Linden, A. Diffusion Kurtosis Imaging and High-Resolution MRI Demonstrate Structural Aberrations of Caudate Putamen and Amygdala after Chronic Mild Stress. PLoS ONE; 2014; 9, e95077. [DOI: https://dx.doi.org/10.1371/journal.pone.0095077]

18. Lemos, J.C.; Wanat, M.J.; Smith, J.S.; Reyes, B.A.S.; Hollon, N.G.; Van Bockstaele, E.J.; Chavkin, C.; Phillips, P.E.M. Severe stress switches CRF action in the nucleus accumbens from appetitive to aversive. Nature; 2012; 490, pp. 402-406. [DOI: https://dx.doi.org/10.1038/nature11436]

19. Nacher, J.; Pham, K.; Gil-Fernandez, V.; McEwen, B.S. Chronic restraint stress and chronic corticosterone treatment modulate differentially the expression of molecules related to structural plasticity in the adult rat piriform cortex. Neuroscience; 2004; 126, pp. 503-509. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2004.03.038]

20. Dunham, S.J.B.; Ellis, J.F.; Li, B.; Sweedler, J.V. Mass Spectrometry Imaging of Complex Microbial Communities. Acc. Chem. Res.; 2017; 50, pp. 96-104. [DOI: https://dx.doi.org/10.1021/acs.accounts.6b00503]

21. Islam, A.; Takeyama, E.; Mamun, M.A.; Sato, T.; Horikawa, M.; Takahashi, Y.; Kikushima, K.; Setou, M. Green Nut Oil or DHA Supplementation Restored Decreased Distribution Levels of DHA Containing Phosphatidylcholines in the Brain of a Mouse Model of Dementia. Metabolites; 2020; 10, 153. [DOI: https://dx.doi.org/10.3390/metabo10040153]

22. Severiano, D.L.R.; Oliveira-Lima, O.C.; Vasconcelos, G.A.; Lemes Marques, B.; Almeida de Carvalho, G.; Freitas, E.M.M.; Xavier, C.H.; Gomez, M.V.; Pinheiro, A.C.O.; Gomez, R.S. et al. Cerebral Lipid Dynamics in Chronic Cerebral Hypoperfusion Model by DESI-MS Imaging. Neuroscience; 2020; 426, pp. 1-12. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2019.11.014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31785353]

23. Mamun, M.A.; Islam, A.; Hasan, M.M.; Waliullah, A.S.M.; Tamannaa, Z.; Huu Chi, D.; Sato, T.; Kahyo, T.; Kikushima, K.; Takahashi, Y. et al. The human vermilion surface contains a rich amount of cholesterol sulfate than the skin. J. Dermatol. Sci.; 2021; 103, pp. 143-150. [DOI: https://dx.doi.org/10.1016/j.jdermsci.2021.07.008] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34334257]

24. Islam, A.; Takeyama, E.; Nabi, M.M.; Zhai, Q.; Fukushima, M.; Watanabe, N.; Al Mamun, M.; Kikushima, K.; Kahyo, T.; Setou, M. Stress upregulates 2-arachidonoylglycerol levels in the hypothalamus, midbrain, and hindbrain, and it is sustained by green nut oil supplementation in SAMP8 mice revealed by DESI-MSI. Biochem. Biophys. Res. Commun.; 2022; 609, pp. 9-14. [DOI: https://dx.doi.org/10.1016/j.bbrc.2022.04.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35413542]

25. Egana-Huguet, J.; Bonilla-Del Rio, I.; Gomez-Urquijo, S.M.; Mimenza, A.; Saumell-Esnaola, M.; Borrega-Roman, L.; del Cano, G.G.; Salles, J.; Puente, N.; Gerrikagoitia, I. et al. The Absence of the Transient Receptor Potential Vanilloid 1 Directly Impacts on the Expression and Localization of the Endocannabinoid System in the Mouse Hippocampus. Front. Neuroanat.; 2021; 15, 645940. [DOI: https://dx.doi.org/10.3389/fnana.2021.645940]

26. Patel, S.; Kingsley, P.J.; Mackie, K.; Marnett, L.J.; Winder, D.G. Repeated Homotypic Stress Elevates 2-Arachidonoylglycerol Levels and Enhances Short-Term Endocannabinoid Signaling at Inhibitory Synapses in Basolateral Amygdala. Neuropsychopharmacology; 2009; 34, pp. 2699-2709. [DOI: https://dx.doi.org/10.1038/npp.2009.101]

27. McDevitt, D.S.; McKendrick, G.; Graziane, N.M. Anterior cingulate cortex is necessary for spontaneous opioid withdrawal and withdrawal-induced hyperalgesia in male mice. Neuropsychopharmacology; 2021; 46, pp. 1990-1999. [DOI: https://dx.doi.org/10.1038/s41386-021-01118-y]

28. Li, X.Y.; Chen, T.; Descalzi, G.; Koga, K.; Qiu, S.; Zhuo, M. Characterization of neuronal intrinsic properties and synaptic transmission in layer I of anterior cingulate cortex from adult mice. Mol. Pain; 2012; 8, 53. [DOI: https://dx.doi.org/10.1186/1744-8069-8-53] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22818293]

29. Parsons, L.H.; Hurd, Y.L. Endocannabinoid signalling in reward and addiction. Nat. Rev. Neurosci.; 2015; 16, pp. 579-594. [DOI: https://dx.doi.org/10.1038/nrn4004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26373473]

30. Berrendero, F.; Sepe, N.; Ramos, J.A.; Di Marzo, V.; Fernandez-Ruiz, J.J. Analysis of cannabinoid receptor binding and mRNA expression and endogenous cannabinoid contents in the developing rat brain during late gestation and early postnatal period. Synapse; 1999; 33, pp. 181-191. [DOI: https://dx.doi.org/10.1002/(SICI)1098-2396(19990901)33:3<181::AID-SYN3>3.0.CO;2-R]

31. Wittmann, G.; Deli, L.; Kallo, I.; Hrabovszky, E.; Watanabe, M.; Liposits, Z.; Fekete, C. Distribution of type 1 cannabinoid receptor (CB1)-immunoreactive axons in the mouse hypothalamus. J. Comp. Neurol.; 2007; 503, pp. 270-279. [DOI: https://dx.doi.org/10.1002/cne.21383]

32. Santori, A.; Morena, M.; Hill, M.N.; Campolongo, P. Hippocampal 2-Arachidonoyl Glycerol Signaling Regulates Time-of-Day- and Stress-Dependent Effects on Rat Short-Term Memory. Int. J. Mol. Sci.; 2020; 21, 7316. [DOI: https://dx.doi.org/10.3390/ijms21197316] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33023013]

33. Di Marzo, V.; Goparaju, S.K.; Wang, L.; Liu, J.; Batkai, S.; Jarai, Z.; Fezza, F.; Miura, G.I.; Palmiter, R.D.; Sugiura, T. et al. Leptin-regulated endocannabinoids are involved in maintaining food intake. Nature; 2001; 410, pp. 822-825. [DOI: https://dx.doi.org/10.1038/35071088] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11298451]

34. Hohmann, A.G.; Suplita, R.L.; Bolton, N.M.; Neely, M.H.; Fegley, D.; Mangieri, R.; Krey, J.F.; Walker, J.M.; Holmes, P.V.; Crystal, J.D. et al. An endocannabinoid mechanism for stress-induced analgesia. Nature; 2005; 435, pp. 1108-1112. [DOI: https://dx.doi.org/10.1038/nature03658] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15973410]

35. Wilson, R.I.; Nicoll, R.A. Endogenous cannabinoids mediate retrograde signalling at hippocampal synapses. Nature; 2001; 410, pp. 588-592. [DOI: https://dx.doi.org/10.1038/35069076] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11279497]

36. Suarez, J.; Ortiz, O.; Puente, N.; Bermudez-Silva, F.J.; Blanco, E.; Fernandez-Llebrez, P.; Grandes, P.; de Fonseca, F.R.; Moratalla, R. Distribution of diacylglycerol lipase alpha, an endocannabinoid synthesizing enzyme, in the rat forebrain. Neuroscience; 2011; 192, pp. 112-131. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2011.06.062] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21756982]

37. Sticht, M.A.; Lau, D.J.; Keenan, C.M.; Cavin, J.B.; Morena, M.; Vemuri, V.K.; Makriyannis, A.; Cravatt, B.F.; Sharkey, K.A.; Hill, M.N. Endocannabinoid regulation of homeostatic feeding and stress-induced alterations in food intake in male rats. Br. J. Pharmacol.; 2019; 176, pp. 1524-1540. [DOI: https://dx.doi.org/10.1111/bph.14453] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30051485]

38. Tsou, K.; Brown, S.; Sañudo-Peña, M.C.; Mackie, K.; Walker, J.M. Immunohistochemical distribution of cannabinoid CB1 receptors in the rat central nervous system. Neuroscience; 1998; 83, pp. 393-411. [DOI: https://dx.doi.org/10.1016/S0306-4522(97)00436-3]

39. Moldrich, G.; Wenger, T. Localization of the CB1 cannabinoid receptor in the rat brain. An immunohistochemical study. Peptides; 2000; 21, pp. 1735-1742. [DOI: https://dx.doi.org/10.1016/S0196-9781(00)00324-7]

40. Marsicano, G.; Lutz, B. Expression of the cannabinoid receptor CB1 in distinct neuronal subpopulations in the adult mouse forebrain. Eur. J. Neurosci.; 1999; 11, pp. 4213-4225. [DOI: https://dx.doi.org/10.1046/j.1460-9568.1999.00847.x]

41. Longaretti, A.; Forastieri, C.; Gabaglio, M.; Rubino, T.; Battaglioli, E.; Rusconi, F. Termination of acute stress response by the endocannabinoid system is regulated through lysine-specific demethylase 1-mediated transcriptional repression of 2-AG hydrolases ABHD6 and MAGL. J. Neurochem.; 2020; 155, pp. 98-110. [DOI: https://dx.doi.org/10.1111/jnc.15000]

42. Bedse, G.; Hill, M.N.; Patel, S. 2-Arachidonoylglycerol Modulation of Anxiety and Stress Adaptation: From Grass Roots to Novel Therapeutics. Biol. Psychiatry; 2020; 88, pp. 520-530. [DOI: https://dx.doi.org/10.1016/j.biopsych.2020.01.015] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32197779]

43. Mecha, M.; Carrillo-Salinas, F.J.; Feliu, A.; Mestre, L.; Guaza, C. Microglia activation states and cannabinoid system: Therapeutic implications. Pharmacol. Ther.; 2016; 166, pp. 40-55. [DOI: https://dx.doi.org/10.1016/j.pharmthera.2016.06.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27373505]

44. Nomura, D.K.; Morrison, B.E.; Blankman, J.L.; Long, J.Z.; Kinsey, S.G.; Marcondes, M.C.G.; Ward, A.M.; Hahn, Y.K.; Lichtman, A.H.; Conti, B. et al. Endocannabinoid Hydrolysis Generates Brain Prostaglandins That Promote Neuroinflammation. Science; 2011; 334, pp. 809-813. [DOI: https://dx.doi.org/10.1126/science.1209200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22021672]

45. Selvaraj, P.; Tanaka, M.; Wen, J.; Zhang, Y.M. The Novel Monoacylglycerol Lipase Inhibitor MJN110 Suppresses Neuroinflammation, Normalizes Synaptic Composition and Improves Behavioral Performance in the Repetitive Traumatic Brain Injury Mouse Model. Cells; 2021; 10, 3454. [DOI: https://dx.doi.org/10.3390/cells10123454]

46. Kirkham, T.C.; Williams, C.M.; Fezza, F.; Di Marzo, V. Endocannabinoid levels in rat limbic forebrain and hypothalamus in relation to fasting, feeding and satiation: Stimulation of eating by 2-arachidonoyl glycerol. Br. J. Pharmacol.; 2002; 136, pp. 550-557. [DOI: https://dx.doi.org/10.1038/sj.bjp.0704767]

47. Corcoran, L.; Mattimoe, D.; Roche, M.; David, P. Attenuation of fear-conditioned analgesia in rats by monoacylglycerol lipase inhibition in the anterior cingulate cortex: Potential role for CB2 receptors. Br. J. Pharmacol.; 2020; 177, pp. 2240-2255. [DOI: https://dx.doi.org/10.1111/bph.14976]

48. Okine, B.N.; Mc Laughlin, G.; Gaspar, J.C.; Harhen, B.; Roche, M.; Finn, D.P. Antinociceptive Effects of the GPR55 Antagonist CID16020046 Injected into the Rat Anterior Cingulate Cortex. Neuroscience; 2020; 443, pp. 19-29. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2020.07.013]

49. Viganò, D.; Grazia Cascio, M.; Rubino, T.; Fezza, F.; Vaccani, A.; Di Marzo, V.; Parolaro, D. Chronic morphine modulates the contents of the endocannabinoid, 2-arachidonoyl glycerol, in rat brain. Neuropsychopharmacology; 2003; 28, pp. 1160-1167. [DOI: https://dx.doi.org/10.1038/sj.npp.1300117]

50. Sustkova-Fiserova, M.; Charalambous, C.; Havlickova, T.; Lapka, M.; Jerabek, P.; Puskina, N.; Syslova, K. Alterations in Rat Accumbens Endocannabinoid and GABA Content during Fentanyl Treatment: The Role of Ghrelin. Int. J. Mol. Sci.; 2017; 18, 2486. [DOI: https://dx.doi.org/10.3390/ijms18112486]

51. Danan, D.; Todder, D.; Zohar, J.; Cohen, H. Is PTSD-Phenotype Associated with HPA-Axis Sensitivity?: The Endocannabinoid System in Modulating Stress Response in Rats. Int. J. Mol. Sci.; 2021; 22, 6416. [DOI: https://dx.doi.org/10.3390/ijms22126416] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34203952]

52. Aoki, J.; Isokawa, M.; Toyoda, M. Space and Time Coherent Mapping for Subcellular Resolution of Imaging Mass Spectrometry. Cells; 2022; 11, 3382. [DOI: https://dx.doi.org/10.3390/cells11213382] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36359779]

53. Sugiyama, E.; Guerrini, M.M.; Honda, K.; Hattori, Y.; Abe, M.; Kallback, P.; Andren, P.E.; Tanaka, K.F.; Setou, M.; Fagarasani, S. et al. Detection of a High-Turnover Serotonin Circuit in the Mouse Brain Using Mass Spectrometry Imaging. Iscience; 2019; 20, pp. 359-372. [DOI: https://dx.doi.org/10.1016/j.isci.2019.09.036] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31614319]