Cerebrospinal Fluid Metabolome in Mood Disorders-Remission State has a Unique Metabolic Profile

doi:10.1038/srep00667

Sci Rep. 2012; 2: 667.

Published online 2012 September 19. doi: 10.1038/srep00667

PMCID: PMC3446657

Cerebrospinal Fluid Metabolome in Mood Disorders-Remission State has a Unique Metabolic Profile

Rima Kaddurah-Daouk,^a¹ Peixiong Yuan,^2,⁹ Stephen H. Boyle,^1,⁹ Wayne Matson,³ Zhi Wang,⁴ Zhao Bang Zeng,⁴ Hongjie Zhu,⁴ George G. Dougherty,^5,⁶ Jeffrey K. Yao,^5,⁶ Guang Chen,² Xavier Guitart,² Paul J. Carlson,² Alexander Neumeister,^2,¹⁰ Carlos Zarate,² Ranga R. Krishnan,^1,⁷ Husseini K. Manji,^2,⁹ and Wayne Drevets^2,^8,⁹

¹Department of Psychiatry and Behavioral Sciences, Duke University, Durham, NC

²NIH/National Institute of Mental Health Intramural Research Program, Bethesda, MD

³Department of Systems Biochemistry, Bedford VA Medical Center, Bedford, MA

⁴Department of Statistics & Bioinformatics Research Center, North Carolina State University, Raleigh, NC

⁵VA Pittsburgh Healthcare System, Pittsburgh, PA

⁶Department of Psychiatry, University of Pittsburgh Medical Center, Western Psychiatric Institute & Clinic, Pittsburgh, PA

⁷Duke-NUS Graduate Medical School, Singapore

⁸Laureate Institute for Brain Research, and Department of Psychiatry, Oklahoma University, Tulsa, OK

⁹Contributed equally as second authors or as senior authors.

¹⁰Current address: New York University, Department of Psychiatry, New York, NY.

^aEmail: Kaddu001/at/mc.duke.edu

Received October 7, 2011; Accepted August 23, 2012.

This work is licensed under a Creative Commons Attribution-NonCommercial-ShareALike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/

Abstract

Targeted metabolomics provides an approach to quantify metabolites involved in specific molecular pathways. We applied an electrochemistry-based, targeted metabolomics platform to define changes in tryptophan, tyrosine, purine and related pathways in the depressed and remitted phases of major depressive disorder (MDD). Biochemical profiles in the cerebrospinal fluid of unmedicated depressed (n = 14; dMDD) or remitted MDD subjects (n = 14; rMDD) were compared against those in healthy controls (n = 18; HC). The rMDD group showed differences in tryptophan and tyrosine metabolism relative to the other groups. The rMDD group also had higher methionine levels and larger methionine-to-glutathione ratios than the other groups, implicating methylation and oxidative stress pathways. The dMDD sample showed nonsignificant differences in the same direction in several of the metabolic branches assessed. The reductions in metabolites associated with tryptophan and tyrosine pathways in rMDD may relate to the vulnerability this population shows for developing depressive symptoms under tryptophan or catecholamine depletion.

The biochemical underpinnings of major depressive disorder (MDD) and of this condition's phenotypically distinct illness phases, namely depression versus remission, have only been partly characterized. While multiple neurotransmitter systems appear dysregulated in MDD, the greatest amount of empirical research has focused upon the role of impairments in central monoaminergic function¹^,2^,3^,4. A major limitation of previous studies, however, is that they have not examined multiple compounds within a pathway or multiple interconnected pathways simultaneously.

The application of novel metabolomics technology enables mapping of biochemical pathways implicated in disease⁵^,6^,7^,8^,9^,10. Through metabolomic profiling researchers recently have identified initial metabolic signatures for a variety of central nervous system (CNS) disorders, including depression, schizophrenia, and Parkinson's disease¹¹^,12^,13^,14^,15^,16^,17^,18. For example, preliminary metabolomic studies using plasma of MDD¹¹^,12 patients revealed that the depressed state is associated with changes in amino acids, neurotransmitters including gamma-aminobutyric acid (GABA), dicarboxylic fatty acids and lipids. Nevertheless, no previous metabolite profiling study has been performed on the cerebrospinal fluid (CSF), the biofluid believed to be most closely linked to brain function¹⁹.

Previous studies examined CSF concentrations of the primary metabolites of serotonin, dopamine and norepinephrine (i.e., 5-hydroxyindoleacetic acid [5-HIAA], homovanillic acid [HVA] and 4-hydroxy-3-methoxyphenyl glycol [MHPG], respectively) in currently depressed subjects. These studies found lower concentrations of HVA, and no significant difference in MHPG, in depressed subjects with melancholic features (i.e., a subgroup characterized by more prominent abnormalities of weight/appetite, thought, movement, and incapacity to experience pleasure) relative to controls²⁰, although variable results were obtained in samples selected more generally according to the MDD criteria alone²¹. Similarly, abnormal reductions in the CSF 5-HIAA concentrations were identified in some but not most currently-depressed samples²²^,23, although in the positive studies this abnormality sometimes was associated more specifically with clinical manifestations of aggression or suicidal behavior²⁴^,25. These previous studies did not, however, assess the expanded metabolomic profile where many monoaminergic neurotransmitter metabolites are measured quantitatively and simultaneously.

Moreover, few studies assessed monoaminergic neurotransmitter metabolites in the remitted phase, and the extant data from these studies remain inconclusive. Preliminary observations suggested that 5-HIAA levels may be abnormally reduced in a subgroup of individuals in the remitted phase of MDD²¹^,26^,27. In contrast, a previous study that assessed CSF HVA levels in remitted MDD subjects²⁷ did not detect a significant difference relative to currently depressed subjects or healthy controls .

Here we applied a targeted electrochemistry-based metabolomics platform to interrogate perturbations in the neurotransmitter pathways involving serotonin, dopamine and norepinephrine, and in related pathways for methionine and purine in both the symptomatic and recovered phases of MDD. A total of 26 metabolites were measured quantitatively and simultanousely that cover the conversion of tryptophan to 5-HIAA and kynurenine (KYN), and the conversion of tyrosine (TYR) to catecholamines and HVA. Our aim was to characterize a possible altered biochemical profile in the CSF of participants with MDD studied while either depressed or remitted.

Results

Participant demographics

The baseline demographic characteristics of the groups are shown in Table 1. The subject groups did not differ significantly in gender, age or body mass index. There were a greater proportion of African-American individuals in the dMDD group.

