Next Article in Journal
Chloride Channel Family in the Euhalophyte Suaeda altissima (L.) Pall: Cloning of Novel Members SaCLCa2 and SaCLCc2, General Characterization of the Family
Next Article in Special Issue
Deciphering the Effect of Different Genetic Variants on Hippocampal Subfield Volumes in the General Population
Previous Article in Journal
Circulating Histones to Detect and Monitor the Progression of Cancer
Previous Article in Special Issue
Macrophage Migration Inhibitory Factor in Major Depressive Disorder: A Multilevel Pilot Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relationship between Urinary Metabolomic Profiles and Depressive Episode in Antarctica

1
Department of Gastrointestinal Surgery, Tokyo Medical University, Tokyo 160-0023, Japan
2
National Institute of Polar Research, Tokyo 190-8518, Japan
3
Department of Breast Oncology and Surgery, Tokyo Medical University, Tokyo 160-0023, Japan
4
Institute of Medical Sciences, Tokyo Medical University, Tokyo 160-0022, Japan
5
Institute for Advanced Biosciences, Keio University, Yamagata 997-0052, Japan
6
Department of Psychiatry, Tokyo Medical University, Tokyo 160-0023, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(2), 943; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24020943
Submission received: 10 December 2022 / Revised: 29 December 2022 / Accepted: 2 January 2023 / Published: 4 January 2023
(This article belongs to the Special Issue Molecular Research on Depression)

Abstract

:
Antarctic expeditions have a high risk of participant depression owing to long stays and isolated environments. By quantifying the stress state and changes in biomolecules over time before the onset of depressive symptoms, predictive markers of depression can be explored. Here, we evaluated the psychological changes in 30 participants in the Japanese Antarctic Research Expedition using the Patient Health Questionnaire-9 (PHQ-9). Urinary samples were collected every three months for a year, and comprehensive urinary metabolomic profiles were quantified using liquid chromatography time-of-flight mass spectrometry. Five participants showed major depressive episodes (PHQ-9 ≥ 10) at 12 months. The urinary metabolites between these participants and the 25 unaffected participants were compared at individual metabolite and pathway levels. The individual comparisons showed the most significant differences at 12 months in 14 metabolites, including ornithine and beta-alanine. Data from shorter stays showed less significant differences. In contrast, pathway and enrichment analyses showed the most significant difference at three months and a less significant difference at longer stays. These time transitions of urinary metabolites could help in the development of urinary biomarkers to detect subjects with depressive episodes at an early stage.

1. Introduction

People are resistant and adaptable to various environmental conditions. In particular, expeditioners in the Antarctic are exposed to considerable stress through isolation and extreme physical conditions [1], leading to unique psychological effects. In addition to the cold, danger, and hardship, the difficulties associated with individual adjustment to the group constitutes a major stressor [2], with members of Antarctic expeditions undergoing psychological changes resulting from exposure to long periods of isolation and confinement [3]. The psychological effects of polar expeditions in patients under abnormal light, cardiopulmonary, and cold conditions have been frequently investigated [3] since the first studies were performed in Polar Regions [2]. Such patients experienced alterations in peripheral circulation along with hypothermia and frostbite, in additional to immune system and hormonal abnormalities. Their symptoms included somatic and gastrointestinal complaints, fatigue, headaches, disturbed sleep, impaired cognition, depressed mood, anger and irritability, anxiety, and interpersonal tension and conflict [3]. Approximately 5% of the South Pole expedition participants fulfilled the criteria for psychiatric disorders as defined in the Diagnostic and Statistical Manual of Mental Disorders IV (DSM-IV) [4,5]. Therefore, the appropriate assessment of psychological status is necessary to prevent these issues.
The identification of psychological risk factors that have a negative effect on mood and behavior is desirable for designing efficient preventive approaches to reduce symptoms and generate satisfactory adaptation [6]. Seasonal patterns in depressive symptoms, such as winter-over syndrome, polar T3 syndrome, and subsyndromal seasonal affective disorder, have been identified [7,8,9]. The psychological profiles of the Japanese Antarctic Research Expedition (JARE) team over three consecutive years showed psychological changes at the end of winter, such as an increase in the planning trial score and a decrease in the body coldness score, despite showing high robustness against stress [10]. Psychological changes in the JARE team evaluated using a psychological scale based on an international comparison project of psychology called the Polar Psychology Project indicated that morale was high at the beginning of the winter season, fell into a midwinter slump during the polar night, and then recovered [11]. Together, these research reports reveal that not only assessing psychological changes but also appropriate early care is necessary to reduce the effects of psychological risks during Antarctic expeditions.
Molecular biomarkers have been identified to predict and quantitatively assess altered mental and immune status. A long-term stay in a confined environment increases stress; therefore, the development of objective molecular markers to assess mental health is required. An analysis of the changes in already established stress biomarkers, such as cortisol and dehydroepiandrosterone (DHEA), in researchers at Antarctic stations revealed that DHEA levels could predict mental health status [12]. The profiling of thyroid-stimulating hormone (TSH), melatonin, and cortisol during the Antarctic expedition suggested that different luminosity and temperature conditions increase psychological stress [13]. The profiling of blood cytokine and interleukin levels showed reduced T-cell function [14]. Altered cytokine production in the blood has also been observed [15,16]. Moreover, total serum thyroxine levels are significantly associated with tension anxiety, depression, and total mood disturbances [17]. Nevertheless, the evaluation of novel biomarkers to determine the association with their mental health in the extreme conditions associated with Antarctic expeditions remains necessary.
Recent omics technologies, which enable the simultaneous identification and quantification of thousands of molecules, have been used to explore new biomarkers associated with mental health. Among molecular profiling approaches, the large-scale analysis of metabolites via metabolomics has been recently used for biomarker discovery in the psychiatry field [18]. Owing to the large diversity of the chemical features of metabolites, a single approach cannot profile all metabolites; thus, the combination of individual separation systems prior to detection has been applied. For example, gas chromatography–mass spectrometry (GC/MS) was used to profile urinary metabolites to predict postpartum [19] and poststroke depression [20]. A combination of GC/MS and nuclear magnetic resonance (NMR) has been used to explore urinary metabolite biomarkers for major depressive disorder (MDD) and bipolar disorder (BD) [21,22,23]. Capillary electrophoresis-MS and liquid-chromatography-MS (LC-MS) were used to profile plasma hydrophilic metabolites in patients with MDD [24,25]. However, these studies were designed as cross-sectional studies, whereas a longitudinal study is necessary to evaluate the trend of changes in biomarker levels accompanied by mental status in order to realize the early identification of mental changes.
In this study, we aimed to elucidate the urinary metabolite biomarkers that have potential to discriminate subjects with depressive episodes from those without these symptoms. Here, we investigated the time course of urinary metabolites in the Japanese Antarctic Expedition team and evaluated their association with depressive episodes. Samples were collected every three months over a one-year period in Antarctica. Comprehensive metabolomic analyses were conducted with LC-MS. The individual metabolites and pathway-level changes associated with the transition of depressive episodes were identified.

