Urinary volatilome analysis in a mouse model of anxiety and depression

Psychiatric disorders including depression and anxiety comprise a broad range of conditions with different symptoms. We have developed a mouse model of depression/anxiety in mice deficient in the St3gal4 gene. In this study, we performed a comparative analysis of urinary volatile organic compounds (VOCs) in St3gal4-deficient (St3gal4-KO) and wild-type mice using gas chromatography-mass spectrometry, and we screened 18 putative VOCs. Principal component analysis (PCA) based on these VOCs identified a major group of 11 VOCs, from which two groups were clarified by hierarchical clustering analysis. One group including six VOCs (pentanoic acid, 4-methyl-, ethyl ester; 3-heptanone, 6-methyl; benzaldehyde; 5,9-undecadien-2-ol, 6,10-dimethyl; and unknown compounds RI1291 and RI1237) was correlated with the startle response (r = 0.620), which is related to an unconscious defensive response. The other group including two VOCs (beta-farnesene and alpha-farnesene) comprised pheromones which increased in KO mice. Next, male mice underwent a social behavior test with female mice in the estrus stage, showing reduced access of KO male mice to female mice. Comparative analysis of urinary VOCs before and after encounters revealed that the six VOCs were not changed by these encounters. However, in WT mice, the two farnesenes increased after the encounters, reaching the level observed in KO mice, which was not altered following the encounter. Taken together, these results indicated that St3gal4 was involved in modulating urinary VOCs. Moreover, VOC clusters discovered by comparison of St3gal4-KO mice with WT mice were correlated with differential emotional behaviors.


