Susceptibility of Human Metabolic Phenotypes to Dietary ... - CiteSeerX

microbial mammalian cometabolite) was higher in the vegetarian than meat diet samples, signaling an alteration of the bacterial composition or .... of...
1 downloads 13 Views 318KB Size
Susceptibility of Human Metabolic Phenotypes to Dietary Modulation Cinzia Stella,† Bridgette Beckwith-Hall,† Olivier Cloarec,† Elaine Holmes,† John C. Lindon,† Jonathan Powell,§ Frans van der Ouderaa,§ Sheila Bingham,‡ Amanda J. Cross,| and Jeremy K. Nicholson*,† Biomolecular Medicine, Division of Surgery, Oncology, Reproductive Biology and Anaesthetics (SORA), Faculty of Medicine, Imperial College London, South Kensington, London SW7 2AZ, United Kingdom, Unilever Corporate Research, Colworth House, Sharnbrook, Bedford MK44 1LQ, United Kingdom, MRC Dunn Human Nutrition Unit, Wellcome Trust/MRC Building, Hills Road, Addenbrookes Hospital, Cambridge CB2 2XY, United Kingdom, and Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, Maryland 20852 Received May 31, 2006

Dietary composition has been shown to influence metabolism and to impact on the prevalence and risk for certain diseases, but hitherto, there have been no systematic studies on the effects of dietary modulation of human metabolic phenotype (metabotype). Here, we have applied 1H NMR spectroscopy in combination with multivariate statistical analysis to characterize the effects of three diets: “vegetarian”, “low meat”, and “high meat” on the metabotype signature of human participants. Twelve healthy male participants (age range of 25-74 years) consumed each of these diets, in a randomized order, for continuous 15-day-periods with an intervening washout period between each diet of 7 days duration. Each participant provided three consecutive 24-hour urine collections on days 13, 14, and 15 of each dietary period, and 1H NMR spectra were acquired on all samples. Pattern recognition analysis allowed differentiation of the characteristic metabolic signatures of the diets with creatine, carnitine, acetylcarnitine, and trimethylamine-N-oxide (TMAO) being elevated in the high-meat consumption period. Application of orthogonal projection to latent structure discriminant analysis (O-PLS-DA) allowed the low-meat diet and vegetarian diet signatures to be characterized, and p-hydroxyphenylacetate (a microbial mammalian cometabolite) was higher in the vegetarian than meat diet samples, signaling an alteration of the bacterial composition or metabolism in response to diet. This work shows the potential for the routine use of metabonomics in nutritional and epidemiological studies, in characterizing and predicting the metabolic effects and the influence of diet on human metabotypes. Keywords: metabonomics • protein • diet • vegetarian • human metabolism • NMR • pattern recognition • metabotype • phenotype

Introduction A western lifestyle, characterized in part by a diet high in fat, red meat, and calories and low in fruits, vegetables, and fiber, has been associated with several chronic diseases including hypertension, diabetes, and coronary heart disease, as well as cancer.1-3 In particular, red meat has been associated with a variety of cancers, such as colorectal adenomas and carcinomas; furthermore, vegetarians have a lower level of glucose, insulin, and insulin resistance (IR), as well as a more favorable lipoprotein profile (increased HDL cholesterol levels and apoA1/apoB ratio in conjunction with reduced LDL levels).4-10 * To whom correspondence should be addressed. E-mail: j.nicholson@ imperial.ac.uk; phone, +44 207 594 3225; fax, +44 207 594 3226. † Unilever Corporate Research. § Biomolecular Medicine. ‡ MRC Dunn Human Nutrition Unit. | National Cancer Institute.

2780

Journal of Proteome Research 2006, 5, 2780-2788

Published on Web 09/13/2006

Although the long-term effects of different diets have been investigated, there is little understanding of the mechanisms of nutritional protection against disease. The evaluation of dietary intake is an essential part of many nutritional studies, and the most common methods of dietary assessment generally include quantitative food diaries and dietary recall or measurement of parameters such as blood glucose, IR, and so forth.11-13 However, these parameters do not tend to generate mechanistic information concerning dietary benefits, and dietary recall can be prone to lack of compliance.14,15 Recent proteomic studies have shown that short chain fatty acids (such as butyrate) generated from microbial digestion of dietary fiber may have important protective functions against gut cancers.16 Chung et al. have shown that N-butyrate specifically arrests the growth of colorectal cancer cells and induces apoptosis, suggesting new butyrate-induced protein-mediated mechanisms for these processes. 10.1021/pr060265y CCC: $33.50

 2006 American Chemical Society

Susceptibility of Human Metabotypes to Dietary Modulation

research articles