Table 1

Clinical and Demographic Characteristics of currently depressed subjects with major depressive disorder (MDD), MDD subjects in full remission, and healthy controls

Metabolic signatures for dMDD and rMDD - Analysis within pathways

Using an electrochemistry based metabolomics platform we quantified 26 metabolites within the tryptophan, tyrosine, purine, and methionine pathway (Supplemental Tables 1,2 and Figures 1a–d). We compared levels and ratios of metabolites within pathways and across pathways between the dMDD, rMDD and HC groups using univariate and multivariate approaches.

Figure 1

(a) Tryptophan pathway metabolites quantitated by the LCECA platform. (b)Tyrosine pathway metabolites quantitated by the LCECA platform. (c) Purine pathway metabolites quantitated by the LCECA platform. (d) Methionine pathway metabolites quantitated by (more ...)

Tryptophan pathway

Analyses of the tryptophan pathway metabolites (N = 4, Table 2) and ratios (N = 4) (Supplemental Tables 1, 2 for names of metabolites and abbreviations) showed significant main effects of group for the 5-HIAA concentration (p = .01, Figure 2a) and the 5-HIAA/TRYP (p = .008, Figure 2b.) and 5-HIAA/KYN (p = .0057) ratios. Follow-up tests (Bonferroni-adjusted alpha = .017 for significance) revealed that the rMDD group had significantly lower levels of 5-HIAA than both the dMDD group (p = .007) and the HC group (p = .014) and significantly smaller 5-HIAA/TRYP and 5-HIAA/KYN ratios than both the dMDD group (ps = .006 and .002, respectively) and the HC group (ps = .006 and .007, respectively).

Table 2

Kruskal-Wallis One-Way ANOVA comparing remitted MDD, depressed MDD, and Healthy Control samples

Figure 2

(a) 5HIAA levels in healthy control, depressed MDD, and remitted MDD samples. (b) 5-HIAA/TRP ratios in healthy control, depressed MDD, and remitted MDD samples. (c) HVA levels in healthy control, depressed MDD, and remitted MDD samples. (d) TYR/4HPLA (more ...)

Tyrosine pathway

Analyses of the tyrosine pathway metabolites (N = 8, Table 2) and ratios (N = 10) (Supplemental Tables 1, 2) showed significant main effects of group for HVA concentrations (p = 0.02, See Figure 2c.) and the TYR/4HPLA ratio (p<0.02, See Figure 2d.). The test for HVA/MHPG (p = 0.01, Bonferroni alpha = 0.008) was not considered significant using the Bonferroni adjustment, but was considered significant when correction for pFDR was applied. Follow-up tests (Bonferroni alpha = .017 for significance) revealed that the rMDD group had significantly lower levels of HVA than the HC group (p = 0.007). Examination of the ratios revealed that the rMDD group had a significantly larger TYR/4HPLA ratio than the HC group (p = 0.008). The rMDD group also had a significantly smaller HVA/MHPG ratio (p = 0.003) than the HC group. There was no significant difference between the HC and dMDD groups.

Examination of the unadjusted p-values associated with the Kruskal-Wallis tests revealed a non significant trend toward a group difference in the HVA/TYR ratio (p = .024). The rMDD group had a smaller HVA/TYR ratio than the HC group.

Purine pathway

None of the tests of the purine pathway metabolites (N = 7, Table 2) or the within pathway ratios (N = 10, Supplemental Tables 1, 2) identified significant differences. The hypoxanthine and xanthine concentrations showed a non significant trend toward being lower in the rMDD group than the HC group.

Methionine and related pathways

Analyses of the metabolites (N = 5, Table 2) and ratios (N = 2) of these pathways (Supplemental Tables 1, 2) revealed significant main effects of group for the methionine concentration (MET; p = 0.02, Figure 3a.) and the glutathione (GSH)/MET ratio (p = 0.003, Bonferroni alpha = 0.008, Figure 3b.). Follow-up tests (Bonferroni alpha = 0.017 for significance) revealed that the rMDD group had significantly higher levels of MET than the HC group (p = 0.007) and a trend toward having higher levels of MET (p = 0.03) than the dMDD group. The rMDD group also had significantly smaller GSH/MET ratios than the dMDD (p = 0.001) and HC (p = 0.002) groups. The GSH/MET effect was due to the rMDD group having both higher MET levels and lower GSH levels than the other two groups. There was no significant difference between the HC and dMDD groups.

Figure 3

(a) Methionine levels in healthy control, depressed MDD, and remitted MDD samples. (b) GSH/MET ratios in healthy control, depressed MDD, and remitted MDD samples. (c) HVA/5HIAA ratios in healthy control, depressed MDD, and remitted MDD samples. (d) XAN/HVA (more ...)

Metabolic Signatures for dMDD and rMDD - Ratios across pathways

Kruskal-Wallis tests on the selected across pathway ratios (N = 5, Table 2) showed significant effects for HVA/5HIAA (p = 0.04, Figure 3c.), XAN/HVA (p = 0.03, Figure 3b.), and 5HIAA/XAN (p = 0.007) ratios. Follow-up tests (Bonferroni alpha = 0.017 for significance) revealed that the rMDD group had significantly larger XAN/HVA ratios than the HC group (p = 0.008) and a significantly larger HVA/5HIAA ratio than the dMDD group (p = 0.014). The rMDD group also had a significantly smaller 5HIAA/XAN ratio than both the HC (p = 0.014) and dMDD (p = 0.003) groups.

Associations of metabolites with disease severity indices

There were several significant correlations between levels of metabolites (N = 26) and indices of disease severity in the dMDD group. TPOL was negatively associated with HAM-D-17 scores (r = −.61, p = .036) indicating that lower levels of this compound are associated with higher levels of melancholic symptoms. TYRA was negatively associated with the HAM-D-A1-8 scores (r = −.58, p = .049) indicating that patients with lower levels of TYRA have more severe symptoms of atypical depression (a phenotype characterized by increases in appetite/weight and sleep, and preservation of mood reactivity in response to pleasurable events). HAM-A scores were positively associated with 5HIAA (r = .60, p = .023), HX (r = .55, p = .04), and MHPG (r = .65, p = .01) and negatively associated with TRYP (r = −.59, p = .026). This suggests that dMDD patients with more severe symptoms of anxiety have higher levels of 5-HIAA, HX, and MHPG and lower levels of TRYP.