2. Results

Thirty participants were enrolled in the study (Table 1). Figure 1 shows the time-series changes in the Patient Health Questionnaire-9 (PHQ-9) scores. Five participants showed a high PHQ-9 score (≥10) at the end of the 12-month stay, and were classified as the depressive episode group. The remaining 25 participants were included in the nondepressive episode group. At baseline (0-month stay), the PHQ-9 values were not significantly different between the two groups. However, after the third month of stay, a significant difference was observed between the groups (p < 0.05, Mann–Whitney U-test). Subsequently, the PHQ-9 values of the depressive episode group increased monotonically and were significantly higher at 9 and 12 months (p < 0.0001).
Hydrophilic urinary metabolites were identified and quantified using LC-MS-based metabolomic analysis. In total, 79 metabolites were frequently detected (≥95% of all samples), and 77 metabolites were used for subsequent analyses after excluding two insufficiently separated metabolites in LC analysis. The deviation of the peak areas among the quality control samples was 9.97 ± 7.09%. Score plots of the principal component analysis were used to reveal the overall similarity and discrepancy of the quantified metabolite patterns among the samples (Figure A1). The score plots for the samples collected at baseline were higher than those of the other samples. This indicates a large change in the metabolite concentration patterns of the baseline samples. However, these changes are considered to be influenced by differences in the sampling conditions and storage, as only the baseline samples were collected in Japan. Therefore, baseline data were excluded from subsequent analyses, and time-series data for the remaining one-year stays in Antarctica were compared.
Metabolites with significant differences (p < 0.05, Mann–Whitney test) of at least one point among the four data points (3- to 12-month stays) between the depressive and nondepressive episode groups were depicted in a heatmap (Figure 2). The time-course data are shown in Figure 3. Most metabolites exhibited significant differences in the data of the 12-month stay, including 2-aminoadipate, 3-methylhistidine, 4-guanidinobutanoate, alpha-aminobutyrate, beta-alanine, carnosine, cystathionine, cytosine, N-acetylhistamine, ornithine, and pyridoxamine. All of these metabolite concentrations were significantly lower in the depressive episode group compared with those of the nondepressive episode group at the 12-month stay (Figure 3a). Among these, only 4-guanidinobutanoate showed a significant difference in the false discovery rate (FDR)-adjusted p < 0.05. In addition, 2-aminoadipate levels were significantly lower in both the 9- and 12-month stays (Figure 3b). Guanine, putrescine, and indoleacetate showed significant differences at different time points during the 12-month stay.
Pathway and enrichment analyses were also performed between the depressive episode and nondepressive episode groups to evaluate the altered pathways. The results of the pathway analysis are shown in Figure 4. Each plot represents a single pathway. A higher Y-value (−log10(P)) indicates a smaller p-value, which means that more alternations were present at the pathway level. The higher X-axis value (pathway impact) indicates that the alteration of the mapped metabolites affects the pathway more widely. Therefore, pathways with higher values on both the X- and Y-axes were more completely altered. The three-month-stay data included five pathways that were significantly altered (p < 0.05): (1) nitrogen metabolites; (2) pyrimidine metabolism; (3) glyoxylate and dicarboxylate; (4) D-glutamine and D-glutamate metabolism; and (5) alanine, aspartate, and glutamate (Figure 4a). Pathways (1) and (2) showed smaller pathway impact values, whereas pathways (4) and (5) had relatively higher pathway impact values. Thiamine metabolism was significantly different at the six-month stay; however, the pathway impact was small (Figure 4b). No significant difference in the data was observed at the nine-month stay (Figure 4c). Propanoate metabolism was significantly different in the data at the 12-month stay; however, the pathway impact was small (Figure 4d). Detailed data are shown in Figure A2, Figure A3 and Figure A4.
An enrichment analysis was also conducted to evaluate the differences between the depressive and nondepressive episode groups (Figure 5). Five pathways at the three-month stay showed significant differences: (1) D-glutamine and D-glutamate metabolism; (2) nitrogen metabolism; (3) pyrimidine metabolism; (4) glyoxylate and dicarboxylate metabolism; and (4) alanine, aspartate, and glutamate metabolism (Figure 5a). Thiamine metabolism differed significantly at six months (Figure 5b). There was no significant difference at nine months (Figure 5c). Two pathways, (1) propanoate metabolism and (2) pyrimidine metabolism, showed significant differences at 12 months (Figure 5d). The detailed data are shown in Figure A5, Figure A6 and Figure A7.

3. Discussion