In this study, we have evaluated the variation in metabolic phenotype (metabotype) of human participants following different diets over time using high-resolution 1H NMR spectroscopy to generate information on the subtle diet-induced changes in the metabotype. The metabotype is defined as a “probabilistic multiparametric description of an organism in a given physiological state based on analysis of its cell types, biofluids or tissues”.17 Any alteration of the physiological status can disrupt homeostasis, resulting in perturbations of the levels of endogenous biochemicals that are involved in key metabolic processes in the cells and tissues of an organism.18 To maintain homeostasis, and to adjust for the changes in tissue biochemistry, the compositions of body fluids are altered accordingly and lead to different and dynamic “metabotypes”, that may not be observable from gross examination of the phenotype. Thus, monitoring perturbations in biofluid composition may yield valuable information regarding underlying molecular mechanisms. The quantitative measurement of the dynamic multiparametric metabolic response of living systems to pathophysiological stimuli or genetic modification is termed metabonomics,19 and largely relies on chemometric analysis of highresolution spectral data. The use of 1H NMR spectroscopy to study the low MW (molecular weight) composition and molecular interactions in biological fluids, cells, and tissue extracts is well-established.20-22 1H NMR spectra of biological materials can be collected quickly and efficiently, are rich in molecular structural and latent biological information, and allow a large range of low MW metabolites to be viewed simultaneously as a ‘metabolic fingerprint’.23-25 These properties make NMR uniquely suited for the noninvasive study of human populations.26 To simplify NMR data with a view to facilitating interpretation, automated data reduction followed by chemometric analysis has been extensively applied to NMR spectral data sets. Physiological variation per se can be of interest, particularly when considering health and lifestyle issues. However, the analysis of data from human studies is challenging due to the genetic and environmental diversity of human populations and the fact that protocols are generally not as tightly regulated as those for traditional toxicology screening in experimental animal models. To identify meaningful patterns in data where the primary aim is to characterize the functional consequences of disease or drug toxicity, elucidation of the effects of nutrition on the urinary biochemical profile must be established. Previous 1H NMR metabonomic investigations on human urine samples obtained under differing physiological states showed partial classification on the basis of participant,26-29 thereby highlighting inter-individual variability. In this study, chemometric methodologies, including PCA (principal component analysis) and O-PLS (orthogonal projection on latent structure) analysis, have been applied to the classification of urine samples obtained from healthy human individuals in relation to their dietary intake, using samples obtained from a controlled ‘metabolic suite’ study. In this metabolic study, healthy participants consumed a low-meat diet, a diet rich in red meat protein, and a diet with equivalent protein levels to the high red meat diet but from nonmeat sources in three randomly assigned 15-day dietary periods.

approved by the Cambridge Local Ethics Committee. Studies with volunteers were carried out in a metabolic suite in two groups of six, where all food was provided and carefully controlled, and where all specimens could be collected and processed immediately. All food items were accurately weighed to the nearest gram, and main courses were prepared in advance and deep-frozen until required. Deionized water was used throughout for cooking and drinking. Three diets (i) lowmeat diet (60 g/day), (ii) high red meat diet (420 g/day), and (iii) vegetarian diet (420 g/day from nonmeat sources) were randomly assigned in a crossover design in order to avoid an order effect. The 60-g meat diet contained 65-g/day of protein, and the high-meat and vegetarian diets contained 143-150g/day of protein. Each dietary period lasted 15 days, and three 24-h urine collections were made after the volunteers had consumed each diet for 10 days. The participants acted as their own control, and the diets were isoenergetic and kept constant in fat (30% total energy) by exchanging protein for carbohydrate. Basal diets contained 10 MJ per day, and individual energy requirements (to maintain constant body weight) were achieved using 1 MJ increments of bread, butter, and marmalade. A 3 day rotating menu was used, but the items were similar throughout to minimize day-to-day variation. Breakfast was composed of 40 g of cereal (cornflakes on the high-meat diet and Weetabix on the vegetarian diet), 120 g of bread (white bread on the high-meat diet and wholemeal bread on the vegetarian diet), 20 g of butter, 40 g of marmalade, and a daily allowance of 300 g of semi-skimmed milk. Lunch comprised of 100 g of roast beef given as a sandwich (120 g of white bread with 25 g of pickle, 20 g of butter, 50 g of cucumber, and 60 g of tomato), together with 150 g of apple. Wholemeal bread and egg or cheese (40 g) were substituted for the beef on the vegetarian diet. The red meat for the evening meal was cooked in sauces given in the form of 320 g (cooked weight) of minced beef pie, sweet and sour pork, and cottage pie. Vegetable and pasta bake, egg and chips, and vegetable and lentil bake were given for the vegetarian evening meals. Canned fruit (100 g) and ice cream were given for dessert. Volunteers were allowed to choose how much tea and coffee they consumed per day, but amounts were constant from day to day.

Experimental Section Study Protocol. Twelve healthy, Caucasian, nonsmoking men, of normal body wight range (60-91 kg, mean of 76 kg) aged 25-74 were enrolled in this study.30 The protocol was

Urine Sample Preparation. The 24 h urine collections (from 7 a.m. to 7 a.m.) were made into 2 L containers containing 2 g of boric acid/L as a preservative.31 To monitor participant and dietary compliance, para-amino benzoic acid (PABA) tablets were given with each meal on the last 3 days of the dietary period, and the PABA-check method was employed to ensure a complete urine collection.32 Every volunteer provided three consecutive 24-h urine collections on days 13, 14, and 15 of each dietary period. To reduce the pH range of the samples, aliquots of urine samples (600 µL) were mixed with 300 µL of 100 mM phosphate buffer solution (made up in H2O, pH 7.4). The samples were allowed to stand for 10 min prior to centrifugation at 8000 rpm for 5 min in order to remove insoluble material. Aliquots of the supernatants (600 µL) were mixed with 50 µL of TSP/D2O solution (3-trimethylsilyl [2,2,3,32H4] propionate in deuterium oxide, 1 mM final concentration). TSP was added to act an internal chemical shift reference (δ1H 0.0), and D2O was added to provide a lock signal for the NMR spectrometer. Order effects in the statistical analysis were avoided due to the randomized crossover design used. 600 MHz 1H NMR Spectroscopic Analysis of Urine Samples. H NMR measurements were made on a Bruker AMX600 spectrometer operating at 600.13 MHz 1H frequency at ambient