Discussion

We compared metabolites from the tryptophan, tryrosine, purine and methionine pathways between rMDD and dMDD subjects and healthy controls using the LCECA platform. This approach detects a subset of the metabolome, namely compounds amenable to oxidation-reduction, with greater sensitivity for studying monoamine neurotransmitter pathways than can be achieved using other technologies. Since cellular metabolism reflects the integrated interconversion of metabolic substrates through enzyme-catalyzed biochemical reactions, an alteration in the enzyme function at one reaction can affect the fluxes of one or several subsequent reactions²⁸. In this CSF metabolite pathway analysis we observed that the levels of several metabolites within the investigated pathways were altered in unmedicated-remitted individuals with MDD. In contrast, unmedicated-depressed individuals did not exhibit statistically significant perturbations within the same metabolic branches, consistent with the results of most previous studies on 5-HIAA, HVA, and purine metabolites in depressed individuals selected according to similar diagnostic criteria.

The counterintuitive findings that the CSF 5-HIAA levels and the 5-HIAA/TRYP ratio were significantly lower in the rMDD versus the dMDD subjects conceivably may reflect greater biochemical heterogeneity within the dMDD group. One source of this heterogeneity may arise from variation in symptom severity, as suggested by the associations observed between the tryptophan, tyrosine, and purine pathway metabolites and depression or anxiety ratings. For example, Nordin²⁹ found that the relationship between 5-HIAA concentrations and depressive symptoms in MDD was curvilinear, such that individuals with the lowest and highest 5-HIAA levels had the lowest depression ratings. Our data appear compatible with this finding, since most of the currently depressed subjects we studied had depression ratings in the moderate severity range. Moreover, the observation that anxiety ratings correlated positively with 5-HIAA levels but negatively with TRYP levels suggests that the greater anxiety symptoms associated with active depressive illness may contribute to the higher 5-HIAA levels and the higher 5-HIAA/ TRYP ratio in the currently ill versus the currently remitted MDD samples. In preclinical studies of anxiety, rats exposed to stressors or threats show adaptive increases in serotonin turnover in several brain structures³⁰, which putatively would increase the CSF 5-HIAA levels. Thus the anxiety symptoms associated with depressive episodes may result in an elevation of 5-HT turnover, potentially obscuring a basal deficit in 5-HIAA levels that is evident during remission.

Another potential explanation for the more prominent difference in the remitted sample is that the entrance criteria used to select the rMDD subjects yielded a more homogenous subphenotype of MDD. The rMDD subjects were characterized by manifesting an episodic course (i.e., episodes of illness separated by epochs of remission) and the ability to remain in remission while unmedicated. In contrast, while the depressed sample included some subjects with an episodic course, most of the dMDD subjects manifested chronic illness. Thus, one interpretation of the marked reduction in 5-HIAA and HVA levels in the rMDD subjects is that the entrance criteria for this sample identifies a subphenotype that is more likely to maintain euthymia because of having a lower dopamine and serotonin turnover rate, potentially due to genetic variation that results in slower metabolism of these neurotransmitters (so that their availability within the synapse is prolonged). Compatible with such a hypothesis, drugs that achieve such a condition through monoamine oxidase inhibition exert potent antidepressant effects. Thus the persistently remitted group may manifest biological differences that confer resilience and prevent illness chronicity in the course of MDD. Nevertheless, the cross-sectional design of our study did not allow us to distinguish whether these CSF abnormalities reflect trait-like abnormalities or compensatory changes arising in response to past illness.

In either case the findings in remitted MDD subjects appear compatible with other evidence that, despite symptom remission, such individuals show behavioral and biochemical deficits relative to never-depressed controls. Clinically, unmedicated-remitted MDD subjects show a greater risk of depressive relapse, higher rates of suicide attempts and ideation, lower pleasure ratings and impaired socio-occupational functioning relative to controls³¹. Such individuals also show neurobiological abnormalities that in some cases resemble³², but in other cases contrast with, those shown by their currently-depressed counterparts. For example, under acute tryptophan depletion or catecholamine depletion, unmedicated rMDD subjects redevelop depressive symptoms, neuropsychological deficits, and neurophysiological abnormalities characteristically seen in the depressed phase³³^,34^,35, while unmedicated-depressed MDD subjects show no exacerbation of depressive symptoms under these challenges³⁶^,37. These findings suggest that alterations in the tryptophan or tyrosine pathways limit the capacity to remain in remission in rMDD subjects, consistent with their abnormal reductions in CSF 5-HIAA and HVA levels shown herein.

The cross-sectional design also did not allow us to exclude the possibility that persistent effects of previous treatment influenced our findings, since only a minority of the subjects was treatment-naïve. Nevertheless, while several studies reported decreased CSF concentrations of 5-HIAA following acute or subacute administration of selective serotonin reuptake inhibitors (SSRIs)³⁸^,39^,40^,41^,42, after chronic SSRI treatment the CSF 5-HIAA levels return to the pre-treatment baseline⁴³. Similarly, HVA levels are not changed significantly by chronic SSRI treatment⁴³. Thus our findings more likely relate to the pathophysiology of MDD, and as such may underlie predisposing factors or adaptive responses to illness, as opposed to homeostatic mechanisms arising in response to past treatment⁴⁴. For example, the abnormal reduction in CSF 5-HIAA levels in the rMDD sample is consistent with evidence that MDD is associated with neuropathological involvement of brainstem serotonergic nuclei⁴⁵, and that maintenance treatment with drugs that increase 5-HT transmission reduce the risk of illness recurrence in rMDD patients⁴⁶^,47.

Our findings in the rMDD sample appear compatible with those of previous studies that measured CSF 5-HIAA in isolation from the remainder of the tryptophan pathway. van Praag and de Hann⁴⁸ reported that 33 of 54 depressed individuals with low post probenecid accumulation of CSF 5-HIAA also showed this low accumulation after recovery. Furthermore, Vestergaard et al.²⁷ reported that CSF 5-HIAA and HVA levels were significantly lower following recovery than during depression in 16 mood disordered subjects, although the differences in the 5-HIAA levels with respect to healthy controls did not reach significance in either the depressed or the recovered phases. Our data thus were the first to document that remitted MDD subjects show significantly reduced mean 5-HIAA levels relative to controls in the absence of probenicid challenge.

