- Open Access
Metabolome in schizophrenia and other psychotic disorders: a general population-based study
Genome Medicine volume 3, Article number: 19 (2011)
Persons with schizophrenia and other psychotic disorders have a high prevalence of obesity, impaired glucose tolerance, and lipid abnormalities, particularly hypertriglyceridemia and low high-density lipoprotein. More detailed molecular information on the metabolic abnormalities may reveal clues about the pathophysiology of these changes, as well as about disease specificity.
We applied comprehensive metabolomics in serum samples from a general population-based study in Finland. The study included all persons with DSM-IV primary psychotic disorder (schizophrenia, n = 45; other non-affective psychosis (ONAP), n = 57; affective psychosis, n = 37) and controls matched by age, sex, and region of residence. Two analytical platforms for metabolomics were applied to all serum samples: a global lipidomics platform based on ultra-performance liquid chromatography coupled to mass spectrometry, which covers molecular lipids such as phospholipids and neutral lipids; and a platform for small polar metabolites based on two-dimensional gas chromatography coupled to time-of-flight mass spectrometry (GC × GC-TOFMS).
Compared with their matched controls, persons with schizophrenia had significantly higher metabolite levels in six lipid clusters containing mainly saturated triglycerides, and in two small-molecule clusters containing, among other metabolites, (1) branched chain amino acids, phenylalanine and tyrosine, and (2) proline, glutamic, lactic and pyruvic acids. Among these, serum glutamic acid was elevated in all psychoses (P = 0.0020) compared to controls, while proline upregulation (P = 0.000023) was specific to schizophrenia. After adjusting for medication and metabolic comorbidity in linear mixed models, schizophrenia remained independently associated with higher levels in seven of these eight clusters (P < 0.05 in each cluster). The metabolic abnormalities were less pronounced in persons with ONAP or affective psychosis.
Our findings suggest that specific metabolic abnormalities related to glucoregulatory processes and proline metabolism are specifically associated with schizophrenia and reflect two different disease-related pathways. Metabolomics, which is sensitive to both genetic and environmental variation, may become a powerful tool in psychiatric research to investigate disease susceptibility, clinical course, and treatment response.
Psychotic disorders are among the most severe and impairing medical diseases . Schizophrenia is the most common of them, with a lifetime prevalence of 1% in a general population . The current view is that schizophrenia is a developmental disorder caused by a combination of genetic vulnerability, early environmental insults, subtle developmental and cognitive impairments, and later influences such as social adversity and drug abuse , with heritability of about 80% [4, 5]. The Diagnostic and Statistical Manual of Mental Disorders (DSM)-IV divides primary psychotic disorders into nine different diagnoses based on symptom patterns, clinical course and outcome, although it is unclear whether this has any etiological justification. Nevertheless, while there is overlap in genetic vulnerability between different psychotic disorders, like schizophrenia and bipolar I disorder, they also have non-shared genetic and environmental risk factors [5, 6]. Given the multi-factorial complexity of psychotic disorders , identification of molecular markers sensitive to the underlying pathogenic factors of specific diseases would be of high relevance, not only to assist in their early detection and diagnosis, but also to subsequently facilitate disease monitoring and treatment responses.
Metabolomics is a discipline dedicated to the global study of small molecules (that is, metabolites) in cells, tissues, and biofluids. Concentration changes of specific groups of circulating metabolites may be sensitive to pathogenically relevant factors, such as genetic variation, diet, age, or gut microbiota [8–12]. Over the past years, technologies have been developed that allow comprehensive and quantitative investigation of a multitude of different metabolites . The study of high-dimensional chemical signatures as obtained by metabolomics may therefore be a powerful tool for characterization of complex phenotypes affected by both genetic and environmental factors . Previous metabolomic studies in schizophrenia and related psychoses have highlighted the importance of glucoregulatory processes [15, 16] and tryptophan metabolism  in psychosis, and lipidomics approaches have identified specific drug-response profiles for three commonly used atypical antipsychotics . However, no metabolomics studies have so far been conducted to discriminate between different groups of psychotic disorders.
Here we sought to determine the serum metabolic profiles associated with different psychotic disorders, clustered into three main categories: schizophrenia, affective psychoses, and other non-affective psychoses (ONAP). A metabolomics approach with broad analytical coverage was applied to serum samples from a well characterized population cohort . We investigated dependencies of the three different diagnostic groups on specific metabolic profiles in the context of metabolic comorbidity, antipsychotic medication as well as other lifestyle variables.
Materials and methods
The subjects are from the Health 2000 survey, which is based on a nationally representative sample of 8,028 people aged 30 years or over from Finland . A two-stage stratified cluster sampling procedure was used. The field work took place between September 2000 and June 2001, and consisted of a home interview and a health examination at the local health center, or a condensed interview and health examination of non-respondents at home. In addition, register information was gathered on the whole sample. The Health 2000 study and the accompanying Psychoses in Finland study were approved by the Ethics Committees of the National Public Health Institute (since 2009 the National Institute for Health and Welfare) and the Hospital District of Helsinki and Uusimaa, and participants gave written informed consent . The response rate in the survey, 93%, was exceptionally high compared with other recent surveys.
In the Psychoses in Finland study, we screened people with possible psychotic disorders from the Health 2000 study sample and interviewed them using the Research Version of the Structured Clinical Interview for DSM-IV (SCID-I) . People were invited to participate in the SCID interview if they reported having been diagnosed with a psychotic disorder, received a diagnosis of a possible or definite psychotic disorder from the physician conducting the health examination, or reported possible psychotic or manic symptoms in the Composite International Diagnostic Interview  conducted as part of the health examination. A register-based screen was also used, including hospital treatment for a diagnosis of any psychotic disorder, reimbursement for antipsychotic medication, receipt of a disability pension because of a psychotic disorder, or use of mood-stabilizing medication without a diagnosis of any relevant medical condition, such as epilepsy .
Of the screen-positive people, 63.4% participated in the SCID interview. We diagnosed those who did not participate in the interview using hospital and outpatient case notes from psychiatric and primary care units. Case notes for those who participated in the interview were also collected. Final DSM-IV-based diagnoses were made by JS, JP, and SIS using all available information. Kappa values between the raters ranged from 0.74 to 0.97 for different psychotic disorders .
In this study, lifetime diagnoses of psychotic disorders were grouped into schizophrenia, ONAPs (schizophreniform disorder, schizoaffective disorder, delusional disorder, brief psychotic disorder, psychotic disorder not otherwise specified), and affective psychosis (major depressive disorder with psychotic features and bipolar I disorder). The final study sample comprised 45 subjects with schizophrenia (19 men), 57 with ONAP (20 men), and 37 with affective psychosis (23 men) for whom serum samples were available. There were more women than men in the schizophrenia and ONAP groups, which reflects the gender distribution in the Finnish general population aged 30 years and over and the higher prevalence of schizoaffective disorder in women than in men . An equal number of controls, matched for age, sex, and region of residence, was selected for each group (Table 1). Most of the antipsychotics used by patients were first-generation antipsychotics (Table 1). A total of 12 subjects in the sample used second-generation antipsychotics, of whom 7 used risperidone, 4 clozapine, and one olanzapine. There were 54 subjects who used first-generation antipsychotics, of which the most commonly used were perphenazine (22 users) and thioridazine (16 users).
Participants were asked to fast a minimum of 4 hours before the examination. Subjects with antidiabetic medication were allowed to take their medication and meals at the time they would usually take them (the number of such subjects was three in the schizophrenia group and two in their controls, six in the ONAP group and two in their controls, none in the affective psychosis group and one in their controls). Blood samples were taken at the beginning of the health examination or home health examination. Serum samples were separated, aliquoted and subsequently stored at -70°C (-94°F).
Total, high-density lipoprotein (HDL), and low-density lipoprotein (LDL) cholesterol, triglycerides and glucose were measured with an AU400 analyzer (Olympus, Japan). The inter-assay coefficient of variation for glucose (Olympus System reagent, O'Callaghan's Mills, Co. Clare, Ireland), triglycerides (Olympus System reagent), total cholesterol (Olympus System reagent), HDL cholesterol (HDL-C Plus, Roche Diagnostics, Mannheim, Germany), and LDL cholesterol (LDL-C Plus, Roche Diagnostics) was 2.3%, 3.2%, 2.2%, 5.3%, and 5.7%, respectively. Serum insulin concentrations were determined with an IMx analyzer (Abbott Laboratories, Abbott Park, IL, USA) by microparticle enzyme immunoassay. C-reactive protein (CRP) was determined using an ultra-sensitive immunoturbidometric test (Orion Diagnostica, Espoo, Finland) on an Optima analyzer (Thermo Electron Corporation, Vantaa, Finland). The inter-assay coefficient of variation of both insulin and CRP assays was 4.5%. The cotinine concentration was determined from serum using a radioimmunoassay methodology (Nicotine Metabolite Double Antibody kit, Diagnostic Products Corporation, Los Angeles, CA, USA). The inter-assay coefficient of variation was 12.3%.
Blood pressure was measured after a 5-minute rest twice from the right upper arm with the person sitting. Values used here are average values from the measurements. Weight was measured during bioimpedance measurement. Waist circumference was measured while standing, midway between the lowest rib and the iliac crest, after a modest expiration .
Type 2 diabetes was diagnosed according to the World Health Organization 1999 criteria , combining information from several sources: self-reported diagnosis of type 2 diabetes that was further confirmed in the clinical examination; antidiabetic medication use based on self-report or health care registers; or fasting plasma glucose ≥126 mg/dl (7.0 mmol/l) or nonfasting glucose ≥200 mg/dl (11.1 mmol/l) . Metabolic syndrome was diagnosed using the National Cholesterol Education Program's Adult Treatment Panel III (ATPIII) criteria .
The quantity of alcohol consumption was investigated by asking the respondents to report their average weekly consumption during the past month, separately for each type of alcoholic beverage. The answers were converted into grams of alcohol per week. Daily smoking was self-reported and was defined as having smoked at least 100 cigarettes, having smoked for at least 1 year, and having smoked during the day of the interview or the day before. Standard, validated diet-related questions were used to assess the habitual use of vegetable oils versus butter, use of and fat content in milk products, and daily use of raw vegetables .
The Beck Depression Inventory (BDI-21)  was used to assess current depressive symptoms.
Lipidomic analysis by ultra-performance liquid chromatography coupled to mass spectrometry
EDTA-blood samples (10 ml) were centrifuged at 3,200 rpm (1600 G) for 15 minutes at room temperature within 2 hours of blood sampling. Serum was separated and stored at -80°C. For lipidomics profiling, 10 μl aliquots of serum were used. The samples were mixed with 10 μl of 0.9% (0.15 M) sodium chloride in Eppendorf tubes, spiked with a standard mixture consisting of 10 lipids (0.2 μg/sample; PC(17:0/0:0), PC(17:0/17:0), PE(17:0/17:0), PG(17:0/17:0), Cer(d18:1/17:0), PS(17:0/17:0), PA(17:0/17:0), MG(17:0/0:0/0:0)[rac], DG(17:0/17:0/0:0)[rac], TG(17:0/17:0/17:0), where PC is phosphatidylcholine, PE is phosphatidylethanolamine, PG is phosphatidylglycerol, Cer is ceramide, PS is phosphatidylserine, PA is phosphatidic acid, MG is monoglyceride, DG is diglyceride, and TG is triglyceride) and extracted with 100 μl of chloroform/methanol (2:1). After vortexing (2 minutes) and standing (1 hour) the tubes were centrifuged at 10,000 rpm (7826 G) for 3 minutes and 60 μl of the lower organic phase was separated and spiked with a standard mixture containing three labeled lipids (0.1 μg/sample; PC(16:0/0:0-D3), PC(16:0/16:0-D6), TG(16:0/16:0/16:0-13C3)).
Lipid extracts were analyzed in a randomized order on a Waters Q-Tof Premier mass spectrometer combined with an Acquity UltraPerformance LC™ system (UPLC) (Waters Corporation, Milford, MA, USA). The column (at 50°C) was an Acquity UPLC™ BEH C18 1 × 50 mm with 1.7 μm particles. The solvent system included ultrapure water (1% 1 M NH4Ac, 0.1% HCOOH) and liquid chromatography/mass spectrometry (MS) grade acetonitrile/isopropanol (5:2, 1% 1 M NH4Ac, 0.1% HCOOH). The gradient started from 65% A/35% B, reached 100% B in 6 minutes and remained there for the next 7 minutes. There was a 5-minute re-equilibration step before the next run. The flow rate was 0.200 ml/minute and the injected amount 1.0 μl (Acquity Sample Organizer). Reserpine was used as the lock spray reference compound. The lipid profiling was carried out using ESI+ mode and the data were collected at a mass range of m/z 300 to 1,200 with a scan duration of 0.2 s.
Metabolomic analysis by two-dimensional gas chromatography coupled to time-of-flight MS
Each serum sample (30 μl) was spiked with an internal standard (7 μl 258 ppm labeled palmitic acid) and the mixture was then extracted with 400 μl of methanol. Labeled d-valine (10 μl, 37 ppm) was added to the extracts as a derivatization standard. After centrifugation the supernatant was evaporated to dryness and the original metabolites were then converted into their trimethylsilyl (TMS) and methoxime derivative(s) by two-step derivatization. First, 25 μl methoxyamine hydrochloride (MOX) reagent was added to the residue and the mixture was incubated for 60 minutes at 45°C. Next, 25 μl N-methy-N-(trimethylsilyl) trifluoroacetamide was added and the mixture was incubated for 60 minutes at 45°C. The derivatized samples were diluted 1:1 with hexane. Finally, a retention index standard mixture (n-alkanes) and an injection standard (4,4'-dibromooctafluorobiphenyl), both in pyridine, were added to the mixture.
For the analysis, a Leco Pegasus 4D GC × GC-TOFMS (two-dimensional gas chromatography coupled to time-of-flight MS) instrument (Leco Corp., St Joseph, MI, USA) equipped with a cryogenic modulator was used. The GC part of the instrument was an Agilent 6890N gas chromatograph (Agilent Technologies, Palo Alto, CA, USA) equipped with a split/splitless injector. For the injection, a pulsed splitless injection (0.5 μl) at 240°C was used, with pulse pressure of 55 psig for 1 minute. The first-dimension chromatographic column was a 10-m RTX-5 capillary column with an internal diameter of 0.18 mm and a stationary-phase film thickness of 0.20 μm, and the second-dimension chromatographic column was a 1.5-m BPX-50 capillary column with an internal diameter of 100 μm and a film thickness of 0.1 μm. A diphenyltetramethyldisilyl deactivated retention gap (3 m × 0.53 mm internal diameter) was used in the front of the first column. High-purity helium was used as the carrier gas at a constant pressure mode (39.6 psig). A 5-s separation time was used in the second dimension. The MS spectra was measured at 45 to 700 amu with 100 spectra per second. Pulsed splitless injection 0.5 μl at 240°C was used. The temperature program was as follows: the first-dimension column oven ramp began at 40°C with a 2-minute hold, after which the temperature was programmed to 295°C at a rate of 7°C/minute and then held at this temperature for 3 minutes; the second-dimension column temperature was maintained 20°C higher than the corresponding first-dimension column. The programming rate and hold times were the same for both columns.
The data were scaled into zero mean and unit variance to obtain metabolite profiles comparable to each other. Bayesian model-based clustering was applied on the scaled data to group lipids with similar profiles across all samples. The analyses were performed using the MCLUST  method, implemented in R  as package 'mclust'. In MCLUST the observed data are viewed as a mixture of several clusters and each cluster comes from a unique probability density function. The number of clusters in the mixture, together with the cluster-specific parameters that constrain the probability distributions, will define a model that can then be compared to others. The clustering process selects the optimal model and determines the data partition accordingly. The number of clusters ranging from 4 to 15 and all available model families were considered in our study. Models were compared using the Bayesian information criterion, which is an approximation of the marginal likelihood. The best model is the one that gives the largest marginal likelihood of data, that is, the highest Bayesian information criterion value.
Descriptive statistical analyses and linear mixed models
Differences between each diagnostic group and their matched controls in metabolic comorbidity, lifestyle-related factors, mood, and glucose and lipid measurements were compared using the χ2 test for categorical variables and Mann-Whitney U test for continuous variables. One-way analysis of variance (ANOVA), implemented in Matlab (MathWorks, Natick, MA, USA), was applied to compare the average metabolite profiles in each metabolite cluster. Individual metabolite levels were visualized using the beanplots , implemented in the 'beanplot' R package . Beanplot provides information on the mean metabolite level within each group, the density of the data-point distribution, as well as shows individual data points. The independent effects of diagnostic categories, current antipsychotic medication use, metabolic comorbidity (that is, type 2 diabetes, metabolic syndrome, and obesity (body mass index ≥30)), diet (use of vegetable oil versus butter, use of milk and cheese with high fat content, daily use of vegetables), and duration of fasting were analyzed using linear mixed models  that took the matching of case-control pairs into account. Because the matching was based on both sex and age, these were not included in the models as independent variables. This analysis was performed using PROC MIXED in SAS statistical software, version 9.1.3 (Cary, NC, USA). Logarithm transformations were applied to the metabolomics cluster values to improve normality.
Partial correlation network analysis
Construction of the dependency network for selected variables was performed using undirected Gaussian graphical Markov networks that represent q-order partial correlations between variables, implemented in the R package 'qpgraph'  from the Bioconductor project . In these networks missing edges denote zero partial correlations between pairs of variables, and thus imply the conditional independence relationships in the Gaussian case.
Structure learning of the Gaussian graphs corresponds to a statistical test such as t-test for the hypothesis that a given q-order partial correlation is zero. If all of such hypotheses of zero q-order partial correlations are rejected, then the two variables are joined by an edge. In practice, we tested the hypothesis by default with four equidistant q-values along the (1, 52) interval, namely q = 1, 13, 26 and 38. For each of the q-values, the test was repeated for each pair of variables by sampling 500 elements randomly selected from the subsets of the data that contain q variables. A missing edge is identified if the proportion of such tests where the null-hypothesis is not rejected - for example, the average non-rejection rate of the hypothesis - is above a certain threshold. A small average non-rejection rate therefore implies strong evidence of dependence. The resulting graph can thus be obtained by removing all the missing edges from the complete graph. Unlike Pearson correlation coefficients, use of partial correlation adjusts for the confounding effects and thus removes spurious associations to a large extent. The network was visualized using Cytoscape  and yED graphical editor .
A logistic regression model implemented in R was applied to discriminate the 45 schizophrenia patients from the 94 other participants diagnosed with psychoses using four selected metabolic markers. In order to assess the best marker combination, 10,000 cross-validation runs were performed. In each run, 93 and 46 samples were selected at random as the training and test sets, respectively, and the best marker combination in the logistic regression model was selected using a stepwise algorithm using Akaike's information criterion . The best model was then applied to the test set samples to calculate their predicted classes. The optimal marker combinations in each of the cross-validation runs, receiver operating characteristic (ROC) curves with area under the curve (AUC) statistics, odds-ratios and relative risks were recorded.
Two analytical platforms for metabolomics were applied to all serum samples: a global lipidomics platform based on UPLC-MS, which covers molecular lipids such as phospholipids, sphingolipids, and neutral lipids; and a platform for small polar metabolites based on GC × GC-TOFMS covers small molecules such as amino acids, free fatty acids, keto-acids, various other organic acids, sterols, and sugars. Both platforms were recently described and applied in a large prospective study in type 1 diabetes . The final dataset from each platform consisted of a list of metabolite peaks (identified or unidentified) and their concentrations, calculated using the platform-specific methods, across all samples. All metabolite peaks were included in the data analyses, including the unidentified ones. We reasoned that inclusion of complete data as obtained from the platform best represents the global metabolome, and the unidentified peaks may still be followed-up later on with de novo identification using additional experiments if deemed of interest.
Associations of global metabolome with psychotic disorders
A total of 360 molecular lipids and 201 metabolites were measured, of which 170 and 155 were identified, respectively. Due to a high degree of co-regulation among the metabolites , one cannot assume that all the 562 measured metabolites are independent. The global metabolome was therefore first surveyed by clustering the data into a subset of clusters using the Bayesian model-based clustering . Lipidomic platform data were decomposed into 13 clusters (LC1 to LC13) and the metabolomic data into 8 clusters (MC1 to MC8). Descriptions of each cluster and representative metabolites are provided in Table 2. As expected, the division of clusters to a large extent follows different metabolite functional or structural groups.
As shown in Figure 1, several of the clusters had different average metabolite profiles across the four diagnostic groups, with the control groups pooled into one in this part of the analysis. The average profiles of the lipid clusters LC4 to LC9, which predominantly contained TGs, were most elevated in the schizophrenia group, although the ONAP and affective psychosis groups also tended to have higher TGs compared to controls. The differences were most pronounced for TGs containing more saturated fatty acids, while the cluster containing TGs with polyunsaturated fatty acids (LC10) did not differ between the groups. Two small-molecule clusters were upregulated in schizophrenia, MC3 and MC5, containing branched chain amino acids (BCAAs) and other amino acids, including proline, phenylalanine and glutamic acid. A cluster containing various sugar molecules, MC1, displayed a similar pattern to those of MC3 and MC5, but at a marginal significance level. Cluster MC2, which contained ketone bodies, keto-acids as well as specific free fatty acids, had a distinct pattern that separated the (high level) ONAP and the (low level) affective psychosis groups.
Metabolic comorbidity, antipsychotic medication use, and other lifestyle
It is known that psychoses are associated with metabolic comorbidities  and that the lipid profiles as measured by lipidomics in schizophrenic patients are greatly affected by the use of specific antipsychotic medication . In order to assess the disease-specificity of the observed metabolic changes, the linear mixed effects models were applied on individual metabolite clusters, which included the three diagnostic categories, metabolic comorbidity, current antipsychotic medication, and diet as well as fasting time as explanatory variables (Table 2).
The TG-containing lipid clusters (LC4 to LC10) all associated with metabolic comorbidity, but most of them were also independently and positively associated with schizophrenia. Diet-related factors also affected most of them. Surprisingly, none of the lipid clusters associated with antipsychotic medication use after taking diagnoses, metabolic comorbidity and diet into account. Metabolite cluster MC5 was positively associated with both schizophrenia and metabolic comorbidity, while one (MC3) was associated only with schizophrenia. The only cluster associated with psychoses other than schizophrenia was MC2, which was negatively associated with schizophrenia and affective psychosis. One cluster, MC4, containing various organic acids and energy metabolites, was specifically negatively associated with antipsychotic use.
The observed associations of lipid and metabolic clusters with schizophrenia remained significant in most clusters if patients diagnosed with type 2 diabetes and their controls were excluded from the analysis (Additional file 1).
The linear mixed model analysis suggests that the dependencies of different metabolite classes and related metabolic phenotypes among themselves and with the specific diagnostic groups are likely complex. We hypothesized that a network approach may help elucidate these dependencies to a greater depth. In addition to diagnostic groups, which included also type 2 diabetes (non-insulin-dependent diabetes mellitus (NIDDM)) and the metabolite clusters, we selected 27 other environmental and phenotypic variables related to antipsychotic medication use, diet and lifestyle, metabolic phenotypes (for example, body mass index, insulin, glucose, HDL-cholesterol, total TG), and other biochemical measures, such as CRP and gamma-glutamyltransferase (GGT). The undirected Gaussian graphical Markov model was applied to estimate partial correlations between the variables (Figure 2).
In addition to variables related to antipsychotic use, schizophrenia was associated with two metabolic variables, lipid cluster LC9 and fasting serum insulin (Insulin in Figure 2). Insulin was further associated with related metabolic variables such as homeostatic model assessment (HOMA in Figure 2) index and glucose, while LC9 was associated with other TG-containing clusters as well as with total triglycerides. Both insulin and LC9 were associated with metabolite cluster MC5, which was directly linked to MC3. Neither the ONAP nor the affective psychosis group was directly associated with the specific metabolic clusters. ONAP was associated with the inflammatory marker CRP and with depressive symptoms. Affective psychosis was directly associated with the liver marker gamma-glutamyltransferase, which not surprisingly was associated with alcohol use.
Feasibility of metabolic profile in assisting schizophrenia diagnosis
We reasoned that due to their independent association with schizophrenia, insulin as well as specific other metabolite clusters reflect the disease process itself, and may thus help discriminate schizophrenia from other psychoses. To assess the feasibility of diagnosis, we selected insulin as well as the top-ranking metabolites from three clusters of most interest based on the network structure in Figure 2: triglyceride TG(18:1/18:0/18:1) (LC9), isoleucine (MC3), and proline (MC5). Only the three psychotic groups were included in the analysis, without the controls, and the comparisons were made between the schizophrenia versus the pooled ONAP and affective psychosis groups.
The best model derived from logistic regression analysis was obtained by combining proline and TG(18:1/18:0/18:1). This combination was selected in 53% of 10,000 cross-validation runs. Other strongly performing models were proline alone (25%) and combined insulin and proline (13%). Figure 3 shows the summary of the combined proline and TG(18:1/18:0/18:1) diagnostic model, based on independently tested data taken from 2,000 samplings.
Our findings, based on a highly phenotypically detailed general population sample of different psychoses, independently associate specific metabolic phenotypes, as measured by metabolomics, with schizophrenia. It is known that schizophrenia is associated with elevated fasting total triglycerides and insulin resistance , but this metabolic abnormality has usually been attributed to antipsychotic drug-specific side effects . The strongest association with schizophrenia based on network analysis as well as linear mixed models was with the lipid cluster LC9, which contains saturated and longer chain triglycerides. In a recent lipidomic study of different lipoprotein fractions in subjects with varying degrees of insulin resistance, we found that the lipids found in LC9 are abundant in liver-produced very low density lipoprotein particles and are associated with insulin resistance . In agreement with this, schizophrenia patients in the present study were insulin resistant and had elevated fasting serum insulin levels. Together, our data indicate that schizophrenia, independent of antipsychotic medication and metabolic comorbidity, is characterized by insulin resistance, and consequently enhanced hepatic very low density lipoprotein production  and thus elevated serum concentrations of specific triglycerides. This is consistent with findings from an earlier study that demonstrated that antipsychotic medication-naïve patients with schizophrenia display hepatic insulin resistance independent of intra-abdominal fat mass or other known factors associated with hepatic insulin resistance .
The possible pathogenic relevance of our findings is supported by recent studies showing that abnormal insulin secretion and response [44–47] and abnormal glucose tolerance and risk of diabetes  are found already in drug-naïve first-episode patients with schizophrenia. In line with this, the insulinotropic [49, 50] BCAAs from the metabolic cluster MC3 were also elevated and specifically associated with schizophrenia. In the context of psychoses, BCAAs are not only important due to their role in stimulating insulin secretion, but also since they compete with aromatic amino acids for transport across the blood-brain barrier . Their increase may thus lead to concentration decreases of neurotransmitters derived from the aromatic amino acids in the brain, specifically catecholamines from tyrosine and phenylalanine (MC3) and serotonin from tryptophan (MC5). However, the effect of BCAA-induced dopamine or serotonin depletion in the brain on schizophrenia-related cognitive performance is currently controversial [52, 53]. Another potential mechanism linking schizophrenia and long-term hyperinsulinemia is dysregulation of insulin-receptor-mediated signaling, which has a role in learning and memory as well as in regionally specific glucose metabolism in the brain .
The metabolic cluster MC5, which included proline and glutamate, was strongly associated with schizophrenia. Glutamate has been hypothesized to play an important role in schizophrenia . Our data show that serum glutamate is elevated in all psychoses compared to controls (Figure 1), supporting the view that glutamate-related metabolic abnormalities may reflect a common pathway across different psychoses . However, one should also note that the dependency of glutamate concentrations between the brain and blood is weak and complex due to restricted and tightly controlled passage of glutamate across the blood-brain barrier .
Upregulation of serum proline was specific to schizophrenia. There is evidence from genetics that polymorphisms in the PRODH gene, encoding proline oxidase, which is located at 22q11, are associated with schizophrenia risk [58, 59] and that the related hyperprolinemia negatively associates with cognitive performance . In particular, functional variants in the PRODH gene that result in reduction of proline oxidase activity and hyperprolinemia are associated with increased risk of schizophrenia and changes in fronto-striatal structure and function [59, 61]. Interestingly, schizophrenia is linked to the same copy number variants spanning the 22q11 region including PRODH as autism and other childhood developmental disorders, whereas bipolar disorder is not . Furthermore, recent functional studies suggest that microdeletions on human chromosome 22 (22q11.2) lead to impaired long-range synchrony of neuronal activity and may thus be an important component of the pathophysiology of schizophrenia .
Having a population-based sample with carefully matched controls was a definite strength of the study. The Psychoses in Finland study has been characterized as 'arguably the most thorough study ever undertaken on the prevalence of psychotic disorders' . In addition to the careful screening and assessment of psychotic disorders, the assessment of health and lifestyle in the Health 2000 survey was comprehensive. Diabetes and metabolic syndrome had been carefully diagnosed [24, 39] and their effects had been controlled for in the analyses. Notably, most antipsychotics used by patients were first-generation antipsychotics, which are less associated with diabetes compared to second-generation antipsychotics . However, the sample was relatively old and the mean duration of illness among subjects with psychotic disorders had been long. Although we controlled for the effects of current lifestyle, all the long-term effects of antipsychotic medication and lifestyle-related factors, like smoking, nutrition and exercise, may not have been captured. Nevertheless, studies on drug-naïve first-episode patients already find impaired glucose tolerance, elevated insulin and metabolic abnormalities [43–48] that are not related to poor health habits . Longitudinal research in prodromal and early psychosis is needed to further elucidate the role of the identified metabolomic changes in psychotic disorders.
Our study suggests that proline-related metabolic abnormalities and insulin secretion-related changes (BCAAs, insulin, triglycerides) reflect two different disease-related pathways. This is further supported by the fact that the best candidate diagnostic model separating schizophrenia from other psychoses is obtained by combining the selected metabolites from each of the two pathways. We believe metabolomics, which is sensitive to both genetic and environmental variation, will be a powerful tool to further investigate susceptibility to psychotic disorders, their clinical course, and treatment responses.
branched chain amino acid
Beck Depression Inventory
Diagnostic and Statistical Manual of Mental Disorders
- GC × GC-TOFMS:
two-dimensional gas chromatography coupled to time-of-flight mass spectrometry
other non-affective psychoses
phosphatidylcholine, PE: phosphatidylethanolamine
Structured Clinical Interview for DSM-IV
ultra-performance liquid chromatography.
Üstün TB, Rehm J, Chatterji S, Saxena S, Trotter R, Room R, Bickenbach J: Multiple-informant ranking of the disabling effects of different health conditions in 14 countries. Lancet. 1999, 354: 111-115.
Perälä J, Suvisaari J, Saarni SI, Kuoppasalmi K, Isometsä E, Pirkola S, Partonen T, Tuulio-Henriksson A, Hintikka J, Kieseppa T, Harkanen T, Koskinen S, Lonnqvist J: Lifetime prevalence of psychotic and bipolar I disorders in a general population. Arch Gen Psychiatry. 2007, 64: 19-28.
Murray RM, Lappin J, Di Forti M: Schizophrenia: from developmental deviance to dopamine dysregulation. Eur Neuropsychopharmacol. 2008, 18: S129-S134. 10.1016/j.euroneuro.2008.04.002.
Sullivan PF, Kendler KS, Neale MC: Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch Gen Psychiatry. 2003, 60: 1187-1192. 10.1001/archpsyc.60.12.1187.
Lichtenstein P, Yip BH, Björk C, Pawitan Y, Cannon TD, Sullivan PF, Hultman CM: Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet. 2009, 373: 234-239. 10.1016/S0140-6736(09)60072-6.
Goldstein JM, Buka SL, Seidman LJ, Tsuang MT: Specificity of familial transmission of schizophrenia psychosis spectrum and affective psychoses in the New England family study's high-risk design. Arch Gen Psychiatry. 2010, 67: 458-467. 10.1001/archgenpsychiatry.2010.38.
van Os J, Kapur S: Schizophrenia. Lancet. 2009, 374: 635-645. 10.1016/S0140-6736(09)60995-8.
Illig T, Gieger C, Zhai G, Romisch-Margl W, Wang-Sattler R, Prehn C, Altmaier E, Kastenmuller G, Kato BS, Mewes HW, Meitinger T, de Angelis MH, Kronenberg F, Soranzo N, Wichmann HE, Spector TD, Adamski J, Suhre K: A genome-wide perspective of genetic variation in human metabolism. Nat Genet. 2010, 42: 137-141. 10.1038/ng.507.
Lenz EM, Bright J, Wilson ID, Hughes A, Morrisson J, Lindberg H, Lockton A: Metabonomics, dietary influences and cultural differences: a 1H NMR-based study of urine samples obtained from healthy British and Swedish subjects. J Pharm Biomed Anal. 2004, 36: 841-849. 10.1016/j.jpba.2004.08.002.
Nikkilä J, Sysi-Aho M, Ermolov A, Seppänen-Laakso T, Simell O, Kaski S, Oresic M: Gender dependent progression of systemic metabolic states in early childhood. Mol Syst Biol. 2008, 4: e197-
Velagapudi VR, Hezaveh R, Reigstad CS, Gopalacharyulu PV, Yetukuri L, Islam S, Felin J, Perkins R, Boren J, Oresic M, Backhed F: The gut microbiota modulates host energy and lipid metabolism in mice. J Lipid Res. 2010, 51: 1101-1112. 10.1194/jlr.M002774.
Oresic M, Simell S, Sysi-Aho M, Näntö-Salonen K, Seppänen-Laakso T, Parikka V, Katajamaa M, Hekkala A, Mattila I, Keskinen P, Yetukuri L, Reinikainen A, Lähde J, Suortti T, Hakalax J, Simell T, Hyöty H, Veijola R, Ilonen J, Lahesmaa R, Knip M, Simell O: Dysregulation of lipid and amino acid metabolism precedes islet autoimmunity in children who later progress to type 1 diabetes. J Exp Med. 2008, 205: 2975-2984. 10.1084/jem.20081800.
van der Greef J, Stroobant P, Heijden Rvd: The role of analytical sciences in medical systems biology. Curr Opin Chem Biol. 2004, 8: 559-565. 10.1016/j.cbpa.2004.08.013.
Oresic M: Metabolomics, a novel tool for studies of nutrition, metabolism and lipid dysfunction. Nutr Metab Cardiovasc Dis. 2009, 19: 816-824. 10.1016/j.numecd.2009.04.018.
Holmes E, Tsang TM, Huang JT-J, Leweke FM, Koethe D, Gerth CW, Nolden BM, Gross S, Schreiber D, Nicholson JK, Bahn S: Metabolic profiling of CSF: evidence that early intervention may impact on disease progression and outcome in schizophrenia. PLoS Med. 2006, 3: e327-10.1371/journal.pmed.0030327.
Schwarz E, Prabakaran S, Whitfield P, Major H, Leweke FM, Koethe D, McKenna P, Bahn S: High throughput lipidomic profiling of schizophrenia and bipolar disorder brain tissue reveals alterations of free fatty acids, phosphatidylcholines, and ceramides. J Proteome Res. 2008, 7: 4266-4277. 10.1021/pr800188y.
Yao JK, Dougherty GG, Reddy RD, Keshavan MS, Montrose DM, Matson WR, Rozen S, Krishnan RR, McEvoy J, Kaddurah-Daouk R: Altered interactions of tryptophan metabolites in first-episode neuroleptic-naive patients with schizophrenia. Mol Psychiatry. 2009, 15: 938-953. 10.1038/mp.2009.33.
Kaddurah-Daouk R, McEvoy J, Baillie RA, Lee D, Yao JK, Doraiswamy PM, Krishnan KRR: Metabolomic mapping of atypical antipsychotic effects in schizophrenia. Mol Psychiatry. 2007, 12: 934-945. 10.1038/sj.mp.4002000.
Health and functional capacity in Finland: baseline results of the Health 2000 Health Examination Survey (Publication B12/2004). [http://www.ktl.fi/terveys2000/index.uk.html]
First MB, Anthony JC, Tepper S, Dryman A: Structured Clinical Interview for DSM-IV Axis I Disorders, Research Version, Nonpatient Edition (SCID-I/NP). 1997, New York: Biometrics Research, New York State Psychiatric Institute
Wittchen HU, Pfister H: DIA-X-Interviews: Manual für screening-verfahren und Interview; Interviewheft Längsschnittuntersuchung (DIA-X-Lifetime); Ergänzungsheft (DIAX-Lifetime); Interviewheft Querschnittuntersuchung (DIA-X-12 Monate); Ergänzungsheft (DIA-X-12 Monate); PC-Programm zur Durchführung des Interviews (Längs- und Querschnittuntersuchung). 1997, Frankfurt, Germany: Swets and Zeitlinger
Methodology report: Health 2000 Survey. [http://www.terveys2000.fi/doc/methodologyrep.pdf]
Definition, diagnosis and classification of diabetes mellitus and its complications (WHO/NCD/NCS/99.2). [http://whqlibdoc.who.int/hq/1999/WHO_NCD_NCS_99.2.pdf]
Suvisaari J, Perälä J, Saarni S, Härkänen T, Pirkola S, Joukamaa M, Koskinen S, Lönnqvist J, Reunanen A: Type 2 diabetes among persons with schizophrenia and other psychotic disorders in a general population survey. Eur Arch Psychiatry Clin Neurosci. 2008, 258: 129-136. 10.1007/s00406-007-0762-y.
Expert Panel on Detection, Evaluation and Treatment of High Blood Cholesterol in Adults: Executive summary of the third report of the National Cholesterol Education Program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (Adult Treatment Panel III). JAMA. 2001, 285: 2486-2497. 10.1001/jama.285.19.2486.
Beck AT, Ward CH, Mendelson M, Mock J, Erbaugh J: An inventory for measuring depression. Arch Gen Psychiatry. 1961, 4: 561-571.
Pluskal T, Castillo S, Villar-Briones A, Oresic M: MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics. 2010, 11: 395-10.1186/1471-2105-11-395.
Yetukuri L, Katajamaa M, Medina-Gomez G, Seppanen-Laakso T, Vidal-Puig A, Oresic M: Bioinformatics strategies for lipidomics analysis: characterization of obesity related hepatic steatosis. BMC Syst Biol. 2007, 1: 12-10.1186/1752-0509-1-12.
Fraley C, Raftery AE: Model-based methods of classification: Using the mclust software in chemometrics. J Stat Soft. 2007, 18: 1-13.
The R Project for Statistical Computing. [http://www.r-project.org/]
Kampstra P: Beanplot: a boxplot alternative for visual comparison of distributions. J Stat Soft. 2008, 28: 1-9.
McCulloch CE, Searle SR, Neuhaus JM: Generalized, Linear, and Mixed Models. 2008, New York: John Wiley and Sons, 2
Castelo R, Roverato A: Reverse engineering molecular regulatory networks from microarray data with qp-graphs. J Comput Biol. 2009, 16: 213-227. 10.1089/cmb.2008.08TT.
Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, Christmas R, Avila-Campilo I, Creech M, Gross B, Hanspers K, Isserlin R, Kelley R, Killcoyne S, Lotia S, Maere S, Morris J, Ono K, Pavlovic V, Pico AR, Vailaya A, Wang PL, Adler A, Conklin BR, Hood L, Kuiper M, Sander C, Schmulevich I, Schwikowski B, Warner GJ, et al: Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007, 2: 2366-2382. 10.1038/nprot.2007.324.
Brohee S, Faust K, Lima-Mendez G, Vanderstocken G, van Helden J: Network Analysis Tools: from biological networks to clusters and pathways. Nat Protoc. 2008, 3: 1616-1629. 10.1038/nprot.2008.100.
Yamashita T, Yamashita K, Kamimura R: A stepwise AIC method for variable selection in linear regression. Commun Stat Theory Methods. 2007, 36: 2395-2403. 10.1080/03610920701215639.
Steuer R, Kurths J, Fiehn O, Weckwerth W: Observing and interpreting correlations in metabolomic networks. Bioinformatics. 2003, 19: 1019-1026. 10.1093/bioinformatics/btg120.
Suvisaari JM, Saarni SI, Perälä J, Suvisaari JV, Härkänen T, Lönnqvist J, Reunanen A: Metabolic syndrome among persons with schizophrenia and other psychotic disorders in a general population survey. J Clin Psychiatry. 2007, 68: 1045-1055. 10.4088/JCP.v68n0711.
Meyer JM, Davis VG, Goff DC, McEvoy JP, Nasrallah HA, Davis SM, Rosenheck RA, Daumit GL, Hsiao J, Swartz MS, Stroup TS, Lieberman JA: Change in metabolic syndrome parameters with antipsychotic treatment in the CATIE Schizophrenia Trial: Prospective data from phase 1. Schizophr Res. 2008, 101: 273-286. 10.1016/j.schres.2007.12.487.
Kotronen A, Velagapudi VR, Yetukuri L, Westerbacka J, Bergholm R, Ekroos K, Makkonen J, Taskinen MR, Oresic M, Yki-Järvinen H: Saturated fatty acids containing triacylglycerols are better markers of insulin resistance than total serum triacylglycerol concentrations. Diabetologia. 2009, 52: 684-690. 10.1007/s00125-009-1282-2.
Kotronen A, Yki-Jarvinen H: Fatty liver: a novel component of the metabolic syndrome. Arterioscler Thromb Vasc Biol. 2008, 28: 27-38. 10.1161/ATVBAHA.107.147538.
van Nimwegen LJ, Storosum JG, Blumer RM, Allick G, Venema HW, de Haan L, Becker H, van Amelsvoort T, Ackermans MT, Fliers E, Serlie MJ, Sauerwein HP: Hepatic insulin resistance in antipsychotic naive schizophrenic patients: stable isotope studies of glucose metabolism. J Clin Endocrinol Metab. 2008, 93: 572-577. 10.1210/jc.2007-1167.
Fernandez-Egea E, Bernardo M, Donner T, Conget I, Parellada E, Justicia A, Esmatjes E, Garcia-Rizo C, Kirkpatrick B: Metabolic profile of antipsychotic-naive individuals with non-affective psychosis. Br J Psychiatry. 2009, 194: 434-438. 10.1192/bjp.bp.108.052605.
Steiner J, Walter M, Guest P, Myint AM, Schiltz K, Panteli B, Brauner M, Bernstein HG, Gos T, Herberth M, Schroeter ML, Schwarz MJ, Westphal S, Bahn S, Bogerts B: Elevated S100B levels in schizophrenia are associated with insulin resistance. Mol Psychiatry. 2010, 15: 3-4. 10.1038/mp.2009.87.
Guest PC, Wang L, Harris LW, Burling K, Levin Y, Ernst A, Wayland MT, Umrania Y, Herberth M, Koethe D, van Beveren JM, Rothermundt M, McAllister G, Leweke FM, Steiner J, Bahn S: Increased levels of circulating insulin-related peptides in first-onset, antipsychotic naive schizophrenia patients. Mol Psychiatry. 2010, 15: 118-119. 10.1038/mp.2009.81.
Guest PC, Schwarz E, Krishnamurthy D, Harris LW, Leweke FM, Rothermundt M, van Beveren NJ, Spain M, Barnes A, Steiner J, Rahmoune H, Bahn S: Altered levels of circulating insulin and other neuroendocrine hormones associated with the onset of schizophrenia. Psychoneuroendocrinology. 2011
Kirkpatrick B, Miller BJ, Garcia-Rizo C, Fernandez-Egea E, Bernardo M: Is abnormal glucose tolerance in antipsychotic-naive patients with nonaffective psychosis confounded by poor health habits?. Schizophr Bull. 2010.
Newsholme P, Bender K, Kiely A, Brennan L: Amino acid metabolism, insulin secretion and diabetes. Biochem Soc Trans. 2007, 35: 1180-1186. 10.1042/BST0351180.
Newgard CB, An J, Bain JR, Muehlbauer MJ, Stevens RD, Lien LF, Haqq AM, Shah SH, Arlotto M, Slentz CA, Rochon J, Gallup D, Ilkayeva O, Wenner BR, Yancy WS, Eisenson H, Musante G, Surwit RS, Millington DS, Butler MD, Svetkey LP: A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab. 2009, 9: 311-326. 10.1016/j.cmet.2009.02.002.
Fernstrom JD: Branched-chain amino acids and brain function. J Nutr. 2005, 135: 1539S-1546.
Neuhaus AH, Goldberg TE, Hassoun Y, Bates JA, Nassauer KW, Sevy S, Opgen-Rhein C, Malhotra AK: Acute dopamine depletion with branched chain amino acids decreases auditory top-down event-related potentials in healthy subjects. Schizophr Res. 2009, 111: 167-173. 10.1016/j.schres.2009.03.023.
Leung S, Croft R, Guille V, Scholes K, O'Neill B, Phan K, Nathan P: Acute dopamine and/or serotonin depletion does not modulate mismatch negativity (MMN) in healthy human participants. Psychopharmacology. 2010, 208: 233-244. 10.1007/s00213-009-1723-0.
Craft S, Stennis Watson G: Insulin and neurodegenerative disease: shared and specific mechanisms. Lancet Neurol. 2004, 3: 169-178. 10.1016/S1474-4422(04)00681-7.
Coyle J: Glutamate and schizophrenia: beyond the dopamine hypothesis. Cell Mol Neurobiol. 2006, 26: 365-384. 10.1007/s10571-006-9062-8.
Cherlyn SY, Woon PS, Liu JJ, Ong WY, Tsai GC, Sim K: Genetic association studies of glutamate, GABA and related genes in schizophrenia and bipolar disorder: a decade of advance. Neurosci Biobehav Rev. 2010, 34: 958-977. 10.1016/j.neubiorev.2010.01.002.
Smith QR: Transport of glutamate and other amino acids at the blood-brain barrier. J Nutr. 2000, 130: 1016S-1022S.
Liu H, Heath SC, Sobin C, Roos JL, Galke BL, Blundell ML, Lenane M, Robertson B, Wijsman EM, Rapoport JL, Gogos JA, Karayiorgou M: Genetic variation at the 22q11 PRODH2/DGCR6 locus presents an unusual pattern and increases susceptibility to schizophrenia. Proc Natl Acad Sci USA. 2002, 99: 3717-3722. 10.1073/pnas.042700699.
Kempf L, Nicodemus KK, Kolachana B, Vakkalanka R, Verchinski BA, Egan MF, Straub RE, Mattay VA, Callicott JH, Weinberger DR, Meyer-Lindenberg A: Functional polymorphisms in PRODH are associated with risk and protection for schizophrenia and fronto-striatal structure and function. PLoS Genet. 2008, 4: e1000252-10.1371/journal.pgen.1000252.
Raux G, Bumsel E, Hecketsweiler B, van Amelsvoort T, Zinkstok J, Manouvrier-Hanu S, Fantini C, Breviere G-MM, Di Rosa G, Pustorino G, Vogels A, Swillen A, Legallic S, Bou J, Opolczynski G, Drouin-Garraud V, Lemarchand M, Philip N, Gerard-Desplanches A, Carlier M, Philippe A, Nolen MC, Heron D, Sarda P, Lacombe D, Coizet C, Alembik Y, Layet V, Afenjar A, Hannequin D, et al: Involvement of hyperprolinemia in cognitive and psychiatric features of the 22q11 deletion syndrome. Hum Mol Genet. 2007, 16: 83-91. 10.1093/hmg/ddl443.
Bender HU, Almashanu S, Steel G, Hu CA, Lin WW, Willis A, Pulver A, Valle D: Functional consequences of PRODH missense mutations. Am J Hum Genet. 2005, 76: 409-420. 10.1086/428142.
Guilmatre A, Dubourg C, Mosca AL, Legallic S, Goldenberg A, Drouin-Garraud V, Layet V, Rosier A, Briault S, Bonnet-Brilhault F, Laumonnier F, Odent S, Le Vacon G, Joly-Helas G, David V, Bendavid C, Pinoit JM, Henry C, Impallomeni C, Germano E, Tortorella G, Di Rosa G, Barthelemy C, Andres C, Faivre L, Frebourg T, Saugier Veber P, Campion D: Recurrent rearrangements in synaptic and neurodevelopmental genes and shared biologic pathways in schizophrenia, autism, and mental retardation. Arch Gen Psychiatry. 2009, 66: 947-956. 10.1001/archgenpsychiatry.2009.80.
Sigurdsson T, Stark KL, Karayiorgou M, Gogos JA, Gordon JA: Impaired hippocampal-prefrontal synchrony in a genetic mouse model of schizophrenia. Nature. 2010, 464: 763-767. 10.1038/nature08855.
McGrath JJ: The surprisingly rich contours of schizophrenia epidemiology. Arch Gen Psychiatry. 2007, 64: 14-16. 10.1001/archpsyc.64.1.14.
Smith M, Hopkins D, Peveler RC, Holt RIG, Woodward M, Ismail K: First- v. second-generation antipsychotics and risk for diabetes in schizophrenia: systematic review and meta-analysis. Br J Psychiatry. 2008, 192: 406-411. 10.1192/bjp.bp.107.037184.
We thank Ulla Lahtinen, Anna-Liisa Ruskeepää, and Sandra Castillo for their help in metabolomic analysis and data processing. This work was supported in part by the EU-funded projects ETHERPATHS (FP7-KBBE-222639 to MO) and TORNADO (FP7-KBBE-222720 to MO), and by the NARSAD Maltz Investigator Award and the Academy of Finland grant (129434; to JS).
The authors declare that they have no competing interests.
JT, MSA, MO, and JS performed the statistical analysis. TSL, IM, and TH carried out metabolomic analyses. JL and JP participated in the study design. SES, SIS, JP, and JS researched the primary clinical data. MO and JS conceived of the study, and participated in its design and coordination and drafted the manuscript. All authors read and approved the final manuscript.