Biomarkers of Inflammation, Immunosuppression and Stress Are Revealed by Metabolomic Profiling of Tuberculosis Patients

Abstract

Although tuberculosis (TB) causes more deaths than any other pathogen, most infected individuals harbor the pathogen without signs of disease. We explored the metabolome of >400 small molecules in serum of uninfected individuals, latently infected healthy individuals and patients with active TB. We identified changes in amino acid, lipid and nucleotide metabolism pathways, providing evidence for anti-inflammatory metabolomic changes in TB. Metabolic profiles indicate increased activity of indoleamine 2,3 dioxygenase 1 (IDO1), decreased phospholipase activity, increased abundance of adenosine metabolism products, as well as indicators of fibrotic lesions in active disease as compared to latent infection. Consistent with our predictions, we experimentally demonstrate TB-induced IDO1 activity. Furthermore, we demonstrate a link between metabolic profiles and cytokine signaling. Finally, we show that 20 metabolites are sufficient for robust discrimination of TB patients from healthy individuals. Our results provide specific insights into the biology of TB and pave the way for the rational development of metabolic biomarkers for TB.

About the Authors:

January Weiner 3rd

Contributed equally to this work with: January Weiner 3rd, Shreemanta K. Parida

* E-mail: [email protected] (SHEK); [email protected] (JW)

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Shreemanta K. Parida

Contributed equally to this work with: January Weiner 3rd, Shreemanta K. Parida

Current address: Vaccine Grand Challenge Program, Department of Biotechnology, Ministry of Science and Technology, Government of India, New Delhi, India

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Jeroen Maertzdorf

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Gillian F. Black

Current address: The Sustainable Livelihoods Foundation, Cape Town, South Africa

Affiliation: Department of Biomedical Sciences, University of Stellenbosch, Cape Town, South Africa

Dirk Repsilber

Affiliation: Biomathematics/Bioinformatics Group, Genetics and Biometry, Leibniz Institute for Farm Animal Biology, FBN, Dummerstorf, Germany

Anna Telaar

Affiliation: Biomathematics/Bioinformatics Group, Genetics and Biometry, Leibniz Institute for Farm Animal Biology, FBN, Dummerstorf, Germany

Robert P. Mohney

Affiliation: Metabolon, Inc., Durham, North Carolina, United States of America

Cordelia Arndt-Sullivan

Current address: Department of Biology, Humboldt University, Berlin, Germany

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Christian A. Ganoza

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Kellen C. Faé

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Gerhard Walzl

Affiliation: Department of Biomedical Sciences, University of Stellenbosch, Cape Town, South Africa

Stefan H. E. Kaufmann

* E-mail: [email protected] (SHEK); [email protected] (JW)

Affiliation: Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany

Introduction

Metabolomics provides unprecedented insights into the biology of an organism in the state of disease. In particular, metabolites can be viewed as a close recapitulation of the disease phenotype, closer to ongoing pathogenesis in an organism than changes in gene expression [1]. Consequently, metabolomics has deepened our understanding of biological mechanisms involved in several noninfectious diseases and provided a platform for the identification of new biomarkers [2]. Metabolic signatures have proven their value in several diseases, such as Alzheimer’s disease [3], Parkinson’s disease [4], myocardial ischemia [5], hypertension [6], cancer [7] and diabetes [8]–[11]. In contrast, fewer studies have specifically addressed the metabolomic alterations that occur in infectious diseases [12]–[16].

Tuberculosis (TB) is caused by the intracellular bacterial pathogen Mycobacterium tuberculosis [17], [18]. It is characterized by chronic latency without clinical symptoms in the vast majority of infected individuals. Of the more than 2 billion people infected with M. tuberculosis globally, less than one-tenth is likely to develop active TB during their lifetime and in the majority of cases infection remains in an asymptomatic stage of latency. Even though TB can be cured, drug treatment is long-lasting and complex, and requires at least three drugs given over a minimum of 6 months. In many regions, notably in developing countries with high TB incidences, diagnosis is neither sensitive nor specific, and an estimated 40% of TB patients fail to be correctly diagnosed [19]. The tuberculin skin test (TST) diagnoses immunologic sensitization to mycobacteria, indicating exposure, and hence cannot be used to distinguish infected healthy individuals from patients with active TB. More recently, gene expression profiling has shown promise in differentiating between infection states in TB [20]–[22].

We investigated the feasibility of identifying small molecule biochemical profiles in serum for gaining novel biological insights into the mechanisms underlying TB. The aim of this study is to use metabolites as biomarkers of TB infection and disease status, and to elucidate their role in protection and pathogenesis in TB. We acquired information on 428 distinct small molecules, including amino acids, short peptides, carbohydrates, fatty acids, nucleotides and cofactors, and compared the characteristic metabolomic profiles of three groups: (i) healthy M. tuberculosis-uninfected controls (TST^–), (ii) latently M. tuberculosis-infected healthy individuals (TST⁺), and (iii) patients with active TB (TB^active; see Table S1).

We defined several remarkable differences between the metabolic profiles of healthy subjects and TB patients, including lower relative abundances of amino acids, medium-chain fatty acids and lysophosphatidylcholines (LPCs), higher relative abundances of fibrinopeptides and adenosine degradation products inosine, hypoxanthine and ribose, as well other compounds, such as bile acids and uremic toxins. The observed differences and correlations between compound abundances allowed us to construct functional clusters and infer their biological significance in TB disease progression. In addition, we could identify a serum metabolomic biomarker set that distinguishes active TB from the two groups of healthy subjects, which is a necessary first step towards the use of metabolic biomarkers in TB diagnosis. We also found that distinct cytokine abundances in several instances correlated with the relative abundances of metabolites, indicating a functional link between these systems. Finally, we studied the role of tryptophan degradation by indoleamine 2,3 dioxygenase 1 (IDO1), a well-known immunosuppressive mechanism. We revealed that IDO1 expression was upregulated in pulmonary lesions from mice suffering from experimental TB. Moreover, in vitro infection of macrophages and dendritic cells (DCs) with M. tuberculosis mirrored findings observed in serum.

Results

Test Groups Differ in the Relative Abundance of Several Metabolites

Differences in the abundance of small molecules in serum between the three study groups were determined by applying a t-test for each of the three possible comparisons (TST^– vs. TST⁺, TST^– vs. TB^active and TST⁺ vs. TB^active), corrected for multiple testing at a significance threshold of 0.05. Results revealed significant differences for 176 compounds between the TB^active group and the two healthy groups (TST^– and TST⁺; see Table S2 and Figure 1).

[Figure omitted. See PDF.]

Figure 1. Examples of metabolite patterns in tuberculosis patients (TB^active), healthy uninfected (TST^–) and latently infected (TST⁺) individuals.

(A) Changes in relative abundance of three exemplary small metabolites: fibrinopeptide A, inosine and kynurenine. Red line indicates sample mean, blue line indicates sample median. Stars indicate significant differences between profiles (result from t-test corrected for multiple testing; *, p<0.05; **, p<0.01; ***, p<0.001) (B) Correlation between abundance of inosine and abundance of hypoxanthine and choline. Colors and symbols denote study groups: grey squares, TST^–; green circles, TST⁺; red triangles, TB^active. Spearman correlation coefficients (rho) and corresponding p-values are given.

https://doi.org/10.1371/journal.pone.0040221.g001

Among various other compounds, we found that several amino acids, such as histidine, cysteine, glutamine, tryptophan, citrulline and creatine, were at lower levels in the TB^active group compared to the two control groups, and only some amino acids were increased in abundance (including kynurenine, phenylalanine and pyroglutamine) in the TB^active group. Other identified compounds with markedly differentiated abundance between TST⁺ and TB^active groups were sialic acid (N-acetylneuraminate), 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF), inosine, xanthine, hypoxanthine, fibrinopeptide A, glucose, and mannose (Figure 2 and Table S2). Intriguingly, we also identified six metabolites that differed significantly between the two healthy control groups, TST^– and TST⁺ (inosine, hypoxanthine, glycylvaline, 5-oxoproline and two unidentified compounds; see Table S2, A and Discussion).

[Figure omitted. See PDF.]

Figure 2. Heat map showing fold changes of small metabolic compounds in the three study groups, TB patients, healthy uninfected and latently infected individuals.

Fold changes are relative to the average abundance in the TST^– group. Red indicates relative abundance higher than average in the TST^– group; blue indicates relative abundance lower than average in the TST^– group. Horizontal axis: samples belonging to the three study groups; vertical axis: top 50 compounds selected by variable importance in RF analysis, including compounds that could not be identified, but were strong predictors of sample status. Color bars above the heat map denote study groups: grey, TST^–; green, TST⁺; red, TB^active. See also Figure S2.

https://doi.org/10.1371/journal.pone.0040221.g002

Differences in Small Metabolites can be Used as Specific and Sensitive Biosignatures

To determine whether differences in metabolite abundance between any two of the three study groups could be efficiently exploited for building a sensitive biosignature of TB status, we applied random forests (RF) classification, a supervised machine learning algorithm [23]. The TB^active group could be sensitively and specifically distinguished from the other two groups, TST^– and TST⁺ with a cross-validation error rate of 3.3%. Similar results were obtained for classification between the TB^active group and each of the healthy groups separately. Classification between TST^– and TST⁺ groups was error prone (i.e., overall cross-validation error rate of 41.3% and 32.6% for the two groups, respectively), due to highly similar metabolite profiles in these two groups of healthy individuals. Additional analyses revealed that less than 20 compounds sufficed for robust classification of TB^active with a cross-validation error rate of approximately 3% (Figures S1 and S2; Tables 1, 2 and S3). For illustration purposes, 10 samples from the TST⁺ and TB^active groups were randomly assigned to a test set, and the remaining training set was used to repeat all calculations. As expected from error rates calculated by an external cross-validation procedure, no misclassifications occurred in the resulting data set (Figure S2).

[Figure omitted. See PDF.]

Table 1. Confusion matrices showing the results of external cross-validation for random forests models of the discrimination procedure among the tuberculosis patients (TB^active) and healthy uninfected (TST^–) and latently infected (TST⁺) individuals using the full set of metabolites.