Introduction
Depression and anxiety are common mental disorders that co-occur frequently [1], and they have been identified as a leading cause of burden globally [2]. Depression is a psychiatric disorder that is characterized by a depressed mood and/or loss of interest or pleasure and causes severe symptoms that affect daily activities, such as sleeping, eating, or working. Anxiety is characterized by feelings of anxiety and fear [3]. In order to help patients who suffer from mental illness, it is important to deepen our understanding of the molecular mechanisms underlying psychiatric disorders and to investigate the potential mechanisms of action of drugs in patients with these disorders. Animal models, particularly rodent models, are valuable tools in this regard [4].
We previously reported that ST3 beta-galactoside alpha-2,3-sialyltransferase IV (St3gal4)deficient mice showed symptoms of anxiety and depression [5]. St3gal4 is one of the six a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 sialyltransferases (St3gal1-6) that generate Sia-2,3Gal linkages on the ends of glycoproteins. The correlation between enzymes involved in synthesis of sialoglycan and psychiatric disorders have been reported in humans. In humans, the GM3 synthase ST3GAL5 is associated with severe infantile-onset seizures, developmental delay, and blindness [6][7][8]. Moreover, mutations in ST3GAL3 result in nonsyndromic intellectual disability [9] and West syndrome [10]. Genome-wide association studies have demonstrated associations between ST3GAL4 gene variants or single-nucleotide polymorphisms and lipid traits, including high ApoB levels, total cholesterol, triglyceride levels, and coronary artery disease [11,12]. In mouse, St3gal4 expression is region specific, showing abundant expression in the thalamic sensory relay nuclei, including the medial and lateral geniculate, inter geniculate, and ventroposterior nuclei [5,13]. Additionally, St3gal4 is upregulated in neurons located within the neural circuits that exhibit propagation of kindling stimulation [13]. In contrast, St3gal4-knockout (KO) mice exhibit no epileptic seizures [5]. Furthermore, emotional behavioral tests, including open-field tests, forced swim tests, and fear conditioning tests, showed depression-and anxiety-like symptoms in these mice [5]. We found that kindling stimulation increased growth hormone (GH) levels in neuronal cells and upregulated St3gal4 [14]. In contrast, St3gal4-KO mice showed downregulation of GH and insulin-like growth factor 1 (Igf1) mRNA in the cerebral cortex [5]. Because GH and IGF1 are the main hormones involved in lipid metabolism, we investigated the relationships between St3gal4 and lipid metabolism using different types of oils to determine the effects of oils on the emotional behaviors of KO mice [15]. Changes in emotional behaviors generated by consumption of a diet containing oil differed between St3gal4-KO mice and their wild-type (WT) littermates. When St3gal4-deficient mice were fed pellets containing saturated fatty acids, the mice showed decreased tone fear, reaching the level of that in WT mice. Thus, these data suggested that St3gal4 may be involved in both lipid metabolism and neuropsychiatric disorders. The metabolic rate of mice is roughly 7 times higher than that of humans [16], and glucose metabolism per gram of brain tissue in mice is 2.9 times higher than that in humans [17]. These suggest that the higher metabolic rate in mice increases the loading to neural functions. The use of depression and anxiety model mice offers an advantage in studying the effects of metabolism on psychiatric conditions. Hence, we decided to investigate the relationship between metabolic changes that occur in St3gal4-KO mice.
Metabolomics is the systematic characterization of the repertoire of small molecules (metabolites) in a tissue or organism and has been widely used recently owing to methodological advances [18][19][20][21]. Studies have shown that the metabolic profiles of biological fluids reflect differences in the endocrine state, age, and genetic makeup of the donors [22][23][24][25]. Urinary profiles have been well characterized for the diagnosis of disease [26,27]. Urine, an excreted body fluid, is readily available and easily collected using noninvasive methods. Urine is rich in chemical information and is a potent source of pheromones [28,29]. Among the several thousand components in urine [30], volatile urinary compounds represent nearly all chemical classes, including aldehydes, alcohols, ketone, and hydrocarbons [31]. Alterations in urine volatiles have been used to detect urinary tract infections as well as bladder, prostate, and other cancers [32].
Previously, we demonstrated that differential peripheral lipid metabolism was involved in determining emotional and cognitive behaviors. Although various small molecules are produced as a result of lipid metabolism, the relationships between their production and emotion have not been examined. In this study, we aimed to investigate metabolic changes related to metabolism in St3gal4-KO mice by directly evaluating volatile organic compounds (VOCs) as end-metabolites using solid-phase micro-extraction (SPME) to extract urinary VOCs [33] followed by gas chromatography-mass spectrometry (GC-MS) to screen for differential expression of urinary VOCs between St3gal4-KO and WT mice. Overall, our results provided insights into the roles of VOCs in neuropsychiatry and a deeper understanding of emotional responses.

Ethics statement
All animals were treated in accordance with the Guidelines for Proper Conduct of Animal Experiments published by the Science Council of Japan (2006). The protocol was approved by the Committee on the Ethics of Animal Experiments of the Kyoto Sangyo University (approval no. 2017-08, 2018-07).

Animals
St3gal4-deficient mice (C57BL/6-St3gal4<tm1.1Bsi>, RBRC02286) were generated previously [5]. Briefly, St3gal4-KO and WT mice carried a C57Bl/6J background and were obtained by genome editing with MS12 embryonic stem cells (derived from C57Bl/6J mice) via electroporation and backcrossing with C57Bl/6J mice. The experiments were conducted using male mice of the 13 th to 20 th generations. Mice were housed in a temperature-controlled room (22-24˚C) with a 12-h light/dark cycle and were provided free access to food (MF; Oriental Yeast Co., Ltd, Tokyo, Japan.) and water. Four mice were maintained in open-top plastic cages (225 mm [width] × 338 mm [length] × 140 mm [height]). The cage tops were covered with stainless steel wire grid lids, and the cage floors were covered with paper-chips (Eco-chip; CLEA Japan, Inc., Tokyo, Japan), which were changed twice weekly. At the completion of the study mice were euthanized with isoflurane (Pfizer, NY, NY, USA).