1

Journal of Proteome Research • Vol. 5, No. 10, 2006 2781

research articles probe temperature (298 ( 1 K). For each sample, 96 free induction decays (FIDs) were collected into 65 536 data points using a spectral width of 12 019 Hz and an acquisition time of 2.72 s. An additional relaxation delay (RD) of 2.0 s was added between pulses to allow T1 relaxation. 1H NMR spectra were acquired with suppression of the water resonance using a standard presaturation pulse sequence comprising (RD-90°t1-90°-tm-90°-acquire FID (Bruker Biospin Germany) where t1 is a short delay of typically 3 µs and tm is a mixing time (150 ms). Irradiation of the water signal was achieved during tm and RD. The summed FIDs were multiplied by an exponential weighting function corresponding to a line broadening of 0.3 Hz prior to Fourier transformation. Preprocessing and 1H NMR Urine Spectra. To use as much relevant data as possible, full resolution 1H NMR spectra in the chemical shift range δ 0.0-10.0 were used for pattern recognition. The effect of variations in the presaturation of the water signal was removed by zeroing the intensity values in the region δ 4.5-4.9. All spectra were phase-corrected and referenced to the TSP resonance set at δ 0.0. A baseline correction was also applied to each spectrum using a polynomial curve fit, carried out using in-house software (Matlab Version 7.0, The Mathworks, MA) after importation of 42 000 points/spectrum, using spline cubic interpolation in order to give to all the spectra the same chemical shift scale.33 Finally, to remove concentration differences between dilute and concentrated urine samples, using full resolution data sets, all spectra were normalized to unit total integral. For some analyses, the aromatic region of the NMR spectra was omitted in order to investigate whether the presence of peaks from PABA metabolites interfered in the sample clustering in the PCA analysis. Chemometric Analysis. 1. Principal Component Analysis (PCA). PCA is a projection method used for overviewing and explaining ‘clustering’ and trends within multivariate data and to detect outlier samples. The data are reduced to a few latent variables (or principal components) describing maximum variation within the data. The principal components (PCs) consist of a set of scores, which highlight clustering, trends, or outliers, in the observation (or sample) direction in the data matrix. A corresponding set of ‘loadings’ focuses on the variation in the variable direction, describing the influence from the variables on the scores. The nonsystematic part of variation not explained by the models forms the residuals.34 2. O-PLS (Orthogonal Projection to Latent Structures) Discriminant Analysis. While PCA explains the maximum variance in the data set, PLS defines the maximum separation between defined class samples in the data, and for this reason, class membership of each sample must be known. Once a PLS model is calculated and validated, it can be used for prediction of class membership for unknown samples. O-PLS is an extension of the PLS chemometric method for multivariate data that models sources of variation from X (NMR spectra) that is not related to Y (class, i.e., dietary regimen).35 One special case of PLS is O-PLS-Discriminant Analysis (O-PLS-DA) where PLS is used to find the relationship between the descriptors in X and the class identity of the samples, described by a dummy matrix containing the class information.36 These sources of variation may disturb the multivariate modeling, cause imprecise predictions for new samples, and also affect the robustness of the method over time. For the removal of such undesirable systematic variation in the data, different types of preprocessing methods prior to PCA or PLS have been commonly used, such as Orthogonal Signal Correction (OSC).37 However, the advan2782

Journal of Proteome Research • Vol. 5, No. 10, 2006

Stella et al.

tage of O-PLS compared to OSC is that the orthogonal variation partitioning is implemented in the algorithm and avoids the two step process, OSC followed by PLS, making the validation of the model easier and more robust.35 The interpretation is also improved because any “structured noise” is modeled separately from the data variation common to X (the NMR spectra) and Y (the class variable). The method used in this study to interpret the O-PLS coefficients (or loadings relating to the relative influence of a particular signal or metabolite on class discrimination) is to calculate a model based on unit variance scaled data, and then, in a first step, back-transform the predictive O-PLS coefficients by multiplying all values by their respective standard deviation. In a second step, the back-transformed loading is plotted using a color corresponding to the coefficient value in the model that represents the correlation of the variable X (NMR data) with Y (class or dietary intervention in this case) for each data point. The scale range of colors is set using the maximum and the minimum of the unit variance scaled model weight. The interpretation of the loadings is therefore straightforward for the spectroscopist because the resulting plot provides a loading with the same shape as that of a spectrum (based on the covariance matrix), but on the same plot, the important variables for the discrimination between the classes (calculated from the correlation matrix) are highlighted by the color code.

Results and Discussion Typical examples of 1H NMR spectra of urine samples obtained from a participant representing each of the dietary regimes are depicted in Figure 1. Resonance assignments were performed according to the existing literature.38 The most marked differences, observable on the 1H NMR spectra between the three diets, included elevated urinary levels of creatinine, taurine, carnitine, trimethylamine-N-oxide (TMAO), and methylhistidine. These perturbations were consistent among the three samples for each participant and relate to the high-meat group. The spectral changes differentiating low-meat and vegetarian groups were, on the other hand, more subtle. The split observed in the NMR signal of citrate (Figure 1) is due to complexation with boric acid. The effects of adding boric acid to biological samples as a preservative has been studied previously and found to have little or no influence on the pattern recognition analysis.39 PABA has long been used as a marker of completeness of 24 h urine samples in nutritional epidemiological studies and prescribed diet regimes.40,41 The conventional PABA recovery cutoff for completeness is 85% or above.31,41 Colorimetric or HPLC methods are used for measuring PABA and its conjugated metabolites in human urine after oral administration. In the current study, dietary compliance was found to be acceptable for all participants on all diets. The chemometric analysis has been performed both with and without inclusion of the aromatic region (δ 10-4.3) of NMR spectra in order to evaluate whether PABA metabolites interfere with the pattern recognition. Results obtained without the inclusion of the aromatic region were not different from those obtained with the full-spectrum range, indicating that the presence of PABA did not influence the pattern recognition (data not shown). Initial analysis of the data set was conducted with no knowledge as to the dietary status of the samples from each individual using PC analysis (with unit variance scaling) with and without taking account of the aromatic region. Interindividual variation dominated the PCA analysis (Figure 2A)

Susceptibility of Human Metabotypes to Dietary Modulation

research articles

Figure 1. 1H NMR spectra of urine from a single individual on each of the three diets (HM ) high meat; LM ) low meat, and VE ) vegetarian): only significant metabolites have been labeled for each of the three urinary metabolic profiles.