The significantly smaller 5HIAA/TRP and 5HIAA/KYN ratios in the rMDD group versus the dMDD and HC groups are noteworthy, in light of evidence that the pathophysiology of MDD involves activation of pro-inflammatory pathways⁴⁹. Activation of the KYN pathway is initiated by pro-inflammatory cytokines via induction of the enzyme, indoleamine 2,3-dioxygenase (IDO), which metabolizes TRP into KYN. Activation of the KYN pathway is hypothesized to contribute to reductions in 5-HT synthesis in MDD and in individuals who develop depressive symptoms during cytokine administration⁵⁰^,51. Potentially compatible with this hypothesis, the reduced 5HIAA/KYN ratio in the rMDD sample suggests that more TRP is being shunted toward the KYN pathway than toward 5-HT synthesis. These findings thus appear compatible with evidence that a pro-inflammatory state exists in remitted MDD subjects, despite their recovery from depressive symptoms⁵².

Within the tyrosine pathway, the rMDD group showed significantly higher TYR/4HPLA ratios, lower TYRA levels, and smaller HVA/MHPG and HVA/TYR ratios versus the HC group. This combination of findings suggests that of the TYR routed toward catecholamine synthesis, relatively more is being routed toward norepinephrine utilization than toward dopamine utilization (figure 1b). Given the role of dopamine in reward processing, these data together with the abnormal reduction in HVA levels in the rMDD group appear compatible with evidence that hypohedonia (diminished ability to enjoy pleasurable activities) is a trait-like marker of MDD⁵³. In addition, the quartile data (table 2) suggested that a subset of the currently depressed subjects also manifested reduced HVA levels, consistent with previous reports that subjects with the melancholic subtype of MDD show reductions in the CSF HVA concentration (see Introduction).

Our data further suggest that a balance between dopaminergic and serotonergic function may be relevant to the switch between depression and remission. Previously Roy et al.⁵⁴ reported that the CSF HVA-to-5-HIAA ratio was lower in depressed MDD subjects than in controls. In our study the significant difference in the HVA/5HIAA ratio across groups (table 2) was attributable to a lower ratio in the dMDD subjects versus the HC and rMDD subjects (consistent with Roy et al.⁵⁴). These data imply that the ratio of dopamine turnover to serotonin turnover is lower in dMDD than in rMDD cases, raising the possibility that imbalances between 5-HT and DA neurotransmission contribute to the development of depressive symptoms. Notably the literature characterizing interactions between dopaminergic and serotonergic neurotransmission has emphasized functionally antagonistic relationships between these systems⁵⁵.

The MET concentration was higher and the GSH/MET ratio was lower in the rMDD group than in the HC and dMDD groups. These data conceivably reflect a reduction in the conversion rate of MET to GSH, possibly implicating a defect in 1-carbon metabolism in MDD. Methionine plays a critical role in protein synthesis, methylation, and polyamine biosynthesis, and its derivative S-adenosyl methionine (SAM) serves as a methyl donor for several processes including DNA methylation. A role for DNA methylation in mood regulation is suggested by the antidepressant effect of S-adenosyl methionine⁵⁶^,57^,58, the post mortem evidence that DNA methylation of genes involved in catecholamine metabolism are altered in mood disorders⁵⁹^,60^,61, and preclinical evidence that epigenetic processes underlie both the resilience against and the vulnerability toward the development of depressive behaviors under stress⁶². Our data raise the possibility that higher MET concentrations support the maintenance of symptom remission in MDD⁵⁸. Conversely, the reduction in GSH levels observed in the remitted sample conceivably may confer the vulnerability of the rMDD population to illness recurrence. Glutathione is the major antioxidant in the brain, and thus plays a key role in defending against oxidative damage. Glutathione is thought to play an important role in the neuroprotective effects of mood stabilizing drugs⁶³, given the accumulating evidence for oxidative stress mechanisms as common pathophysiological pathways in mood disorders⁶⁴^,65^,66^,67. We previously reported¹⁰ that glycine levels differed significantly between remitters and non-remitters to SSRI treatment, and proposed that glycine's link to folate and methionine cycles and methylation processes might play a role in recovery from the depressed state.

Although none of the purine metabolites showed significant effects on their own, some ratios involving the purine metabolites and the metabolites within other pathways differed across groups (i.e., XAN/HVA and 5-HIAA/XAN). Previous studies reported strong positive correlations between CSF levels of xanthine and hypoxanthine with corresponding levels of HVA and 5HIAA, despite finding no significant difference in the mean concentrations of these purine metabolites between depressives and controls⁶⁸^,69. The strong correlation in the lumbar CSF pool between HVA and XAN across the purine and tyrosine pathways has been observed in previous studies⁷⁰. This cross pathway correlation appears to extend the relationship between XAN and tyrosine/tryptophan evident in the brain and the ventricular CSF to the lumbar CSF, as the end product metabolites of tyrosine and tryptophan resulting from spinal cord metabolism. The purine and tyrosine pathways also share the pterin compounds as common cofactors, potentially contributing both to the similarities in the direction of correlation and to the differences in the ratio obtained across the rMDD group and the other two groups.

This study provides new insights about the metabolic state associated with remission from MDD, and distinguishes this state from the normative metabolic state. The relatively small sample size in our study limited the statistical sensitivity to detect other potentially meaningful differences between groups. The lack of difference between the dMDD and HC groups may reflect both the limited statistical power and the greater biochemical heterogeneity seen among chronically depressed subjects. Moreover, greater coverage of the metabolome using complementary platforms may reveal perturbations in dMDD that were not detected in our study. Finally, correlations between biochemical changes in the CSF with those observed in plasma might enable the connection of central and peripheral changes.

In conclusion, the correlation of metabolites within and across the tryptophan, tyrosine and methionine pathways revealed differences in remitted MDD subjects which suggest alterations in the regulatory feedback within these pathways. Our data thus demonstrate that metabolomics approaches hold the potential to elucidate failures in the regulation of specific biochemical pathways and networks in mood disorders. The differentiation of whether these abnormalities reflect consequences of prior depressive episodes or vulnerabilities toward the development of future mood episodes awaits longitudinal studies in at-risk samples that assess the function of these pathways before and following illness-onset. Furthermore, the biochemical heterogeneity extant within the depressed MDD population suggests that future studies should characterize sufficiently large samples that subphenotypes can be identified which manifest more homogenous metabolomic signatures (potentially analogous to those found in the rMDD sample).

Methods

Participants met DSM-IV-TR criteria⁷¹ for MDD in full remission (rMDD; n = 14) or MDD in a current major depressive episode (dMDD; n = 14). The control group was psychiatrically healthy (n = 18). Diagnosis was established using the Structured Clinical Interview for DSM-IV⁷² and an unstructured interview with a psychiatrist. To assess general health, physical examination, urine drug screening, and laboratory testing of electrolytes, blood count, thyroid function, hepatic and renal function, viral titers, and pregnancy status were performed. Volunteers were excluded if they had major medical or neurological disorders (including hypertension or endocrine disorders), past head injury, pregnancy, electrolyte disturbance, anemia, positive drug or viral (HIV, hepatitis) screen, current serious suicidal ideation or recent suicidal behavior, delusions or hallucinations, or met DSM-IV-TR criteria for substance abuse within one year or substance dependence within the lifetime. The rMDD sample additionally met the entrance criterion of having been in full remission (DSM-IV-TR) and free of exposure to psychotropic medications >3 months. The dMDD and HC subjects were excluded if they had been exposed to psychotropic or other medications likely to alter monoamine neurochemistry within 3 weeks. In no case was a medication discontinued for the purpose of this study; volunteers either were naïve to psychotropic medications or had discontinued medication for other reasons. The HC subjects met the additional exclusion criteria of having had no lifetime history of a psychiatric disorder and no first-degree relative with a mood or anxiety disorder, as established using the SCID and the Family Interview for Genetic Studies (FIGS)⁷³. Females had not received hormonal treatments other than oral contraceptives. Subjects gave written informed consent, as approved by the NIH Combined Neuroscience IRB.

The 17 item Hamilton Rating Scale for Depression (HAM-D) and the Montgomery-Asberg Depression Rating Scale were administered to provide an overall depression severity index. The 8 item “atypical feature” component from the 25 item HAM-D reflecting reversed neurovegetative changes (i.e., increased sleep, appetite and weight) also was administered. The Hamilton Anxiety Rating Scale (HAM-A) was used to rate anxiety symptom severity.

CSF collection

Baseline cerebrospinal fluid (CSF) samples were collected via lumbar puncture (LP) into the subarachnoid space between the L3-4 or L4-5 interspaces. During the 3 days prior to LP, subjects were requested to observe a monoamine free diet. All subjects were admitted to the NIMH Inpatient Unit at 7 pm the day prior to LP, and LP was performed at 9 am after overnight stay.

Metabolomics analysis: liquid chromatography electrochemical array platform

Metabolic profiling of the samples was performed using LCECA coulometric array detection, which has been previously described in other studies¹⁷^,18^,19^,74^,75^,76^,77 (See Supplement for a detailed description). This platform quantified 26 known compounds, primarily within the tyrosine, tryptophan, purine and methionine pathways. Levels of anthranilic acid and guanine were below the detectable limit of the platform in over 50% of the samples, so they were not included in the analysis. Of the remaining 24 compounds (See Supplementary Table 1), tryptophol (TPOL), L-DOPA (LD), and 7-methylxanthine (7MXAN) had fewer than 20% such samples and all other metabolites had no censored values.

Statistical analysis

Descriptive statistics were computed for each group. The raw data were viewed by quantile-quantile normal and chi-square plots, and by variable-pair scatterplots, to assess normality and nonlinear relationships. As most analytes were not approximately normally distributed, nonparametric Kruskal Wallis tests were used to test for differences across groups. A significant Kruskal-Wallis test was followed by post-hoc comparisons. When ties were present, and notably when values were censored below their minimum value, the normal approximations for the above Kruskal Wallis test statistics⁷⁸ were used in the tests. The first set of analyses focused on a series of compounds and ratios of compounds for which there was a clear basis from previous literature to be involved in the pathophysiology of MDD. These included the following compounds: tryptophan (TRYP), 5-HIAA, KYN, TYR, LD, HVA, MHPG, and methionine (MET). We also examined the following ratios: HVA/MHPG, 5HIAA/TRYP, xanthine(XAN)/HVA, HVA/5-HIAA, 5-HIAA/XAN, 5HIAA/HX, TYR/4HPLA, TYR/TRYP. Ratios of compounds inform us about the relative effectiveness of the enzymes in regulating the pathway or pathways involved. Thus, the analysis of ratios can help identify pathways that may play a role in MDD. For compounds and ratios reported to be abnormal in MDD in previous studies, we did not apply a correction for the number of Kruskal Wallis tests performed. When making post-hoc tests, a Bonferroni corrected p-value for three comparisons (.05/3 = 0.017) was used in judging significance. The imbalance in Race across the groups could affect the results for an analyte if that analyte were correlated with Race. We addressed this possibility by testing each analyte for a nonzero Kendall's tau correlation with Race, rejecting the null for p<0.05. Kendall's tau is known to handle ties well, even up to dichotomous variables. In such cases, we used McSweeney and Porter's rank ANCOVA⁷⁹ (e.g., see Barrett, 2011) with 3 groups and Race the covariate, to test for group differences, and the McSweeney p-value is presented in lieu of the Kruskal-Wallis p-value (Table 2).

The second set of analyses was exploratory and involved the 38 remaining compounds and ratios within pathways (see Table 2). The metabolites and ratios of metabolites were grouped according to their biochemical roles in four pathways (tryptophan, tyrosine, purine, methionine). Control of Type I error for this multiple testing process was effected by applying the Bonferroni correction to alpha = .05 for the (exploratory) remainder of analytes and ratios within each pathway, separately for each pathway. Here again, Kruskal Wallis or McSweeney tests of group differences were done.

Correlations of metabolites with disease severity

Spearman correlations were used to examine associations of metabolite levels to disease severity indicators in the dMDD group. Disease severity was measured with the HAM-D-17, HAM-D-A8, and HAM-A scales. The HAM-D-17 and HAM-D-A8 scales provide ratings that emphasize the melancholic and atypical features of depression, respectively. The HAM-A scale rates the severity of anxiety symptoms.

Author Contributions

R.K.D., was involved in study design, data interpretation and drafting of the manuscript. P.Y. was involved in the design of the study and data interpretation. S.H.B. was involved in data interpretation and drafting of the manuscript. W.M. profiled the CSF samples and contributed to the biochemical interpretation of the data. Z.W., Z.B.Z. and H.Z. played a role in data analysis and interpretation. G.G.D., conducted data analysis and provided input on interpretation of the data . J.Y. was involved in the biochemical interpretation of the data. G.C. was involved with research design, data acquisition, and data interpretation. X.G. was involved with research design, data acquisition, and data interpretation. P.C. oversaw the storage and organization of the cerebrospinal fluid samples and participated in preparing the data for publication. A.N. recruited the subject participants and obtained the cerebrospinal fluid samples and participated in the study design. C.Z. was involved in subject recruitment and evaluation. R.R.K. was involved in data interpretation. H.K.M. was involved with research design, data acquisition, and data interpretation. W.D. contributed to the study design, the subject clinical assessment, the interpretation of the results, and the drafting of the manuscript. He also oversaw the laboratory in which the research samples were collected and obtained funding for this part of the study. All authors reviewed the manuscript.

Supplementary Material

Supplementary Information

Supplementary Tables

Click here to view.^{(95K, pdf)}

Acknowledgments

This research was supported in part by funding from NIH R24 GM078233, “The Metabolomics Research Network” (R.K.-D.) and the Department of Veterans Affairs Senior Research Scientist Award (JKY). Earle Bain, M.D. performed many of the lumbar punctures. Michele Drevets, Ruth Tinsley and other clinical staff at Mood and Anxiety Disorder Program of NIMH assisted with patient recruitment and evaluation, and contributed to collecting and storing the CSF samples. The contents of this article do not represent the views of the Department of Veterans Affairs or the United States Government.

Footnotes

Rima Kaddurah-Daouk is equity holder in Metabolon Inc., a biotechnology company in the metabolomics domain, and also an inventor on patents in the metabolomics field. She has received funding or consultancy fees for BMS, Pfizer Inc., AstraZeneca and Lundbeck. George G. Dougherty owns shares of Pfizer and Merck. Carlos Zarate is listed as a co-inventor on a patent for the use of ketamine in major depression. Dr. Zarate has assigned his patent rights on ketamine to the U.S. government. Ranga R Krishnan has received consultancy fees from Amgen, Bristol-Myers Squibb, CeNeRx, Corcept, GlaxoSmithKline, Johnson & Johnson, Lundbeck, Merck, Organon, Pfizer, Sepracor, and Wyeth. Wayne Drevets has received consultancy fees from Pfizer, Rules Based Medicine, Eisai, Abbott, and Johnson & Johnson, and has a use patent filed on scopolamine in the treatment of depression. Peixiong Yuan, Stephen H. Boyle, Wayne Matson, Zhi Wang, Zhao-bang Zeng, Hongjie Zhu, Jeffrey Yao, Guang Chen, Xavier Guitart, Paul J. Carlson, Alexander Neumeister, and Husseini K. Manji have no relevant financial interests to disclose.

References