In this study, the relationship between major depressive episodes and urinary metabolites were investigated during the long-term stays of the 59th JARE in an isolated environment. Of the thirty members, five showed a high PHQ-9 score (≥10) during the 12-month stay. We conducted LC-MS-based metabolomics and compared the urinary metabolites between participants with and without depressive episodes. Changes in metabolites after three-, six-, and nine-month stays were also analyzed to obtain time-series data. These data were analyzed at both the molecular and pathway levels.
A comparison of individual metabolites in the 12-month stay data revealed 11 significantly decreased metabolites. Here, we evaluated the consistency and discrepancy between these metabolites and data reported in other mental illness studies. Among these 11 metabolites, 4-guanidinobutanoate showed the most significant decrease in the depressive episode group (fold change in the depressive episode group/nondepressive episode group, FC = 0.37, p = 0.00041). Altered urinary 4-guanidinobutanoate levels have been reported in patients with liver cirrhosis [26]; however, no study has reported its relationship with mental illnesses. This metabolite binds to water and is converted to 4-aminobutanoate (GABA) and urea by the metabolic enzyme guanidinobutyrase. Accordingly, consistent with our 4-guanidinobutanoate data, GABA (FC = 0.64, p = 0.065) and urea (FC = 0.82, p = 0.12) showed a decreasing trend in the depressive episode group, although it did not decrease significantly.
A significant decrease in 3-methylhistidine levels was observed in the depressive episode group (FC = 0.69, p = 0.037). This metabolite is produced by the post-translational methylation of histidine molecules. Serum concentrations of 3-methylhistidine were significantly reduced in depression model rats [27]. Urinary 3-methylhistidine levels in rats with chronic-unpredictable mild stress-induced depression were also significantly decreased [28]. These results are consistent with our data. However, the urinary concentration of this metabolite in patients with postpartum depression was elevated [19], which is in contrast to our data. The urinary concentration of this metabolite is also elevated under fatigue conditions [29]. The concentration of this metabolite may therefore be affected by both mental illness and fatigue conditions.
In the present study, 2-aminoadipate showed a significant decrease in the depressive episode group at both 9 and 12 months. This intermediate metabolite in lysine metabolism has neuroexcitatory actions and regulates kynurenic acid metabolism in the brain [25]. Ornithine levels were also significantly decreased (FC = 0.22, p = 0.049). Ornithine is a major intermediate metabolite synthesized from arginine in the arginine–citrulline cycle. Arginine–proline metabolism has been reported to be altered in both patients with MDD and BD [30]. Differences in amino acids have been frequently reported in comparisons between patients with MDD and healthy controls, but no significant differences were found in this study. For example, alanine and glycine levels were elevated in MDD; however, our data showed a decreasing trend without significance (FC = 0.66, p = 0.23, and FC = 0.73, p = 0.25) [31].
Pathway and enrichment analyses were also performed. In the comparison of the 12-month-stay data, propanoate metabolism was listed as a pathway with significant differences and changes in both analyses (Figure 4d and Figure 5d). This pathway has also been reported as a potential therapeutic target for BD [32]. However, our pathway analysis showed that the pathway impact was low and only β-alanine was mapped to this pathway (Figure A4). Pyrimidine metabolism also showed significant differences (Figure 4d and Figure 5d). Both pathways include β-alanine, an inhibitory neurogenic amino acid. Notably, this metabolite is altered in urinary samples from patients with MDD and BD [33].
A comparison of the 3-month-stay data showed five significantly altered pathways (Figure 4a). These pathways were also ranked high by enrichment analysis (Figure 5a). These pathways included glutamine and glutamate. D-glutamine and D-glutamate metabolism describes the glutamate–glutamine cycle in glial and neuronal cells, with glutamate serving as an excitatory neurotransmitter in the central nervous system [33]. Altered (1) D-glutamine and D-glutamate metabolism and (2) nitrogen metabolism in urine samples have been frequently reported in patients with MDD and BD [31]. In our study, the changes in these pathways were observed only in the 3-month-stay data, and these alterations tended to be smaller in the six-, nine-, and twelve-month data.
No metabolomic study has analyzed the relationship between urinary metabolites and mental health in JARE. However, various teams have investigated the relationship between mental illnesses and metabolic changes in the Antarctic area [34,35]. For example, metabolism related to bone-mass density is disturbed by vitamin D deficiency caused by limited ultraviolet light exposure [36,37]. The season-dependent alteration of constituents of stress-related metabolism, such as cortisol and immunoglobulins, was also observed [38]. The seasonality of plasma insulin and growth hormones has also been reported [39]. The metabolic pathways related to these studies include fat-soluble intermediate metabolites. Alternatively, the current study analyzed only hydrophilic metabolites, which may explain the few overlaps in altered metabolites between our data and the above-mentioned studies [34,35,36,37,38,39].
This study has several limitations. First, the depression screening tool, PHQ-9, was used as an index to divide the participants into depressive and nondepressive episode groups. Second, changes in urinary metabolites over time were examined to investigate changes under long-term isolated stressful conditions. A problem with using metabolites as markers is that many factors may cause unexpected changes in metabolism, such as the diet and lifestyle of individuals. In this study, the cohorts were designed to reduce these effects; participants ate the same diet and stayed in the same environment. However, we did not control for urine sampling timing, presumably because of diurnal variation. Third, the association between urinary catecholamines and depression has been frequently reported [40,41]. Our pretreatment and measurement method could detect dopamine, a catecholamine, in the standard mixture. However, the peak area of dopamine in the urinary samples was below the limits of detection of our method (Figure A8). A similar phenomenon was observed with serotonin. The protocol used herein was able to evaluate a wide range of hydrophilic metabolites within the linearity limits. However, different sample treatments and/or measurement methods should be applied to quantify these catecholamines. Finally, whereas tryptophan, a precursor of these catecholamines [42], was detected at levels above the detection limits, no significant difference was identified between the depressive and nondepressive episode groups.

4. Materials and Methods

4.1. Subjects and Sample Collections