with some participants showing a much greater metabolic variation than others. For example, samples from participant 1 show a tighter clustering than the samples from participant 6. However, intra-individual variation was dominated by dietary exposure, and analysis of samples from each participant revealed three clusters (Figure 2B) relating to each of the three diets. With reference to the loadings (variable coefficients) for the spectral descriptors, it could easily be established which cluster related to the high-meat diet. The most marked differences relating to the high-meat diet included elevated urinary levels of creatinine, creatine, TMAO, taurine, and 1- and 3- methylhistidine. These perturbations were consistent among the samples identified as high-meat, and metabolites such as carnitine and 1- and 3-methylhistidine have been cited in the literature as markers of meat ingestion.42 In one epidemiological study, monitoring of urinary methylhistidine excretion as a marker of meat consumption was used to successfully predict vegetarian status in a study of 126 individuals based on food frequency questionnaires.43

For the two remaining diets (low-meat and vegetarian), the spectral changes between the two were more subtle, and this was confirmed with PCA analysis (data not shown): clear, but incomplete separation of the high-meat diet was observed, while the low-meat and vegetarian diet participants were partially comapping. O-PLS-DA analysis was applied to the data matrix to remove variation in X (NMR spectra) that was unrelated (orthogonal) to Y (diet). In this study, removal of nonrelated noise resulted in a much improved separation of the low-meat diet and the vegetarian diet, and all three dietary regimes formed distinct clusters as shown in Figure 3A. The application of O-PLS-DA produced a regression model (Figure 3B) describing the maximum separation between predefined classes, which focused on the actual class discriminating variation in the data. Since the principle aim of the study was to explore the metabolic consequences of a diet with high-meat protein content, analysis was focused on the differentiation between the high-meat and vegetarian diets which contained the same amount of protein (143-149-g protein/day). The O-PLS-DA loadings for the pairwise discrimination of the highJournal of Proteome Research • Vol. 5, No. 10, 2006 2783

research articles

Stella et al.

Figure 3. (A) O-PLS-DA score plot analysis and (B) summary of the O-PLS-DA models coefficients for the discrimination between HM, LM, and VE classes (Q2Yhat, cross-validated prediction of Y; R2Yhat, predicted variation of Y; R2X, modeled variation of X; R2Xcorr, modeled joint X-Y covariation of Y). Figure 2. (A) PC analysis performed on all data and color-coded according to participant. (B) PC analysis performed on one participant and color-coded according to the diet.

meat and vegetarian groups are presented in Figures 4 and 5, respectively. The color scale corresponds to the UV model variable weights and allows the identification of significant class-discriminating metabolites (Table 1), responsible for the clustering patterns. The O-PLS-DA coefficients of the high-meat group (Figure 4) show that the urinary concentration of creatine and creatinine (Figure 4) are much greater than those from the lowmeat (figure not shown) and vegetarian groups (Figure 5). Creatine is a metabolite that is produced naturally by the human body, and it is found mainly in the red muscle tissue, but it is also present in other tissues with high and fluctuating energy demands such as brain, retina, and testes.44 Normally, creatine, which functions to increase the availability of cellular ATP (adenosine triphosphate), is acquired through regular dietary intake of products such as meat and fish, which are high in creatine/creatinine, so that urinary levels increase in response to meat consumption. Although creatinine excretion is related to body weight, creatinine excretion therefore is not constant as it varies from day to day according to the amount of meat consumed.45 The O-PLS-DA coefficients of the carnitine peak and acetylcarnitine indicate that the highest concentrations of these metabolites are found in the high-meat group (Figure 4 and Table 1), and these metabolites are thus responsible for the clustering separation of this group from the low-meat and vegetarian groups. Carnitine, a trimethylated amino acid, is a 2784

Journal of Proteome Research • Vol. 5, No. 10, 2006

cofactor required for transformation of free long-chain fatty acids into acylcarnitines and for their subsequent transport into the mitochondrial matrix, where they undergo β-oxidation for cellular energy production. Mitochondrial fatty acid oxidation is the primary energy source in heart and skeletal muscle. It is assumed that a diet adequate in protein will supply enough exogenous, and promote any additional endogeneous, synthesis needed to supply an individual’s requirements for L-carnitine. However, since L-carnitine is found primarily in animal proteins, with red meat regarded as the richest source, it is theoretically possible to consume a high protein diet consisting of beans, legumes, or egg whites and still promote a relative deficiency. Since a vegetarian diet is very low in exogenous carnitine, and is potentially low in some of the substrates required for carnitine synthesis, there is some concern that following a strict vegetarian diet might produce a carnitine deficiency in some individuals.46 Since adequate dietary lysine is required as a substrate for carnitine synthesis, a deficiency of this amino acid or the other cofactorssiron, vitamin C, B6, niacinsmight also compromise carnitine status. The urinary concentration of trimethylamine-N-oxide (TMAO) and taurine after ingestion of a high-meat diet is higher in comparison to the low-meat and vegetarian diets (Figure 4 and Table 1). TMAO has been correlated to meat and fish ingestion, but the formation of TMAO and γ-butyrobetaine from Lcarnitine by enteric bacteria has also been demonstrated in rats and humans.47 While TMA (trimethylamine) and TMAO are generally regarded as nontoxic substances, they are of clinical interest because of their potential to form the carcinogen N-nitrosodimethylamine.48 The excretion of the amino acid taurine, a conditionally essential amino acid playing an im-

Susceptibility of Human Metabotypes to Dietary Modulation

research articles

Figure 4. The loading coefficient plot of the high-meat group corresponding to the scores plot represented in Figure 3A: the signal orientation indicates the higher (positive orientation) or lower (negative orientation) concentration of characteristic metabolites. The color of the signal reflects the significance of the correlation with diet (red indicates highest significance and blue indicates no significance).

Figure 5. The loading coefficient plot of the vegetarian group corresponding to the scores plot represented in Figure 3A: the signal orientation indicates the higher (positive orientation) or lower (negative orientation) concentration of characteristic metabolites. The color of the signal reflects the significance of the correlation with diet (red indicates highest significance and blue indicates no significance). The region between δ 6.85 and 7.20 is expanded to show the other resonances of p-hydroxyphenylacetate.

portant role in many physiological functions, is also increased following a high-meat diet (Figure 4 and Table 1). This metabolite is commonly associated with ingestion of a high

amount of proteins (Atkins diet), and urinary taurine excretion in vegan subjects is about half the value of that in the omnivores.49,50 Journal of Proteome Research • Vol. 5, No. 10, 2006 2785

research articles

Stella et al.

Table 1. List of Most Characteristic Metabolites Found in the High-Meat and Vegetarian Groups Obtained with O-PLS-DA Analysis assigned metabolite

chemical shift

correlation (HM)

correlation (LM)

correlation (VE)

covariance (HM)

covariance (LM)

covariance (VE)

creatine taurine TMAO carnitine acetylcarnitine creatine glutamine p-hydroxyphenylacetatea N-acetylglutamate N,N,N trimethyllysine

3.934 (s) 3.425 (t) 3.269 (s) 3.228 (s) 3.195 (s) 3.041 (s) 2.445 (m) 3.452 (s) 2.035 (s) 3.112 (s)

0.75 0.66 0.71 0.81 0.70 0.72 0.72 -0.44 -0.49 -0.11

-0.38 -0.32 -0.47 -0.49 -0.42 -0.33 -0.32 -0.24 -0.31 0.57

-0.38 -0.35 -0.23 -0.33 -0.28 -0.38 -0.40 0.69 0.82 -0.47

0.45 0.05 0.71 0.17 0.03 0.49 0.01 -0.04 -0.03 -0.02

-0.23 -0.02 -0.47 -0.10 -0.02 -0.23 -0.01 -0.02 -0.02 0.08

-0.22 -0.02 -0.23 -0.07 -0.01 -0.26 -0.01 0.07 0.04 -0.07

a

Doublet at 6.867 and doublet at 7.172.