https://doi.org/10.1371/journal.pone.0040221.t001

[Figure omitted. See PDF.]

Table 2. The top 20 named metabolites identified in the random forests (RF) analysis as most important for discrimination between the TST⁺ and TB^active groups.

https://doi.org/10.1371/journal.pone.0040221.t002

Metabolomic signatures that correctly classified TB^active patients comprised distinct clusters, which provided deeper insights into pathogenic mechanisms of TB. We applied hierarchical clustering with correlations between profiles as a distance measure, and bootstrapping for statistical assessment of the quality of clustering. This method created a tree of metabolites, and at each tree node a p-value was obtained by bootstrapping. Consequently, statistical significance could be evaluated for each potential cluster. Using a p-value threshold of 0.05, we found 18 clusters with at least five compounds. Among these, 12 clusters were of biological interest as they contained at least one compound showing significant differences between the three study groups (Table S2). Similar results were obtained using a substantially different method, sparse PLS discriminant analysis (SPLS-DA) [24].

Biological Significance of Discriminative Clusters

To understand the relationships between and within clusters, as well as to visually explore the underlying biological mechanisms, we constructed a network based on partial correlations between compounds. Group effects were removed such that compounds showing similar differences between groups, but no within-group correlation that would indicate a functional relationship, were not clustered together. For each compound, we fitted a linear model using the three study groups as a predictor variable and then used the residuals from this model to calculate the correlations between compounds. Based on these correlations, we constructed a network graph that can be explored visually (Figures 3 and S3).

[Figure omitted. See PDF.]

Figure 3. Network showing functional relationships between the small metabolic compounds in TB patients, healthy uninfected and latently infected individuals.

Nodes correspond to metabolites; edges correspond to statistically significant correlation between residual small metabolite profiles corrected for study classes. Colors correspond to differences between the TST⁺ and classes (see Figure S3 for additional comparisons). Color intensity indicates significance of difference with darker colors corresponding to more significant differences. Metabolites with adjusted p value >0.05 are not colored. Line thickness corresponds to the absolute Spearman correlation coefficients corrected for groups (see Methods). See also Figure S3.

https://doi.org/10.1371/journal.pone.0040221.g003

In a first step, we validated our approach to identify functionally related compounds. We predicted that tightly coregulated compounds (e.g., a class of fatty acids) would form significant clusters. Consistent with this prediction, numerous groups of related compounds were reliably identified and grouped in tight clusters. For example, all medium-chain fatty acids (MCFAs) formed a uniform group of correlated biochemical compound profiles. Moreover, correct identification of functionally related compounds should allow close correlation of xenobiotics with their metabolites. Consistent with this prediction, degradation products of xenobiotics, including caffeine, stachydrine, salicylate and acetaminophen (Tylenol), were tightly correlated (Figures 3 and S3).

We labeled the clusters by the predominant pathway or type of metabolic compound found within the particular cluster (Table S4). Several of the clusters we identified included metabolic compounds that significantly differed between the study groups. Notably, we determined alterations in clusters of the following: amino acids, medium-chain fatty acids, carnitines, LPCs and fibrinopeptides. We did not find significant differences in the clusters of hormones, xenobiotics or long chain fatty acids.

Abundance of Metabolites is Linked to Immune Response

To determine whether there is a direct functional connection between serum abundance of metabolites and of immune mediators, we analyzed 35 cytokines/chemokines in 99 out of the 136 serum samples (note: sample number was reduced due to limited availability of aliquots). Significantly elevated serum concentrations of granulocyte colony-stimulating factor (G-CSF), interferon-gamma (IFN-γ), interleukin 6 (IL-6), C-X-C motif chemokine 10 (CXCL-10), soluble alpha chain of the IL-2 receptor alpha (sIL2ra) and vascular endothelial growth factor (VEGF), as well as lower concentrations of macrophage-derived cytokine (MDC) were detected in TB^active patients as compared to the two healthy groups (Figure 4). Notably, several of the differentially abundant cytokines in the study anticorrelated significantly with serum metabolites. For example, IL-6, IP-10 and sIL2ra were significantly correlated with the amino acids tryptophan and glutamine (p<10⁻⁵, r < –0.3). In total, out of 15,408 correlations between cytokines/chemokines and metabolic compounds tested, 211 significant correlations (p<0.05) remained after correction for multiple testing. We wanted to determine categories of metabolic compounds that frequently correlated with the abundances of cytokines and chemokines. To this end, we performed an enrichment analysis, which revealed significant enrichment of amino acids (p<10^–6) and carbohydrates (p = 0.004). Therefore, this establishes a link between changes in the metabolic profile and the cytokine response in TB.

[Figure omitted. See PDF.]

Figure 4. Abundance of cytokines and their correlation with selected metabolites in TB patients, healthy uninfected and latently infected individuals.

(A) Strip charts showing abundances of eight cytokines that differed significantly between the study groups. Significance thresholds for a two-tailed t-test corrected for multiple testing: *, p<0.05; **, p<0.01; ***, p<0.001. Blue line indicates group median, red line indicates group mean. (B) Correlations between metabolic compounds and cytokines. Grey squares, healthy controls (TST^–); green circles, latently infected individuals (TST⁺); brown triangles, active TB patients (TB^active). Spearman correlation coefficient (rho) and p-values are given.

https://doi.org/10.1371/journal.pone.0040221.g004

Increased IDO1 Activity and Kynurenine Production

Sera of TB patients were characterized by a significant increase in the abundance of kynurenine, a product of tryptophan metabolism. The enzyme IDO1 catalyzes the degradation of tryptophan to kynurenine, with well-established immunosuppressive effects in mammalian pregnancy [25], [26], tumor resistance [27], [28], chronic infection [29]–[33] and autoimmune diseases [34]. Furthermore, IDO1 is critically involved in CD4 and CD8 effector T cell suppression, as well as in generation and activation of regulatory T cells [35], [36]. Therefore, we decided to investigate the role of IDO and its metabolites in immunity against TB.

Consistent with elevated tryptophan degradation to kynurenine in TB, IDO1 protein was highly expressed in pulmonary lesions in mice suffering from experimental TB (Figure 5A). Human monocyte-derived DCs and macrophages strongly upregulated the expression of IDO1 after infection with M. tuberculosis; live bacteria induced stronger expression compared to heat-killed or irradiated bacteria (Figure 5B and C). In agreement with this finding, the enzymatic activity of IDO1, measured by the conversion of tryptophan to kynurenines by high-performance liquid chromatography (HPLC), was increased in infected host cells (Figure 5D and E). T cells stimulated with DCs pulsed with purified protein derivative (PPD) plus mycobacterial lipid (mannosylated lipoarabinomannan, ManLAM) and treated with the specific IDO1 inhibitor 1-methyl-DL-tryptophan (1-MT-DL) showed elevated proliferation compared to cells treated only with PPD and ManLAM. However, addition of kynurenine reversed this effect (Figure 5F), emphasizing the regulation of T cells by kynurenines.

[Figure omitted. See PDF.]

Figure 5. Demonstration of IDO1 expression and kynurenine production in response to M. tuberculosis infection and regulation of M. tuberculosis-specific T-cell responses by kynurenines.

(A) Immunohistochemistry staining of formalin-fixed, paraffin-embedded tissue of a murine pulmonary TB lesion stained with anti-IDO polyclonal antibody; staining representative of lesions from five animals. Bar is equal to 200 nm. Human monocyte-derived dendritic cells (DCs) (B) and macrophages (C) were infected with M. tuberculosis H37Rv or stimulated with irradiated and heat-killed M. tuberculosis H37Rv for 24 h and indoleamine 2,3 dioxygenase 1 (IDO1) gene expression was measured by qPCR. Mean and standard deviation (SD) of fold-change IDO1 gene expression of one donor representative of four. Line indicates minimal significant fold change threshold equal to 1.5. DCs (D) and macrophages (E) were infected with M. tuberculosis H37Rv and cell culture supernatants were collected at indicated times for measurement of kynurenines by HPLC; kynurenine levels from uninfected controls were subtracted. Means and SD of four donors are depicted (C and D). Star indicates significance (p<0.05) in Friedman test.(F) Human monocyte-derived DCs were pulsed with purified protein derivative (PPD) and mannosylated lipoarabinomannan (ManLAM) and co-cultured for 4 days with autologous CFSE-labeled purified T cells (DC:T cell ratio 1∶20) in the presence or absence of 1-methyl-DL-tryptophan (1-MT-DL) and 3-OH-kynurenine (Kyn). Cell proliferation was assessed by CFSE dilution using flowcytometry. Figure representative of three independent experiments (ANOVA F = 4.1 for CD⁺CD4⁺ cells and F = 3.3 for CD3⁺CD8⁺ cells, p<0.05).

https://doi.org/10.1371/journal.pone.0040221.g005

Discussion

Our analysis of several hundred serum small molecules has provided intriguing insights into the metabolomic status of three groups in M. tuberculosis infection: no M. tuberculosis infection (TST^–); M. tuberculosis latent infection (TST⁺); and active TB disease (TST^active). These metabolomic characteristics could be translated into a custom-made biosignature capable of discriminating healthy individuals with and without M. tuberculosis infection from TB patients.

In two recent studies, tissues from mice infected with M. tuberculosis [16] and lung granulomas from M. tuberculosis-infected guinea pigs [15] were profiled for metabolites using nuclear magnetic resonance (NMR). Although the number of metabolites profiled in these studies was significantly lower, several findings of these studies are congruent with our findings.

We applied a clustering procedure for groups of compounds that were functionally related. In each group, the compounds were significantly correlated, indicating their codependence on a specific metabolic process or systemic response. Consistent with this notion, several compounds expected to be present in a given group were significantly correlated. This procedure allowed us to link additional compounds that were not predicted to form a cluster. Notably, in many cases, differential abundance of serum metabolites between the different study groups could not be directly attributed to specific metabolic pathways. The presence of metabolic compounds in serum is rarely a simple consequence of metabolic processes within cells, but rather the result of a complex interplay between different parameters, such as unique metabolic and secretory features of cells in different organs. In the following, we discuss the most prominent metabolite clusters in the context of their biological relevance and their potential as custom-tailored biosignatures for TB.