The 59th JARE stayed at the Showa Station from March 2017 to January 2018. Urinary samples were collected in Japan before the participants visited Antarctica; this served as the baseline data. Participants were required to fill in the following self-administered questionnaires: EPQ-R neuroticism short version and STAI-Y trait anxiety to evaluate personality characteristics, PHQ-9 to evaluate depression episodes, and the Pittsburgh Sleep Quality Index (PSQI) to evaluate sleep quality. In Antarctica, urinary samples and questionnaires were collected every three months (May, Jun., Sep., and Dec. 2018). At each time point, 5 mL of urinary samples were collected in 1.8 mL Nunc Cryotubes (Thermo Fisher, Tokyo, Japan) and were immediately stored at −80 °C. This study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of Tokyo Medical University (protocol code no. SH3712, date 29 May 2017).

4.2. Sample Preparation for Metabolomic Analysis

A urine sample (10 μL) was mixed with MeOH containing 149.6 mM NH4OH and 1.5 µM internal standards (d8-spermine, d8-spermidine, d6-N1-acetylspermidine, d3-N1-acetylspermine, d6-N1,N8-diacetylspermidine, d6-N1,N12-diacetylspermine, hypoxanthine-13C2,15N, and 1,6-diaminohexsne). Samples were centrifuged at 20,380× g for 10 min at 4 °C (MDX-310, Tomy Seiko, Tokyo, Japan), and 90 μL of the supernatant were transferred to a new tube and were vacuum dried for 1 h at room temperature (VC-96W; TAITEC, Saitama, Japan). The samples were reconstituted with 10 μL of 90% methanol (v/v), and 190 μL of H2O were added to the sample solution. The sample was vortexed and centrifuged at 20,380× g for 10 min at 4 °C. The supernatant (150 μL) was injected into the LC-MS system. For creatinine determination, a portion of the supernatant was diluted 5000 fold with H2O. Firstly, all the samples were measured in four batches. Between every 13 samples, quality control samples were measured and the relative deviation of the peak areas of the detected metabolites were determined. Trimethylaniline N-oxide, creatinine, creatine, carnitine, and urea, exhibited peak sizes larger than those of the other metabolites. Thus, secondly, these metabolites in all samples were measured in a single batch with 500-fold dilution. The measurement orders of the samples were randomized, and the standard mixtures were measured consecutively.

4.3. Instrument Parameters for Metabolomic Analysis

The instrument parameters of the LC-time-of-flight mass spectrometry (LC-TOFMS) system used in this study were as described previously, with slight modifications [43]. The LC system was an Agilent Technologies 1290 Infinity instrument (Santa Clara, CA, USA), consisting of an autosampler, quaternary pump, and column compartment. MS detection was conducted using an Agilent Technologies G6230B TOF mass spectrometer. In the present study, we used 40 °C as the temperature of the LC columns. The samples were analyzed in positive ion mode.

4.4. Data Processing

Agilent MassHunter Qualitative Analysis software (version B.08.00; Agilent Technologies) was used for data processing. Peaks were extracted according to the retention time and m/z values corresponding to those of our standard libraries. Peaks out of the range for a retention time of 0.1 min and m/z of 50 ppm were manually integrated for peak area calculations. The peaks below the detection limit were used as 0. The other peaks were within the upper limit of the linearity range. The peak area (no unit) of each metabolite was divided by the area of the corresponding internal standard, resulting in a relative peak area to eliminate the unexpected bias of the sensitivity drift of MS and reduce the effect of ion suppression depending on the relation times. Subsequently, the relative peak area of each metabolite was divided by that of creatinine for the statistical analyses.

4.5. Data Analysis