A high concentration of N-acetyl-5-hydroxytryptamine was found in the high-meat urinary profile. This compound, also known as N-acetylserotonin, is produced in one of the steps of tryptophan metabolism (tryptophan f 5-hydroxytryptophan f 5-hydroxytryptamine f N-acetyl-5-hydroxytryptamine). The amount of tryptophan (an amino acid necessary for normal growth in infants and for nitrogen balance in adults) in the foods that are eaten has only a small influence upon the amount of tryptophan that enters the brain. Large-molecule amino acids, among them tryptophan, compete with each other to enter “gates” between the circulating blood stream and the relatively confined brain fluids.51 A high-protein meal (e.g., high in meat, dairy products, and eggs) provides many other amino acids that compete with tryptophan for entry into the brain; the end result is less tryptophan passing into the brain, a decrease in the synthesis of serotonin in the brain, which may explain the high amount of N-acetyl-5-hydroxytryptamine in the urine.52 Conversely, a low-protein, carbohydrate-rich diet (full of starches, vegetables, and fruits) results in the highest levels of serotonin in the brain, because fewer large-molecule amino acids are competing with tryptophan to enter the brain.53,54 Another characteristic metabolite of the high-meat group is glutamine, which is not an essential amino acid, since it can be readily produced from glutamate in all tissues, with muscle tissue being the primary source of glutamine in the blood. Although there is no dietary reference intake (for under normal conditions, glutamine is not a necessary dietary constituent), during critical illness, severe trauma, intestinal disease, starvation, total parenteral nutrition (intravenous feeding), wasting (excessive loss of lean body mass), and extreme endurance exercise, the body’s need and consumption of glutamine can exceed the ability of tissues to produce this amino acid. Under such stress conditions, dietary glutamine is beneficial. However, while glutamine is nontoxic and probably harmless, the benefits of consuming dietary supplements containing this amino acid have not been proven.55 Low glutamine levels in plasma may reflect chronic insufficient essential amino acid dietary intake, and this may lead to ammonia detoxification problems.56 Characteristic NMR peaks found for the urines of the vegetarian group can be assigned to p-hydroxyphenylacetate (Figure 5). Tyrosine from dietary protein is the parent compound from which p-cresol, p-hydroxybenzoate, and p-hydroxyphenylacetate are formed (Figure 6A). These compounds are not products of normal human metabolism but are produced by bacteria and protozoa that can populate the gut.57 Both the transamination of tyrosine to form p-hydrohyphenylacetate (HPA) and the decarboxylation to p-cresol are carried out by a limited number of microfloral species including 2786

Journal of Proteome Research • Vol. 5, No. 10, 2006

Figure 6. (A) Tyrosine metabolism; (B) L-carnitine biosynthesis.

Clostridium difficile. Since Proteus vulgaris can do only the first of these steps, HPA will increase in urine if P. vulgaris is the predominant organism. This might indicate that some compounds that appear in urine, such as p-hydroxyphenylacetate, are unique metabolic products of the microbes that inhabit the lumen of gut, which is strictly dependent on the dietary regime.58-60 In addition, p-hydroxyphenylacetic aciduria has been found useful in detecting small bowel disease associated with anaerobic bacterial overgrowth.61 Another compound characteristic of the urinary profile of the vegetarian group is N6,N6,N6-trimethyllysine (TML): the negative loading coefficient indicates a low amount of this metabolite at urinary level following a vegetarian diet. TML is generated within lysosomes by the hydrolysis of proteins containing lysines that are trimethylated at their -amino group by a protein-dependent methyltransferase (Figure 6B). Apart from the dietary intake of carnitine, most eukaryotes are able to synthesize this compound from TML, and this would explain the low levels of TML following a vegetarian diet.

Susceptibility of Human Metabotypes to Dietary Modulation

research articles

One nutrient that is almost totally absent in a vegetarian diet is L-carnitine, which plays a key role within several energy production pathways. If intake of L-carnitine is low, the body must almost entirely rely on the endogenous synthesis to meet the needs, and this could explain the lower level of TML in urine. The changes in the urinary metabolite profile associated with a high-meat content diet in comparison with a vegetarian diet have been characterized using a combination of NMR spectroscopy with chemometric analyses. This work should provide a baseline for extended dietary studies in which associations between perturbations in specific metabolic pathways and disease risk or prevalence may be explored.

