Malnutrition and altered metabolome in patients with schizophrenia
During May ~ Aug 2018, the study enrolled 100 SZ patients from a sanatorium for patients of mental illness in County Gai, Liaoning, China and 52 age- and gender-matched healthy controls who were staff members of the sanatorium or inhabitants of the same town in order to minimize the influence of food sources and dietary habits in gut microbiota (23). Demographic and clinical information, and dietary records of the recent month before enrollment, were collected in detail for each subject (Table 1). Stool and plasma samples of all subjects were collected upon informative consent and cryopreserved for subsequent shotgun sequencing and metabolome analysis.
Table 1
Demographic, clinical characteristics and dietary patterns of all subjects.
Variable | Schizophrenia (n = 100) | healthy controls (n = 53) | p-value |
---|
Demographic characteristics |
Age, means(SD) | 43.11(9.81) | 48.92(12.95) | 0.005 |
Education year, means (SD) | 7.80(2.79) | 9.34(3.29) | 0.003 |
Height, means(SD) | 1.64(0.07) | 1.64(0.08) | 0.645 |
weight, means(SD) | 60.53(7.86) | 65.04(8.20) | 0.001 |
BMI, means (SD) | 22.59(2.44) | 24.11(2.45) | ༜0.001 |
Female, No.(%) | 34(34) | 21(39.6) | 0.490 |
Married, No.(%) | 39(39) | 48(90.6) | ༜0.001 |
Smoking, No.(%) | 0 | 11(20.8) | ༜0.001 |
Drinking, No.(%) | 0 | 10(18.9) | ༜0.001 |
Clinical characteristics |
Fasting plasma glucos, means (SD) | 4.26(0.65) | 4.40(0.35) | 0.252 |
Triglycerides, means (SD) | 1.06(0.83) | 1.07(0.33) | 0.888 |
Cholesterol, means (SD) | 3.80(0.74) | 3.77(0.85) | 0.844 |
Course of disease, means(SD) | 13.35(8.53) | NA | NA |
Monotherapy, No.(%) | 81(81) | NA | NA |
Equivalent dose of chlorpromazine, means(SD) | 246.00(179.51) | NA | NA |
positive syndrome score, means(SD) | 17.45(5.56) | NA | NA |
negative syndrome score, means(SD) | 24.36(7.26) | NA | NA |
General syndrome score, means(SD) | 39.76(5.88) | NA | NA |
PANSS total score, means(SD) | 81.51(10.95) | NA | NA |
Dietary characteristics, means(SD) |
Energy | 1595.03(106.97) | 1817.25(295.00) | ༜0.001 |
Water | 982.15(32.53) | 996.84(122.75) | 0.396 |
Protein | 58.15(8.16) | 97.57(19.67) | ༜0.001 |
Fatty | 82.31(8.52) | 83.95(26.04) | 0.658 |
Carbohydrate | 155.27(9.74) | 167.76(20.86) | ༜0.001 |
Dietary fiber | 12.84(1.61) | 16.37(2.65) | ༜0.001 |
Retinol equivalent | 788.98(84.04) | 841.69(143.98) | 0.017 |
Vitamin B1 | 0.78(0.06) | 1.13(0.27) | ༜0.001 |
Vitamin B2 | 0.96(0.05) | 1.49(0.42) | ༜0.001 |
Vitamin PP | 20.46(2.98) | 21.32(3.98) | 0.173 |
Vitamin E | 18.56(3.95) | 26.87(8.62) | ༜0.001 |
Na | 1260.04(188.79) | 2118.78(449.25) | ༜0.001 |
Ca | 432.72(31.07) | 543.71(134.25) | ༜0.001 |
Fe | 17.75(5.22) | 25.58(5.32) | ༜0.001 |
Vitamin C | 119.39(5.48) | 114.46(8.46) | ༜0.001 |
Cholesterol | 471.40(39.49) | 805.65(277.21) | ༜0.001 |
We first noticed apparent underweight in SZ patients and a significantly reduced average BMI despite adequate food supply for all patients (Fig. 1A), which observation was in concordance with previous reports from other countries (24, 25). The daily total calorie intake calculated from the dietary records in SZ patients was significantly reduced (Fig. 1B), mainly due to the almost halved protein intake (Fig. 1C) and proportion of calories provided by ingested protein (%P, Fig. 1D). This result suggested a general malnutrition status in SZ patients, especially the inadequate ingestion of food proteins.
Non-targeted metabolomics analysis for plasma and fecal samples were first performed to investigate the influence of the aberrant macronutrients on metabolism. Comparison of the plasma metabolism between cases and controls identified little differences related to protein catabolism, and among the 14 metabolites significantly altered (Wilcoxon rank-sum test, adj. p < 0.05, Fig S1A), only urea and α-hydroxyisobutyric acid (elevated in control samples) were products of protein catabolism. In contrast, the fecal metabolome revealed more extensively altered protein metabolism in SZ patients, and 12 of the 35 significantly altered metabolites (Wilcoxon rank-sum test, adj. p < 0.05, Fig. 1E) were products or related derivates of protein catabolism. The most notable alterations were the elevated levels of the neuroactive dopamine and decrease in the neuroleptic 𝛾-aminobutyric acid in SZ patients (Fig. 1F), both of which are metabolites derived from amino acids and important neurotransmitters in the pathogenesis of SZ.
We then performed metabolome-wide association studies (MWASs) using random forest-based machine learning variable selection techniques to identify fecal metabolite features that deviated SZ. The permutation importance of each metabolite showed that three out of the eight top-ranking metabolites favoring SZ were metabolites derived from amino acids, whereas none metabolites favoring control related to protein catabolism (Fig S1B). The apparent differences in fecal metabolism suggested potentially dysregulated protein metabolism in SZ patients, possibly the results of both aberrant protein ingestion and the intestinal microbiome dysbiosis.
To test the alteration in fatty acid metabolism in previous observation (26), we performed absolute quantification of the plasma and fecal levels of medium- and long-chain free fatty acids using GC-MS (gas chromatography-mass spectrometry). Our results confirmed significant deficiency in the plasma levels of some free fatty acids in SZ patients (Wilcoxon rank-sum test, adj. p < 0.05, Fig S2). However, the deficiency seemed not to stem from the intestinal metabolism as no difference in fatty acids was identified in the stool (Wilcoxon rank-sum test, adj. p < 0.05, Fig S2), given the fact that fats from food are primarily absorbed in the jejunum with little proportion delivered to the intestinal microbiota. Thus, the altered fatty acid profile in plasma is less likely associated with the intestinal microbiome dysbiosis.
Shifted gut microbial fermentation from carbohydrates to proteins in schizophrenia
To investigate how gut microbes participate in the metabolic disturbance in SZ, we performed shotgun-sequencing of the metagenome for all participants. In the analysis of differentially represented species (Wilcoxon rank-sum test, p < 0.05, Table S1), we noticed that 4/22of the species enriched in patients with highest significance and fold-change were asaccharolytic, i.e., Fusobacterium mortiferum, Desulfovibrio piger, Phascolarctobacterium succinatutens, and Sutterella wadsworthensis, whereas no asaccharolytic species enriched in controls. In the microbiome-wide association studies (MWASs) using random forest model to select taxonomic features of SZ, and the top three species of highest permutation importance in favor of SZ were all asaccharolytic (Fig S3A). In contrast, no asaccharolytic species were in the list of top permutation importance in favor of controls (Fig S3A). Furthermore, we combed all asaccharolytic species represented in our samples and found that most of them, especially those of higher abundance, were enriched in SZ patients, leading to the total abundance of asaccharolytic species significantly elevated in patients (Fig. 2A).
Comparison in metabolic pathways between cases and controls discovered that the total abundance of carbohydrate catabolism pathways, as well as two major ones, i.e., starch degradation and anaerobic glycolysis, were significantly reduced in patients (Fig S3B). The reduction, when normalized by the proportion of calories provided by ingested carbohydrates in daily calorie intake (%C), was more remarkable (Fig. 2B), indicating the carbohydrate catabolism was hypoactive in SZ. The hypoactivity in microbial carbohydrate metabolism, together with the enrichment of asaccharolytic species in patients, implies reduced carbohydrate supply to intestinal microbes and hypoactive carbohydrate fermentation by them in the intestine of SZ patients.
Further, according to the MEROPS (database of proteolytic enzymes) and CAZymes (Carbohydrate-active enzymes database), we annotated all genes of peptidases and carbohydrate-active enzymes (CAZy), which account for the hydrolyzation of protein and carbohydrates into monomers. The initial comparison of their abundance between patients and controls resulted in no significant difference. However, when considering the great difference between the two groups in the ratio of protein to carbohydrate intake (Fig. 2C), the ratio between peptidase and CAZys was significantly reversed (Fig. 2D). Additionally, when normalized by the proportion of calories provided by ingested protein (%P) or carbohydrate (%C) in daily caloric intake, the abundance of peptidases was significantly increased in SZ (Fig. 2E), and reduced in CAZy (Fig. 2F).
To test whether more proteins were hydrolyzed in patients’ gut, we performed targeted chromatographic assay to quantify the concentrations of all amino acids in the stool. The result exhibited that most amino acids’ fecal concentrations increased in patients, and the total concentration of amino acids was significantly higher than controls, which were more significant when normalized by daily protein intake (Fig. 2G). Notably, we observed significantly elevated fecal concentrations of phenylalanine, tyrosine, and glutamine, which are neurotoxic at high levels in the brain (27) and neural-effective on enteric nerves (28, 29). All the above observations suggested that, in SZ patients, more undigested proteins reach the colon and are hydrolyzed by microbes there, instead of mostly being hydrolyzed and absorbed in the small intestine in normal conditions.
Enhanced microbial amino acid catabolism instead of protein biosynthesis in schizophrenia
Microbes are powerful in catabolizing amino acids into a great variety of derivatives (30, 31), many of which, including amines, NO, indole, kynurenine, quinolone, and so on, are neurologically active or affect human behavior through actions on the immune or endocrine system (32, 33). Key enzymes of amino acid catabolism fall into three categories that participate in the decarboxylation, transamination, and deamination of amino acids, respectively. We then compared their abundance and found that SZ patients harbored more enzymes in all three categories, and the difference became extremely significant when normalized by their protein intake (Fig. 3A).
Decarboxylases are the major enzymes in the generation of various amines and other derivatives. Using targeted chromatography, we quantitated the major decarboxylation derivatives of amino acids in stool and found that all derivatives were elevated in SZ patients, confirming the activated microbial decarboxylation of amino acids (Fig. 3B). Among these derivatives, indole, kynurenine, and IAA (indole-3-acetic acid), all derived from tryptophan, showed a significant increase in SZ, which was confirmed by the significantly elevated abundance of enzymes responsible for converting tryptophan into kynurenine and IAA as well (Fig S4A). As it has been reported that a variety of tryptophan metabolites are neurologically or immunologically active, the significant increase in tryptophan fermentation products in our study is concordant with some previous reports where tryptophan intake was associated with other brain disorders (32, 34), supporting the potential roles of tryptophan and its metabolites in the pathogenesis of SZ.
Transaminases and deaminases, encoded by both host and microbes, account for catabolizing amino acids into ammonium and urea. The urea concentration, when quantitatively measured and normalized by protein intake, showed significantly elevated in stool in SZ patients but no difference in plasma, which indicated activated microbial production of urea in SZ (Fig. 3C). Microbial fermentation is the primary source of short-chain fatty acids (SCFAs), including branched SCFAs (BSCFAs) and straight-chain SCFAs. BSCFAs, including isobutyric and isovaleric acid, are only generated from transamination of branched-chain amino acids, i.e., Ile, Leu, and Val (35), while straight-chain SCFAs derived from both amino acids (i.e., lysine) and carbohydrates (i.e., dietary fibers) (30). We then performed targeted chromatographic analysis to quantify the concentration of all SCFAs in both plasma and stool. The results showed that both fecal and plasma concentrations of BSCFAs significantly increased in patients (Fig. 3D). In contrast, no significant differences were detected in both plasma and fecal concentrations of SCFAs (Fig S4B), supporting the enhanced catabolism of amino acids in the gut of SZ patients.
To define the deviations in the metabolic profiles between SZ and controls, we utilized a recently developed method for metagenome analysis—Quasi-paired cohort (21). In the list of differential metabolic pathways identified by the Quasi-paired cohort (Wilcoxon signed-rank test for paired samples, FDR < 0.001, mean abundance in control > 10− 5, Table S2), among pathways enriched in controls, 8/31(26%) were in the category of amino acid biosynthesis, whereas none in this category overrepresented in SZ. Notably, we also compared the abundance of aminoacyl-tRNA synthases (EC 6.1.1) between the paired SZ-control samples, and 12 out of the 22 synthases were significantly enriched in controls (Wilcoxon signed-rank test for paired samples, FDR < 0.05) with the other ten enzymes showed no difference in-between. Comparing the total abundance of the 38 amino acid biosynthesis pathways and the 22 aminoacyl-tRNA transferases revealed extremely significant deficiency in the functional potential of protein synthesis in SZ patients (Fig. 3E).
We further constructed random forest classifiers based on the abundance of the 38 amino acid biosynthesis pathways and the 22 aminoacyl-tRNA transferases, respectively, and the models performed excellently in discriminating cases and controls and achieved AUC (area under the curve) of 0.85 and 0.91, respectively, when evaluated with ROC (receiver operating characteristic) curve (Fig. 3F). This result indicated that a shift from protein synthesis in normal conditions to protein catabolism was a major deviation in the intestinal metabolism in SZ patients.
Associations between gut protein fermentation and the impairment of social behaviors
Finally, we investigated whether the protein fermentation in SZ patients correlated to the severity of their psychiatric symptoms, which was reassessed by the PANNS scale (20) for all patients. For each patient, the total score of negative symptoms (N), positive symptoms (P), general symptoms (G), and total score of all symptoms (T) were calculated and used as indicators for the severity of the disease.
First, we investigated the correlations of amino acid catabolizing enzymes to clinical scores and observed that most enzymes exhibited positive correlations to N, G, and T, but sight negative correlations to P (Fig. 4A). Supportively, the fecal concentration of amino acids and their derivatives also presented similar correlations to clinical scores (Fig. 4B). The irrelevance of positive symptoms to intestinal protein fermentation is not surprising as many cofactors such as apastia (food refusal) and psychiatric medications might alleviate positive symptoms (36) and conceal the exacerbating effects of protein fermentation on them. These results suggested the association of intestinal microbial protein fermentation to the severity of psychiatric symptoms, especially the negative ones, which representing the unsolved central pathogenesis of SZ.
In concordance with the hypothesis that the enhanced intestinal protein fermentation might impair human behaviors, we even observed positive correlations of daily protein intake to N, G, and T, but not to P (Fig. 4C). In contrast, the daily carbohydrate intake showed no significant correlation to any psychiatric symptoms (Fig S5). Therefore, it is the food supply of excess proteins but not carbohydrates associated with severe symptoms, highlighting the role of macronutrient intake in the pathogenesis of SZ. In this point of view, the reduction in the daily protein intake in SZ patients might protect against harmful protein fermentation products from the gut in case of ordinary food supplies.