The time course of the quantified PHQ-9 values and urinary metabolite concentrations were obtained. Higher PHQ-9 scores (≥10) and lower scores (<9) were used to classify participants into the major depressive episode and non-major-depressive episode groups, respectively [44]. The difference between these groups was assessed using the Mann–Whitney U test for quantitative values. The p-values of each metabolite were adjusted using FDR (Benjamini–Hochberg method) to control for the α error caused by multiple independent tests. Pathway and enrichment analyses were conducted using MetaboAnalyst software (ver 5, https://www.metaboanalyst.ca/, accessed on 1 June 2022) [45]. Quantile normalization was performed for each sample. The metabolite concentration was normalized using log transformation and autoscaling options. The Homo sapiens pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.genome.jp/kegg/kegg1.html, accessed on 1 June 2022) were used for both analyses. A principal component analysis was conducted using all samples.
Statistical analyses and data visualizations were conducted using JMP Pro (ver 16.0.0; SAS Institute, Cary, NC, USA) and GraphPad Prism (ver 9.2.0, GraphPad Software, San Diego, CA, USA). Pathway and enrichment analyses were conducted using MeV TM (ver 4.7.4, Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA, http://mev.tm4.org/, accessed on 1 June 2022) [46] and MetaboAnalyst, respectively.

5. Conclusions

The psychological changes in 30 participants in the JARE were assessed using PHQ-9 levels. Urinary samples were collected every three months over one year, and metabolomic analyses were conducted. Five participants showed major depressive episodes (PHQ-9 ≥ 10) during the 12-month stay. In the comparison of the data at 12 months, various metabolites showed significant differences; 2-aminoadipate decreased significantly in the depressive episode group at the 9- and 12-month periods, suggesting that such markers may lead to the early detection of depression.

Author Contributions

Conceptualization, K.K., T.I. (Takashi Ishikawa), S.I. and T.I. (Takeshi Inoue); methodology, M.S. and T.I. (Takeshi Inoue); formal analysis, M.S. and T.I. (Takeshi Inoue); investigation, K.K. and T.I. (Takeshi Inoue); writing—original draft preparation, M.S.; writing—review and editing, K.K.; visualization, M.S.; supervision, T.I. (Takeshi Inoue); project administration, K.K. and T.I. (Takeshi Inoue); funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JSPS KAKENHI (grant number: JP20H05743).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Tokyo Medical University (protocol code no. SH3712, date 29 May 2017).

Informed Consent Statement

Written informed consent was obtained from the patients to publish this paper.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the participants of the JARE.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Score plots of the principal component (PC) analysis. The analytical results using (a) the peak area divided by the peak area of the internal standard and (b) the peak areas divided by the peak area of creatinine. PC1 and PC2 indicate the first and second PC, respectively. Contribution ratios (%) are shown in parentheses.
Figure A1. Score plots of the principal component (PC) analysis. The analytical results using (a) the peak area divided by the peak area of the internal standard and (b) the peak areas divided by the peak area of creatinine. PC1 and PC2 indicate the first and second PC, respectively. Contribution ratios (%) are shown in parentheses.
Ijms 24 00943 g0a1
Figure A2. Detailed results of the pathway analysis in Figure 4a. Nodes with mapped data are colored. In each box plot, C and P indicate the data without and with depressive episodes, respectively. ** p < 0.01.
Figure A2. Detailed results of the pathway analysis in Figure 4a. Nodes with mapped data are colored. In each box plot, C and P indicate the data without and with depressive episodes, respectively. ** p < 0.01.
Ijms 24 00943 g0a2
Figure A3. Detailed results of the pathway analysis in Figure 4b. Nodes with mapped data are colored. In each box plot, C and P indicate the data without and with depressive episodes, respectively. * p < 0.05.
Figure A3. Detailed results of the pathway analysis in Figure 4b. Nodes with mapped data are colored. In each box plot, C and P indicate the data without and with depressive episodes, respectively. * p < 0.05.
Ijms 24 00943 g0a3
Figure A4. Detailed results of the pathway analysis in Figure 4d. Nodes with mapped data are colored. In each box plot, C and P indicate the data without and with depressive episodes, respectively. * p < 0.05.
Figure A4. Detailed results of the pathway analysis in Figure 4d. Nodes with mapped data are colored. In each box plot, C and P indicate the data without and with depressive episodes, respectively. * p < 0.05.
Ijms 24 00943 g0a4
Figure A5. Detailed results of the enrichment analysis in Figure 5a. In each box plot, C and P indicate the data without and with depressive episodes, respectively.
Figure A5. Detailed results of the enrichment analysis in Figure 5a. In each box plot, C and P indicate the data without and with depressive episodes, respectively.
Ijms 24 00943 g0a5
Figure A6. Detailed results of the enrichment analysis in Figure 5b. In each box plot, C and P indicate the data without and with depressive episodes, respectively.
Figure A6. Detailed results of the enrichment analysis in Figure 5b. In each box plot, C and P indicate the data without and with depressive episodes, respectively.
Ijms 24 00943 g0a6
Figure A7. Detailed results of the enrichment analysis in Figure 5c. In each box plot, C and P indicate the data without and with depressive episodes, respectively.
Figure A7. Detailed results of the enrichment analysis in Figure 5c. In each box plot, C and P indicate the data without and with depressive episodes, respectively.
Ijms 24 00943 g0a7
Figure A8. Representative extracted ion chromatographs of dopamine (m/z = 154.0863, positive mode). (a) Standard mixture. (b) Urinary sample.
Figure A8. Representative extracted ion chromatographs of dopamine (m/z = 154.0863, positive mode). (a) Standard mixture. (b) Urinary sample.
Ijms 24 00943 g0a8

References

  1. Rothblum, E.D. Psychological factors in the antarctic. J. Psychol. 1990, 124, 253–273. [Google Scholar] [CrossRef] [PubMed]
  2. Mullin, C.S., Jr. Some psychological aspects of isolated Antarctic living. Am. J. Psychiatry 1960, 117, 323–325. [Google Scholar] [CrossRef] [PubMed]
  3. Palinkas, L.A.; Suedfeld, P. Psychological effects of polar expeditions. Lancet 2008, 371, 153–163. [Google Scholar] [CrossRef] [PubMed]
  4. Palinkas, L.A.; Glogower, F.; Dembert, M.; Hansen, K.; Smullen, R. Incidence of psychiatric disorders after extended residence in Antarctica. Int. J. Circumpolar. Health 2004, 63, 157–168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Norman, J.N. A comparison of the patterns of illness and injury occurring on offshore structures in the northern North Sea and the stations of the British Antarctic Survey. Arctic Med. Res. 1991, 719–721. [Google Scholar]
  6. Zimmer, M.; Cabral, J.C.C.R.; Borges, F.C.; Côco, K.G.; Hameister, B.d.R. Psychological changes arising from an Antarctic stay: Systematic overview. Estud. Psicol. (Camp.) 2013, 30, 415–423. [Google Scholar] [CrossRef]
  7. Palinkas, L.A.; Cravalho, M.; Browner, D. Seasonal variation of depressive symptoms in Antarctica. Acta Psychiatr. Scand. 1995, 91, 423–429. [Google Scholar] [CrossRef]
  8. Palinkas, L.A. Effects of physical and social environments on the health and well-being of Antarctic winter-over personnel. Environ. Behav. 1991, 23, 782–799. [Google Scholar] [CrossRef]
  9. Reed, H.L.; Reedy, K.R.; Palinkas, L.A.; Van Do, N.; Finney, N.S.; Case, H.S.; LeMar, H.J.; Wright, J.; Thomas, J. Impairment in cognitive and exercise performance during prolonged antarctic residence: Effect of thyroxine supplementation in the polar triiodothyronine syndrome. J. Clin. Endocrinol. Metab. 2001, 86, 110–116. [Google Scholar] [CrossRef]
  10. Weiss, K.; Suedfeld, P.; Steel, G.D.; Tanaka, M. Psychological adjustment during three Japanese Antarctic research expeditions. Environ. Behav. 2000, 32, 142–156. [Google Scholar] [CrossRef]
  11. Ikegawa, M.; Kimura, M.; Makita, K.; Itokawa, Y. Psychological studies of a Japanese winter-over group at Asuka Station, Antarctica. Aviat Space Environ. Med. 1998, 69, 452–460. [Google Scholar] [PubMed]
  12. Alfano, C.A.; Bower, J.L.; Connaboy, C.; Agha, N.H.; Baker, F.L.; Smith, K.A.; So, C.J.; Simpson, R.J. Mental health, physical symptoms and biomarkers of stress during prolonged exposure to Antarctica’s extreme environment. Acta Astronaut. 2021, 181, 405–413. [Google Scholar] [CrossRef]
  13. Moraes, M.M.; Bruzzi, R.S.; Martins, Y.A.T.; Mendes, T.T.; Maluf, C.B.; Ladeira, R.V.P.; Núñez-Espinosa, C.; Soares, D.D.; Wanner, S.P.; Arantes, R.M.E. Hormonal, autonomic cardiac and mood states changes during an Antarctic expedition: From ship travel to camping in Snow Island. Physiol. Behav. 2020, 224, 113069. [Google Scholar] [CrossRef] [PubMed]
  14. Tingate, T.R.; Lugg, D.J.; Muller, H.K.; Stowe, R.P.; Pierson, D.L. Antarctic isolation: Immune and viral studies. Immunol. Cell Biol. 1997, 75, 275–283. [Google Scholar] [CrossRef]
  15. Shearer, W.T.; Lee, B.N.; Cron, S.G.; Rosenblatt, H.M.; Smith, E.O.; Lugg, D.J.; Nickolls, P.M.; Sharp, R.M.; Rollings, K.; Reuben, J.M. Suppression of human anti-inflammatory plasma cytokines IL-10 and IL-1RA with elevation of proinflammatory cytokine IFN-gamma during the isolation of the Antarctic winter. J. Allergy Clin. Immunol. 2002, 109, 854–857. [Google Scholar] [CrossRef]
  16. Shirai, T.; Magara, K.K.; Motohashi, S.; Yamashita, M.; Kimura, M.; Suwazomo, Y.; Nogawa, K.; Kuriyama, T.; Taniguchi, M.; Nakayama, T. TH1-biased immunity induced by exposure to Antarctic winter. J. Allergy Clin. Immunol. 2003, 111, 1353–1360. [Google Scholar] [CrossRef]
  17. Xu, C.; Zhu, G.; Xue, Q.; Zhang, S.; Du, G.; Xi, Y.; Palinkas, L.A. Effect of the Antarctic environment on hormone levels and mood of Chinese expeditioners. Int. J. Circumpolar. Health 2003, 62, 255–267. [Google Scholar] [CrossRef]
  18. Shih, P.B. Metabolomics Biomarkers for Precision Psychiatry. Adv. Exp. Med. Biol. 2019, 1161, 101–113. [Google Scholar] [CrossRef]
  19. Lin, L.; Chen, X.M.; Liu, R.H. Novel urinary metabolite signature for diagnosing postpartum depression. Neuropsychiatr. Dis. Treat. 2017, 13, 1263–1270. [Google Scholar] [CrossRef] [Green Version]
  20. Xie, J.; Han, Y.; Hong, Y.; Li, W.W.; Pei, Q.; Zhou, X.; Zhang, B.; Wang, Y. Identification of Potential Metabolite Markers for Middle-Aged Patients with Post-Stroke Depression Using Urine Metabolomics. Neuropsychiatr. Dis. Treat. 2020, 16, 2017–2024. [Google Scholar] [CrossRef]
  21. Chen, J.J.; Xie, J.; Li, W.W.; Bai, S.J.; Wang, W.; Zheng, P.; Xie, P. Age-specific urinary metabolite signatures and functions in patients with major depressive disorder. Aging (Albany N. Y.) 2019, 11, 6626–6637. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, J.J.; Zhou, C.J.; Liu, Z.; Fu, Y.Y.; Zheng, P.; Yang, D.Y.; Li, Q.; Mu, J.; Wei, Y.D.; Zhou, J.J.; et al. Divergent Urinary Metabolic Phenotypes between Major Depressive Disorder and Bipolar Disorder Identified by a Combined GC-MS and NMR Spectroscopic Metabonomic Approach. J. Proteome Res. 2015, 14, 3382–3389. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, J.J.; Liu, Z.; Fan, S.H.; Yang, D.Y.; Zheng, P.; Shao, W.H.; Qi, Z.G.; Xu, X.J.; Li, Q.; Mu, J.; et al. Combined application of NMR- and GC-MS-based metabonomics yields a superior urinary biomarker panel for bipolar disorder. Sci. Rep. 2014, 4, 5855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kawamura, N.; Shinoda, K.; Sato, H.; Sasaki, K.; Suzuki, M.; Yamaki, K.; Fujimori, T.; Yamamoto, H.; Osei-Hyiaman, D.; Ohashi, Y. Plasma metabolome analysis of patients with major depressive disorder. Psychiatry Clin. Neurosci. 2018, 72, 349–361. [Google Scholar] [CrossRef] [Green Version]
  25. Homorogan, C.; Nitusca, D.; Enatescu, V.; Schubart, P.; Moraru, C.; Socaciu, C.; Marian, C. Untargeted Plasma Metabolomic Profiling in Patients with Major Depressive Disorder Using Ultra-High Performance Liquid Chromatography Coupled with Mass Spectrometry. Metabolites 2021, 11, 466. [Google Scholar] [CrossRef]
  26. Marescau, B.; De Deyn, P.P.; Holvoet, J.; Possemiers, I.; Nagels, G.; Saxena, V.; Mahler, C. Guanidino compounds in serum and urine of cirrhotic patients. Metabolism 1995, 44, 584–588. [Google Scholar] [CrossRef]
  27. Heath, A.; Lindberg, D.R.; Makowiecki, K.; Gray, A.; Asp, A.J.; Rodger, J.; Choi, D.-S.; Croarkin, P.E. Medium-and high-intensity rTMS reduces psychomotor agitation with distinct neurobiologic mechanisms. Transl Psychiatry 2018, 8, 126. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, C.; Mo, Y.Y.; Feng, S.S.; Meng, M.W.; Chen, S.Y.; Huang, H.M.; Ling, X.; Song, H.; Liang, Y.H.; Ou, S.F.; et al. Urinary metabonomics study of anti-depressive mechanisms of Millettia speciosa Champ on rats with chronic unpredictable mild stress-induced depression. J. Pharm. Biomed. Anal. 2021, 205, 114338. [Google Scholar] [CrossRef]
  29. Manabe, S.; Kurroda, I.; Okada, K.; Morishima, M.; Okamoto, M.; Harada, N.; Takahashi, A.; Sakai, K.; Nakaya, Y. Decreased blood levels of lactic acid and urinary excretion of 3-methylhistidine after exercise by chronic taurine treatment in rats. J. Nutr. Sci. Vitaminol. (Tokyo) 2003, 49, 375–380. [Google Scholar] [CrossRef] [Green Version]
  30. Cao, B.; Deng, R.; Wang, D.; Li, L.; Ren, Z.; Xu, L.; Gao, X. Association between arginine catabolism and major depressive disorder: A protocol for the systematic review and meta-analysis of metabolic pathway. Medicine (Baltimore) 2020, 99, e21068. [Google Scholar] [CrossRef]
  31. MacDonald, K.; Krishnan, A.; Cervenka, E.; Hu, G.; Guadagno, E.; Trakadis, Y. Biomarkers for major depressive and bipolar disorders using metabolomics: A systematic review. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2019, 180, 122–137. [Google Scholar] [CrossRef]
  32. Zarate, C.A., Jr.; Singh, J.; Manji, H.K. Cellular plasticity cascades: Targets for the development of novel therapeutics for bipolar disorder. Biol. Psychiatry 2006, 59, 1006–1020. [Google Scholar] [CrossRef] [PubMed]
  33. Auer, D.P.; Pütz, B.; Kraft, E.; Lipinski, B.; Schill, J.; Holsboer, F. Reduced glutamate in the anterior cingulate cortex in depression: An in vivo proton magnetic resonance spectroscopy study. Biol. Psychiatry 2000, 47, 305–313. [Google Scholar] [CrossRef] [PubMed]
  34. Nirwan, M. Human psychophysiology in Antarctica. Sri. Ramachandra J. Health Sci. 2022, 2, 12–18. [Google Scholar] [CrossRef]
  35. Spinelli, E.; Werner, J.J. Human adaptative behavior to Antarctic conditions: A review of physiological aspects. WIREs Mech. Dis. 2022, 14, e1556. [Google Scholar] [CrossRef]
  36. Iuliano-Burns, S.; Wang, X.F.; Ayton, J.; Jones, G.; Seeman, E. Skeletal and hormonal responses to sunlight deprivation in Antarctic expeditioners. Osteoporos. Int. 2009, 20, 1523–1528. [Google Scholar] [CrossRef]
  37. Smith, S.M.; Gardner, K.K.; Locke, J.; Zwart, S.R. Vitamin D supplementation during Antarctic winter. Am. J. Clin. Nutr. 2009, 89, 1092–1098. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Yadav, A.P.; Mishra, K.P.; Ganju, L.; Singh, S.B. Wintering in Antarctica: Impact on immune response of Indian expeditioners. Neuroimmunomodulation 2012, 19, 327–333. [Google Scholar] [CrossRef] [PubMed]
  39. Malhotra, A.S.; Pal, K.; Prasad, R.; Bajaj, A.C.; Kumar, R.; Sawhney, R.C. Plasma insulin and growth hormone during antarctic residence. Jpn. J. Physiol. 1998, 48, 167–169. [Google Scholar] [CrossRef] [Green Version]
  40. Homandberg, L.K.; Fuller-Rowell, T.E. Experiences of discrimination and urinary catecholamine concentrations: Longitudinal associations in a college student sample. Ann. Behav. Med. 2020, 54, 843–852. [Google Scholar] [CrossRef]
  41. Kokka, I.; Chrousos, G.P.; Darviri, C.; Bacopoulou, F. Measuring Adolescent Chronic Stress: A Review of Established Biomarkers and Psychometric Instruments. Horm. Res. Paediatr. 2022. [Google Scholar] [CrossRef]
  42. Gelenberg, A.J.; Gibson, C.J. Tyrosine for the treatment of depression. Nutr. Health 1984, 3, 163–173. [Google Scholar] [CrossRef] [PubMed]
  43. Tomita, A.; Mori, M.; Hiwatari, K.; Yamaguchi, E.; Itoi, T.; Sunamura, M.; Soga, T.; Tomita, M.; Sugimoto, M. Effect of storage conditions on salivary polyamines quantified via liquid chromatography-mass spectrometry. Sci. Rep. 2018, 8, 12075. [Google Scholar] [CrossRef] [PubMed]
  44. Inoue, T.; Tanaka, T.; Nakagawa, S.; Nakato, Y.; Kameyama, R.; Boku, S.; Toda, H.; Kurita, T.; Koyama, T. Utility and limitations of PHQ-9 in a clinic specializing in psychiatric care. BMC Psychiatry 2012, 12, 73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef] [PubMed]
  46. Saeed, A.I.; Sharov, V.; White, J.; Li, J.; Liang, W.; Bhagabati, N.; Braisted, J.; Klapa, M.; Currier, T.; Thiagarajan, M.; et al. TM4: A free, open-source system for microarray data management and analysis. Biotechniques 2003, 34, 374–378. [Google Scholar] [CrossRef]
Figure 1. Average and standard deviation of time-series changes in PHQ-9. The solid line indicates the depressive episode group (n = 5), and the dotted line indicates the nondepressive episode group (n = 25) with PHQ-9 reaching ≥10 at the 12-month stay. Data between the two groups at each time point were compared using Mann–Whitney tests. * p < 0.05, ** p < 0.01, **** p < 0.0001.
Figure 1. Average and standard deviation of time-series changes in PHQ-9. The solid line indicates the depressive episode group (n = 5), and the dotted line indicates the nondepressive episode group (n = 25) with PHQ-9 reaching ≥10 at the 12-month stay. Data between the two groups at each time point were compared using Mann–Whitney tests. * p < 0.05, ** p < 0.01, **** p < 0.0001.
Ijms 24 00943 g001
Figure 2. Heatmap of time-series changes in metabolites in the depressive episode and nondepressive episode groups. Mean values were calculated within each group and divided by the mean value for each metabolite to determine the appropriate color. Red and blue indicate values higher and lower than the average, respectively, and white indicates the average. The order of metabolites was clustered using Pearson’s correlation coefficient, and the metabolites with high similarity of change patterns in time series are arranged close to each other. In the graph on the right, black squares indicate data showing a significant difference between the two groups using a Mann–Whitney test (p < 0.05).
Figure 2. Heatmap of time-series changes in metabolites in the depressive episode and nondepressive episode groups. Mean values were calculated within each group and divided by the mean value for each metabolite to determine the appropriate color. Red and blue indicate values higher and lower than the average, respectively, and white indicates the average. The order of metabolites was clustered using Pearson’s correlation coefficient, and the metabolites with high similarity of change patterns in time series are arranged close to each other. In the graph on the right, black squares indicate data showing a significant difference between the two groups using a Mann–Whitney test (p < 0.05).
Ijms 24 00943 g002
Figure 3. Time-series changes in metabolites that differed between the depressive episode and nondepressive episode groups. Metabolites showing significant differences at the (a) 12-month and (b) shorter stays are shown. Mean values and standard deviations are shown as plots and error bars, respectively. The X-axis indicates the length of stay (months) and the Y-axis indicates the metabolite concentration divided by the creatinine concentration (unitless). * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 3. Time-series changes in metabolites that differed between the depressive episode and nondepressive episode groups. Metabolites showing significant differences at the (a) 12-month and (b) shorter stays are shown. Mean values and standard deviations are shown as plots and error bars, respectively. The X-axis indicates the length of stay (months) and the Y-axis indicates the metabolite concentration divided by the creatinine concentration (unitless). * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 24 00943 g003
Figure 4. Pathway analysis results between depressive and nondepressive episode groups. Analysis results for (a) 3-, (b) 6-, (c) 9-, and (d) 12-month stays are shown. Metabolism is omitted from metabolite names. Underlined pathways indicate significant differences (p < 0.05).
Figure 4. Pathway analysis results between depressive and nondepressive episode groups. Analysis results for (a) 3-, (b) 6-, (c) 9-, and (d) 12-month stays are shown. Metabolism is omitted from metabolite names. Underlined pathways indicate significant differences (p < 0.05).
Ijms 24 00943 g004
Figure 5. Enrichment analysis between depressive and nondepressive episode groups. In each panel, pathways are aligned according to p-values, and the length of each bar represents the enrichment ratio. Analysis results for (a) 3-, (b) 6-, (c) 9-, and (d) 12-month stays are shown.
Figure 5. Enrichment analysis between depressive and nondepressive episode groups. In each panel, pathways are aligned according to p-values, and the length of each bar represents the enrichment ratio. Analysis results for (a) 3-, (b) 6-, (c) 9-, and (d) 12-month stays are shown.
Ijms 24 00943 g005
Table 1. Subject characteristics.
Table 1. Subject characteristics.
CharacteristicAverage/n/hStandard Deviation
Age40.2 9.2
SexM 29, F 1
Education Years15.32.2
Subjective Social Status 15.51.2
Marital StatusMarried 18, Unmarried 11
Presence of Physical DiseaseY 1, N 29
Presence of Psychiatric DiseaseY 0, N 30
History of Psychiatric DiseaseY 2, N 28
Family history of Psychiatric DiseaseY 2, N 28
Alcohol DrinkingY 24, N 6
Smoking Y 6, N 24
Overtime work (h/month)
   ≤2019
   ≤304
   ≤404
   ≤502
Parental Bonding Instrument
   Paternal Care21.76.7
   Paternal Overprotection10.85.7
   Maternal Care26.45.4
   Maternal Overprotection11.56.2
PHQ-9 total scores 233.5
   State anxiety of STAI-Y41.79.2
   Trait anxiety of STAI-Y40.19.6
   Neuroticism (0−12)3.23
   Pittsburgh Sleep Quality Index Global Scores4.51.9
1 1 and 10 are the lowest and the highest, respectively. 2 Depressive symptoms.
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.

Share and Cite

MDPI and ACS Style

Kasuya, K.; Imura, S.; Ishikawa, T.; Sugimoto, M.; Inoue, T. Relationship between Urinary Metabolomic Profiles and Depressive Episode in Antarctica. Int. J. Mol. Sci. 2023, 24, 943. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24020943

AMA Style

Kasuya K, Imura S, Ishikawa T, Sugimoto M, Inoue T. Relationship between Urinary Metabolomic Profiles and Depressive Episode in Antarctica. International Journal of Molecular Sciences. 2023; 24(2):943. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24020943

Chicago/Turabian Style

Kasuya, Kazuhiko, Satoshi Imura, Takashi Ishikawa, Masahiro Sugimoto, and Takeshi Inoue. 2023. "Relationship between Urinary Metabolomic Profiles and Depressive Episode in Antarctica" International Journal of Molecular Sciences 24, no. 2: 943. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24020943

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop