Blood and urinary metabolomic evidence validating traditional Chinese medicine diagnostic classification of major depressive disorder (original) (raw)

Chinese Medicine volume 13, Article number: 53 (2018)Cite this article

Abstract

Background

Major depressive disorder (MDD) is a highly heterogeneous disease. Further classification may characterize its heterogeneity. The purpose of this study was to examine whether metabolomic variables could differentiate traditional Chinese medicine (TCM) diagnostic subtypes of MDD.

Methods

Fifty medication-free patients who were experiencing a recurrent depressive episode were classified into Liver Qi Stagnation (LQS, n = 30) and Heart and Spleen Deficiency (HSD, n = 20) subtypes according to TCM diagnosis. Healthy volunteers (n = 28) were included as controls. Gas chromatography-mass spectrometry (GC–MS) was used to examine serum and urinary metabolomic profiles.

Results

Twenty-eight metabolites were identified for good separations between TCM subtypes and healthy controls in serum samples. Both TCM subtypes had similar profiles in proteinogenic branched-chain amino acids (BCAAs) (valine, leucine, and isoleucine) and energy metabolism-related metabolites that were differentiated from healthy controls. The LQS subtype additionally differed from healthy controls in multiple amino acid metabolites that are involved in biosynthesis of monoamine and amino acid neurotransmitters, including phenylalanine, 3-hydroxybutric acid, _o_-tyrosine, glycine, l-tryptophan, and _N_-acetyl-l-aspartic acid. Threonic acid, methionine, stearic acid, and isobutyric acid are differentially associated with the two subtypes.

Conclusions

While both TCM subtypes are associated with aberrant BCAA and energy metabolism, the LQS subtype may represent an MDD subpopulation characterized by abnormalities in the biosynthesis of monoamine and amino acid neurotransmitters and closer associations with stress-related pathophysiology. The metabolites differentially associated with the two subtypes are promising biomarkers for predicting TCM subtype-specific antidepressant response [registered at http://www.clinicaltrials.gov (NCT02346682) on January 27, 2015].

Background

Major depressive disorder (MDD) is a highly heterogeneous mental illness with a wide range of clinical manifestations and multi-system etiopathogenesis. The core symptoms of low mood and loss of interest and pleasure are often accompanied by somatic and other psychiatric symptoms, such as pain, sleep disturbance, anxiety, and neurocognitive dysfunction [1]. Although a great attempt has been made to characterize its heterogeneity by subdividing MDD based on its clinical features, severity, and polygenic features, convincing evidence for the existence of depressive symptom dimensions and symptomatic subtypes is lacking, mainly due to symptomatic diversity and the absence of patterns [2,3,4,5]. Therefore, a novel classification system that could improve clinical applicability and more precisely identify pathological heterogeneity of MDD is highly desirable.

As an ancient medical practice, traditional Chinese medicine (TCM) has an established distinct diagnostic system. One key approach to diagnosis is the differentiation of etiopathological pattern which is generally defined by a comprehensive analysis of clinical symptoms and signs collected through inspection, auscultation, olfaction, interrogation, and palpation of the pulses [6]. Over the past decade, several large-scale studies have revealed 12 different TCM patterns of MDD on the basis of the modern psychiatric diagnostic instruments and analysis tools, such as latent tree model analysis [7,8,9,10,11,12]. Despite disagreements on some patterns, a consensus on clinical diagnostic criteria for the two most common and opposing patterns, Liver Qi Stagnation (LQS) and Heart and Spleen Deficiency (HSD), which account for approximately 2/3 of depressed patients, have been reached as shown in Table 1 [7, 9, 10].

Table 1 Clinical manifestations and diagnostic criteria of TCM-based subtypes of MDD

Full size table

Differential profiles of functional brain connectivity has been observed between “deficiency” and “excessive” patterns of MDD [13]. One recent study has further revealed differences between patients with liver-depression and spleen-deficiency syndrome and healthy volunteers in multiple serum metabolic pathways associated with amino acid and energy metabolism [14]. On the other hand, aberrant amino acid metabolism that is involved in protein synthesis, energy metabolism, monoamine and amino acid neurotransmitter biosynthesis has been suggested in the pathogenesis of MDD and antidepressant response [15]. Differential metabolomic profiles are also associated with the severity of MDD and whether the patients had early life stress and suicidal ideation [16,17,18,19,20]. These studies have led to the hypothesis that certain metabolic pathways particularly associated with protein synthesis, energy metabolism, and biosynthesis of monoamine and amino acid neurotransmitters could serve valid biomarkers for differentiating TCM diagnostic subtypes of MDD.

To test this hypothesis, this study was designed to examine whether blood and urinary metabolomic profiles were different between the LQS and HSD subtypes of MDD with the inclusion of healthy volunteers as controls.

Materials and methods

Setting and participants

This study was conducted at Tongde Hospital of Zhejiang Province in Hangzhou, China between April 2015 and March 2017. The study protocol was approved by the Medical Ethical Committee [equivalent to the Institutional Review Board (IRB)] of Tongde Hospital of Zhejiang Province on March 11, 2014 (The approved document is attached in Additional file 1) and registered at http://www.clinicaltrials.gov (NCT02346682) on January 27, 2015 (https://clinicaltrials.gov/ct2/show/NCT02346682). All participants gave voluntary, written, informed consent before entering the study. We reported this study according to the Minimum Standards of Reporting Checklist (see Additional file 2).

Participants were sought from out-inpatients. Subjects were eligible for this study if they: (a) were aged 18–65 years; (b) were currently experiencing a recurrent, moderate or severe depressive episode according to the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5), as evidenced by a score of at least 21 on the 24-item Hamilton Rating Scale for Depression (HAMD-24) [21]; (c) met the diagnostic criteria of the LQS or HSD subtype as defined in Table 1 with representative tongue manifestations as shown in Fig. 1; and (d) had received no treatment with antidepressants and other psychotropic drugs in the previous 3 months. Tongue diagnosis is a critical diagnostic approach in discriminating different TCM patterns according to manifestations of the tongue coating and body [22].

Fig. 1

figure 1

Representative tongue pictures taken from a healthy volunteer (a), depressed patients with Liver Qi Stagnation (b) and Heart Spleen Deficiency (c) subtypes

Full size image

Patients were excluded from this study if they: (a) were experiencing their first episode of depression; (b) had serious comorbid cardiac, hepatic, or renal conditions; (c) had a history of brain injury or surgery; (d) had a history of manic, hypomanic, or mixed episodes; (e) had received investigational drug treatment within the previous 6 months; (f) had experienced alcohol or drug abuse within the previous 12 months; or (g) were pregnant or breastfeeding.

A group of healthy volunteers who had no personal or family history of significant mental and physical illness was recruited from Tongde Hospital and local community to serve as controls.

Screening, assessment, and TCM diagnosis

The screening of patients and healthy volunteers was done by a psychiatrist and a TCM practitioner. The severity of the depression symptoms of the patients and healthy controls was assessed using HAMD-24 [21]. The assessment was conducted by a trained rater (L.Y.L., C.Q.Z.). TCM subtypes diagnoses were made by at least two senior TCM practitioners (Y.P.L., C.Q.Z.) and a third TCM practitioner (L.Y.L.) was involved in the diagnosis process if the first two practitioners could not reach an agreement.

To ensure consistency of assessment and diagnosis, a manual was provided and a training workshop was carried out on videotaped patients with different TCM subtypes before the study was initiated. An inter-rater reliability coefficient of > 0.85 was achieved after the completion of the training workshop.

Collection and preparation of serum and urine samples

Blood and urine sample for each were collected within 2 days after the completion of clinical assessment. Following the sample collection, patients immediately received antidepressant and other psychotropic treatment at psychiatrists’ discretion based on their clinical condition, with routine subsequent monitoring.

In the morning (09:00–10:00) after an overnight fasting, 10 ml of midstream urine were collected; 15 ml of blood were drawn and sera were immediately separated. Samples were further prepared for metabolomic measurement as described previously [17, 23]. Briefly, 500 μl of urine or sera was mixed with 50 μl heptanoic acid and 10 mg/ml para-chlorophenylalanine as the internal standard solution. 500 μl methanol was added to precipitate protein, mixed for 5 min, and then centrifuged at 12,000 rpm/min for 15 min at 4 °C. A 500-μl aliquot of supernatant was transferred to a clean Eppendorf tube and dried under a low temperature vacuum drier. The residue was derivatized by adding 100 μl methoxyamine hydrochloride (15 mg/ml in pyridine) at 30 °C for 90 min. For each sample, 100 μl BSTFA (1% TMCS) was added and the mixture was heated at 70 °C for 60 min. The derivative was cooled and filtered in 0.45 μm membrane prior to GC–MS analysis. Quality control samples and the reference standard were processed as done for experimental samples.

GC–MS acquisition

One μl aliquot of derivatized sample was injected into a Varian 450-GC/240-MS equipped with 19091 N-113 capillary column (30 m × 0.32 mm × 0.25 μm, Agilent J & W Scientific, USA) at a split ratio of 10:1. Helium was used as the carrier gas with a constant flow rate of 1 ml/min. The initial temperature was set at 70 °C for 4 min, elevated to 300 °C at a rate of 8°C/min, and then maintained for 3 min. Temperature for the injector, transfer line, and ion source was set at 280 °C, 250 °C, and 220 °C, respectively. The mass range (50–800 m/z) in a full-scan mode for electron impact ionization (1.0 kV) was applied. The solvent delay time was set to 6 min.

Statistical analysis

There were no studies detecting metabolomic effects in TCM subtypes of MDD. Sample size estimation was based on one previous study that has revealed differential metabolomic profiles between MDD patients with and without early life stress [17]. As the LQS and HSD subtypes are appear to be differentially associated with stress-related (reactive) and endogenous (melancholic) depression [2, 9], we assumed that metabolomic differences between the two TCM subtypes was similar to that of MDD patients with and without early life stress [17]. The study has shown an averaged 36% difference in plasma level of major metabolites between MDD patients with and without early life stress with an averaged standard deviation of 41% [17]. A sample size of 22 each group would be sufficient to yield an 80% power at a statistical level of 0.05.

For baseline data, one-way analysis of variance (ANOVA) was used to detect differences in continuous variables among healthy controls and the two TCM subtypes. Student’s _t_-test was used to detect differences in continuous variables between the two TCM subtypes. Categorical baseline variables were analyzed using Chi square (χ2) test.

For metabolomic data, the pretreatment process was performed, including novel nonlinear retention time alignment, baseline filtration, peak identification, matching, and integration. The resulting data matrix consisting of variables, sample code, and peak area was further processed using Microsoft Excel program. The original spectral data obtained from GC-MS spectroscopy were scaled to unit variance (z), which was calculated from the formula z = (x – y)/s, where x, y, and s represent the level of the particular metabolite in one subject, the mean level and the standard deviation of this metabolite across all subjects, respectively.

Nonparametric Mann–Whitney U test was used to detect significantly differential metabolites. All metabolites were determined by standard samples and/or a similarity of > 70% that was obtained by comparing with the mass spectral database of the US National Institute for Standards and Technology (NIST).

Metabolites were determined using variable importance in the projection (VIP) which value was defined as > 1 and _t_-test was set at a level of 95%. The principal component analysis (PCA) was used to discriminate metabolic patterns among the three groups. The orthogonal projection to latent structures-discriminant analysis (OPLS-DA) model was further constructed to identify meaningful metabolites that could differentiate between the three groups using SIMCA-P 13.0 program. A regression method was applied to establish the optimal discriminant model with controlling gender as a confounding factor. The quality of the model was tested with cross-validation and R2X, R2Y, and Q2 values were obtained. OPLS-DA models which all R2X, R2Y, and Q2 values were ≥ 0.5 were acceptable. Hierarchical clustering patterns were further established with heatmap analysis using R i386 3.3.0 software for visual verification of OPLS-DA models. The impact pathway was determined using MetaboAnalyst 3.0 (http://www.metaboanalyst.ca), a web-based tool for pathway analysis, with the criteria as _p_-value < 0.05, false discovery rate (FDR) < 0.05, and impact value > 0. Visualization of metabolomic correlation network was drawn with CytoScape 3.3.0. The receiver operating characteristic (ROC) curve analysis was conducted to determine the optimal metabolite combination patterns that could well dichotomize the subtypes and healthy controls at acceptable sensitivity and specificity (defined as greater than 80% for both).

Statistical analysis was conducted using SPSS 18.0 (SPSS Inc., Chicago, IL, USA) and, unless otherwise indicated, significance level was set at a two-tailed P < 0.05.

Results

Participant characteristics

Flowchart of screening and recruitment is shown in Fig. 2. Of 243 depressed patients screened, 189 met the diagnostic criteria of LQS (n = 102) and SHD (n = 87) subtype; 27.0% (51/189) rejected to participate the study; 30 LQS subtypes and 20 SHD subtypes who met the inclusion criteria but did not meet any exclusion criteria were recruited. Twenty-eight healthy volunteers were recruited to serve as controls. The proportion of female participants of the LQS subtype was significantly higher than that of the HSD subtype (χ2 = 6.696, p = 0.035) (Table 2). Other demographic and clinical variables were not different among the three groups.

Fig. 2

figure 2

Flowchart of screening and recruitment. LQS Liver Qi Stagnation, HSD Heart and Spleen Deficiency

Full size image

Table 2 Demographic and clinical characteristics of participants

Full size table

Clustering patterns of differential metabolite profiles

A total of 28 metabolites were identified to well differentiate the two TCM subtypes from healthy controls (Table 3); 8 of them were present in both serum and urine samples and/or in the two subtypes. The OPLS-DA models revealed that serum metabolite profiles of both subtypes were clearly separated from healthy controls with R2X values of ≥ 0.842, R2Y values of ≥ 0.670, and Q2 values of ≥ 0.502 (Fig. 3a, b). Urinary metabolite profiles of the two subtypes were not well separated from healthy controls, with considerable overlaps and low R2Y and Q2 values (Fig. 2d, e). The two subtypes showed large overlaps in serum and urinary metabolite profiles with low R2Y and Q2 values (Fig. 2c, f). Hierarchical clustering analysis displayed remarkable differences in mentalism between the three groups, visually validating OPLS-DA models in serum and urine samples (Figs. 4 and 5).

Table 3 Differential metabolites identified from MDD patients with TCM subtypes and healthy controls a

Full size table

Fig. 3

figure 3

Clustering analysis of serum (a, b, c) and urine (d, e, f) metabolomic profiles. OPLS-DA models were built in patients with Liver Qi Stagnation (LQS) and Heart and Spleen Deficiency (HSD) subtypes of MDD and healthy controls (HC). Acceptable criteria for the models were defined as all R2X, R2Y, and Q2 values of ≥ 0.5. Comparisons were conducted between HC and LQS (a, d), HC and HSD (b, e), and LQS and HSD (c, f)

Full size image

Fig. 4

figure 4

Heat maps generated from hierarchical Pearson clustering show metabolites in serum samples obtained from Liver Qi Stagnation (LQS) and Heart and Spleen Deficiency (HSD) subtypes of MDD and healthy controls (HC)

Full size image

Fig. 5

figure 5

Heat maps generated from hierarchical Pearson clustering show metabolites in urine samples obtained from Liver Qi Stagnation (LQS) and Heart and Spleen Deficiency (HSD) subtypes of MDD and healthy controls (HC)

Full size image

Differential metabolites of TCM subtypes

Between LQS subtype and healthy controls

24 metabolites were valid in discriminating the LQS subtype from healthy controls (Fig. 6, Table 3). In serum samples, levels of l-valine, l-phenylalanine, l-lysine, l-proline, l-leucine, and isobutyric acid were markedly higher than those of healthy controls and positively correlated with the LQS subtype. l-lactic acid level of the LQS subtype was greater than that of healthy subjects with negative correlation with the subtype. Levels of glycine and threonic acid were significantly lower than those of healthy controls and negatively correlated with the LQS subtype.

Fig. 6

figure 6

Correlation network analysis of 28 differential metabolites among Liver Qi Stagnation (LQS) and Heart and Spleen Deficiency (HSD) subtypes of MDD versus healthy controls. Green and red lines with ovals indicates negative and positive correlation, respectively. The ovals with half green and half red colors indicate mixed correlations. The oval size represents p value compared with healthy controls. Isobutyric acid (dark brown oral) has a significant difference between the two subtypes. B blood samples, U urine samples

Full size image

In urine samples, levels of l-sorbose, diacetyl, _N_-acetyl-d-glucosamine, l-methionine, 3-hydroxybutyric acid, pyruvic acid, stearic acid, trans-aconitic acid, glycine, _N_-acetyl-l-aspartic acid, and threonic acid were significantly higher than those of healthy subjects and positively correlated with the LQS subtype. Levels of l-tryptophan, α-lactose, l-lactic acid, palmitic acid, indoxyl sulfate, and o-tyrosine were significantly greater than those of healthy subjects but negatively correlated with the LQS subtype. Levels of _p_-hydroxyphenylacetic acid and citric acid were significantly lower higher than those of healthy subjects and negatively correlated with the LQS subtype.

Between HSD subtype and healthy controls

13 metabolites were identified to be effective in distinguishing the HSD subtype from healthy controls (Fig. 6, Table 3). In serum samples, levels of l-valine, l-leucine, and l-lysine were remarkably higher than those of healthy subjects and positively correlated with the HSD subtype. Levels of l-threonine, d-chiro-Inositol, l-methionine, and l-lactic acid were significantly lower than those of healthy subjects and negatively correlated with the HSD subtype.

In urine samples, levels of l-Isoleucine, isobutyric acid, threonic acid, and _N_-acetyl-d-glucosamine were considerably higher than those of healthy subjects and positively correlated with the HSD subtype. Level of stearic acid was also significantly higher than that of healthy controls but negatively correlated with the HSD subtype. Level of α-lactose was significantly lower than that of healthy controls and negatively correlated with the HSD subtype.

Between LQS and HSD subtypes

10 metabolites showed correlations with both subtypes (Fig. 6). They are urine α-lactose, _N_-acetyl-d-glucosamine, and stearic acid; serum l-lactic acid, l-lysine, l-valine, and l-leucine; and isobutyric acid, threonic acid, and l-methionine in both serum and urine samples. Most metabolite correlates with the two subtypes were consistent. Serum threonic acid correlate was negative with the LQS subtype but positive with the HSD subtype. Urine stearic acid correlate was positive with the LQS subtype but negative with the HSD subtype. l-Methionine correlate was positive with the LQS subtype in urine samples but negative with the HSD subtype in serum samples. The LQS subtype had a markedly lower urinary level of isobutyric acid than the HSD subtype (P = 0.047).

Metabolite-associated biochemical pathways

The impact pathway analysis revealed that the differential metabolites of the LQS subtype were significantly associated with aminoacyl-tRNA biosynthesis; valine, leucine and isoleucine biosynthesis; butanoate metabolism; phenylalanine metabolism; and glycine, serine and threonine metabolism (Fig. 7a). The differential metabolites of the HSD subtype were heavily involved in aminoacyl-tRNA biosynthesis; valine, leucine and isoleucine biosynthesis; and valine, leucine and isoleucine degradation (Fig. 7b).

Fig. 7

figure 7

Biochemical pathway analysis shows biological impacts in differentiating Liver Qi Stagnation (LQS, a) and Heart and Spleen Deficiency (HSD, b) subtypes of MDD from healthy controls. X- and Y-axis indicates the magnitude of the impact and p value compared with healthy controls, respectively. Aminoacyl-tRNA biosynthesis and valine, leucine and isoleucine biosynthesis pathways are most closely associated with the two subtypes. The HSD subtype was additionally associated with valine, leucine and isoleucine degradation

Full size image

The optimal metabolite profiles for separating TCM subtypes and healthy controls

ROC analysis showed that the optimal metabolite profile with a combination of five serum metabolites (leucine, valine, phenylalanine, threonic acid, and glycine) well separated the LQS subtype and healthy controls with a sensitivity of 0.87 and a specificity of 0.80 (Fig. 8a). A combination of five serum metabolites (leucine, valine, threonine, methionine, and inositol) optimally dichotomized the HSD subtype and healthy controls with a sensitivity of 0.95 and a specificity of 0.88 (Fig. 8b). There were no optimal urinary metabolite profiles that could separate TCM subtypes and healthy controls at the acceptable sensitivity and specificity (data not shown).

Fig. 8

figure 8

Receiver operating curve (ROC) analysis show the optimal serum metabolite panels that well differentiate Liver Qi Stagnation (LQS, a) and Heart and Spleen Deficiency (HSD, b) subtype from healthy controls

Full size image

Discussion

This study represents the first attempt to introduce TCM-based diagnostic subtypes into the modern classification of depressive disorders. It not only could improve personalized treatment of depression, but also promote the integration of TCM theoretical concepts into the modern biomedical knowledge system. The sample recruited in this study was representative of MDD patients seen in China [7, 9, 11, 16].

The purpose of this study was to explore blood and urinary metabolomic evidence that could validate TCM diagnostic classification of MDD. A total of 28 meaningful metabolites were identified as potential biomarkers for differentiating MDD and healthy subjects. These metabolites mainly included amino acid metabolites involved in the synthesis of monoamine and amino acid neurotransmitters (l-tryptophan, phenylalanine, _o_-tyrosine, glycine, 3-hydroxybutric acid, and _N_-acetyl-l-aspartic acid), proteinogenic branched-chain amino acids (BCAAs) (valine, leucine, and isoleucine), oxaloacetate-derived amino acids (l-lysine, l-threonine and l-methionine) which are mainly involved in the biosynthesis of proteins and energy production via gluconeogenesis, and energy metabolism-related metabolites (l-lactic acid, palmitic acid, citric acid, stearic acid, pyruvic acid, and α-lactose) [24].

Clustering analyses revealed that the two TCM subtypes of MDD were well separated from healthy controls on serum metabolic profiles, but not on the urinary profiles. Likewise, ROC analysis displayed that the two optimal metabolite combination patterns for dichotomizing the TCM subtypes of MDD and healthy subjects were identified from serum samples, but not from urine samples. Similar results also have been observed in previous studies that attempted to identify blood and urine metabolic biomarkers for the severity of MDD [16, 20]. Blood metabolites seem to be more sensitive than urinary metabolites in distinguishing MDD from healthy subjects. This may be largely because blood metabolites can be maintained at a relatively constant state compared to urinary metabolites [25].

Although clustering analyses showed large overlaps in serum and urinary metabolite patterns between the two TCM subtypes of MDD, the LQS subtype had strong mixed correlations with multiple amino acid metabolites involved in the biosynthesis of monoamine and amino acid neurotransmitters, including phenylalanine, 3-hydroxybutric acid, _o_-tyrosine, glycine, l-tryptophan, and _N_-acetyl-l-aspartic acid; but the HSD subtype did not. ROC analysis further indicated that phenylalanine and glycine were among the optimal metabolite combination pattern for discriminating the LQS subtype and healthy controls. A preliminary study also has observed differences in multiple brain monoamine neurotransmission among TCM subtypes of MDD [26]. According to the presence and absence of stress prior to its onset, MDD can be classified as reactive and endogenous (melancholic) subtypes [2]. Stress-associated depression has been suggested to be more closely associated with dysfunction of brain catecholamine and amino acid neurotransmitters, in particular norepinephrine, glutamate and γ-aminobutyric acid systems [27, 28]. Clinically, the LQS subtype patients often exhibit stress-related comorbid symptoms and signs, such as nervousness, irritability, agitation, and a wiry pulse [2, 10]. It thus appears that the LQS subtype may represent a subpopulation of MDD characterized by abnormal biosynthesis of monoamine and amino acid neurotransmitters and closer associations with stress-related pathophysiology.

Moreover, both subtypes had opposing correlations with serum threonic acid, serum/urine l-methionine, and urine stearic acid. Serum threonic acid and l-methionine were also among the optimal metabolite combination patterns for respectively discriminating the LQS and HSD subtype from healthy controls. Threonic acid is a metabolite of ascorbic acid (vitamin C) and methionine is an essential amino acid in humans [24]. Dietary treatment with threonic acid- and methionine-derived supplements (e.g., magnesium-L-threonate (MgT) and S-adenosyl-methionine (SAMe)) had the benefits in improving multiple brain functions, in particular memory, depression, and anxiety [29,30,31]. Stearic acid is a saturated fatty acid which plasma concentrations have been found to be abnormally altered in patients with winter episode of MDD [32] and in the superior temporal gyrus and orbitofrontal cortex of bipolar patients [33, 34]. We further found that the LQS subtype had a markedly lower urinary level of isobutyric acid, a short-chain fatty acid than the HSD subtype. Plasma level of β-amino-isobutyric acid of autistic children has been shown to be significantly higher than that of healthy children [35]. These results suggest that threonic acid, methionine, stearic acid, and isobutyric acid are differentially associated with the two TCM subtypes and may be promising biomarkers for TCM classification of MDD.

However, the present study revealed that both subtypes had large overlaps in serum and urinary metabolite panels and shared several similar profiles, in particular proteinogenic BCAAs (valine and leucine) and energy metabolism-related metabolites (α-lactose and l-lactic acid), manifesting as consistent positive correlations with serum levels of valine and leucine, and negative correlations with levels of urine α-lactose and serum l-lactic acid. Biochemical pathway analysis displayed that both subtypes had consistent and closer links with valine, leucine, and isoleucine biosynthesis and aminoacyl-tRNA biosynthesis. ROC analysis also revealed that valine and leucine were involved in the two optimal metabolite combination patterns identified for discriminating the two subtypes and healthy controls. BCAAs and energy metabolism dysfunction is related to the pathogenesis of multiple psychiatric disorders and a common comorbidity in depressive disorders [36, 37]. The BCAAs play a key role in maintaining the supply of the neurotransmitters glutamate and γ-aminobutyric acid (GABA) [38]. Our study suggests that aberrant BCAA and energy metabolism seems to be distinct biochemical profiles in the pathophysiology of MDD and the BCAA-glutamate/GABA metabolic cycle in the brain may be potential biomarkers and therapeutic targets for MDD [36, 38, 39].

Several limitations of this study should be noted. First, there are 12 subtypes according to TCM classification of MDD [7,8,9,10,11], but we only examined the two most common and opposing subtypes of MDD. Whether other subtypes also have differential profiles needs further investigation. Second, the TCM diagnostic procedure was conducted mainly based on empirical practice and subjective assessments. Objective and precise TCM diagnostic approaches should be sought. Third, the sample size was relatively small with medication-free participant and gender imbalance. This may limit the applicability of the findings of this study. Whether similar results could be achieved in first-episode and treatment-resistant depressed populations needs further determination. Finally, we did not examine the differential response of the two subtypes to antidepressant treatment. Whether metabolomic biomarkers identified are effective in predicting TCM subtype-specific treatment outcomes deserves further examination.

Conclusions

While the two subtypes are associated with aberrant BCAA and energy metabolism, the LQS subtype may represent a subpopulation of MDD characterized by abnormal biosynthesis of monoamine and amino acid neurotransmitters and closer associations with stress-related pathophysiology. Threonic acid, methionine, stearic acid, and isobutyric acid are differentially associated with the two subtypes and may be promising biomarkers for TCM classification of MDD. This study suggests that the TCM diagnostic subtypes could serve a valid classification of MDD.

Abbreviations

ANOVA:

analysis of variance

BCAAs:

proteinogenic branched-chain amino acids

DSM-5:

Statistical Manual of Mental Disorders, Fifth Edition

GABA:

γ-aminobutyric acid

GC–MS:

gas chromatography-mass spectrometry

HAMD-24:

24-item Hamilton Rating Scale for Depression

HSD:

Heart and Spleen Deficiency

LQS:

Liver Qi Stagnation

MDD:

major depressive disorder

OPLS-DA:

latent structures-discriminant analysis

PCA:

principal component analysis

ROC:

receiver operating characteristic

TCM:

traditional Chinese medicine

VIP:

variable importance in the projection

References

  1. Kennedy SH. Core symptoms of major depressive disorder: relevance to diagnosis and treatment. Dialogues Clin Neurosci. 2008;10:271–7.
    PubMed PubMed Central Google Scholar
  2. Malki K, Keers R, Tosto MG, et al. The endogenous and reactive depression subtypes revisited: integrative animal and human studies implicate multiple distinct molecular mechanisms underlying major depressive disorder. BMC Med. 2014;12:73.
    Article Google Scholar
  3. van Loo HM, de Jonge P, Romeijn JW, Kessler RC, Schoevers RA. Data-driven subtypes of major depressive disorder: a systematic review. BMC Med. 2012;10:156.
    Article Google Scholar
  4. Ilgen MA, Downing K, Zivin K, Hoggatt KJ, Kim HM, Ganoczy D, Austin KL, McCarthy JF, Patel JM, Valenstein M. Exploratory data mining analysis identifying subgroups of patients with depression who are at high risk for suicide. J Clin Psychiatry. 2009;70:1495–500.
    Article Google Scholar
  5. Milaneschi Y, Lamers F, Peyrot WJ, Abdellaoui A, Willemsen G, Hottenga JJ, Jansen R, Mbarek H, Dehghan A, Lu C, CHARGE inflammation working group, Boomsma DI, Penninx BW. Polygenic dissection of major depression clinical heterogeneity. Mol Psychiatry. 2016;21:516–22.
    Article CAS Google Scholar
  6. Cheung F. Modern TCM: enter the clinic. Nature. 2011;480:S94–5.
    Article CAS Google Scholar
  7. Chen Z, Hu S, Zhang H, Zhang H, Zheng L, Li X. Study on criteria for commonly seen TCM syndrome of depression. JTCM. 2005;46:47–9 (Article in Chinese).
    Google Scholar
  8. Guo R, Yu M, Wang JL, et al. Establishing a TCM pattern element identification scale for depression. J Beijing Univ TCM. 2015;38:561–5 (Article in Chinese).
    Google Scholar
  9. Hu SY, Wang SE, Zhang CH. Relationship between symptom stratification and syndrome differentiation of traditional Chinese medicine for depressive episode. Zhong Xi Yi Jie He Xue Bao. 2011;9:933–6 (Article in Chinese).
    Article Google Scholar
  10. China Academy of Chinese Medical Sciences and China Association for Acupuncture and Moxibustion. Evidence-based guidelines of clinical practice in Chinese medicine: acupuncture. Beijing: China Traditional Chinese Medicine Publishing House; 2011. p. 107.
    Google Scholar
  11. Zhou L, Chen WK, Mei XY. Study on the characteristics of population distribution of TCM syndromes and its related factors in patients of depression. Zhongguo Zhong Xi Yi Jie He Za Zhi. 2006;26:106–9 (Article in Chinese).
    PubMed Google Scholar
  12. Yeung WF, Chung KF, Zhang NL, Zhang SP, Yung KP, Chen PX, Ho YY. Identification of Chinese medicine syndromes in persistent insomnia associated with major depressive disorder: a latent tree analysis. Chin Med. 2016;11:4.
    Article Google Scholar
  13. Zhang YF, Han Y, Wang YZ, et al. Characterization of resting-state fMRI-derived functional connectivity in patients with deficiency versus excess patterns of major depression. Complement Ther Med. 2015;23:7–13.
    Article Google Scholar
  14. Fu S, Xie ZY, Yu CJ, Gao XY. Methods to establishment of serum induced syndrome of liver-depression and spleen-deficiency cell model. Zhongguo Zhong Yao Za Zhi. 2018;43:2999–3005.
    PubMed Google Scholar
  15. Peng GJ, Tian JS, Gao XX, Zhou YZ, Qin XM. Research on the pathological mechanism and drug treatment mechanism of depression. Curr Neuropharmacol. 2015;13:514–23.
    Article CAS Google Scholar
  16. Chen JJ, Zhou CJ, Zheng P, et al. Differential urinary metabolites related with the severity of major depressive disorder. Behav Brain Res. 2017;332:280–7.
    Article CAS Google Scholar
  17. Ding X, Yang S, Li W, et al. The potential biomarker panels for identification of Major Depressive Disorder (MDD) patients with and without early life stress (ELS) by metabonomic analysis. PLoS ONE. 2014;9:e97479.
    Article Google Scholar
  18. Liu X, Zheng P, Zhao X, et al. Discovery and validation of plasma biomarkers for major depressive disorder classification based on liquid chromatography-mass spectrometry. J Proteome Res. 2015;14:2322–30.
    Article CAS Google Scholar
  19. Liu Y, Yieh L, Yang T, et al. Metabolomic biosignature differentiates melancholic depressive patients from healthy controls. BMC Genomics. 2016;17:669.
    Article Google Scholar
  20. Setoyama D, Kato TA, Hashimoto R, et al. Plasma metabolites predict severity of depression and suicidal ideation in psychiatric patients—a multicenter pilot analysis. PLoS ONE. 2016;11:e0165267.
    Article Google Scholar
  21. Hamilton M. A rating scale for depression. J Neurol Neurosurg Psychiatry. 1960;23:56–62.
    Article CAS Google Scholar
  22. Ye J, Cai X, Cao P. Problems and prospects of current studies on the microecology of tongue coating. Chin Med. 2014;9:9.
    Article Google Scholar
  23. Zheng P, Chen JJ, Huang T, et al. A novel urinary metabolite signature for diagnosing major depressive disorder. J Proteome Res. 2013;12:5904–11.
    Article CAS Google Scholar
  24. Nelson DL, Cox MM. Principles of Biochemistry. 4th ed. New York: W. H. Freeman; 2005.
    Google Scholar
  25. Barrios C, Spector TD, Menni C. Blood, urine and faecal metabolite profiles in the study of adult renal disease. Arch Biochem Biophys. 2016;589:81–92.
    Article CAS Google Scholar
  26. Kim KN, Li YH, Xiang TY. TCM syndrome differentiation of depression and the clinical S-Et analysis. Chin J Exp Tradit Med Formulae. 2010;16:215–8.
    Google Scholar
  27. Gao SF, Bao AM. Corticotropin-releasing hormone, glutamate, and γ-aminobutyric acid in depression. Neuroscientist. 2011;17:124–44.
    Article CAS Google Scholar
  28. Moret C, Briley M. The importance of norepinephrine in depression. Neuropsychiatr Dis Treat. 2011;7(Suppl 1):9–13.
    CAS PubMed PubMed Central Google Scholar
  29. Huang Y, Huang X, Zhang L, et al. Magnesium boosts the memory restorative effect of environmental enrichment in Alzheimer’s disease mice. CNS Neurosci Ther. 2018;24:70–9.
    Article CAS Google Scholar
  30. Sharma A, Gerbarg P, Bottiglieri T, et al. S-Adenosylmethionine (SAMe) for neuropsychiatric disorders: a clinician-oriented review of research. J Clin Psychiatry. 2017;78:e656–67.
    Article Google Scholar
  31. Slutsky I, Abumaria N, Wu LJ, et al. Enhancement of learning and memory by elevating brain magnesium. Neuron. 2010;65:165–77.
    Article CAS Google Scholar
  32. Otoki Y, Hennebelle M, Levitt AJ, et al. Plasma phosphatidylethanolamine and triacylglycerol fatty acid concentrations are altered in major depressive disorder patients with seasonal pattern. Lipids. 2017;52:559–71.
    Article CAS Google Scholar
  33. McNamara RK, Jandacek R, Rider T, et al. Deficits in docosahexaenoic acid and associated elevations in the metabolism of arachidonic acid and saturated fatty acids in the postmortem orbitofrontal cortex of patients with bipolar disorder. Psychiatry Res. 2008;160:285–99.
    Article CAS Google Scholar
  34. McNamara RK, Rider T, Jandacek R, Tso P. Abnormal fatty acid pattern in the superior temporal gyrus distinguishes bipolar disorder from major depression and schizophrenia and resembles multiple sclerosis. Psychiatry Res. 2014;215:560–7.
    Article CAS Google Scholar
  35. Zaki MM, Abdel-Al H, Al-Sawi M. Assessment of plasma amino acid profile in autism using cation-exchange chromatography with postcolumn derivatization by ninhydrin. Turk J Med Sci. 2017;47:260–7.
    Article Google Scholar
  36. Baranyi A, Amouzadeh-Ghadikolai O, von Lewinski D, et al. Branched-chain amino acids as new biomarkers of major depression—a novel neurobiology of mood disorder. PLoS ONE. 2016;11:e0160542.
    Article Google Scholar
  37. Zuccoli GS, Saia-Cereda VM, Nascimento JM, Martins-de-Souza D. The energy metabolism dysfunction in psychiatric disorders postmortem brains: focus on proteomic evidence. Front Neurosci. 2017;11:493.
    Article Google Scholar
  38. Cole JT, Sweatt AJ, Hutson SM. Expression of mitochondrial branched-chain aminotransferase and α-keto-acid dehydrogenase in rat brain: implications for neurotransmitter metabolism. Front Neuroanat. 2012;6:18.
    Article CAS Google Scholar
  39. Mathews DC, Henter ID, Zarate CA. Targeting the glutamatergic system to treat major depressive disorder: rationale and progress to date. Drugs. 2012;72:1313–33.
    Article CAS Google Scholar

Download references

Authors’ contributions

LYL and ZZJ were involved in conception and design of the study, data analysis, and preparation of the manuscript. HJZ, JBP, and WQL were responsible for metabolomic experiments and data analysis. LYLuo, CQZ, YPL, PRW, YYZ, and CYY recruited subjects, conducted clinical assessment, and collected clinical data. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The materials and data of this study are available to other researchers upon request.

All authors have provided consent for publication in the journal of Chinese Medicine.

This study was approved by the Medical Ethical Committee of Tongde Hospital and registered at http://www.clinicaltrials.gov (NCT02346682) before the recruitment started. All participants gave voluntary, written, informed consent before entering the study.

Funding

This study was supported by the National Natural Science Foundation of China (81403502, L.Y.L.) and Research Grant Council of Hong Kong (17115017, 17124119, Z.J.Z.).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

  1. Department of Psychiatry, Tongde Hospital of Zhejiang Province, Hangzhou, 310012, Zhejiang, China
    Lan-Ying Liu, Li-Yuan Luo, Chun-Qin Zhu, Pei-Rong Wang, Yuan-Yuan Zhang & Chun-Yu Yang
  2. Zhejiang Academy of Traditional Chinese Medicine, Hangzhou, 310007, Zhejiang, China
    Hong-Jian Zhang, Jin-Bao Pu & Wei-Qing Liang
  3. Department of Internal Chinese Medicine, Tongde Hospital of Zhejiang Province, Hangzhou, 310012, Zhejiang, China
    Ya-Ping Li
  4. School of Chinese Medicine, LKS Faculty of Medicine, The University of Hong Kong, 10 Sassoon Road, Pokfulam, Hong Kong, China
    Zhang-Jin Zhang

Authors

  1. Lan-Ying Liu
    You can also search for this author inPubMed Google Scholar
  2. Hong-Jian Zhang
    You can also search for this author inPubMed Google Scholar
  3. Li-Yuan Luo
    You can also search for this author inPubMed Google Scholar
  4. Jin-Bao Pu
    You can also search for this author inPubMed Google Scholar
  5. Wei-Qing Liang
    You can also search for this author inPubMed Google Scholar
  6. Chun-Qin Zhu
    You can also search for this author inPubMed Google Scholar
  7. Ya-Ping Li
    You can also search for this author inPubMed Google Scholar
  8. Pei-Rong Wang
    You can also search for this author inPubMed Google Scholar
  9. Yuan-Yuan Zhang
    You can also search for this author inPubMed Google Scholar
  10. Chun-Yu Yang
    You can also search for this author inPubMed Google Scholar
  11. Zhang-Jin Zhang
    You can also search for this author inPubMed Google Scholar

Corresponding author

Correspondence toZhang-Jin Zhang.

Additional files

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Liu, LY., Zhang, HJ., Luo, LY. et al. Blood and urinary metabolomic evidence validating traditional Chinese medicine diagnostic classification of major depressive disorder.Chin Med 13, 53 (2018). https://doi.org/10.1186/s13020-018-0211-z

Download citation

Keywords