Tryptophan and Kynurenine

In the TB^active group, but not in the TST^−/+ groups, kynurenine and its derivative, 3-hydroxykynurenine were elevated. The activity of IDO1, an enzyme that converts tryptophan to kynurenine, is increased under various inflammatory conditions resulting in heightened abundance of kynurenines [37]. Therefore, increased serum concentrations of kynurenines suggest upregulation or increased IDO1 activity in the TB^active group.

We tested the hypothesis that increased IDO1 activity is a direct consequence of stimulation of immune cells by M. tuberculosis and that this increased enzymatic activity leads to an increased level of kynurenines. We identified IDO1 in situ in granulomatous lesions obtained from the lungs of mice suffering from experimental TB. Furthermore, we observed elevated levels of IDO gene expression and enzymatic activity in human DCs and macrophages upon infection with M. tuberculosis (Figure 5B–E). Finally, we demonstrate human T cell regulation by kynurenines in TB (Figure 5F). These findings support the immune regulatory function of the active catabolites generated by the increased degradation of tryptophan during TB, suggesting that the tryptophan degradation pathway controlled by IDO1 is involved in the pathogenesis of this disease.

Increased abundance of kynurenine in serum from the TB^active group was significantly correlated with similarly increased abundance of the immunosuppressive stress hormone cortisol (p<10^–5, r = 0.38). Glucocorticoids can activate IDO1 via glucocorticoid-induced tumor necrosis factor receptor (GITR) and protect against allergic bronchopulmonary aspergillosis [38]. Hence, the tryptophan catabolism pathway could be a mechanism by which corticoids exert immunosuppression [38]. Independently, in a recent publication, Suzuki et al. [39] found that serum IDO activity and kynurenine abundance is higher in TB patients.

Fibrinopeptide A

Three forms of fibrinopeptide A were significantly correlated and elevated in both TST⁺ and TB^active groups with highest abundance in TB patients. Tuberculous granulomas are walled off from the surrounding tissue by a fibrinous wall. Fibrin deposition in granulomas resulting in release of fibrinopeptides could explain the high serum abundance of fibrinopeptides in TB patients and latently infected individuals. Interestingly, in female members of the study groups, another short peptide corresponding to a fragment of the C3 complement component was correlated with fibrinopeptide A and found at significantly higher serum concentrations in TB^active than in TST^– and TST⁺ groups.

LPCs

One of the most prominent clusters in our analysis consisted of different LPCs, which were found at lower levels in the TB^active group but did not significantly differ between TST^– and TST⁺ groups. A major source of LPCs is the proinflammatory phospholipase A2, as plasma concentrations of LPCs significantly correlate with the activity of this enzyme [40]. Lower levels of LPCs have been observed in some types of cancer [41], [42], as well as in sepsis patients in whom lower LPC concentrations correlate with mortality [43]. Lower abundance of LPCs in TB patients could be mechanistically related to the finding that M. tuberculosis can induce macrophage apoptosis by inhibition of phospholipase A2 [44]. In contrast, Somashekar et al. [15] found a signature for lipolysis in the lung granulomas from M. tuberculosis-infected guinea pigs.

Hypoxia Cluster

Inosine, hypoxanthine, ribose and xanthine formed one cluster, in which the compounds were present in high abundance in several individuals belonging to the TST^– and TB^active groups, and uniformly low in abundance in the TST⁺ group (Figure 1B). These compounds are part of the inosine breakdown pathway in which inosine is processed to hypoxanthine and ribose by the 5′ nucleosidase, and hypoxanthine is subsequently oxidized to xanthine by xanthine oxidase. Inosine, hypoxanthine and ribose have been defined as biomarkers of hypoxia, hypoxemia and ischemic brain injury [45], [46]. Inosine is a natural analog of adenosine and binds adenosine receptors A2A and A3 [47], [48], thus decelerating inflammation.

In tuberculous granulomas, hypoxia prevails and is thought to contribute to the containment of M. tuberculosis. An abundance of inosine catabolites in serum from the TB^active group is consistent with this notion. Reasons for higher abundance of hypoxia biomarkers in serum of the TST^– group as compared to the TST⁺ group remain to be clarified. A high prevalence of current and previous smoking was previously reported in the study community [49]. Uniformly low serum concentrations of breakdown products of inosine in the TST⁺ group argues against a major role of hypoxia in latent infection. Rather, since individuals with higher abundance of inosine are absent from the TST⁺ group, but present in the TST^– group, it is tempting to speculate that the occurrence of hypoxia is coincident with progression from latent infection to active TB disease. Notably, elevated serum levels of hypoxanthine have been also reported in severe childhood pneumonia [12]. Also, a hypoxia signature has been identified in lung granulomas of guinea pigs infected with M. tuberculosis [15] as well as in the organs of M. tuberculosis-infected mice [16].

Several other metabolites were correlated within the inosine pathway cluster. Methionine was present in serum at low abundance in the TB^active group. The biochemical derivative of methionine, N-acetyl-L-methionine (NAM), was detected in more than one-third of both of the TB^active and TST^– groups, but was undetectable in the TST⁺ group. The levels of methionine and NAM were inversely correlated. Although the physiological role of NAM during infection remains unclear [50], a relevant enzyme, methionine adenosyltransferase, is known to be regulated by hypoxia [51]. Choline abundance correlates significantly with abundance of inosine and xanthine, and is known to be a reliable biomarker in ischemia and acute coronary syndrome, released in the course of ischemic membrane damage and the activity of phopholipase D [37]. Lactate was found at higher levels in the TB^active group, coinciding with a previous finding in mice and possibly related to increased glycolysis in granulomatous inflammation [16].

MCFAs

Eight to twelve carbon-containing MCFAs are absorbed into the bloodstream and primarily metabolized in the liver. MCFA serum concentrations were significantly reduced in the TB^active group. Medium-chain triglycerides inhibit free radical formation and tumor necrosis factor-alpha (TNF-α) production [52]. Hence, our findings point to a role of MCFAs in protection against collateral damage caused by free radicals and TNF, which are essential elements of protective immunity against M. tuberculosis.

Complement System Peptide

The peptide C3f is a breakdown product of the C3 component of the complement system, which is generated by the combined actions of factors I and H [53]. C3, a central molecule of the complement cascade, triggers a variety of important immunologic sequellae [54]. Elevated levels of C3f have been correlated with coronary disease and vascular events in women [55]. In our study, the peptide C3f was elevated in female members of the TB^active group.

Uremic Toxins

Significantly reduced serum concentrations in the TB^active group of CMPF, 3-indoxyl sulfate (IS), hippurate and indole acetate, suggest uremic cytotoxic activity. Although these compounds have not been analyzed in the context of TB, uremic toxins have been related to vitamin D metabolism [56]. Vitamin D is involved in macrophage activation and has been described as critical for host defense against TB [57]. Reduced serum concentrations in TB patients could be explained by accumulation of respective metabolites in granulomas or by the assumption that vitamin D deficiency is directly related to elevated susceptibility in TB [58], [59]. Moreover, uremic toxins can modulate apoptosis and promote neutrophil clearance [60].

Stachydrine

Stachydrine and putative demethylated stachydrine (X-11513) were correlated in serum abundance. Their average abundance was only insignificantly elevated in the TB^active group. In-depth analysis of the data revealed marked differences amongst participants in this group with eight outliers showing highly elevated serum concentrations of these compounds, compared with only three outliers in the TST⁺ group and four outliers in the TST^– group.

Stachydrine and homostachydrine are plant alkaloids found, for example, in alfalfa (lucern, Medicago sativa). Stachydrine is also found in Capparis tomentosa (Woolly caper-bush), which is used as a traditional remedy in South Africa [61] for treatment of cough and chest pain. These findings emphasize the broad scope of metabolome analyses, which not only encompasses pathogen, host and microbiome metabolism, but also environmental factors such as drug intake.

Cytokines Linked to Metabolic Changes

The serum levels of several cytokines that were significantly upregulated in samples from TB patients were significantly anticorrelated with the relative serum concentrations of a number of amino acids (see Figure 4). The relative abundance of these amino acids was significantly reduced in the TB^active group compared with TST⁺ and TST^– groups, whereas cytokines were found at elevated concentrations. However, also within the TB^active group these two classes of molecules were negatively correlated. This indicates that in patients with the strongest inflammation, amino acid abundance was diminished and IL-6 and other proinflammatory cytokines were elevated (see Figure 4).

Strongest negative correlations between amino acids and cytokines were observed for the amino acid glutamine. Glutamine is the most abundant amino acid in the body, but can become limiting under stress conditions [62], and consequently impair immune functions [63]. It is tempting to speculate that the observed increase in cytokine serum concentrations accompanied by decreased serum abundance of amino acids are indicators for disease progression in which the organism is gradually depleted of resources, resulting in cachexia. Indeed, the ratio glutamine/glutamate indicates a risk for loss of mass [64]. We found that the glutamine/glutamate ratio in serum of TB^active was significantly lower than in TST^– and TST⁺ groups (p<2×10^–5).

The concentrations of G-CSF were significantly anticorrelated with the relative abundances of gamma-glutamyl glutamine and gamma-glutamyl leucine (p<10^–3, r < –0.55 in the TB^active group). This could indicate an effect on gamma-glutamyl transpeptidase, known to be elevated in patients treated with G-CSF [65], [66].

Specificity of the Biosignature

Since the experimental groups tested did not include patients with diseases other than TB, we cannot ascertain the specificity of the resulting biosignature. It can be expected that several of the observed changes are the result of a general inflammatory process rather than a specific response to TB. Nonetheless, it is tempting to speculate that three clusters of compounds included a specific profile. Both the inosine cluster and fibrinopeptide cluster revealed significant differences between the two healthy groups (TST^– and TST⁺), indicating that these effects are specific. Furthermore, we demonstrated that M. tuberculosis directly stimulates kynurenine production at the site of infection. We assume that this effect includes a specific pattern because it represents a direct response to the pathogen, although it has been described in other diseases, as well. In addition, TB is linked to general wasting (cachexia), which has been generally observed in TB and AIDS and less so in other infectious diseases. Therefore, there is reason to speculate that definition of a specific metabolite biosignature of TB is feasible.

Concluding Remarks

TB is a chronic disease known to cause profound metabolic changes [67] as indicated by its various designations as wasting disease, consumption or phthisis in pre-antibiotic times. Metabolomics not only allows robust differentiation of patients with active TB from healthy individuals (be they M. tuberculosis-infected or not), it also reveals remarkable differences in the metabolite profile between individual TB cases (see Figure 1A). In contrast, healthy individuals from both TST^– and TST⁺ groups were more homogeneous for the most discriminative metabolites.

Our study aimed at a better understanding of the biological processes operative during infection and disease in TB, and we identified several biological mechanisms that apparently play a major role. In addition, our study provides a first step towards the development of metabolomic-based diagnosis of TB. Moreover, understanding host metabolic processes during TB could ultimately inform individualized adjunctive therapy in combination with chemotherapy, notably, during prolonged treatment of patients with multidrug-resistant or extensively drug-resistant TB. Clearly, further development of a tailored biosignature suitable for TB diagnosis will require verification in independent cohorts and comparison with other pulmonary infectious diseases. Yet, our current analyses of several hundred metabolites has led to the definition of a signature comprising fewer than 20 metabolites for reliable discrimination between active TB and latently infected or uninfected healthy individuals, emphasizing the power of metabolomic profiling in TB. It is tempting to speculate that one day, metabolic signatures may be harnessed for the development of a simple and robust dipstick test for point-of-care diagnosis of TB.

Methods

Ethics Statement

Human serum samples: the study was approved by the ethical committee of the University of Stellenbosch (Stellenbosch, South Africa) and written informed consent was obtained from all study participants. Human cell samples were obtained in accordance with the local ethical committee (Charite Ethikkommission, Berlin, Germany, EA1/200/08). Animal experiments were performed at the Max Planck Institute for Infection Biology in Berlin, Germany, and experimental protocols were approved by the State Office of Health and Social Affairs (Landesamt für Gesundheit und Soziales), Berlin, Germany.

Specimens

Serum samples collected at the University of Stellenbosch, in Stellenbosch, South Africa comprised three groups: 46 healthy controls with no clinical signs of TB infection (TST^–), 46 latently infected subjects with no clinical signs of TB, confirmed by TST (TST⁺) and 44 patients with clinical signs of pulmonary TB (TB^active; see Table S1). TB patients had the following symptoms: cough for more than 2 weeks and at least two additional symptoms from the following: hemoptysis (coughing up blood), breathing difficulty, fever, night sweats, weight loss, chest pain, fatigue. All patients had two positive (≥1+) sputum smear stains for acid-fast bacilli (Auramine stain or Ziehl-Neelsen). In 24 out of the 34 patients a Mycobacteria Growth Indicator Tube (MGIT) culture was performed, was positive and was speciated for M. tuberculosis. All subjects were HIV^– and bacille Calmette–Guérin (BCG)-vaccinated as infants, in accordance with the national vaccination program. Subjects were balanced between the groups with respect to age and gender. Clinically relevant information on patients is available upon request. Samples were stored at –80°C. At the time of sample collection, none of the TB patients was a recipient of TB treatment.

Sample Preparation

Samples were stored at –70°C until processed. Sample preparation was carried out as described previously [68] at Metabolon, Inc. Briefly, recovery standards were added prior to the first step in the extraction process for quality control purposes. To remove protein, dissociate small molecules bound to protein or trapped in the precipitated protein matrix, and to recover chemically diverse metabolites, proteins were precipitated with methanol under vigorous shaking for 2 min (Glen Mills Genogrinder 2000) followed by centrifugation. The resulting extract was divided into four fractions: one for analysis by ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS; positive mode), one for analysis by UPLC-MS/MS (negative mode), one for analysis by gas chromatography–mass spectrometry (GC-MS), and one sample was reserved for backup.

Three types of controls were analyzed in concert with the experimental samples: samples generated from a pool of human plasma (extensively characterized by Metabolon, Inc.) served as technical replicate throughout the data set; extracted water samples served as process blanks; and a cocktail of standards spiked into every analyzed sample allowed instrument performance monitoring. Instrument variability was determined by calculating the median relative standard deviation (RSD) for the standards that were added to each sample prior to injection into the mass spectrometers (median RSD = 10%; n = 30 standards). Overall process variability was determined by calculating the median RSD for all endogenous metabolites (i.e., non-instrument standards) present in 100% of the pooled human plasma samples (median RSD = 15%; n = 239 metabolites). Experimental samples and controls were randomized across the platform run.

Mass Spectrometry Analysis

Non-targeted MS analysis was performed at Metabolon, Inc. Extracts were subjected to either GC-MS [68] or UPLC-MS/MS [69]. The chromatography was standardized and once the method was validated, no further changes were made. As part of Metabolon’s general practice, all columns were purchased from a single manufacturer’s lot at the outset of experiments. All solvents were similarly purchased in bulk from a single manufacturer’s lot in sufficient quantity to complete all related experiments. For each sample, vacuum-dried samples were dissolved in injection solvent containing eight or more injection standards at fixed concentrations, depending on the platform. The internal standards were used both to assure injection and chromatographic consistency. Instruments were tuned and calibrated for mass resolution and mass accuracy daily.

The UPLC-MS/MS platform utilized a Waters Acquity UPLC and a ThermoFisher LTQ mass spectrometer, which included an electrospray ionization source and a linear ion-trap mass analyzer. The instrumentation was set to monitor for positive ions in acidic extracts or negative ions in basic extracts through independent injections. The instrument was set to scan 99–1000 m/z and alternated between MS and MS/MS scans. The scan speed was approximately six scans per s (three MS and three MS/MS scans). MS/MS scans were collected using dynamic exclusion, a process in which after an MS/MS scan of a specific m/z has been obtained, then that m/z is placed on a temporary MS/MS exclude list for a user-set period of time to allow greater MS/MS coverage of ions present in the MS scan because the instrument will not trigger an MS/MS scan of the same ion repeatedly. Extracts were loaded onto columns (Waters UPLC BEH C18-2.1×100 mm, 1.7 µm) and gradient-eluted with water and 95% methanol containing 0.1% formic acid (acidic extracts) or 6.5 mM ammonium bicarbonate (basic extracts). Columns were washed and reconditioned after every injection.

The samples destined for analysis by GC-MS were dried under vacuum desiccation for a minimum of 18 h prior to being derivatized under dried nitrogen using bistrimethyl-silyltrifluoroacetamide. Derivatized samples were separated on a 5% phenyldimethyl silicone column with helium as carrier gas and a temperature ramp from 60° to 340°C within a 17-min period. All samples were analyzed on a Thermo-Finnigan Trace DSQ MS operated at unit mass resolving power with electron impact ionization and a 50–750 atomic mass unit scan range.

Compound Identification, Quantification, and Data Curation

Metabolites were identified by automated comparison of the ion features in the experimental samples to a reference library of chemical standard entries that included retention time, molecular weight (m/z), preferred adducts, and in-source fragments as well as associated MS spectra and curated by visual inspection for quality control using software developed at Metabolon [70]. Identification of known chemical entities is based on comparison to metabolomic library entries of purified standards. Over 2,400 commercially available purified standard compounds have been acquired and registered into LIMS for distribution to both the LC/MS and GC/MS platforms for determination of their detectable characteristics. An additional 5,300 mass spectral entries have been created for structurally unnamed biochemicals, which have been identified by virtue of their recurrent nature (both chromatographic and mass spectral). These compounds have the potential to be identified by future acquisition of a matching purified standard or by classical structural analysis. Peaks were quantified using area under the curve.

Multiplex Cytokine/chemokine Analysis

Cytokine and chemokine levels in serum were measured using a premixed 42-plex MILLIPLEX MAP Human Cytokine/Chemokine kit (Millipore GmbH, Germany) according to manufacturer’s instructions. A total of 25 µl of serum was analyzed for the presence of the following analytes: epidermal growth factor (EGF), eotaxin, fibroblast growth factor (FGF)-2, Flt3-L, fractalkine, G-CSF, granulocyte macrophage colony-stimulating factor (GM-CSF), GRO, IFN-α2, IFN-γ, IL-1a, IL-1b, IL-1ra, IL-2, IL-3, IL-4, IL-5, IL-6, IL-7, IL-8, IL-9, IL-10, IL-12p40, IL-12p70, IL-13, IL-15, IL-17, interferon-inducible protein (IP)10, monocyte chemotactic protein (MCP)-1, MCP-3, MDC, macrophage inflammatory protein 1-(MIP1)a, MIP1b, platelet-derived growth factor (PDGF)-AA, PDGF-AB/BB, RANTES, sCD40L, sIL2ra, transforming growth factor-alpha (TGF)-α, TNF-α, TNF-β and VEGF. Cytokine and chemokine levels were measured on a Luminex® 100/200TM system (Luminex Corporation, Austin, TX, USA).

Experimental TB Infection and Histology

Animals (C57BL/6 mice, 8–12 weeks of age) were infected using a Glas-Col inhalation exposure system with M. tuberculosis H37Rv (∼250 colony forming units per mouse). The initial challenge dose was verified 24 h p.i. by plating complete lung homogenates onto Middlebrook 7H11 agar plates. For immunohistochemistry, fixed lung sections (2 µm) were de-waxed, rehydrated and subjected to antigen retrieval, blocking and exposed to goat polyclonal IDO antibody (LifeSpan Biosciences), followed by rabbit anti-goat alkaline phosphatase-conjugated secondary antibody. Alkaline phosphatase activity was visualized using new fuchsin (Dako). Sections were counterstained with hemalaun (Merck).

Ex-vivo Human Cells

Human DCs and macrophages were generated from buffy coats obtained from healthy volunteers recruited at the German Red Cross (Deutsches Rotes Kreuz). Peripheral blood mononuclear cells were recovered from interphase of Ficoll gradient according to standard protocols. To obtain untouched monocytes, cells were subjected to negative MACS selection using the Monocyte Isolation Kit II (Miltenyi Biotec). Flow-through cells were collected and grown in complete RPMI1640 medium containing 800 U/ml hGM-CSF and 500 U/ml hIL-4 (Strathmann) for 7 days to obtain immature DCs. Macrophages were generated from CD14⁺ cells by adherence to plastic culture flasks. For T cell assays, cells from healthy PPD⁺ donors were used to obtain DCs and T cells. DCs were generated as previously described and CD3⁺CD4⁺ and CD3⁺CD8⁺ T cells were obtained by negative MACS selection according to manufacturer’s instructions (Pan T cell Isolation Kit, Miltenyi).

Kynurenine Detection

Cell culture supernatants from infected and stimulated cells (DCs and macrophages) were deproteinized and kynurenines and tryptophan concentrations were measured by high-performance liquid chromatography using a Waters reverse phase C-18 column with a photodiode array detector (Waters 996). IDO1 enzymatic activity was assessed as a direct function of kynurenine production [71].

T Cell Assay

Monocyte-derived DCs treated with 1-methyl-tryptophan-DL (1-MT-DL) or with 1-MT-DL plus 3-OH-kynurenine (Kyn) and pulsed with PPD (10 µg/ml) and ManLAM (15 µg/ml) were cocultured with autologous T cells at a ratio of 1∶20 (10⁴ DCs and 2×10⁵ T cells in triplicates). T cells were labeled with carboxyfluorescein diacetate succinimidyl ester (CFSE, CellTrace CFSE Cell Proliferation Kit, Molecular Probes) at 2.5 uM according to manufacturer’s instructions. After 4 days of coculture, T cell proliferation was assessed by flow cytometry using anti-CD3, -CD4 and -CD8 monoclonal antibodies (all from BD) in combination with CFSE. Lymphocytes were gated based on forward and side scatter, and propidium iodide-positive cells were excluded from the analysis. The reduction in CFSE fluorescence intensity was analyzed on CD3⁺CD4⁺ and CD3⁺CD8⁺ gated cells. The ratio of proliferating to nonproliferating T cells was calculated and normalized to untreated controls. Data were analyzed by Anova with Dunnet’s multiple comparison test.

Quantitative RT-PCR

RNA was isolated with TRIzol (Invitrogen) and High Pure RNA Isolation Kit (Roche). cDNA was synthesized using the SuperScript III First-Strand Synthesis System and random hexamer primers (Invitrogen). RT-PCR was performed in triplicates using the ABIPRISM SDS 7900 system (Applied Biosystems). Fragments were amplified using the SYBR Green I Reaction Mix (Applied Biosystems) and specific primers for IDO1 (Hs_IDO1_1_SG QuantiTect Primer Assay, Qiagen) and normalized to the expression of human acidic ribosomal protein (HuPO) housekeeping gene (F-5′-GCTTCCTGGAGGGTGTCC-3′; R-5′-GGACTCGTTTGTACCCGTTG-3′). Threshold cycle values CT for each PCR product were determined using ABIPRISM SDS7900 software and ΔCT values were calculated. Fold differences (FD) in expression levels between IDO1 and HUPO were calculated according to the formula FD = 2^ΔCT.

Statistical Analysis

Relative abundances of biochemical compounds obtained by MS were normalized using the aquantile normalisation method from the R package limma [72] and used as an entry point for further statistical analyses. For a general comparison, biochemical compounds with altered levels in the different groups (TST^–, TST⁺ and TB^active) were obtained and compared using t-test and Wilcoxon sum rank test (see Table S2). Significant p values were adjusted for multiple testing using the Benjamini and Hochberg correction [73]. Also, several groups of biochemically similar compounds were tested in an ANOVA. Second, a RF classifier was used to identify biomarkers sufficiently distinct between groups. Functional linkage between different compounds was obtained by first fitting a simple linear model for each compound and subsequently calculating the correlation coefficient between the models’ residuals applying a bootstrapping procedure. In the applied linear model Y = aX+b, where Y equals group classification (TST^–, TST⁺ and TB^active) and X is the vector of the relative levels of the compounds. R scripts used to generate the data are available upon request. Significantly correlated compound profiles were used to infer clusters of functionally related small biochemical compounds. To explore the link between metabolic and cytokine profiles, we calculated Spearman correlation coefficients between the analyzed cytokine profiles and small metabolic compounds, and corrected for multiple testing. Next, to determine which categories of metabolic compounds correlated with cytokine and chemokine profiles, we performed an enrichment analysis using a hypergeometric test for the main categories of metabolic compounds (amino acids, carbohydrates, cofactors and vitamins, energy metabolism, lipids, nucleotides, peptides and xenobiotics).

Classification Analysis

The predictive value of small metabolic compounds as biomarkers in TB diagnosis was estimated using a supervised machine learning algorithm [23], as implemented in the R package “randomForest” (RF) [74]. The run parameters selected for sensitivity were mtry = 50 and ntree = 1,000. The classifiers were obtained for all three experimental groups as one data set, and in pairwise comparisons between each of the two groups as another data set. The findings were independently confirmed using SPLS-DA [24]. To estimate the number of variables needed for a robust classification, we applied an external leave-one-out (LOO) cross-validation. In each LOO round, one sample was removed from the data set. An RF model was constructed based on all variables, and the variables were sorted based on the importance from that model. Next, a series of secondary RF models were constructed based on the top 2, 3, 4, etc. variables from the importance list. Each of these secondary models was then validated against the LOO sample. The results were summarized for each number of variables tested, and plotted on Figure S1.

Clustering

Biochemical compounds were clustered using normalized relative compound levels and averaged hierarchical clustering, correlation-based distance and with bootstrapping implemented in the R package pvclust [75]. Clusters obtained for 10,000 bootstraps with approximately unbiased (AU) bootstrap confidence values >0.95 were considered significant.

Data Availability

All data, cluster assignments, scripts, and additional plots are available upon request.

List of Abbreviations

1-MT-DL, 1-methyl-DL-tryptophan

AU, approximately unbiased

BCG, bacille Calmette-Guérin

CFSE, carboxyfluorescein diacetate succinimidyl ester

CMPF, 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid

CXCL-10, C-X-C motif chemokine 10

DC, dendritic cell

EGF, epidermal growth factor

FGF, fibroblast growth factor

GC–MS, gas chromatography–mass spectrometry

G-CSF, granulocyte colony-stimulating factor

GITR, glucocorticoid-induced tumor necrosis factor receptor

GM-CSF, granulocyte macrophage colony-stimulating factor

HPLC, high-performance liquid chromatography

HuPO, human acidic ribosomal protein

IDO1, indoleamine 2,3 dioxygenase 1

IFN-γ, interferon-gamma

IL, interleukin

IP, interferon-inducible protein

LC–MSMS, liquid chromatography-tandem mass spectrometry

LOO, leave-one-out cross-validation

LPCs, lysophosphatidylcholines

ManLAM, mannosylated lipoarabinomannan

MCFA, medium-chain fatty acid

MCP, monocyte chemotactic protein

MDC, macrophage-derived cytokine

MGIT, Mycobacteria Growth Indicator Tube

MIP, macrophage inflammatory protein

NAM, N-acetyl-L-methionine

NMR, nuclear magnetic resonance

PDGF, platelet-derived growth factor

PPD, purified protein derivative

RF, random forests

sIL2ra, soluble alpha chain of the IL-2 receptor

SPLS-DA, sparse PLS discriminant analysis

TB, tuberculosis

TB^active, patients with active TB

TGF, transforming growth factor

TNF, tumor necrosis factor

TST, tuberculin skin test

TST^–, healthy M. tuberculosis-uninfected controls

TST⁺, latently M. tuberculosis-infected healthy individuals

VEGF, vascular endothelial growth factor

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Fifteen to twenty metabolites suffice to distinguish between TB patients and latently infected groups. Figure shows the decrease of the average classification error rate as a function of the number of different small metabolic compounds chosen for the classification. Error bars denote standard error of the mean of 50 re-sampling procedures.

https://doi.org/10.1371/journal.pone.0040221.s001

(TIF)

Figure S2.

Heatmap showing relative levels of small metabolic compounds in TB patients (TB^active) and latently infected individuals (TST⁺). For purposes of illustration, 10 profiles from each group were randomly assigned to a test set, and calculations were repeated for the remaining training set. Left, training-set levels; only 50 compounds selected by variable importance from the RF model were chosen. Right, test-set levels; all test set samples were correctly assigned to a study class.

https://doi.org/10.1371/journal.pone.0040221.s002

(TIF)

Figure S3.

Network showing functional relationships between the small metabolic compounds in TB patients, healthy uninfected and latently infected individuals. Nodes correspond to metabolites; edges correspond to statistically significant correlation between residual small metabolite profiles corrected for study classes. Colors correspond to differences between the TST^– and TB^active (A), TST⁺ and TB^active (B) or TST^– and TST⁺ classes (C). Color intensity indicates significance of difference with darker colors corresponding to more significant differences. Metabolites with adjusted p value >0.05 are not colored. Line widths correspond to the absolute Spearman correlation coefficients corrected for groups (see “Methods”). Figure S3.A is the same as Figure 3 in the Manuscript and has been included here for completeness.

https://doi.org/10.1371/journal.pone.0040221.s003

(TIF)

Table S1.

Demographic characteristics of study subjects.

https://doi.org/10.1371/journal.pone.0040221.s004

(DOCX)

Table S2.

Significant differences between serum concentrations of small metabolic compounds between each two groups. A), differences between TST^- and TST⁺; B), differences between TST^- and TB^active; C) differences between TST⁺ and TB^active. p denotes the p value from Wilcoxon rank sum test; q is the p value corrected for multiple testing using false discovery rate correction (Benjamini 1995). Only compounds for which the q value was smaller than 0.01 are shown.

https://doi.org/10.1371/journal.pone.0040221.s005

(DOCX)

Table S3.

Top 50 compounds sorted by importance from the random forests analysis of biochemical compound discriminatory power between TST⁺ and TB^active. Compounds with names starting with an “X” could not be uniquely identified in the analysis. The columns TST⁻, TST⁺ and TB^active give the average relative abundance levels in the corresponding study groups. Top twenty compounds are sufficient to sensitively discriminate latent infection from active TB. See text for details.

https://doi.org/10.1371/journal.pone.0040221.s006

(DOCX)

Table S4.

Significant clusters of small metabolic compounds containing at least five different members.

https://doi.org/10.1371/journal.pone.0040221.s007

(DOCX)

Acknowledgments

We wish to thank Anca Dorhoi and Silke Bandermann for support, and Mary Louise Grossman for her commitment in editing and correcting this manuscript.

Author Contributions

Conceived and designed the experiments: SHEK SKP JW RM. Performed the experiments: RM CAG KCF JM. Analyzed the data: JW CAS RM DR AT. Contributed reagents/materials/analysis tools: GW GFB. Wrote the paper: JW SHEK.

Citation: Weiner J 3rd, Parida SK, Maertzdorf J, Black GF, Repsilber D, Telaar A, et al. (2012) Biomarkers of Inflammation, Immunosuppression and Stress Are Revealed by Metabolomic Profiling of Tuberculosis Patients. PLoS ONE 7(7): e40221. https://doi.org/10.1371/journal.pone.0040221

References

1. Illig T, Gieger C, Zhai G, Römisch-Margl W, Wang-Sattler R, et al. (2009) A genome-wide perspective of genetic variation in human metabolism. Nat Genet 42: 137–141.T. IlligC. GiegerG. ZhaiW. Römisch-MarglR. Wang-Sattler2009A genome-wide perspective of genetic variation in human metabolism.Nat Genet42137141

2. Kaddurah-Daouk R, Kristal BS, Weinshilboum RM (2008) Metabolomics: a global biochemical approach to drug response and disease. Pharmacol Toxicol 48: 653–683.R. Kaddurah-DaoukBS KristalRM Weinshilboum2008Metabolomics: a global biochemical approach to drug response and disease.Pharmacol Toxicol48653683

3. Han X, Holtzman M, McKeel DW Jr, Kelley J, Morris JC (2002) Substantial sulfatide deficiency and ceramide elevation in very early Alzheimer’s disease: potential role in disease pathogenesis. J Neurochem 82: 809–818. 0997 [pii]. X. HanM. HoltzmanDW McKeel JrJ. KelleyJC Morris2002Substantial sulfatide deficiency and ceramide elevation in very early Alzheimer’s disease: potential role in disease pathogenesis. J Neurochem 82: 809–818.0997 [pii]

4. Bogdanov M, Matson WR, Wang L, Matson T, Saunders-Pullman R, et al. (2008) Metabolomic profiling to develop blood biomarkers for Parkinson’s disease. Brain 131: 389–396.M. BogdanovWR MatsonL. WangT. MatsonR. Saunders-Pullman2008Metabolomic profiling to develop blood biomarkers for Parkinson’s disease.Brain131389396

5. Sabatine MS, Liu E, Morrow DA, Heller E, McCarroll R, et al. (2005) Metabolomic identification of novel biomarkers of myocardial ischemia. Circulation 112: 3868–3875.MS SabatineE. LiuDA MorrowE. HellerR. McCarroll2005Metabolomic identification of novel biomarkers of myocardial ischemia.Circulation11238683875

6. Brindle JT, Nicholson JK, Schofield PM, Grainger DJ, Holmes E (2003) Application of chemometrics to 1 H NMR spectroscopic data to investigate a relationship between human serum metabolic profiles and hypertension. Analyst 128: 32–36.JT BrindleJK NicholsonPM SchofieldDJ GraingerE. Holmes2003Application of chemometrics to 1 H NMR spectroscopic data to investigate a relationship between human serum metabolic profiles and hypertension.Analyst1283236

7. Yang J, Xu G, Zheng Y, Kong H, Pang T, et al. (2004) Diagnosis of liver cancer using HPLC-based metabonomics avoiding false-positive result from hepatitis and hepatocirrhosis diseases. J Chromatogr B 813: 59–65.J. YangG. XuY. ZhengH. KongT. Pang2004Diagnosis of liver cancer using HPLC-based metabonomics avoiding false-positive result from hepatitis and hepatocirrhosis diseases.J Chromatogr B8135965

8. Wang C, Kong H, Guan Y, Yang J, Gu J, et al. (2005) Plasma phospholipid metabolic profiling and biomarkers of type 2 diabetes mellitus based on high-performance liquid chromatography/electrospray mass spectrometry and multivariate statistical analysis. Anal Chem 77: 4108–4116.C. WangH. KongY. GuanJ. YangJ. Gu2005Plasma phospholipid metabolic profiling and biomarkers of type 2 diabetes mellitus based on high-performance liquid chromatography/electrospray mass spectrometry and multivariate statistical analysis.Anal Chem7741084116

9. Yuan K, Kong H, Guan Y, Yang J, Xu G (2007) A GC-based metabonomics investigation of type 2 diabetes by organic acids metabolic profile. J Chromatogr B 850: 236–240.K. YuanH. KongY. GuanJ. YangG. Xu2007A GC-based metabonomics investigation of type 2 diabetes by organic acids metabolic profile.J Chromatogr B850236240

10. Salek RM, Maguire ML, Bentley E, Rubtsov DV, Hough T, et al. (2007) A metabolomic comparison of urinary changes in type 2 diabetes in mouse, rat, and human. Physiological genomics 29: 99–108.RM SalekML MaguireE. BentleyDV RubtsovT. Hough2007A metabolomic comparison of urinary changes in type 2 diabetes in mouse, rat, and human.Physiological genomics2999108

11. Suhre K, Wallaschofski H, Raffler J, Friedrich N, Haring R, et al. (2011) A genome-wide association study of metabolic traits in human urine. Nat Genet 43: 565–569.K. SuhreH. WallaschofskiJ. RafflerN. FriedrichR. Haring2011A genome-wide association study of metabolic traits in human urine.Nat Genet43565569

12. Laiakis EC, Morris GAJ, Fornace AJ, Howie SRC (2010) Metabolomic analysis in severe childhood pneumonia in the Gambia, West Africa: findings from a pilot study. PloS one 5: e12655.EC LaiakisGAJ MorrisAJ FornaceSRC Howie2010Metabolomic analysis in severe childhood pneumonia in the Gambia, West Africa: findings from a pilot study.PloS one5e12655

13. Wang Y, Utzinger J, Saric J, Li JV, Burckhardt J, et al. (2008) Global metabolic responses of mice to Trypanosoma brucei brucei infection. Proc Natl Acad Sci USA 105: 6127–6132.Y. WangJ. UtzingerJ. SaricJV LiJ. Burckhardt2008Global metabolic responses of mice to Trypanosoma brucei brucei infection.Proc Natl Acad Sci USA10561276132

14. Slupsky CM, Rankin KN, Fu H, Chang D, Rowe BH, et al. (2009) Pneumococcal pneumonia: potential for diagnosis through a urinary metabolic profile. J Proteome Res 8: 5550–5558.CM SlupskyKN RankinH. FuD. ChangBH Rowe2009Pneumococcal pneumonia: potential for diagnosis through a urinary metabolic profile.J Proteome Res855505558

15. Somashekar BS, Amin AG, Rithner CD, Troudt JL, Basaraba R, et al. (2011) Metabolic profiling of lung granuloma in Mycobacterium tuberculosis infected guinea pigs: ex vivo 1 H magic angle spinning NMR studies. J Proteome Res 10: 4186–4195.BS SomashekarAG AminCD RithnerJL TroudtR. Basaraba2011Metabolic profiling of lung granuloma in Mycobacterium tuberculosis infected guinea pigs: ex vivo 1 H magic angle spinning NMR studies.J Proteome Res1041864195

16. Shin JH, Yang JY, Jeon BY, Yoon YJ, Cho SN, et al. (2011) 1H-NMR-based metabolomic profiling in mice infected with Mycobacterium tuberculosis. J Proteome Res 10: 2238–2247.JH ShinJY YangBY JeonYJ YoonSN Cho20111H-NMR-based metabolomic profiling in mice infected with Mycobacterium tuberculosis.J Proteome Res1022382247

17. Lawn SD, Zumla AI (2011) Tuberculosis. Lancet 378: 57–72.SD LawnAI Zumla2011Tuberculosis.Lancet3785772

18. Philips JA, Ernst JD (2012) Tuberculosis Pathogenesis and Immunity. Ann Rev Pathol Mech Dis 7: 353–384.JA PhilipsJD Ernst2012Tuberculosis Pathogenesis and Immunity.Ann Rev Pathol Mech Dis7353384

19. Lonnroth K, Castro KG, Chakaya JM, Chauhan LS, Floyd K, et al. (2010) Tuberculosis control and elimination 2010–50: cure, care, and social development. Lancet 375: 1814–1829.K. LonnrothKG CastroJM ChakayaLS ChauhanK. Floyd2010Tuberculosis control and elimination 2010–50: cure, care, and social development.Lancet37518141829

20. Jacobsen M, Repsilber D, Gutschmidt A, Neher A, Feldmann K, et al. (2007) Candidate biomarkers for discrimination between infection and disease caused by Mycobacterium tuberculosis. J Mol Med 85: 613–621.M. JacobsenD. RepsilberA. GutschmidtA. NeherK. Feldmann2007Candidate biomarkers for discrimination between infection and disease caused by Mycobacterium tuberculosis.J Mol Med85613621

21. Berry MPR, Graham CM, McNab FW, Xu Z, Bloch SAA, et al. (2010) An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature 466: 973–977.MPR BerryCM GrahamFW McNabZ. XuSAA Bloch2010An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis.Nature466973977

22. Maertzdorf J, Repsilber D, Parida SK, Stanley K, Roberts T, et al. (2011) Human gene expression profiles of susceptibility and resistance in tuberculosis. Genes Immun 12: 15–22.J. MaertzdorfD. RepsilberSK ParidaK. StanleyT. Roberts2011Human gene expression profiles of susceptibility and resistance in tuberculosis.Genes Immun121522

23. Breiman L (2001) Random forests. Mach Learn 45: 5–32.L. Breiman2001Random forests.Mach Learn45532

24. Lê Cao K-A, González I, Déjean S (2009) integrOmics: an R package to unravel relationships between two omics datasets. Bioinformatics 25: 2855.K-A Lê CaoI. GonzálezS. Déjean2009integrOmics: an R package to unravel relationships between two omics datasets.Bioinformatics252855

25. Mellor AL, Munn DH (2004) IDO expression by dendritic cells: tolerance and tryptophan catabolism. Nature Rev Immunol 4: 762–774.AL MellorDH Munn2004IDO expression by dendritic cells: tolerance and tryptophan catabolism.Nature Rev Immunol4762774

26. Munn DH, Zhou M, Attwood JT, Bondarev I, Conway SJ, et al. (1998) Prevention of allogeneic fetal rejection by tryptophan catabolism. Science 281: 1191.DH MunnM. ZhouJT AttwoodI. BondarevSJ Conway1998Prevention of allogeneic fetal rejection by tryptophan catabolism.Science2811191

27. Munn DH, Sharma MD, Hou D, Baban B, Lee JR, et al. (2004) Expression of indoleamine 2, 3-dioxygenase by plasmacytoid dendritic cells in tumor-draining lymph nodes. J Clin Invest 114: 280–290.DH MunnMD SharmaD. HouB. BabanJR Lee2004Expression of indoleamine 2, 3-dioxygenase by plasmacytoid dendritic cells in tumor-draining lymph nodes.J Clin Invest114280290

28. Katz JB, Muller AJ, Prendergast GC (2008) Indoleamine 2, 3-dioxygenase in T-cell tolerance and tumoral immune escape. Immunol Rev 222: 206–221.JB KatzAJ MullerGC Prendergast2008Indoleamine 2, 3-dioxygenase in T-cell tolerance and tumoral immune escape.Immunol Rev222206221

29. Popov A, Abdullah Z, Wickenhauser C, Saric T, Driesen J, et al. (2006) Indoleamine 2, 3-dioxygenase-expressing dendritic cells form suppurative granulomas following Listeria monocytogenes infection. J Clin Invest 116: 3160.A. PopovZ. AbdullahC. WickenhauserT. SaricJ. Driesen2006Indoleamine 2, 3-dioxygenase-expressing dendritic cells form suppurative granulomas following Listeria monocytogenes infection.J Clin Invest1163160

30. Boasso A, Herbeuval JP, Hardy AW, Anderson SA, Dolan MJ, et al. (2007) HIV inhibits CD4+ T-cell proliferation by inducing indoleamine 2, 3-dioxygenase in plasmacytoid dendritic cells. Blood 109: 3351.A. BoassoJP HerbeuvalAW HardySA AndersonMJ Dolan2007HIV inhibits CD4+ T-cell proliferation by inducing indoleamine 2, 3-dioxygenase in plasmacytoid dendritic cells.Blood1093351

31. Romani L, Fallarino F, De Luca A, Montagnoli C, D’Angelo C, et al. (2008) Defective tryptophan catabolism underlies inflammation in mouse chronic granulomatous disease. Nature 451: 211–215.L. RomaniF. FallarinoA. De LucaC. MontagnoliC. D’Angelo2008Defective tryptophan catabolism underlies inflammation in mouse chronic granulomatous disease.Nature451211215

32. Favre D, Mold J, Hunt PW, Kanwar B, Loke P, et al. (2010) Tryptophan catabolism by indoleamine 2, 3-dioxygenase 1 alters the balance of TH17 to regulatory T cells in HIV disease. Science Trans Med 2: 32ra36.D. FavreJ. MoldPW HuntB. KanwarP. Loke2010Tryptophan catabolism by indoleamine 2, 3-dioxygenase 1 alters the balance of TH17 to regulatory T cells in HIV disease.Science Trans Med232ra36

33. Makala LHC, Baban B, Lemos H, El-Awady AR, Chandler PR, et al. (2011) Leishmania major attenuates host immunity by stimulating local indoleamine 2, 3-dioxygenase expression. J Infect Dis 203: 715.LHC MakalaB. BabanH. LemosAR El-AwadyPR Chandler2011Leishmania major attenuates host immunity by stimulating local indoleamine 2, 3-dioxygenase expression.J Infect Dis203715

34. Platten M, Ho PP, Youssef S, Fontoura P, Garren H, et al. (2005) Treatment of autoimmune neuroinflammation with a synthetic tryptophan metabolite. Science 310: 850.M. PlattenPP HoS. YoussefP. FontouraH. Garren2005Treatment of autoimmune neuroinflammation with a synthetic tryptophan metabolite.Science310850

35. Sharma MD, Hou DY, Liu Y, Koni PA, Metz R, et al. (2009) Indoleamine 2, 3-dioxygenase controls conversion of Foxp3+ Tregs to TH17-like cells in tumor-draining lymph nodes. Blood 113: 6102.MD SharmaDY HouY. LiuPA KoniR. Metz2009Indoleamine 2, 3-dioxygenase controls conversion of Foxp3+ Tregs to TH17-like cells in tumor-draining lymph nodes.Blood1136102

36. Sharma MD, Hou DY, Baban B, Koni PA, He Y, et al. (2010) Reprogrammed Foxp3+ regulatory T cells provide essential help to support cross-presentation and CD8+ T cell priming in naive mice. Immunity 33: 942–954.MD SharmaDY HouB. BabanPA KoniY. He2010Reprogrammed Foxp3+ regulatory T cells provide essential help to support cross-presentation and CD8+ T cell priming in naive mice.Immunity33942954

37. Müller A, Heseler K, Schmidt SK, Spekker K, MacKenzie CR, et al. (2009) The missing link between indoleamine 2, 3-dioxygenase mediated antibacterial and immunoregulatory effects. J Cell Mol Med 13: 1125–1135.A. MüllerK. HeselerSK SchmidtK. SpekkerCR MacKenzie2009The missing link between indoleamine 2, 3-dioxygenase mediated antibacterial and immunoregulatory effects.J Cell Mol Med1311251135

38. Grohmann U, Volpi C, Fallarino F, Bozza S, Bianchi R, et al. (2007) Reverse signaling through GITR ligand enables dexamethasone to activate IDO in allergy. Nature Med 13: 579–586.U. GrohmannC. VolpiF. FallarinoS. BozzaR. Bianchi2007Reverse signaling through GITR ligand enables dexamethasone to activate IDO in allergy.Nature Med13579586

39. Suzuki Y, Suda T, Asada K, Miwa S, Suzuki M, et al. (2012) Serum indoleamine 2, 3-dioxygenase activity predicts prognosis of pulmonary tuberculosis. Clinical and Vaccine Immunology. Y. SuzukiT. SudaK. AsadaS. MiwaM. Suzuki2012Serum indoleamine 2, 3-dioxygenase activity predicts prognosis of pulmonary tuberculosis.Clinical and Vaccine Immunology

40. Mehta D, Gupta S, Gaur SN, Gangal SV, Agrawal KP (1990) Increased leukocyte phospholipase A2 activity and plasma lysophosphatidylcholine levels in asthma and rhinitis and their relationship to airway sensitivity to histamine. Am Rev Resp Dis 142: 157–161.D. MehtaS. GuptaSN GaurSV GangalKP Agrawal1990Increased leukocyte phospholipase A2 activity and plasma lysophosphatidylcholine levels in asthma and rhinitis and their relationship to airway sensitivity to histamine.Am Rev Resp Dis142157161

41. Okita M, Gaudette DC, Mills GB, Holub BJ (1997) Elevated levels and altered fatty acid composition of plasma lysophosphatidlycholine(lysoPC) in ovarian cancer patients. Int J Cancer 71: 31–34.M. OkitaDC GaudetteGB MillsBJ Holub1997Elevated levels and altered fatty acid composition of plasma lysophosphatidlycholine(lysoPC) in ovarian cancer patients.Int J Cancer713134

42. Taylor L, Arends J, Hodina A, Unger C, Massing U (2007) Plasma lyso-phosphatidylcholine concentration is decreased in cancer patients with weight loss and activated inflammatory status. Lipids Health Dis 6: 17.L. TaylorJ. ArendsA. HodinaC. UngerU. Massing2007Plasma lyso-phosphatidylcholine concentration is decreased in cancer patients with weight loss and activated inflammatory status.Lipids Health Dis617

43. Drobnik W, Liebisch G, Audebert FX, Frohlich D, Gluck T, et al. (2003) Plasma ceramide and lysophosphatidylcholine inversely correlate with mortality in sepsis patients. J Lipid Res 44: 754–761.W. DrobnikG. LiebischFX AudebertD. FrohlichT. Gluck2003Plasma ceramide and lysophosphatidylcholine inversely correlate with mortality in sepsis patients.J Lipid Res44754761

44. Duan L, Gan H, Arm J, Remold HG (2001) Cytosolic Phospholipase A2 Participates with TNF-{{alpha}} in the Induction of Apoptosis of Human Macrophages Infected with Mycobacterium tuberculosis H37Ra. J Immunol 166: 7469–7476.L. DuanH. GanJ. ArmHG Remold2001Cytosolic Phospholipase A2 Participates with TNF-{{alpha}} in the Induction of Apoptosis of Human Macrophages Infected with Mycobacterium tuberculosis H37Ra.J Immunol16674697476

45. Lewis GD, Wei RU, Liu E, Yang E, Shi X, et al. (2008) Metabolite profiling of blood from individuals undergoing planned myocardial infarction reveals early markers of myocardial injury. J Clin Invest 118: 3503.GD LewisRU WeiE. LiuE. YangX. Shi2008Metabolite profiling of blood from individuals undergoing planned myocardial infarction reveals early markers of myocardial injury.J Clin Invest1183503

46. Bell MJ, Kochanek PM, Carcillo JA, Mi Z, Schiding JK, et al. (1998) Interstitial adenosine, inosine, and hypoxanthine are increased after experimental traumatic brain injury in the rat. J Neurotraum 15: 163–170.MJ BellPM KochanekJA CarcilloZ. MiJK Schiding1998Interstitial adenosine, inosine, and hypoxanthine are increased after experimental traumatic brain injury in the rat.J Neurotraum15163170

47. Jin X, Shepherd RK, Duling BR, Linden J (1997) Inosine binds to A3 adenosine receptors and stimulates mast cell degranulation. J Clin Invest 100: 2849–2857.X. JinRK ShepherdBR DulingJ. Linden1997Inosine binds to A3 adenosine receptors and stimulates mast cell degranulation.J Clin Invest10028492857

48. Gomez G, Sitkovsky MV (2003) Differential requirement for A2a and A3 adenosine receptors for the protective effect of inosine in vivo. Blood 102: 4472–4478.G. GomezMV Sitkovsky2003Differential requirement for A2a and A3 adenosine receptors for the protective effect of inosine in vivo.Blood10244724478

49. Den Boon S, Van Lill SWP, Borgdorff MW, Verver S, Bateman ED, et al. (2005) Association between smoking and tuberculosis infection: a population survey in a high tuberculosis incidence area. Thorax 60: 555–557.S. Den BoonSWP Van LillMW BorgdorffS. VerverED Bateman2005Association between smoking and tuberculosis infection: a population survey in a high tuberculosis incidence area.Thorax60555557

50. Dever JT, Elfarra AA (2006) In Vivo Metabolism of l-Methionine in Mice: Evidence for Stereoselective Formation of Methionine-d-Sulfoxide and Quantitation of Other Major Metabolites. Drug Metab Dispos 34: 2036–2043.JT DeverAA Elfarra2006In Vivo Metabolism of l-Methionine in Mice: Evidence for Stereoselective Formation of Methionine-d-Sulfoxide and Quantitation of Other Major Metabolites.Drug Metab Dispos3420362043

51. Avila MA, Carretero MV, Rodriguez EN, Mato JM (1998) Regulation by hypoxia of methionine adenosyltransferase activity and gene expression in rat hepatocytes. Gastroenterology 114: 364–371.MA AvilaMV CarreteroEN RodriguezJM Mato1998Regulation by hypoxia of methionine adenosyltransferase activity and gene expression in rat hepatocytes.Gastroenterology114364371

52. Kono H, Enomoto N, Connor HD, Wheeler MD, Bradford BU, et al. (2000) Medium-chain triglycerides inhibit free radical formation and TNF-alpha production in rats given enteral ethanol. Am J Physiol Gastrointest Liver Physiol 278: G467–G476.H. KonoN. EnomotoHD ConnorMD WheelerBU Bradford2000Medium-chain triglycerides inhibit free radical formation and TNF-alpha production in rats given enteral ethanol.Am J Physiol Gastrointest Liver Physiol278G467G476

53. Harrison RA, Farries TC, Northrop FD, Lachmann PJ, Davis AE (1988) Structure of C3f, a small peptide specifically released during inactivation of the third component of complement. Complement 5: 27–32.RA HarrisonTC FarriesFD NorthropPJ LachmannAE Davis1988Structure of C3f, a small peptide specifically released during inactivation of the third component of complement.Complement52732

54. Sakamoto M, Fujisawa Y, Nishioka K (1998) Physiologic role of the complement system in host defense, disease, and malnutrition. Nutrition 14: 391–398.M. SakamotoY. FujisawaK. Nishioka1998Physiologic role of the complement system in host defense, disease, and malnutrition.Nutrition14391398

55. G. SzeplakiZ. ProhaszkaJ. DubaS. Rugonfalvi-KissI. Karadi 2004 Association of high serum concentration of the third component of complement (C3) with pre-existing severe coronary artery disease and new vascular events in women. Atherosclerosis 177: 383–389. S0021–9150(04)00394–6 [pii];10.1016/j.atherosclerosis.2004.07.022. [doi]

56. Chan TYK (2000) Vitamin D Deficiency and Susceptibility to Tuberculosis. Calcified Tissue Int 66: 476–478.TYK Chan2000Vitamin D Deficiency and Susceptibility to Tuberculosis.Calcified Tissue Int66476478

57. Liu PT, Stenger S, Li H, Wenzel L, Tan BH, et al. (2006) Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response. Science 311: 1770.PT LiuS. StengerH. LiL. WenzelBH Tan2006Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response.Science3111770

58. Wilkinson RJ, Llewelyn M, Toossi Z, Patel P, Pasvol G, et al. (2000) Influence of vitamin D deficiency and vitamin D receptor polymorphisms on tuberculosis among Gujarati Asians in west London: a case-control study. Lancet 355: 618–621.RJ WilkinsonM. LlewelynZ. ToossiP. PatelG. Pasvol2000Influence of vitamin D deficiency and vitamin D receptor polymorphisms on tuberculosis among Gujarati Asians in west London: a case-control study.Lancet355618621

59. Ustianowski A, Shaffer R, Collin S, Wilkinson RJ, Davidson RN (2005) Prevalence and associations of vitamin D deficiency in foreign-born persons with tuberculosis in London. J Infection 50: 432–437.A. UstianowskiR. ShafferS. CollinRJ WilkinsonRN Davidson2005Prevalence and associations of vitamin D deficiency in foreign-born persons with tuberculosis in London.J Infection50432437

60. Cohen G, Rudnicki M, Hoorl WH (2001) Uremic toxins modulate the spontaneous apoptotic cell death and essential functions of neutrophils. Kidney Int 59: S48–S52.G. CohenM. RudnickiWH Hoorl2001Uremic toxins modulate the spontaneous apoptotic cell death and essential functions of neutrophils.Kidney Int59S48S52

61. Hutchings A, van Staden J (1994) Plants used for stress-related ailments in traditional Zulu, Xhosa and Sotho medicine. Part 1: Plants used for headaches. J Ethnopharmacol 43: 89–124. 0378–8741(94)90008–6 [pii]. A. HutchingsJ. van Staden1994Plants used for stress-related ailments in traditional Zulu, Xhosa and Sotho medicine. Part 1: Plants used for headaches. J Ethnopharmacol 43: 89–124.0378–8741(94)90008–6 [pii]

62. Chang WK, Yang KD, Shaio MF (1999) Lymphocyte proliferation modulated by glutamine: involved in the endogenous redox reaction. Clin Exp Immunol 117: 482–488.WK ChangKD YangMF Shaio1999Lymphocyte proliferation modulated by glutamine: involved in the endogenous redox reaction.Clin Exp Immunol117482488

63. Parry-Billings M, Calder PC, Newsholme EA, Evans J (1990) Does glutamine contribute to immunosuppression after major burns? Lancet 336: 523–525.M. Parry-BillingsPC CalderEA NewsholmeJ. Evans1990Does glutamine contribute to immunosuppression after major burns?Lancet336523525

64. Kinscherf R, Hack V, Fischbach T, Friedmann B, Weiss C, et al. (1996) Low plasma glutamine in combination with high glutamate levels indicate risk for loss of body cell mass in healthy individuals: the effect of N-acetyl-cysteine. J Mol Med 74: 393–400.R. KinscherfV. HackT. FischbachB. FriedmannC. Weiss1996Low plasma glutamine in combination with high glutamate levels indicate risk for loss of body cell mass in healthy individuals: the effect of N-acetyl-cysteine.J Mol Med74393400

65. Huettmann A, Duehrsen U, Stypmann J, Noppeney R, Nueckel H, et al. (2006) Granulocyte colony-stimulating factor-induced blood stem cell mobilisation in patients with chronic heart failure. Basic Res Cardiol 101: 78–86.A. HuettmannU. DuehrsenJ. StypmannR. NoppeneyH. Nueckel2006Granulocyte colony-stimulating factor-induced blood stem cell mobilisation in patients with chronic heart failure.Basic Res Cardiol1017886

66. To LB, Haylock DN, Simmons PJ, Juttner CA (1997) The biology and clinical uses of blood stem cells. Blood 89: 2233.LB ToDN HaylockPJ SimmonsCA Juttner1997The biology and clinical uses of blood stem cells.Blood892233

67. Schwenk A, Macallan DC (2000) Tuberculosis, malnutrition and wasting. Curr Opin Clin Nutr Metab Care 3: 285.A. SchwenkDC Macallan2000Tuberculosis, malnutrition and wasting.Curr Opin Clin Nutr Metab Care3285

68. Evans AM, DeHaven CD, Barrett T, Mitchell M, Milgram E (2009) Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal Chem 81: 6656–6667.AM EvansCD DeHavenT. BarrettM. MitchellE. Milgram2009Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems.Anal Chem8166566667

69. Sha W, da Costa KA, Fischer LM, Milburn MV, Lawton KA, et al. (2010) Metabolomic profiling can predict which humans will develop liver dysfunction when deprived of dietary choline. FASEB J 24: 2962–2975.W. ShaKA da CostaLM FischerMV MilburnKA Lawton2010Metabolomic profiling can predict which humans will develop liver dysfunction when deprived of dietary choline.FASEB J2429622975

70. DeHaven C, Evans A, Dai H, Lawton K (2010) Organization of GC/MS and LC/MS metabolomics data into chemical libraries. Journal of cheminformatics 2: 9.C. DeHavenA. EvansH. DaiK. Lawton2010Organization of GC/MS and LC/MS metabolomics data into chemical libraries.Journal of cheminformatics29

71. Fallarino F, Grohmann U, Hwang KW, Orabona C, Vacca C, et al. (2003) Modulation of tryptophan catabolism by regulatory T cells. Nat Immunol 4: 1206–1212.F. FallarinoU. GrohmannKW HwangC. OrabonaC. Vacca2003Modulation of tryptophan catabolism by regulatory T cells.Nat Immunol412061212

72. Smyth GK (2005) Limma: linear models for microarray data. In: Bioinformatics and Computational Biology Solutions Using R and Bioconductor. New York: Springer. pp. 397–420.GK Smyth2005Limma: linear models for microarray data.In: Bioinformatics and Computational Biology Solutions Using R and Bioconductor. New York: Springer.397420

73. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B Met 57: 289–300.Y. BenjaminiY. Hochberg1995Controlling the false discovery rate: a practical and powerful approach to multiple testing.J Roy Stat Soc B Met57289300

74. Liaw A, Wiener M (2002) Classification and Regression by randomForest. R news 2: 18–22.A. LiawM. Wiener2002Classification and Regression by randomForest.R news21822

75. Suzuki R, Shimodaira H (2006) Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22: 1540–1542.R. SuzukiH. Shimodaira2006Pvclust: an R package for assessing the uncertainty in hierarchical clustering.Bioinformatics2215401542