(9) Valachovicova, M.; Krajcovicova-Kudlackova, M.; Blaziceck, P.; Babinska, K. Eur. J. Nutr. 2006, 45, 52-54. (10) Szeto, Y.; Kwok, T.; Benzie, I. Nutrition 2004, 20, 863-866. (11) Cade, J.; Thompson, R.; Burley, V.; Warm, D. Public Health Nutr. 2002, 5, 567-587. (12) Lean, M. E. J.; Anderson, A. S.; Morrison, C.; Currall, J. Eur. J. Clin. Nutr. 2003, 57, 667-673. (13) Masson, L. F.; McNeill, G.; Tomany, J. O.; Simpson, J. A.; Peace, H. S.; Wei, L.; Grubb, D. A.; Bolton Smith, C. Public Health Nutr. 2003, 6, 313-321. (14) Roberts, S. B.; Ferland, G.; Young, V. R.; Morrow, F. M.; Heyman, M. B.; Melanson, K.; G. E., D. Am. J. Clin. Nutr. 1991, 54, 774782. (15) Leiba, A.; Vald, A.; Peleg, E.; Shamiss, A.; Gross, A. E. Nutrition 2005, 21, 462-466. (16) Tan, H. T.; Zubaidah, R. M.; Tan, S.; Hooi, S. C.; Chung, M. C. J. Proteome Res. 2006, 5, 1098-1106. (17) Gavaghan, C. L.; Holmes, E.; Lenz, E. M.; Wilson, I. D.; Nicholson, J. K. FEBS Lett. 2000, 484, 169-174. (18) Nicholson, J. K.; Wilson, I. D. Nat. Rev. Drug Discovery 2003, 668676. (19) Nicholson, J. K.; Lindon, J. C.; Holmes, E. Xenobiotica 1999, 29, 1181-1189. (20) Nicholson, J. K.; Wilson, I. D. Prog. NMR Spectrosc. 1989, 21, 449501. (21) Bollard, M. E.; Holmes, E.; Lindon, J. C.; Mitchell, S. C.; Branstetter, D.; Zhang, W.; Nicholson, J. K. Anal. Biochem 2001, 295, 194-202. (22) Nicholls, A. W.; Lindon, E.; Lindon, J. C.; Shockcor, J. P.; Farrant, R. D.; Haselden, J. N.; Damment, S. J. P.; Waterfield, C. J.; Nicholson, J. K. Chem. Res. Toxicol. 2001, 14, 975-987. (23) Spraul, M.; Neidig, P.; Klauck, U.; Kessler, P.; Holmes, E.; Nicholson, J. K.; Sweatman, B. C.; Salman, S. R.; Farrant, R. D.; Rahr, E.; Beddel, C. R.; Lindon, J. C. J. Pharm. Biomed. Anal. 1994, 12, 1215-1225. (24) Holmes, E.; Foxall, P. J.; Spraul, M.; Farrant, R. D.; Nicholson, J. K.; Lindon, J. C. J. Pharm. Biomed. Anal. 1997, 15, 1647-1659. (25) Beckwith-Hall, B., Bridgette M.; Nicholson, J. K.; Nicholls, A. W.; Foxall, P. J.; Lindon, J. C.; Connor, S. C.; Abdi, M.; Connelly, J.; Holmes, E. Chem. Res. Toxicol. 1998, 11, 260-272. (26) Dumas, M. E.; Maibaum, E. C.; Teague, C.; Ueshima, H.; Zhou, B.; Lindon, J. C.; Nicholson, J. K.; Stamler, J.; Elliott, P.; Chan, Q.; Holmes, E. Anal. Chem 2006, 78, 2199-2208. (27) Beckwith-Hall, B., Bridgette M.; Brindle, J. T.; Barton, R. H.; Coen, M.; Holmes, E.; Nicholson, J. K.; Antti, H. Analyst 2002, 127, 12831288. (28) Lenz, E. M.; Bright, J.; Wilson, I. D.; Hughes, A.; Morrisson, J.; Lindberg, H.; Lockton, A. J. Pharm. Biomed. Anal. 2004, 36, 841849. (29) Lenz, E. M.; Bright, J.; Wilson, I. D.; Morgan, S. R.; Nash, A. F. J. Pharm. Biomed. Anal. 2003, 33, 1103-1115. (30) Cross, A. J.; Pollock, J. R. A.; Bingham, S. A. Cancer Res. 2003, 63, 2358-2360. (31) Bingham, S. A.; Cummings, J. Clin. Sci. 1983, 64, 629-635. (32) Jakobsen, J.; Ovesen, L.; Fagt, S.; Pedersen, A. N. Eur. J. Clin. Nutr. 1997, 51, 514-519. (33) Ebbels, T. M. D.; Lindon, J. C.; Nicholson, J. K.; Holmes, E. Patent Number 2001002938200112200, US 2002145245, 2002. (34) Jolliffe, I. T. Principal Component Analysis; Springer: New York, 1986. (35) Trygg, J.; Wold, S. J. Chemom. 2002, 16, 283-293. (36) Sjostrom, M.; Wold, S.; Soderstrom, B. PLS Discriminant Plots; Elsevier Science Publ B. V.: North Holland, Amsterdam, 1986. (37) Wold, S.; Antti, H.; Lindgren, F.; Ohman, J. Chemom. Intell. Lab. Syst. 1998, 44, 175-185. (38) Nicholson, J. K.; Foxall, P. J.; Spraul, M.; Farrant, R. D.; Lindon, J. C. Anal. Chem. 1995, 67, 793-811. (39) Teague, C., Ph.D. Theses, University of London-Imperial College, London, U.K., 2004. (40) Jacobsen, J.; Pedersen, A. N.; Ovesen, L. Eur. J. Clin. Nutr. 2003, 57, 138-142. (41) Leclerq, C.; Maiani, G.; Polito, A.; Ferro-Luzzi, A. Nutrition 1991, 7, 350-354. (42) Elia, M.; Carte, A.; Bacon, S.; Smith, R. Clin. Sci. 1980, 59, 509511. (43) Myint, T.; Fraser, G. E.; Lindsted, K. D.; Knutsen, S. F.; Hubbard, R. W.; Bennett, H. W. Am. J. Epidemiol. 2000, 152, 752-755. (44) Wyss, M.; Kaddurah-Daouk, R. Physiol. Rev. 2000, 80, 1107-1213. (45) Bingham, S. A.; Cummings, J. Hum. Nutr.: Clin. Nutr. 1985, 39c, 343-353.

Conclusions A clear separation of the urinary metabolic profiles according to the type of diet and the use of O-PLS-DA analysis allowed the identification of the characteristic metabolites for each dietary regime. These results indicate the potential of the routine use of metabonomics in nutritional and epidemiological studies, where it is useful to be able to distinguish individuals according to their meat intake. Some of these metabolites, for example, p-hydroxyphenylacetate in the vegetarian group, are not products of human metabolism, but of the gut microflora, indicating an alteration of the bacterial composition or metabolism in response to diet. O-PLS-DA analysis also has the potential to generate new pathway information that can subsequently be tested, for example by using labeled substrates. Human metabotypes are determined by a combination of genetic and environmental factors, with diet representing a source of significant environmental diversity. We have shown conclusively the heavy influence of diet on the short-term metabotype which could provide valuable information in demographic and population studies on dietary variation in man and the degree of response to dietary modulation shown by individuals. This may ultimately give insight into the role of diet and nutrition in the development of human diseases.

Acknowledgment. Unilever Research (Sharnbrook, Bedford MK44 1LQ, U.K.) is acknowledged for the funding of Dr. Cinzia Stella. References (1) Cross, A. J.; Peters, U.; Kirsh, V. A.; Andriole, G. L.; Reding, D.; Hayes, R. B.; Sinha, R. Cancer Res. 2005, 65, 11779-11784. (2) Verhoef, P.; van Vliet, T.; Olthof, M. R.; Katan, M. B. Am. J. Clin. Nutr. 2005, 82, 553-558. (3) Ahmed, F. E. J. Environ. Sci. Health, Part C 2004, C22, 91-147. (4) Park, Y.; Hunter, D. J.; Spiegelman, D.; Bergkvist, L.; Berrino, F.; van den Brandt, P. A.; Buring, J. E.; Colditz, G. A.; Freudenheim, J. L.; Fuchs, C. S.; Giovannucci, E.; Goldbohm, R. A.; Graham, S.; Harnack, L.; Hartman, A. M.; Jacobs, D. R., Jr.; Kato, I.; Krogh, V.; Leitzmann, M. F.; McCullough, M. L.; Miller, A. B.; Pietinen, P.; Rohan, T. E.; Schatzkin, A.; Willett, W. C.; Wolk, A.; ZeleniuchJacquotte, A.; Zhang, S. M.; Smith-Warner, S. A. JAMA, J. Am. Med. Assoc. 2005, 294, 2849-2857. (5) Lin, J.; Zhang, S. M.; Cook, N. R.; Rexrode, K. M.; Liu, S.; Manson, J. E.; Lee, I. M.; Buring, J. E. Cancer, Causes Control 2005, 16, 225-233. (6) Michels, K. B.; Fuchs, C. S.; Giovannucci, E.; Colditz, G. A.; Hunter, D. J.; Stampfer, M. J.; Willett, W. C. Cancer Epidemiol., Biomarkers Prev. 2005, 14, 842-849. (7) Slattery, M. L.; Curtin, K. P.; Edwards, S. L.; Schaffer, D. M. Am. J. Clin. Nutr. 2004, 79, 274-281. (8) French, L.; Kendall, S.; Stephens, M. B. J. Fam. Pract. 2003, 52, 892-893.

Journal of Proteome Research • Vol. 5, No. 10, 2006 2787

research articles (46) Etzioni, A.; Levy, J.; Nitzan, M. Arch. Dis. Child. 1984, 59, 177179. (47) Rebouche, C. J.; Seim, H. Ann. Rev. Nutr. 1998, 18, 39-61. (48) Cross, A. J.; Pollock, J. R. A.; Bingham, S. A. Cancer Res. 2003, 63, 2358-2360. (49) Laidlaw, S. A.; Shultz, T. D.; Cecchino, J. T.; Kopple, J. D. Am. J. Clin. Nutr. 1988, 47, 660-663. (50) Rana, S. K.; Sanders, T. A. B. Br. J. Nutr. 1986, 56, 17-27. (51) Tews, J. K.; Kim, Y. W.; Harper, A. E. J. Nutr. 1979, 109, 304-315. (52) Young, I. M.; Leone, R. M.; Silman, R. E. Biomed. Mass Spectrom. 1985, 12, 319-337. (53) Fernstrom, J. D. Metabolism 1977, 26, 207-223. (54) Fernstrom, J. D.; Wurtman, R. J. Science 1971, 174, 1023-1025. (55) Abcouwer, S. F. Encyclopedia of Dietary Supplements, 2005.

2788

Journal of Proteome Research • Vol. 5, No. 10, 2006

Stella et al. (56) Newsholme, P.; Lima, M. M. R.; Procopio, J.; Pithon-Curi, T. C.; Doi, S. Q.; Bazotte, R. B.; Curi, R. Braz. J. Med. Biol. Res. 2003, 36, 153-163. (57) Selmer, T.; Andrei, P. I. Eur. J. Biochem. 2001, 268, 13631372. (58) Cole, C. B.; Fuller, R. Food Microbiol. 1987, 4, 77-81. (59) Bouhnik, Y.; Raskine, L.; Simoneau, G.; Paineau, D.; Bornet, F. Nutr. J. 2006, 5, 8. (60) Magee, E. A.; Edmond, L. M.; Tasker, S. M.; Kong, S. C.; Curno, R.; Cummings, J. H. Nutr. J. 2005, 4, 7. (61) Chalmers, R. A., Valman, H. B., Liberman, M. M. Clin. Chem. 1979, 25, 1791-1794.

PR060265Y