Delgado P. L. Common pathways of depression and pain. J Clin Psychiatry 65 Suppl 12, 16–19 (2004). [PubMed]
Delgado P. L. & Moreno F. A. Role of norepinephrine in depression. J Clin Psychiatry 61 Suppl 1, 5–12 (2000). [PubMed]
Brambilla P., Perez J., Barale F., Schettini G. & Soares J. C. GABAergic dysfunction in mood disorders. Mol Psychiatry 8, 721–737, 715 (2003). [PubMed]
Sanacora G. & Saricicek A. GABAergic contributions to the pathophysiology of depression and the mechanism of antidepressant action. CNS Neurol Disord Drug Targets 6, 127–140 (2007). [PubMed]
Kaddurah-Daouk R. & Krishnan K. R. Metabolomics: a global biochemical approach to the study of central nervous system diseases. Neuropsychopharmacology 34, 173–186 (2009). [PubMed]
Kaddurah-Daouk R., Kristal B. S. & Weinshilboum R. M. Metabolomics: a global biochemical approach to drug response and disease. Annu Rev Pharmacol Toxicol 48, 653–683 (2008). [PubMed]
Kristal B. S., Shurubor Y. I., Kaddurah-Daouk R. & Matson W. R. Metabolomics in the study of aging and caloric restriction. Methods Mol Biol 371, 393–409 (2007). [PubMed]
Kristal B. S., Shurubor Y. I., Kaddurah-Daouk R. & Matson W. R. High-performance liquid chromatography separations coupled with coulometric electrode array detectors: a unique approach to metabolomics. Methods Mol Biol 358, 159–174 (2007). [PubMed]
Lindon J. C., Holmes E. & Nicholson J. K. Metabonomics in pharmaceutical R&D. Febs J 274, 1140–1151 (2007). [PubMed]
Ji Y. et al. Glycine and a glycine dehydrogenase (GLDC) SNP as citalopram/escitalopram response biomarkers in depression: pharmacometabolomics-informed pharmacogenomics. Clin Pharmacol Ther 89, 97–104, 10.1038/clpt.2010.250 (2011). [PMC free article] [PubMed]
Paige L. A., Mitchell M. W., Krishnan K. R., Kaddurah-Daouk R. & Steffens D. C. A preliminary metabolomic analysis of older adults with and without depression. Int J Geriatr Psychiatry 22, 418–423 (2007). [PubMed]
Steffens D. C. et al. Metabolomic differences in heart failure patients with and without major depression. J Geriatr Psychiatry Neurol 23, 138–146 (2010). [PMC free article] [PubMed]
Kaddurah-Daouk R. et al. Metabolomic mapping of atypical antipsychotic effects in schizophrenia. Mol Psychiatry 12, 934–945 (2007). [PubMed]
Yao J. K. et al. Altered interactions of tryptophan metabolites in first-episode neuroleptic-naive patients with schizophrenia. Mol Psychiatry (2009). [PMC free article] [PubMed]
Wu J., An Y., Yao J., Wang Y. & Tang H. An optimised sample preparation method for NMR-based faecal metabonomic analysis. Analyst 135, 1023–1030 (2010). [PubMed]
Johansen K. K. et al. Metabolomic profiling in LRRK2-related Parkinson's disease. PLoS One 4, e7551 (2009). [PMC free article] [PubMed]
Bogdanov M. et al. Metabolomic profiling to develop blood biomarkers for Parkinson's disease. Brain 131, 389–396 (2008). [PubMed]
Lan M. J. et al. Metabonomic analysis identifies molecular changes associated with the pathophysiology and drug treatment of bipolar disorder. Mol Psychiatry 14, 269–279 (2009). [PubMed]
Wikoff W. R., Pendyala G., Siuzdak G. & Fox H. S. Metabolomic analysis of the cerebrospinal fluid reveals changes in phospholipase expression in the CNS of SIV-infected macaques. J Clin Invest 118, 2661–2669 (2008). [PMC free article] [PubMed]
Asberg M. et al. CSF monoamine metabolites in melancholia. Acta Psychiatr Scand 69, 201–219 (1984). [PubMed]
Traskman-Bendz L., Asberg M., Bertilsson L. & Thoren P. CSF monoamine metabolites of depressed patients during illness and after recovery. Acta Psychiatr Scand 69, 333–342 (1984). [PubMed]
Mayberg H. S. Limbic-cortical dysregulation: a proposed model of depression. J Neuropsychiatry Clin Neurosci 9, 471–481 (1997). [PubMed]
Van Praag H. M. New evidence of serotonin-deficient depressions. Neuropsychobiology 3, 56–63 (1977). [PubMed]
Ferrier I. N. et al. Postmortem neurochemical studies in depression. Ann N Y Acad Sci 487, 128–142 (1986). [PubMed]
McKeith I. G. et al. 5-HT receptor binding in post-mortem brain from patients with affective disorder. J Affect Disord 13, 67–74 (1987). [PubMed]
van Praag H. M. & de Haan S. Central serotonin deficiency--a factor which increases depression vulnerability? Acta Psychiatr Scand Suppl 280, 89–96 (1980). [PubMed]
Vestergaard P. et al. Biogenic amine metabolites in cerebrospinal fluid of patients with affective disorders. Acta Psychiatr Scand 58, 88–96 (1978). [PubMed]
Lee D. S. et al. The implications of human metabolic network topology for disease comorbidity. Proc Natl Acad Sci U S A 105, 9880–9885 (2008). [PubMed]
Nordin C. Relationships between clinical symptoms and monoamine metabolite concentrations in biochemically defined subgroups of depressed patients. Acta Psychiatr Scand 78, 720–729 (1988). [PubMed]
Inoue T., Tsuchiya K. & Koyama T. Regional changes in dopamine and serotonin activation with various intensity of physical and psychological stress in the rat brain. Pharmacol Biochem Behav 49, 911–920 (1994). [PubMed]
Kennedy N. & Paykel E. S. Residual symptoms at remission from depression: impact on long-term outcome. J Affect Disord 80, 135–144 (2004). [PubMed]
Santesso D. L. et al. Enhanced negative feedback responses in remitted depression. Neuroreport 19, 1045–1048 (2008). [PMC free article] [PubMed]
Hasler G. et al. Neural response to catecholamine depletion in unmedicated subjects with major depressive disorder in remission and healthy subjects. Arch Gen Psychiatry 65, 521–531 (2008). [PMC free article] [PubMed]
Hasler G. et al. Reward processing after catecholamine depletion in unmedicated, remitted subjects with major depressive disorder. Biol Psychiatry 66, 201–205 (2009). [PMC free article] [PubMed]
Neumeister A., Young T. & Stastny J. Implications of genetic research on the role of the serotonin in depression: emphasis on the serotonin type 1A receptor and the serotonin transporter. Psychopharmacology (Berl) 174, 512–524 (2004). [PubMed]
Delgado P. L. Depression: the case for a monoamine deficiency. J Clin Psychiatry 61 Suppl 6, 7–11 (2000). [PubMed]
Moreno F. A. et al. CSF neurochemicals during tryptophan depletion in individuals with remitted depression and healthy controls. Eur Neuropsychopharmacol 20, 18–24 (2010). [PMC free article] [PubMed]
Carpenter L. L., Anderson G. M., Siniscalchi J. M., Chappell P. B. & Price L. H. Acute changes in cerebrospinal fluid 5-HIAA following oral paroxetine challenge in healthy humans. Neuropsychopharmacology 28, 339–347 (2003). [PubMed]
Lundmark J., Walinder J., Alling C., Manniche P. M. & Dalgaard L. The effect of paroxetine on cerebrospinal fluid concentrations of neurotransmitter metabolites in depressed patients. Eur Neuropsychopharmacol 4, 1–6 (1994). [PubMed]
Malone K. M. et al. Fenfluramine challenge test as a predictor of outcome in major depression. Psychopharmacol Bull 29, 155–161 (1993). [PubMed]
Martensson B. et al. Effects of clomipramine treatment on cerebrospinal fluid monoamine metabolites and platelet 3H-imipramine binding and serotonin uptake and concentration in major depressive disorder. Acta Psychiatr Scand 83, 125–133 (1991). [PubMed]
Sheline Y., Bardgett M. E. & Csernansky J. G. Correlated reductions in cerebrospinal fluid 5-HIAA and MHPG concentrations after treatment with selective serotonin reuptake inhibitors. J Clin Psychopharmacol 17, 11–14 (1997). [PubMed]
Backman J., Alling C., Alsen M., Regnell G. & Traskman-Bendz L. Changes of cerebrospinal fluid monoamine metabolites during long-term antidepressant treatment. Eur Neuropsychopharmacol 10, 341–349 (2000). [PubMed]
Bhagwagar Z. & Cowen P. J. ‘It's not over when it's over’: persistent neurobiological abnormalities in recovered depressed patients. Psychol Med 38, 307–313 (2008). [PubMed]
Stockmeier C. A. Involvement of serotonin in depression: evidence from postmortem and imaging studies of serotonin receptors and the serotonin transporter. J Psychiatr Res 37, 357–373 (2003). [PubMed]
Melfi C. A. et al. The effects of adherence to antidepressant treatment guidelines on relapse and recurrence of depression. Arch Gen Psychiatry 55, 1128–1132 (1998). [PubMed]
Claxton A. J., Li Z. & McKendrick J. Selective serotonin reuptake inhibitor treatment in the UK: risk of relapse or recurrence of depression. Br J Psychiatry 177, 163–168 (2000). [PubMed]
van Praag H. & de Hann S. Depression vulnerability and 5-hydroxytryptophan prophylaxis. Psychiatry Res 3, 75–83 (1980). [PubMed]
Drexhage RC K. E., Padmos RC,. Heul-Nieuwenhuijzen L. & Beumer W. V. M. Drexhage HA The mononuclear phagocyte system and its cytokine inflammatory networks in schizophrenia and bipolar disorder. Expert Rev Neurother 10. [PubMed]
Wichers M. C. et al. IDO and interferon-alpha-induced depressive symptoms: a shift in hypothesis from tryptophan depletion to neurotoxicity. Mol Psychiatry 10, 538–544 (2005). [PubMed]
Capuron L. et al. Interferon-alpha-induced changes in tryptophan metabolism. relationship to depression and paroxetine treatment. Biol Psychiatry 54, 906–914 (2003). [PubMed]
Kling M. A. et al. Sustained low-grade pro-inflammatory state in unmedicated, remitted women with major depressive disorder as evidenced by elevated serum levels of the acute phase proteins C-reactive protein and serum amyloid A. Biol Psychiatry 62, 309–313 (2007). [PMC free article] [PubMed]
Hasler G., Drevets W. C., Manji H. K. & Charney D. S. Discovering endophenotypes for major depression. Neuropsychopharmacology 29, 1765–1781 (2004). [PubMed]
Roy A. et al. Reduced CSF concentrations of homovanillic acid and homovanillic acid to 5-hydroxyindoleacetic acid ratios in depressed patients: relationship to suicidal behavior and dexamethasone nonsuppression. Am J Psychiatry 143, 1539–1545 (1986). [PubMed]
Cools R., Nakamura K. & Daw N. D. Serotonin and Dopamine: Unifying Affective, Activational, and Decision Functions. Neuropsychopharmacology, | doi:10.1038/npp.2010.121 (2010). [PMC free article] [PubMed]
McGowan P. O. & Kato T. Epigenetics in mood disorders. Environ Health Prev Med 13, 16–24 (2008). [PMC free article] [PubMed]
Cantoni G. L., Mudd S. H. & Andreoli V. Affective disorders and S-adenosylmethionine: a new hypothesis. Trends Neurosci 12, 319–324 (1989). [PubMed]
Kagan B. L., Sultzer D. L., Rosenlicht N. & Gerner R. H. Oral S-adenosylmethionine in depression: a randomized, double-blind, placebo-controlled trial. Am J Psychiatry 147, 591–595 (1990). [PubMed]
Abdolmaleky H. M. et al. Hypomethylation of MB-COMT promoter is a major risk factor for schizophrenia and bipolar disorder. Hum Mol Genet 15, 3132–3145 (2006). [PMC free article] [PubMed]
Abdolmaleky H. M., Smith C. L., Zhou J. R. & Thiagalingam S. Epigenetic alterations of the dopaminergic system in major psychiatric disorders. Methods Mol Biol 448, 187–212 (2008). [PubMed]
Dempster E. L., Mill J., Craig I. W. & Collier D. A. The quantification of COMT mRNA in post mortem cerebellum tissue: diagnosis, genotype, methylation and expression. BMC Med Genet 7, 10 (2006). [PMC free article] [PubMed]
Krishnan V. et al. Molecular adaptations underlying susceptibility and resistance to social defeat in brain reward regions. Cell 131, 391–404 (2007). [PubMed]
Cui J., Shao L., Young L. T. & Wang J. F. Role of glutathione in neuroprotective effects of mood stabilizing drugs lithium and valproate. Neuroscience 144, 1447–1453 (2007). [PubMed]
Frey B. N. et al. Increased oxidative stress in submitochondrial particles after chronic amphetamine exposure. Brain Res 1097, 224–229 (2006). [PubMed]
Ozcan M. E., Gulec M., Ozerol E., Polat R. & Akyol O. Antioxidant enzyme activities and oxidative stress in affective disorders. Int Clin Psychopharmacol 19, 89–95 (2004). [PubMed]
Tsaluchidu S., Cocchi M., Tonello L. & Puri B. K. Fatty acids and oxidative stress in psychiatric disorders. BMC Psychiatry 8 Suppl 1, S5 (2008). [PMC free article] [PubMed]
Gawryluk J. W., Wang J. F., Andreazza A. C., Shao L. & Young L. T. Decreased levels of glutathione, the major brain antioxidant, in post-mortem prefrontal cortex from patients with psychiatric disorders. Int J Neuropsychopharmacol, 1–8 (2011). [PubMed]
Agren H., Niklasson F. & Hallgren R. Brain purinergic activity linked with depressive symptomatology: hypoxanthine and xanthine in CSF of patients with major depressive disorders. Psychiatry Res 9, 179–189 (1983). [PubMed]
Niklasson F., Agren H. & Hallgren R. Purine and monoamine metabolites in cerebrospinal fluid: parallel purinergic and monoaminergic activation in depressive illness? J Neurol Neurosurg Psychiatry 46, 255–260 (1983). [PMC free article] [PubMed]
LeWitt P. A. et al. Markers of dopamine metabolism in Parkinson's disease. The Parkinson Study Group. Neurology 42, 2111–2117 (1992). [PubMed]
American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders, Fourth edition. Washington, DC American Psychiatric Association; 1994.
First M. B., Spitzer R. L., Gibbon M. & William J. B. W. Structured Clinical Interview fpr DSM-IV-TR Axis 1 disorders, Research Version, Patient Edition (SCID-I/P) New York, NY: New York State Psychiatric Institute, Biometric research; 2002. 73 Maxwell, M. (Bethesda, MD, 1992).
Maxwell E. Family Interview for Genetic Studies (FIGS): A Manual For FIGS, Bethesda MD: Clinical Neurogenetics Branch, Intrmural Research Program National Institute of Mental Health; 1992.
Rozen S. et al. Metabolomic analysis and signatures in motor neuron disease. Metabolomics 1, 101–108 (2005). [PMC free article] [PubMed]
Kristal B. S., Vigneau-Callahan K. E. & Matson W. R. Simultaneous analysis of the majority of low-molecular-weight, redox-active compounds from mitochondria. Anal Biochem 263, 18–25 (1998). [PubMed]
Shi W. & Palmer C. P. Effect of pendent group structures on the chemical selectivity and performance of sulfonated copolymers as novel pseudophases in electrokinetic chromatography. Electrophoresis 23, 1285–1295 (2002). [PubMed]
Vigneau-Callahan K. E., Shestopalov A. I., Milbury P. E., Matson W. R. & Kristal B. S. Characterization of diet-dependent metabolic serotypes: analytical and biological variability issues in rats. J Nutr 131, 924S–932S (2001). [PubMed]
Lehmann E. L. Nonparametrics: Statistical Methods Based on Ranks. Oakland, Springer Science & Business Media, LLC., New York, NY. 5, 202–249 (2006).
Barrett T. Computations using analysis of covariance. Computational Statistics 3, 260–268 (2011).

Articles from Scientific Reports are provided here courtesy of
Nature Publishing Group