Extraction of urinary volatile compounds by solid-phase microextraction
Urinary volatile compounds were extracted using headspace (HS)-SPME fibers. We used 2 cm 50/30 μM divinylbenzene/carboxen/polydimethylsiloxane (DVB/CAR/PDMS; Supelco, Bellefonte, PA, USA) as SPME fibers to cover a wide range of molecular weights. The fibers were preconditioned at 250˚C for 30 min and inserted into the headspace of the vial after the vial was equilibrated for 1 min at 45˚C. The volatile compounds in the headspace were extracted using SPME fibers at 45˚C for 60 min as previously described [35,36].

GC-MS
VOCs in middle-aged mice at 20-35 weeks were measured twice (experiment 1 and experiment 2) with a triple quadrupole GC-MS system (TQ-8040; Shimadzu, Kyoto, Japan). The SPME fiber with absorbed volatile compounds was inserted into the GC injection port and desorbed for 3 min at 240˚C. The injection was pulsed splitless for 3 min. To separate the desorbed volatiles, InertCap PureWAX with ProGuard and T.L. column (60 m + 10 m proguard line and 2 m transfer line, 0.25 mm internal diameter, 0.5 μm film thick; GL science, Tokyo, Japan) was used as previously described by Fujita et al. [36]. The oven temperature was programmed as follows: held at 40˚C for 10 min, ramped to 240˚C at 5˚C/min, and held for 10 min. The injector temperature was set at 240˚C. Helium was used as carrier gas at a constant column flow rate of 2 mL/min for 0-3 min and 1 mL/min for 3-60 min. The operating parameters for MS were as follows: ion source temperature, 200˚C; ionizing energy, 70 eV; scanning frequency, 0.3 s/spectrum from 35 to 300 m/z for experiment 1 and 0.2 s/spectrum for experiment 2. For urinary VOC analysis in mouse at 10 to 15 weeks of age, the procedure was the same as in experiment 2 except that the constant linear velocity of helium was changed to 20 cm/s by QP-2010 Ultra (Shimadzu). VOC identification was carried out based on comparisons of retention times with those of standards and with probability matching against a mass spectral library (NIST/EPA/NIH mass spectral library, NIST14; Mass Spectral Library of Drugs, Poisons, Pesticides, Pollutants, and Their Metabolites, Wiley10). The retention index (RI) of VOCs were calculated from retention time of the n-alkane standard and searched the data in NIST Chemistry WebBook [37]. The qualities of commercial texanol were analyzed using 1H NMR spectra of compounds recorded on a JEOL JNM-ECS 400 spectrometer with TMS as the internal standard [38] as texanol included two compounds, texanol and texanol isomer.

Data processing and statistical analyses
The chromatographic peak areas were integrated using GC-MS Solution ver.4.45 software (Shimadzu). Detection and integration of the generated ion peaks from electron ionization was performed using XCMS software package version 1.3.2 (http://masspec.scripps.edu) [39,40] in R version 3.2.3 (http://cran.r-project.org/). XCMS software was used to process peak-matching, non-linear retention time alignment, and quantitation of mass spectral ion intensities. The total ion current (TIC) for each peak was manually inspected to validate the detected peak. Refinement of VOCs was carried out by manual deconvolution based on the retention times and peak shapes of the various extracted ion peaks (m/z).
First, differential urinary VOCs between St3gal4-KO and WT mice at 20-35 weeks of age were screened. VOCs in mice were measured with GC-MS three times. Triplicate urinary VOCs in St3gal4-KO (n = 3 in experiments 1; n = 3 in experiment 2) and WT (n = 3 in experiment 1; n = 3 in experiment 2) mice were examined by XCMS to detect differential ion peaks (Table 1). Finally, the area of the ion peaks was quantified using GC-MS Solution software. To compile each dataset (experiment 1 and experiment 2), we used the average of three absolute ion peak areas to calculate the relative intensity of the average area of KO mice relative to that e) The NIST library identified 4-octen-3-one, 6-ethyl-7-hydroxy (ketone) as a TIC, and the maximum SI with the MS spectrum was 87%. However, RI information was not obtained from the NIST chemistry book.
f) The NIST library identified 2-hexenal, 2-ethyl-(aldehyde) as a TIC, and the maximum SI with the MS spectrum was 87%. However, the retention time was not the same as that of the commercial standard. g) The NIST library identified 5-hexen-2-one, 5-methyl-/3-hexen-2-one, 5-methyl-(ketone) as a TIC, and the maximum SI with the MS spectrum was 90%. However, RI information was not obtained from the NIST chemistry book.
h) The NIST library identified 5-oxohexanenitrile (ketone) as a TIC, and the maximum SI with the MS spectrum was 90%. However, the retention time was not the same as that of the commercial standard. The second hit from the NIST library was 5-hepten-2-one, and the SI with the MS spectrum was 89%. However, RI information was not obtained from the NIST chemistry book.
i) The NIST library identified 3-heptanone, 5-methylene (ketone) as a TIC, and the maximum SI with the MS spectrum was 90%. However, RI information was not obtained from the NIST chemistry book. j) n.d., not determined. Since no ion peaks were detected in the urine of St3gal4-KO mice and WT littermates in experiment 2, the fold-change, SEM, and p value were not calculated. k) The NIST library identified eicosyl acetate (hydrocarbon) as a TIC, and the maximum SI with the MS spectrum was 95%. However, RI information was not obtained from the NIST chemistry book. l) The chemical class is shown as having a hydroxy group with only carbon-to-carbon double bonds as unsaturation, in Pubchem [PubChem Identifier: CID 5370125]. † VOCs were confirmed by identification using commercial standard references. # The NIST library identified 2,2,4-trimethyl-1,3-pentanediol diisobutyrate (TXIB) with 94% SI as a TIC. However, the MS spectrum of compound (16) did not include m/z 159 and m/z 243, which were detected in commercial TXIB by GC-MS. On the other hand, the commercial texanol standard included two peaks; one represented texanol, and the other showed the same retention time as compound (16) with a similar MS spectrum. Based on NMR analysis, the peak in the commercial texanol standard corresponding to compound (16) was determined to be texanol isomer, which we were unable to find in the NIST library.
of the WT littermates (n = 6 pairs). The absolute values of ion peak area from urinary VOCs were compared between young KO and WT mice (male, 10-16 weeks old). Differential levels of VOCs were analyzed using Excel and Prism 6 by Mann Whitney U-tests and Kruskal-Wallis tests. Among multivariate statistical analyses, principal component analysis (PCA) score plots were obtained from relative values of 17 VOCs in six pairs of St3gal4-KO mice and WT littermates using Excel, KyPlot 5.0 (KyensLab Inc., Tokyo, Japan), and IBM SPSS Statistics Version 25 (Armonk, NY, USA). Hierarchical cluster analysis was conducted using absolute values of ion peak areas with IBM SPSS Statistics Version 25. A polynomial regression curve to evaluate whether VOCs were correlated with startle responses was established using a combination of Excel and IBM SPSS Statistics Version 25.

Startle response
The acoustic startle response was tested in mice to observe the sensorimotor fright response using a startle reflex test unit for mice (O'Hara & Co., Ltd., Tokyo, Japan) [41]. A mouse was inserted to an acrylic chamber, the chamber was placed onto a sensor unit, and both the chamber and sensor unit were fixed using a clamp. The speakers conveyed continuous background noise of 70 dB (white-noise) for a 5-min acclimation period and throughout the test session. The mouse was then exposed to a 120-dB noise pulse for 20 ms per 500 ms 10 times. Vibrations caused by the whole-body startle response of the mouse were displayed numerically with the acceleration sensor. The numerical signal was amplified and converted into the voltage, which was displayed as the shape of waves. The peak-to-peak value of this wave was adopted as the strength of the startle response.

Open field test
Open field tests for encounters with females were conducted on St3gal4-KO and WT male mice at 10-16 weeks of age. Each of the animals was acclimatized to being alone in a separate cage in a waiting room for at least 30 min before the behavioral sessions. Behavioral tests were performed between 0800 and 1300 by an experimenter blinded to the mouse genotype. The open field test was performed by placing the mice in a square-shaped arena (490 mm wide × 490 mm long × 195 mm high) that had a ceiling height 125 cm above the field floor. The light intensity of the fluorescent lamp that lit the open field and the acclimatization waiting room with the separate cage was set to 70 lux.
Female mice (ddY, 8-12 weeks old) underwent vaginal smear tests just before the open field tests to determine estrous cycle stages [42]. The smears were stained with Giemsa solution (Wako Pure Chemical Industries, Ltd.) and examined by Eclipse 80i bright-field microscopy (Nikon, Tokyo, Japan).
A male mouse was first placed in the field, and 1 min later, a female mouse in the proestrus period was placed in the field for 30 min. Both mice were then returned to their home cages. Data were collected with a camera placed 470 mm above the field floor using TimeOFCR4 (O'Hara & Co., Ltd.). After recording, the numbers of male accesses per 1 s from male nose to female nose, from male nose to female circumanal region, and from male hand to female body were counted for the first 10 min.

Specific differences in urinary VOCs in St3gal4-KO mice
Urine samples from middle-aged mice at 20-35 weeks old (n = 3, KO, n = 3, WT in experiment 1; n = 3, KO, n = 3, WT in experiment 2) were analyzed using headspace-SPME GC-MS. The representative TIC peaks and list of VOCs estimated by a similarity index above 85% from experiment 1 and experiment 2 are shown in Figs 1 and S1. Forty-two VOCs were detected in experiment 1, and 75 VOCs were detected in experiment 2 (S1 and S2 Tables). These data indicated that the urine of C57Bl/6J mice contained at least 97 VOCs.
To investigate which screened VOCs were related to emotional behaviors, we performed the startle test and social open field test and analyzed urine samples collected before and after these social behavior tests from St3gal4-KO (n = 9) and WT (n = 11) male mice at 10-16 weeks of age (Table 2). In the young male mice, RI2297 (18) was undetectable. The remaining VOCs were used for further analysis. To evaluate kidney glomerular filtration, we measured urinary creatinine, which exhibited minimal differences between St3gal4-KO and WT mice.

Correlations between startle response and production of urinary VOCs
Next, to analyze correlations between urinary VOC quantities and emotional response, we used a startle response test. First, we confirmed there were few differences in auditory activity or startle response between St3gal4-KO and WT mice (p = 0.1882 by Mann-Whitney U-tests in Fig 3A). Urine was collected after the startle response test. Second, to investigate the correlations between urinary VOC quantities and startle response, PCA scores using absolute values of the 'six VOCs' classified in Fig 2B were constructed. Because the percentage of variance contribution in the 'six VOCs' was captured 89.8% by PC1, which was very high, PC1 and startle response were distributed on the XY axis, resulting in a second-order polynomial regression curve with high fitting (r = 0.886; Fig 3B). Thirty-three percent of St3gal4-KO mice produced higher amounts of the 'six VOCs' in the urine and exhibited strong correlations between the production of the 'six VOCs' and startle response, whereas the remaining 67% of KO mice showed low amounts of the 'six VOCs'. There was little correlation between production of the 'six VOCs' and the startle response in WT mice (Fig 3B).

Correlation between access of male mice to female mice and production of urinary VOCs
To analyze another correlation between urinary VOC quantities and emotional responses, we performed social behavior tests (Fig 4A). The access of sexually naïve male mice to female mice at diestrus did not differ between WT and KO mice. The number of accesses of WT male mice to female mice at the proestrus/estrus (P/E) stage was increased 2.0-fold (means) compared with accesses to females at the diestrus stage, while accesses of St3gal4-KO male mice did not change regardless of female stage. This indicates that St3gal4-KO males failed to access females, even when female mice were in the P/E stage. Hence, we analyzed the effects of encounters with female mice at the P/E stage on male urinary VOC quantities to observe the emotional response.
https://doi.org/10.1371/journal.pone.0229269.g002 Table 2. Absolute values of ion peaks in male urinary VOCs before and after encountering female mice at the proestrus stage.  VOCs in male urine samples collected before and after encounters with females at the proestrus/estrus stages were compared between St3gal4-KO and wild-type mice at 10-16 weeks of age. The order number (NO.) is the same as that in Table 1.

No. VOCs
As orders of areas differed among VOCs, the number of significant digits was unified to five columns. # Slight difference in quantity in male urine after encounters with females at the proestrus stage. $ Increase in VOC quantity after encounters with females.
a) The ion peak (m/z) extracted in Table 1 was used to calculate the absolute area of the ion peak. Mann Whitney U-test: 2T, two-tailed; 1T, one-tailed. Two-tailed (one-tailed) p values were used to evaluate significant correlations. c) Urinary creatinine concentrations were determined by LabAssay TM Creatinine colorimetry kits based on the Jaffé method [34].
https://doi.org/10.1371/journal.pone.0229269.t002 rate [43] (S5 Table). The trends by two multiple comparisons tests were similar. We found that there were three patterns of VOCs showing significant correlations with the number of male accesses. The first was RI1101 (3), β-farnesene (13), and α-farnesene (14), which were present at high levels in the urine of KO mice compared to those in WT mice before access (Mann-Whitney U-test, two-tailed p = 0.0043, one-tailed p = 0.0411 and 0.0260, respectively). After encountering female mice, these VOCs increased in WT male mice to the levels observed in KO mice (green in Figs 4D, 4M, 4N and 5 and Table 2), while the levels in KO mice after the encounters were the same as before the encounters. The second group was styrene (9), texanol (15), and texanol isomer (16), the levels of which changed after encounters with female mice in both St3gal4-KO and WT male mice (p = 0.0283, 0.0042, and 0.0154, respectively; blue in Figs  Table 2). Styrene (9) decreased after encounters with female mice in the urine of both St3gal4-KO and wild-type male mice. Texanol (15) increased 2.0-and 2.6-fold after encounters with female mice in the urine of St3gal4-KO and WT mice, respectively, whereas texanol isomer (16) increased 1.4-and 1.7-fold, respectively.
Notably, the third group, termed the 'six VOCs', which were identified by hierarchical cluster analysis, were correlated with the startle response in the urine of St3gal4-KO mice and were higher than those in WT mice, regardless of encounters with female mice (red in Figs 4F, 4G, 4I, 4K, 4L, 4Q and 5 and Table 2). In other words, encounters with female mice did not elevate the amounts of the 'six VOCs'. Furthermore, the absolute values of ion peaks for the 'six VOCs' from WT mice showed low variance (Fig 4F, 4G, 4I, 4K, 4L and 4Q).

Identification of 18 VOCs
In the present study, we sought to determine whether St3gal4-KO mice showed variations in volatile urinary chemicals. These VOCs were extracted at 45˚C for 60 min in the headspace and were concentrated on SPME fiber. Within these conditions we were unable to detect any non-volatile urinary compounds such as the major urinary proteins (MUPs) secrete by mice. These proteins harvest and sequester urinary nonpolar VOCs in hydrophobic pockets [44] and function as transporters and stabilizers of pheromones in rodents [45,46]. However, the ability to detect VOCs bound to MUPs may be limited in the current study and thus, we only effectively identified those that were unbound. Nevertheless, via GC-MS analysis, we detected 97 VOCs from middle-aged St3gal4-KO and WT mice. Among these VOCs, 18 VOCs were screened as significantly changed in St3gal4-KO mice by XCMS analysis.

Characterization of six VOCs
To examine the influence of the 'six VOCs' in the auditory thalamic sensory relay nuclei, where St3gal4 is specifically expressed [13], we investigated the startle response, an unconscious defensive response to threatening stimuli [41] and analyzed urinary VOCs collected after the startle experiment. These 'six VOCs' were elevated in St3gal4-deficient mice, and showed a correlation with the startle response in St3gal4-KO mice. The startle response is related to the auditory pathway, in which neural circuits and sound information is transmitted to the cochlear nucleus (CN). The auditory neurons in the CN ascend to the superior olivary complex (SO), and the neurons of the SO send axons to the inferior colliculus (IC) of the midbrain. The neurons of IC send axons to the medical geniculate body (MGN), and the axons of the MGN send axons to the auditory cortex [47]. The deficiency of St3gal4, which is specifically expressed in the MGN, appeared to influence the expression of the 'six VOCs'. Further, the abundance of urinary VOCs showed variability in each St3gal4-KO mouse and exhibited large variance in Fig 4 in addition to Fig 3A. This high variability implies that St3gal4 deletion leads to dysregulation of the 'six VOCs' metabolism. Additionally, a portion of the St3gal4-KO mice that exhibited a startle response correlated with 'six VOCs' abundance, causing variability in Fig 3A. These suggested that the amounts of VOCs synchronously progress to out of regulation with behaviors in a part of St3gal4-KO mice.
Moreover, a previous study reported that within a bacterial metabolic network, functionality and constraints drive the tradeoff between robustness and fragility with high variability [48]. St3gal4 is a sialyltransferase and in humans SNPs or gene variants of ST3GAL4 affect lipid traits as well as the subsequent development of coronary artery disease [11,12]. Thus, further investigation is required to better understand how lipid metabolism alterations induced by St3gal4 deficiency affect the high variability of metabolic networks.

Classification of VOCs by genotype and social behaviors in young mice
The number of accesses of male KO mice to female mice was unchanged by estrous stage, whereas the number of accesses of male WT mice to female mice in the P/E stage was increased 2.2-fold compared with that in the diestrus stage. We investigated the effects of social behavior in WT and KO mice on urinary VOCs. Encounters with female mice altered the amounts of urinary VOCs in naïve male mice. These VOC quantities effectively separated groups by genotype as well as by those from before and after social behavior by PCA. β- The vertical line shows the number of male accesses at the P/E and diestrus stages per 1 s for 10 min. Different quantities of urinary VOCs were compared between St3gal4-KO and WT mice and between urine collected before and after encounters with female mice (B-Q). Urine was collected for 1 week before and after encounters with female mice. On the Y-axis, red indicates that St3gal4-KO mice exhibited higher amounts of VOCs in the urine than WT mice with low variance; green indicates that there were low amounts of VOCs in WT mice before encounters with female mice; and blue indicates that the concentrations of VOCs were changed in both St3gal4-KO and WT mice after encounters with female mice. Kruskal-Wallis test and Dunn's multiple comparisons tests were used. WT, wild-type; KO, St3gal4-deficient; P/E, proestrus/estrus; D, diestrus stages. https://doi.org/10.1371/journal.pone.0229269.g004 Farnesene and α-farnesene levels in group 2 were higher in KO mice than in WT mice before social interaction, whereas there was little difference after social interaction. These compounds are pheromones that are produced in males and have roles in attraction to females and puberty acceleration [49,50]. Our data showed that the amounts of urinary farnesenes were regulated in WT mice, but not KO mice, as a response to social behaviors. Encounters with female mice increased levels of group 3 compounds texanol and texanol isomer levels, regardless of genotype, while decreasing urinary styrene (Figs 4O and 4P and 5). It is surprising that these nonpheromonal VOCs were influenced by emotional behaviors. The 'six VOCs' were not correlated with the access of sexually naive male mice to female mice at the P/E stage. This suggests that urinary VOCs are differentially generated in response to emotional behaviors. Thus, our findings indicate that end metabolic products reflect emotional consequences.
It remains unclear how the levels of these VOCs are altered in response to emotional behavior. Further studies are also needed to determine how VOCs are produced and are involved in emotional behavior. VOCs including pheromones are likely to influence to emotion via the olfactory system, including the anterior olfactory nucleus, olfactory tubercle, piriform cortex, and cortical amygdala [51]. Table 3 lists the VOCs that were screened in St3gal4-KO mice, along with their chemical classes and whether they have been reported in mice and humans. These chemical classes include ketones, aldehydes, and hydrocarbons. Specifically, the lipid metabolism-related VOCs presented in Tables 1 and 2 included: pentanoic acid, 4-methyl, ethyl ester (5); 3-heptanone, 6-methyl (6); or β-farnesene (13) and α-farnesene (14), which were classified as ester, ketone, or hydrocarbon, respectively. Ketones are generated from the catabolism of medium chain fatty acids. Ketones in urine are likely partially derived from bacterial action in the gut. In fact, Holland et al found that ketones were produced in conventional  Table 1. Names of VOCs written in red, green, and blue are identical to those written on the vertical lines in Fig 4. https://doi.org/10.1371/journal.pone.0229269.g005 Urinary volatiles and emotional behaviors urine yet were decreased in germ free mice [52]. Since urine was not collected from germ free mice in the current study, ketones may have been produced. Hydrocarbons produced via fatty acids are the products of lipid peroxidation. Alpha-farnesene and beta-farnesene are known to be produced from acetyl-CoA via the mevalonate pathway of cholesterol synthesis [53]. Additionally, 5,9-undecadien-2-ol, 6,10-dimethyl has been described in the Pubchem database as unsaturated with only carbon-carbon double bonds and a hydroxy group [PubChem Identifier: CID 5370125, 54]. It is also produced via the mevalonate pathway [55,56]. Therefore, 5,9-undecadien-2-ol, 6,10-dimethyl is included as lipid metabolism-related VOC in Tables 1  and 2. Since VOCs are the end products of endogenous and exogeneous pathways they often serve as biomarkers. In fact, multiple studies have identified different VOCs as markers of increased/dysregulated lipid metabolism [27,57,58]. Here it is inferred that ST3Gal IV influences VOC production by altering lipid metabolism.
Furthermore, identified groups via PCA and hierarchical cluster analysis were divided by their functions. In group 1, the 'six VOCs' included pheromones: 5,9-undecadien-2-ol, 6,10-dimethyl and benzaldehyde are known as insect pheromones [59,60,61]; 3-heptanone 6-methyl, a derivative of 6-hydroxy-6-methyl-3-heptanone which is a known mouse pheromone [62]. The 'six VOCs' in group 1 were found to be correlated with the startle test. Alternatively, in group 2, mice pheromones classified as hydrocarbons were included, namely alphafarnesene and beta-farnesene, which were altered in St3gal4-KO mice and following social behavior. Hence, St3gal4 affected the expression of pheromones that were dependent on differential behaviors. Lastly, in group 3, texanol and texanol isomer were observed, which are common plasticizers in water-based paints and certain plastic materials [63]. Hence, in young mice, social behaviors may have affected the metabolisms of exogenous compounds regardless of St3gal4-KO/WT. Taken together, the proposed metabolic pathways from the VOC chemical classes as well as the observed relationships between VOCs and emotions, suggest an Urinary volatiles and emotional behaviors association between metabolism and emotional behavior. Identification of the biological pathways of these metabolites may provide insights into the relationships between peripheral metabolites and psychiatric effects.  Table. Volatile organic compound (VOC) names, similarity indexes, chemical formulas, CAS nos., and molecular weights, as analyzed by GC-MS using an InertCap PureWAX column under the operating parameters of experiment 1. � Compounds that are also shown in S2 Table are labeled. Operation parameters for the mass spectrometer in experiment 1 were described in the Materials and Methods. (DOCX) S2 Table. VOC names, similarity indexes, chemical formulas, CAS nos., and molecular weights, as analyzed by GC-MS using an InertCap PureWAX column under the operating parameters of experiment 2. � Compounds that are also shown in S1