Next Article in Journal
Consumption of Ultra-Processed Foods Is Associated with Free Sugars Intake in the Canadian Population
Previous Article in Journal
The Paradox of the Mediterranean Diet in Pediatric Age during the COVID-19 Pandemic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Biology of Veganism: Plasma Metabolomics Analysis Reveals Distinct Profiles of Vegans and Non-Vegetarians in the Adventist Health Study-2 Cohort

by
Fayth L. Miles
1,2,3,4,
Michael J. Orlich
1,2,3,
Andrew Mashchak
1,
Paulette D. Chandler
5,
Johanna W. Lampe
6,
Penelope Duerksen-Hughes
4 and
Gary E. Fraser
1,2,7,*
1
Adventist Health Study, Research Affairs, Loma Linda University, Loma Linda, CA 92350, USA
2
Center for Nutrition, Healthy Lifestyle and Disease Prevention, School of Public Health, Loma Linda University, Loma Linda, CA 92350, USA
3
Department of Preventive Medicine, School of Medicine, Loma Linda University, Loma Linda, CA 92350, USA
4
Department of Basic Science, School of Medicine, Loma Linda University, Loma Linda, CA 92350, USA
5
Division of Preventive Medicine, Brigham and Women’s Hospital, Boston, MA 02215, USA
6
Public Health Sciences Division, Cancer Prevention Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA
7
Department of Medicine, School of Medicine, Loma Linda University, Loma Linda, CA 92350, USA
*
Author to whom correspondence should be addressed.
Nutrients 2022, 14(3), 709; https://doi.org/10.3390/nu14030709
Submission received: 7 January 2022 / Revised: 31 January 2022 / Accepted: 31 January 2022 / Published: 8 February 2022
(This article belongs to the Section Nutrition and Public Health)

Abstract

:
It is unclear how vegetarian dietary patterns influence plasma metabolites involved in biological processes regulating chronic diseases. We sought to identify plasma metabolic profiles distinguishing vegans (avoiding meat, eggs, dairy) from non-vegetarians (consuming ≥28 g/day red meat) of the Adventist Health Study-2 cohort using global metabolomics profiling with ultra-performance liquid chromatography mass spectrometry (UPLC-MS/MS). Differences in abundance of metabolites or biochemical subclasses were analyzed using linear regression models, adjusting for surrogate and confounding variables, with cross-validation to simulate results from an independent sample. Random forest was used as a learning tool for classification, and principal component analysis was used to identify clusters of related metabolites. Differences in covariate-adjusted metabolite abundance were identified in over 60% of metabolites (586/930), after adjustment for false discovery. The vast majority of differentially abundant metabolites or metabolite subclasses showed lower abundance in vegans, including xanthine, histidine, branched fatty acids, acetylated peptides, ceramides, and long-chain acylcarnitines, among others. Many of these metabolite subclasses have roles in insulin dysregulation, cardiometabolic phenotypes, and inflammation. Analysis of metabolic profiles in vegans and non-vegetarians revealed vast differences in these two dietary groups, reflecting differences in consumption of animal and plant products. These metabolites serve as biomarkers of food intake, many with potential pathophysiological consequences for cardiometabolic diseases.

1. Introduction

Findings from the Adventist Health Study-2 Cohort have demonstrated that vegetarian dietary patterns have been associated with many positive health outcomes, including lower risk of metabolic syndrome (56%) [1], lower incidence of diabetes (39–62%) [2], lower overall mortality (hazard ratio (HR): 0.88, 95% confidence interval (CI): 0.80–0.97), and cardiovascular disease mortality (HR: 0.71, 95% CI: 0.57–0.90) for males [3], as well as lower risk of gastrointestinal cancers (HR: 0.76, 95% CI: 0.63–0.90) and cancer overall (HR: 0.92, 95% CI: 0.85–0.99) [4,5]. Vegans, particularly have lower risk of female-specific cancers (HR: 0.66, 95% CI: 0.47–0.92) [5], prostate cancer (HR: 0.65; 95% CI: 0.49, 0.85) [6], and cancer overall (HR: 0.84; 95% CI: 0.72–0.99) [5] relative to non-vegetarians. Metabolic pathways associated with such health-promoting dietary patterns, and particularly a vegan diet, remain to be elucidated.
The metabolome may be modulated by select phytochemical-rich plant-based foods. Many studies have shown higher abundance of potentially beneficial metabolites with higher consumption of fruits and vegetables. Metabolic profiles have been compared in subjects adhering to western or prudent dietary patterns in cross-sectional studies, where differences in plasmalogens, phosphatidylcholines, cholesteryl esters, and acylcarnitines, among other compounds, were reported [7,8]. Additionally, a Mediterranean dietary pattern or score has been associated with changes in carbohydrate and lipid metabolites including phosphatidylcholines and acylcarnitines, amino acids, biogenic amines, xenobiotics, and microbial cometabolites [9,10,11]. Metabolic profiles have also been characterized and compared among dietary patterns defined by various healthy eating indices [9,12].
Few studies have compared metabolic profiles between habitual vegetarians, especially vegans, and omnivores. Long-term, habitual dietary patterns may be associated with more stable alterations in the metabolome. In one plasma metabolomics study where omnivores were compared to vegans and other vegetarians, omnivores were found to have highest concentrations of lipid metabolites, acylcarnitines, glycerophospholipids, and sphingolipids, with many of these mediating biological processes relevant to various pathophysiological conditions [13]. Vegans have also been found to have higher abundance of metabolites related to microbial and polyphenol metabolism in spite of minimal differences in gut microbial composition [14].
We previously measured biomarkers of dietary intake in biospecimens from vegetarians and non-vegetarians in the Adventist Health Study-2 (AHS-2) cohort participating in a calibration study (n = 1011), and reported notable differences, including increased abundance of phytochemicals and healthier profiles of fatty acids in the vegetarians [15]. We also recently reported differences in DNA methylation patterns between vegans and non-vegetarians, noting some alterations in genes with relevance to cancer development [16]. It is probable that dietary profiles characterized by high consumption of phytochemical and polyphenolic compounds present in plant-based foods stimulate metabolic pathways associated with disease prevention.
Metabolomics may prove to be very helpful in elucidating diet-mediated biochemical functions by linking eating patterns with the pathophysiology of chronic diseases. In light of our previous findings associating vegetarian dietary patterns with reduced risk of chronic conditions and metabolic diseases, we sought to analyze differences in profiles of metabolic compounds comparing vegans to non-vegetarians in the AHS-2 cohort, using an untargeted metabolomics approach.

2. Materials and Methods

2.1. Study Design

The AHS-2 cohort (established 2002–2007) consists of over 96,000 Seventh Day Adventists age ≥29, with roughly 55% distributed among vegan (8–9%), and lacto-ovo-(29%), pesco-(10%) and semi-non-vegetarian diet groups (5–6%), and the remainder categorized as non-vegetarians (45%). The cohort was established to examine the influence of diet and lifestyle on various health outcomes including cancer, cardiovascular disease, obesity, diabetes, and total mortality. All participants completed a food frequency questionnaire (FFQ) at enrollment (2002–2007), which was used to assign dietary status [17]. The FFQ was a self-administered semi-quantitative instrument collecting basic demographic information including sex, age, and ethnicity, and assessing dietary habits over the past year. The FFQ collected information on frequency of consumption of over 200 foods, including fresh, cooked or canned fruits and vegetables, legumes (lentils, soybeans, and other beans), breads and grains, soy foods/drinks/supplements, dairy, eggs, red meats, processed meats, poultry, and fish consumption, among other items. Frequency categories varied with food type, and portion sizes included three levels: standard, ½ or less, and 1½ or more. Nutrient intake data were collected with the use of the Nutrition Data System for Research software versions 4.06 and 5.03 (The Nutrition Coordinating Center, University of Minnesota, Minneapolis, MN, USA) [18]. Validity of dietary intake has been assessed extensively using 24-h diet recalls and biomarkers [17], and methods developed for dealing with measurement error [17,19,20]. Biennial health and hospitalization history forms (HHF) captured changes in exposures as well as lifetime dietary pattern trends for AHS-2 members [21]. The current study included vegans and non-vegetarians from the AHS-2 cohort who were previously recruited to participate in one of various substudies, where they were asked to attend clinic and provide blood, urine, saliva, and/or adipose samples. These included a calibration study [17], Biological Manifestations of Religion Substudy (BioMRS), which was nested within the Biology, Religion, and Health substudy [22] and other pilot studies (Supplementary Materials, Figure S1). For the calibration and pilot studies, fasting blood samples were collected at field clinics held in church halls, as described previously, where healthy participants were selected randomly by church (Adventist churches were randomly selected from those within the US and Canada, weighted by church size) and then subject within church [23], while BioMRS participants were recruited to clinic sites in Loma Linda, Riverside, and Los Angeles, CA, USA [22]. Anthropometric data were also collected during clinic visits. For the current cross-sectional study, 96 subjects, including vegans and non-vegetarians were selected with stratified random sampling, balancing by sex and race, and excluding individuals with extreme body mass index (BMI) (<14 or >50) and total caloric intake of <500 or >4500 kcal/day. Vegans were defined as never or rarely (less than once per month) consuming eggs, dairy, fish, and other meats, based on responses to the FFQ. Non-vegetarians were defined as consuming non-fish meats at least once a month and any meat (including but not only fish) more than once per week. Vegans were of primary interest for this pilot study because of their complete avoidance of animal products and high consumption of fruits, vegetables, whole grains, legumes, and soy [24]. Non-vegetarians, who consume relatively low amounts of non-fish meats relative to the general population, were selected only if red meat consumption was ≥28 g/day (1 ounce) as an attempt to maximize the contrast between the two groups, although the majority (>70%) consumed more than 56 g/day (2 ounces). Viable plasma samples from a total of 46 non-vegetarians and 47 vegans were included in the current study for metabolomics profiling. The procedures followed were in accordance with the ethical standards of Loma Linda University and approved by the institutional review board for research involving human subjects.

2.2. Metabolomics Profiling

Ninety-three heparin plasma samples were analyzed by Metabolon Inc. (Durham, NC, USA). Samples were profiled using a Global Assay (DiscoveryHD4) (Metabolon Inc., Durham, NC, USA) that provides a comprehensive picture of metabolites within central carbon metabolism and various biochemical classes/pathways including amino acid, nucleotide, carbohydrate, lipid, xenobiotic, microbial, and others. The Metabolon platform yielded 1017 compounds of known identity.
Samples were prepared using the automated MicroLab STAR® system from Hamilton company (Reno, NV, USA), as described previously [25]. After precipitation of proteins, the resulting extracts were analyzed on four independent ultra-high-performance liquid chromatography-tandem mass spectroscopy (UPLC-MS/MS) platforms. A thorough description of the metabolic platform and quality control procedures have been described previously [25]. Various controls were used in tandem with experimental samples, prior to injection into the mass spectrometer. A cocktail of QC standards was added to each sample to monitor instrument performance and aid in chromatographic alignment. Aliquots of a pooled plasma sample (obtained by taking a small amount of each experimental sample) were included as technical replicates. Extracted water samples served as process blanks. Experimental samples were randomized across the platform run with QC samples spaced evenly among the injections. Biochemical compounds were identified by comparison to library entries of purified, authenticated standards or recurrent unknown entities. The median relative standard deviation (RSD) was calculated for the standards added to each sample as a measure of instrument variability, as well as for all endogenous metabolites within the pooled plasma samples (non-instrument standards) as a measure of overall process variability. The median RSD for internal standards and endogenous metabolites was 6% and 14%, respectively. Samples were measured in one batch and randomized by diet group.

2.3. Statistical Analysis

2.3.1. Data Transformation and Linear Regression

Raw metabolite values were median scaled (divided by the median value), and missing metabolite values (below the detection threshold) were imputed with the minimum value for a given metabolite. Data was subsequently log transformed. Metabolites that were below the detection limit for >50% vegans and omnivores were excluded from the analysis, yielding an analytical set of 930 metabolites.
Linear regression models were generated to determine if individual plasma metabolites differed between vegans and non-vegetarians, using smart surrogate variable analysis (SmartSVA). With this approach, regression residuals, obtained from regression of log-transformed metabolite abundance (response) on diet group and covariates—age at blood collection, sex (male vs. female), race (Black vs. White), and BMI (continuous)—were used to obtain surrogate variables representing other unknown, unwanted sources of variation. A linear model was then fitted where the dependent and independent variables were the residuals obtained from regressing metabolite abundance, and dietary pattern, respectively, on surrogate variables and the other covariates, excluding BMI, which was considered as a mediating variable. Models including energy intake, which differed between vegans and non-vegetarians at baseline, did not show appreciable changes in results, so energy intake was not included in final models. The resulting univariate beta coefficients and adjusted predicted means with 95% confidence intervals were then obtained for each metabolite. The adjusted means estimate log (geometric means of untransformed data), and the difference in these between vegans to non-vegetarians, are estimates of log (fold change) metabolite abundance. Linear regression models were also generated without SmartSVA, regressing metabolite abundance on dietary pattern, adjusting only for age, sex, and race.
An analysis of differential abundance of metabolites at the subclass level used the same approach, but operated on the average of the numerators of the individual metabolite t statistics, i.e., regression beta coefficients that contributed to that subclass. Composite t-statistics were produced by dividing by the standard deviations of the averaged numerators taking account of the covariances between subclass members. Each component metabolite had been measured as a multiple of its own median, then log-transformed. The differences between predicted means of subclass members in vegans compared to non-vegetarians were calculated. These estimated logs (the ratio of composite geometric means of untransformed values) were finally exponentiated to produce estimated fold changes. All statistical analyses were performed in R Statistical Software (version 4.0.2; R Core Team, Vienna, Austria).

2.3.2. Adjustment for False Discovery

The partial t-statistics from the regressions described above correspond to differences in metabolite abundance according to dietary pattern. To control for multiple testing, an adapted Storey et al. [26] permutation approach was used. The residualized dietary pattern variables were permuted as a means of defining the null distribution of the t scores for metabolite abundance [27], thereby retaining covariances between residualized metabolite abundances. Estimating the proportion of null metabolites allows an estimate of the false discovery rate (FDR) avoiding the over-conservative Benjamini–Hochberg approach [26] and the consequent selection of metabolites with small FDR. An identical procedure was also used at the subclass level.

2.3.3. Principal Components Analysis

Principal components analysis (PCA) was performed using the FactoMineR package in R on 930 log-transformed metabolites as a dimension reduction approach to identify principal components (PCs) or axes that maximize/explain the variation in metabolite abundance. This was done to test the hypothesis that a vegan or non-vegetarian diet could be defined by select groups of metabolites. A PCA plot was generated by obtaining individual scores (coordinates) of vegans and non-vegetarians for top PCs (PC1 and PC3) in 92 subjects, excluding one outlier. Ten PCs explained 50% of the variance. These were retained for regression analysis (eigenvalues >20), where associations of the 12 PCs with diet group and other dietary covariates (energy-adjusted) of interest were examined. Partial correlation coefficients between PCs and diet group or other variables of interest were also obtained (using ppcor package in R). For the PCA regression, variables with missing values (BMI, n = 2; Kcal, n = 1) were imputed with the mean value. Any missing dietary data was handled using multiple imputation with appropriate standard errors. Metabolites with loadings of 0.5 or higher (representing metabolite correlations with PCs) were identified for 10 components.

2.3.4. Random Forest Analysis

Random forest analysis was also used as a supervised approach for classifying metabolites to identify the most informative metabolites distinguishing vegans from non-vegetarians. All 930 log transformed metabolites were considered in the analysis, with each of 50,000 trees learning from a random sample of fifty percent of all the data (23 vegans and 23 non-vegetarians) without replacement, and the remaining data (representing the out of bag variables) passed down the tree for class prediction to calculate the out of bag (OOB) error. Mean decrease in accuracy was calculated by randomly permuting a variable (metabolite), and subsequently passing the data down the trees for re-assessment of class prediction. Hence, the most influential metabolites were determined after permuting each predictor variable and measuring the change or decrease in predictive accuracy [28].

2.3.5. Bootstrap Regression

Our bootstrap procedure chose 930 subjects with replacement, and then in each of 100 such choices performed the FDR analysis. At the metabolite level, 129 of the 930 metabolites were significant (FDR < 0.05) in at least 90% of all bootstrap samples (44 significant 100% of samples). At the subclass level, there were 18 significant (FDR < 0.05) at least 90% of the time, and 9 significant in 100% of the samples. The metabolites and subclasses are listed in the Supplementary Materials (Tables S16–S19).
While the FDR < 0.05, should be a relatively unbiased estimator of the fact that only 5% of such metabolites will be erroneously selected as significant, it gives no information about truly differential metabolites that by chance had a t-score that just missed the FDR cut-point (Type II error). This was assessed by identifying metabolites significant in >50% of bootstrap samples. The bootstrap result provided added information by tending to identify metabolites that may be part of the 5% of false positives (they would show significance in few samples), while also identifying some metabolites that by chance missed significance in the parent sample, but achieved it in many of the bootstrap samples.

2.3.6. Cross-Validation

A separate consideration is that the regressions run to adjust for confounding are somewhat overoptimistic in the t-scores that they produce for the beta coefficients, as likelihoods to some extent are being maximized to also reflect random idiosyncrasies of our sample. This can be largely overcome by a cross-validation procedure, as follows: Randomly, divide the participant sample to K parts. Excluding one part at a time, develop a regression model using subjects in the remaining K-1 parts, which is then used to predict new predicted dependent values (Y(p)) of all metabolites for subjects in the Kth partition. Finally, all subjects receive such a set of Y(p) values. New improved estimates of the residual regression variances are then Sum(Y−Y(p))2/N for each metabolite, and these are always larger than those estimated from the total sample. Noting that var(beta) for a particular metabolite in the full sample equals residual variance/(N.var(X)), a new somewhat smaller t score (beta/sqrt(Var(beta)) is calculated using the improved estimates of residual variances, and these smaller t-scores are submitted to the FDR procedure.
For this study, K = 10 partitions were used, and the increase in residual variances, and corresponding changes in t-scores was small. Very few metabolites lost significance as compared to the original results (see Supplementary Materials, Tables S14–S15).

3. Results

3.1. Baseline Characteristics

Plasma metabolomics profiling was performed on 93 participants—47 vegans and 46 non-vegetarians of the AHS-2 cohort. Roughly equal numbers of male and female and African American and Caucasian participants were included in this study and balanced among the two diet groups. Hence, there were no significant differences at baseline in race or sex comparing these two groups (Table 1). However, vegans were older than non-vegetarians (mean 66.5 vs. 60.8 years), and BMI was significantly higher in non-vegetarians (31.3 vs. 24.7), consistent with greater mean dietary kcals/day. Intakes of select foods or nutrients differed greatly by diet group, as non-vegetarians had significantly higher intakes of red meat, total meat, poultry, fish, dairy, and saturated fat (p < 0.001), and vegans had higher intakes of fiber, fruit, vegetables, soy, legumes (p < 0.001), and whole grains (p < 0.022). Significant differences in all lifestyle factors were seen when comparing vegans with non-vegetarians, most notably coffee drinking, where 39% of non-vegetarians consumed coffee (once or more per month) compared to 0% of vegans (p < 0.001). A significantly greater proportion of non-vegetarians also had a history of smoking (p = 0.008), alcohol drinking (p = 0.015), and used aspirin or non-steroidal anti-inflammatory drugs (NSAIDS) (p = 0.038), while the number of minutes of exercise per week for vegans was more than twice as long as non-vegetarians (132 vs. 64 min/wk; p = 0.003).

3.2. Linear Regression—Abundance of Metabolites or Biochemical Subclasses in Vegans and Non-Vegetarians

3.2.1. Linear Regression to Analyze Abundance of Individual Metabolites

Linear regression of individual metabolites with SmartSVA yielded a total of 586 differential metabolites after adjustment for false discovery. The top 40 metabolites present at higher abundance in vegans relative to non-vegetarians ordered by fold change are shown in Table 2. Fold changes of metabolites comparing vegans with non-vegetarians were derived from geometric mean ratios, based on adjusted means. Differences up to nearly 7-fold were observed. Several top metabolites with fold changes above 2-fold were lipids and xenobiotics. In addition, among metabolites showing greatest abundance in vegans were compounds within amino acid subclasses—including methionine/cysteine metabolism, urea cycle, tryptophan, and tyrosine metabolism, besides lipid subclasses including dicarboxylic acids and bile acid metabolism, among others (Table 2).
However, the vast majority of differential metabolites differentially abundant at FDR <0.05 (422/586 = 72%) were decreased in vegans. Metabolites showing the most marked decreases in vegans were again prominently xenobiotics, followed by lipids and amino acids (Table 3). This included metabolites from primarily xanthine, histidine, food component, dicarboxylic acid, and drug metabolism, where values were ~3 to ~25-fold higher in the non-vegetarians. Differentially abundant (positively and negatively associated) metabolites belonging to all major classes were identified (Supplementary Materials, Tables S1–S5).
When surrogate variables were not included in models, a smaller number of metabolites (346) were found to be differentially abundant in vegans relative to non-vegetarians (Supplementary Materials, Tables S6–S11). There was, however, considerable overlap between differentially abundant metabolites (FDR < 0.05) when comparing the regression approaches with and without Smart SVA, as the vast majority of differential metabolites (FDR < 0.05) detected in the regression excluding surrogate variables were identified with inclusion of the SmartSVA approach in the analysis. Metabolites showing greatest fold changes were in strong agreement comparing the two approaches (Table 2 and Table 3, Supplementary Materials Tables S6 and S7).

3.2.2. Linear Regression to Analyze Metabolite Subclasses

Regression analysis of metabolite subclasses with at least two component metabolites identified 50 differentially abundant subclasses at FDR < 0.05 (Table 4) of 93 total subclasses. Subclasses with increased abundance in vegans prominently included ketone bodies, followed by vitamin A metabolism, inositol, fatty acid acyl glycine metabolism, lactosylceramides, and benzoate metabolism subclasses. For each of these subclasses, all or the majority of the component metabolites significant at FDR < 0.05 were positively associated with a vegan dietary pattern. The directionality of the significant metabolites is notable considering the analysis included all metabolites in a subclass that were represented on the panel. Hence, differential abundance of subclasses, and directionality, were highly driven by component metabolites which had reached statistical significance in linear regression analysis.
Similar to findings of differential abundance of individual metabolites, the majority of subclasses showing statistical differences between vegans and non-vegetarians were negatively associated with vegans (Table 4). The subclasses with the greatest negative associations in vegans (>1.5-fold change) included xanthine metabolism, drug, branched fatty acid, histidine metabolism, acetylated peptides, ceramides, and dihydroceramides, where all or the vast majority of the statistically significant component metabolites were in decreased abundance. Additionally, all or the vast majority of statistically significant metabolites within many other negatively associated subclasses were also inversely associated with a vegan diet. These included subclasses representing long-chain acyl carnitine metabolism, long-chain saturated fatty acids, lysoplasmalogens, phenylalanines, long-chain monounsaturated fatty acids, monoacylglycerols, and leucine, isoleucine and valine, which are branched chain amino acids (BCAA), among other subclasses. Very similar results were obtained in the analysis excluding surrogate variables, where the top subclasses showing the greatest positive or negative associations with vegans were represented (Supplementary Materials, Table S12) (Names of individual, statistically significant metabolites in each of these differentially abundant subclasses are listed in Supplementary Materials, Table S13).

3.2.3. Cross-Validation and Bootstrapped Regression for Error Analysis

Linear regression with cross-validation was performed to identify possible metabolites or subclasses falsely rejected as null or non-null. Cross validation yielded results very similar to those obtained with the entire sample, with only 15 metabolites originally found to be differential at FDR < 0.05 not showing significance with cross-validation (type I error) (Supplementary Materials, Table S14), and no metabolites showing significance with cross-validation that were not identified in the analysis of the entire sample (type II error). All 50 subclasses found to be differential with analysis of the entire sample were also significantly differential with cross-validation, with potentially three additional subclasses that were not identified in analysis of the entire sample (Supplementary Materials, Table S15). Bootstrapped regression models also very much coincided with results of the non-bootstrapped regression analysis, with no more than three potentially new, non-null metabolites or subclasses identified (Supplementary Materials, Table S16). The majority of bootstrapped samples showed numbers of differential metabolites or metabolite subclasses comparable to those obtained in the non-bootstrapped analysis, and all metabolites or subclasses that were differential in at least 90% of bootstrapped samples were also differential in the non-bootstrapped regression analysis (Supplementary Materials, Tables S17–S19).

3.3. Random Forest and Principal Components Analyses

3.3.1. Random Forest Analysis for Classification by Diet Group

Random forest analysis was used to determine the ability of metabolites to identify the vegan and non-vegetarian dietary classes and to identify metabolites with the greatest predictive accuracy. The highest ranked metabolites with the ability to distinguish vegans from non-vegetarians (with greatest mean decrease accuracy determined by the out-of-bag error/permutation method) included 3-bromo-5-chloro,2-6dihydroxybenzoic acid, 3-methylhistidine, 1-methyl-5-imidazoleacetate, (14 or 15)-methyl palmitate (a17:0 or i17:0), n,n,n-trimethyl-5-aminovalerate, and sphingomyelin (d18:1/17:0, d17:1/18:0, d19:1/16:0) (Figure 1). Classification by diet group showed a predictive accuracy of 92.5%, with a misclassification (out of bag) error of 7.5%. These influential metabolites overlapped with metabolites determined to be highly statistically significant or have large fold changes from the linear regression analysis.

3.3.2. Principal Component Analysis for Dimension Reduction

Principal component analysis was used to collapse data into orthogonal components and generate clusters of correlated metabolites. Nine hundred thirty metabolites were collapsed into 91 PCs explaining 100% of the variance, for 92 subjects (with exclusion of one outlier). Ten PCs with eigenvalues >20 explained 50% of the variance, with eight of these PCs containing metabolites with loadings >0.5. Metabolites driving PC1 included predominantly long-chain acylcarnitines and ceramides, followed by branched fatty acids and long-chain saturated fatty acids (Supplementary Materials, Figure S2 and Table S20). Most prominent subclasses represented in PC3 included long-chain polyunsaturated fatty acids, and dicarboxylic acids, and there was representation of other types of fatty acids (monohydroxy fatty acids, long-chain monounsaturated fatty acids). Lysophospholipids were the predominant subclass represented in PC4 (Supplementary Materials, Table S20). The top four PCs accounted for nearly 1/3 of the variance, with PC1 and PC3 most clearly separating dietary groups. (Figure 2).
Linear regression was performed to examine associations of each of the 10 PCs with various dietary variables of interest, and correlations determined. Significantly associated PCs (p ≤ 0.05) are shown in Table 5. PC1, PC3, and PC4 were associated with a vegan diet (combined r = −0.5, p = 7.9 × 10−7). These PCs were also associated with consumption of red meat, total meat, processed meat, and poultry, along with fiber and saturated fat, with partial correlations ranging from r = +/− 0.41 to 0.61. The first four PCs were highly correlated with consumption of fish (r = 0.52, p = 2.6 × 10−7), and dairy kcals (r = 0.61, p = 6.4 × 10−10). Correlations of red meat, total meat, and processed meat with these PCs were largely attenuated in expanded models including additional dietary covariates. Inclusion of poultry, fish, dairy, and fruit in the model attenuated associations between red meat or total meat and PCs, and an association remained only with PC1 for processed meat. After adjustment for saturated fat and fiber, significant associations remained for red meat and total meat with PC4 (β = 1.7, p = 0.009; β = 1.47, p = 0.008) (Supplementary Materials, Table S21).

4. Discussion

Among participants in the AHS-2 cohort, distinct metabolic profiles for vegans and non-vegetarians were discovered, with over 60% of metabolites being significantly discriminatory after adjustment for false discovery. Clearly, the serological characteristics of vegans and non-vegetarians differ substantially. The most notable metabolites more abundant in vegans belonged to categories mostly related to plant-food intakes. Those more abundant in non-vegetarians included subclasses of lipids and amino acids, which are mostly related to intakes of animal foods, besides xenobiotics reflecting other lifestyle behaviors such as caffeine consumption and medication use.
Prominent among metabolites more abundant in vegans were products of benzoate metabolism derived from polyphenols in plant foods, possibly also reflecting gut microbial activity. For example, the top compound, 4-ethylphenyl sulfate, is generated by the metabolism of soy protein by gut bacteria [29,30,31], and other metabolites may be generated through microbial metabolism of dietary polyphenols (hippurate metabolites, catechol sulfate, and others). Higher abundance of other food component/plant metabolites (glucopyranoside metabolites, stachydrine), besides vitamins in vegans, similarly reflects consumption of fruits, vegetables, and herbs [32,33], which have roles in reducing risk of chronic diseases. These findings in vegans are consistent with those obtained by Wu et al., where increases were found in benzoate metabolism products and polyphenolic plant compounds as well as gut microbial metabolites (hippurate, 4-ethylphenyl sulfate, 4-hydoxyhippurate, catechol sulfate, phenol sulfate). Notably, vegans showed a significant increase (30%) in butyrate, a short-chain fatty acid that is generated with increased fermentation of nondigestible carbohydrates or dietary fiber, and has a role in regulating inflammation and epithelial barrier function.
Other metabolites that were more abundant in vegan serum have relevance to lipid metabolism and metabolic homeostasis. The cause of the observed higher abundance of ketone bodies in vegans is not clear but may be related to the length of the overnight fast, possibly longer in vegans, caloric restriction, or exercise, which was more frequent in vegans. The acyl glycine subclass of lipids also showed higher abundance in vegans. Acyl glycines are metabolites of fatty acids with important roles in lipid signaling, some with anti-inflammatory ability [34]. There is evidence of negative regulation of acyl glycines by branched chain amino acids [35] through lowering of glycine. Hence, the lower abundance of branched chain fatty acids in vegans in our study might explain in part the higher abundance of acyl glycines, and of primary glycine-conjugated bile acids (glycohyocholate, glycochenodeoxycholate, glycohyocholate, glycol-beta-muricholate.).
The biological and pathophysiologic effects of these differences are several. Carotenoids, associated with vitamin A, and other polyphenols as well as microbial metabolites produced during breakdown of dietary fiber, have anti-inflammatory and antioxidant properties. These compounds counteract oxidative stress and support immune function and gut health to prevent cancer, diabetes, cardiovascular and other diseases. This may happen through inhibition of nuclear factor of activated B-cells (NFkB), regulation of inflammatory cytokines through epigenetic modifications, and increased transcription of antioxidant defense and xenobiotic detoxification genes [36,37,38]. Ketone bodies have beneficial roles in energy metabolism and glucose homeostasis, and thus may prevent or counteract inflammation and oxidative stress [39,40]. Glycine levels have also shown inverse associations with cardiometabolic disease phenotypes [41,42,43,44]. Higher abundance of glycine, along with glycine-conjugated bile acids, has been reported in other cross-sectional studies of vegans and vegetarians [13,45,46,47,48].
The vast majority of metabolites showed lower abundance in vegan serum most often reflecting the absence of animal products in the diet, or the much lower intakes of caffeine and use of medications. Many metabolites within the xanthine metabolism subclass, significantly different in vegans and non-vegetarians, are metabolites of caffeine—theophylline, paraxanthine, 5-acetylamino-6-amino-3-methyluracil, and several others, and would reflect the greater coffee consumption (or perhaps medication use) in non-vegetarians. Other xenobiotic metabolites present at significantly lower abundance in vegans included chemical and drug metabolites, which can be explained by the increased use of acetaminophen and other NSAIDs in non-vegetarians. The large (3-fold) decrease in the acetylated peptides (phenylacetylglutamine, phenylacetyl carnitine) in vegans is likely a consequence of metabolism of phenylalanine (converted to phenylacetate), along with glutamine or L-carnitine, which may increase with higher intake of animal meats, protein, or certain pharmaceutics (NSAIDs). More notable was the high abundance of a large number of histidine metabolites and branched chain amino acids—isoleucine, leucine, and valine—in non-vegetarians, likely reflecting dietary consumption of meat and animal products [49,50,51]. 1- and 3-methylhistidine, for example, are biomarkers of skeletal muscle protein breakdown [51,52]. Histidine and other imidazole-containing compounds may also be related to use of pharmaceuticals [49].
The markedly lower abundance of various lipid metabolites in vegans such as ceramides and dihydroceramides, long-chain acylcarnitines, long-chain saturated fatty acids, monoacylglycerols, and branched chain fatty acids may reflect reduced intake of saturated fats and lipids/sphingolipids derived from animal sources in the diet. Only lactosylceramides were increased in vegans, possibly reflecting consumption of the precursor glucosylceramide from dietary plant sources (i.e., soy, wheat) [52].
Our findings are largely consistent with previous studies examining metabolic profiles associated with vegetarian/vegan, and plant-based diets. Decreases in sphingolipids and some acylcarnitines have been observed in vegans relative to omnivores in a cross-sectional study of individuals in the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort following habitual dietary patterns [13]. In another study of individuals following a vegan diet for at least six months, decreases in phospholipids, saturated fatty acids, 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF), and methylhistidine, but increases in plant-derived and microbial metabolites, were observed [14]. Further, the Mediterranean dietary pattern has been associated with alterations in acylcarnitines and phospholipids [10]. Additionally, ceramides, plasmalogens and acylcarnitines were positively associated with a Western dietary pattern in a cross-sectional study examining associations of plasma metabolites with Western or Prudent dietary patterns in the Women’s Health Initiative cohort [8].
There are important biological and pathophysiologic implications of these differences. While histidine may have some beneficial health effects (i.e., antioxidant activity), histidine metabolites and related imidazole derivatives may be associated with impaired insulin signaling, type 2 diabetes, and kidney disease [53,54]. Several cross-sectional studies have shown positive associations of branched chain amino acids with cardiovascular disease risk, metabolic dysregulation, impaired glucose signaling and insulin resistance [55,56,57], and a number of studies have shown positive associations with type 2 diabetes [58,59,60]. Ceramides are critical lipid signaling molecules, and have various and complex biological roles. Both ceramides and sphingolipids may play a role in insulin resistance [61,62,63]. Ceramides have been found to accumulate in tissues of obese individuals, and have been associated with inflammation [64], and also have roles in the regulation of apoptosis and development of neurological disorders [65,66,67]. Lactosylceramides may have a role in promoting innate immunity [68], besides other biological functions. Phenylacetylglutamine, among others of these metabolites, may reflect activity of gut microbiota, and has been associated with cardiovascular disease risk [69,70,71]. Acylcarnitines are formed in the mitochondria during beta oxidation of fatty acids, but increases in plasma may reflect metabolic disorders [72], as these compounds have been associated with insulin resistance, diabetes, and increased cardiovascular disease risk. Similarly, long-chain saturated fatty acids are implicated in increased risk of obesity and cardiometabolic diseases [73]. Mono-acylglycerols are converted to triacylglycerols, which are associated with heart disease, obesity, and metabolic syndrome [74]. Branched chain fatty acids, on the other hand, which are derived largely from dairy and meat products (though also synthesized by gut bacteria from branched chain amino acids in herbivores [75]), may favorably influence metabolic health. Accordingly, they have potentially beneficial effects on insulin sensitivity and weight management, and may attenuate inflammation [76,77,78], but further human studies are needed.
Vegetarians and particularly vegans have notably higher consumption of fruits, vegetables, legumes, whole grains, soy foods, and nuts, and markedly reduced or absent intake of animal products, as demonstrated previously in the AHS-2 cohort [25]. Additionally, compared to non-vegetarians, they have more favorable fatty acid profiles, including lower saturated fatty acids, and higher total omega-3, along with higher levels of phytochemicals such as carotenoids, enterolactone, and isoflavones in plasma or urine [15], but lower inflammatory cytokines [79]. Importantly, vegetarians have shown significantly reduced risks of diabetes, hypertension, cardiovascular disease, select cancers, and all-cause mortality relative to non-vegetarians [1,2,3,80]. These differences in disease and biomarker profiles between vegans and non-vegetarians coincide with the differences in plasma metabolites identified in the present study.
Strengths of the current study include well-defined diet groups (vegans and non-vegetarians) reflecting habitual dietary patterns [21], and the inclusion of higher meat consuming non-vegetarians. An additional strength is the analytical rigor applied with the use of multiple regression and other approaches, including adjustment for surrogate variables to address additional confounding and unwanted variation, besides classification and dimensionality reduction approaches. Further, adjustment for false discovery with use of the adapted Storey et al. [26] method allowed for an accurate estimation of the proportion of t-statistics from a non-null distribution, which provides an improvement in power compared to other approaches. AHS-2 participants have provided extensive dietary data, besides demographic and medical data, strengthening the analysis of the current study. Limitations of the study are the single measurement of plasma metabolites, and the somewhat limited sample size, although there was sufficient power to detect a large number of statistically significant differences between vegans and non-vegetarians. Estimates of results from an independent sample are obtained by cross-validation. There are notable lifestyle differences in this Adventist cohort compared to other vegan and non-vegetarian individuals. Seventh-day Adventists place an emphasis on health and wellness, and there are religious guidelines prohibiting certain lifestyle behaviors (smoking, alcohol drinking, biblically unclean meats), but no prohibitions for clean meats in general or dairy. These healthy practices, which limit confounding by these factors, may possibly also limit generalizability of results. The overall reduced consumption of non-fish, and particularly red meats among AHS-2 non-vegetarians relative to the general population, was a limitation, although this was partially overcome by selection of subjects with higher meat consumption. As many diet-derived metabolites are converted by gut bacteria, it remains to be understood how differences in these metabolites might be linked with alterations in composition of the gut microbiome.

5. Conclusions

In conclusion, in this study we report marked differences in metabolic profiles between vegans and non-vegetarians. Our results suggest that multiple potentially bioactive metabolites are increased by consumption of plant-based foods, and may lower the risk of metabolic diseases through anti-inflammatory mechanisms. On the other hand, diets high in animal products may lead to increases in various amino acids and lipid species (acyl carnitines, saturated fatty acids, ceramides, branched chain amino acids) that promote chronic diseases by increasing inflammation and insulin dysregulation, so disrupting metabolic homeostasis. The exact roles or physiological functions of other differentially abundant metabolites are not clear. It may be that some differentially abundant metabolites in vegans and non-vegetarians serve only as markers of different foods or eating patterns, while others also have important pathophysiological consequences. This study helps lay the foundation for a deeper understanding of the relationship of diet-associated metabolites to the pathophysiology of chronic diseases.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu14030709/s1; Table S1. Amino acid metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression models with SmartSVA approach; Table S2: Lipid metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression models with SmartSVA approach; Table S3: Carbohydrate, cofactor/vitamin, and energy metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression models with SmartSVA approach; Table S4: Nucleotides, partially characterized molecules, and peptides associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression models with SmartSVA approach; Table S5: Xenobiotic metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression models with SmartSVA approach; Table S6: Top 40 metabolites positively associated with a vegan dietary pattern at FDR < 0.05 in linear regression analysis without SmartSVA; Table S7: Top 40 metabolites inversely associated with a vegan dietary pattern at FDR < 0.05 in linear regression analysis without SmartSVA; Table S8: Amino acid metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression analysis without SmartSVA; Table S9: Lipid metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression analysis without SmartSVA; Table S10: Carbohydrate, cofactors/vitamins, energy, nucleotide, partially characterized, and peptide metabolites associated with a vegan dietary pattern in linear regression analysis without SmartSVA; Table S11: Xenobiotic metabolites associated with a vegan (relative to non-vegetarian) dietary pattern at FDR < 0.05 in linear regression analysis without SmartSVA; Table S12: Metabolite subclasses associated with diet group (vegan vs. non-vegetarian) at FDR < 0.05 without SmartSVA; Table S13: Component metabolites of each subclass associated with a vegan (relative to non-vegetarian dietary pattern at FDR < 0.05; Table S14: Metabolites differentially abundant (FDR < 0.05) in linear regression analysis with full sample that were nondifferential post cross-validation; Table S15: Subclasses differential after cross-validation that were not differential in regression analysis with entire sample; Table S16: Metabolites or metabolite subclasses showing differential abundance between vegans and non-vegetarians (FDR < 0.05) in > 50% of bootstrapped linear regression analyses, but not differential in non-bootstrapped regression analysis; Table S17: Numbers of differential metabolites or metabolite subclasses (at FDR < 0.05) in regression analysis with bootstrap sampling; Table S18: List of 129 metabolites showing differential abundance (FDR < 0.05) in at least 90% of bootstrap regressions; Table S19: Metabolite subclasses showing differential abundance (FDR < 0.05) in at least 90% of bootstrap regressions; Table S20: Top components from principal components analysis and most influential metabolites; Table S21:Adjusted linear regression predicting red meat, processed, and total meat consumption from top principal com-ponents from principal components analysis, with adjustment for additional potential dietary confounders; Figure S1. Study design for individuals in the Metabolomics Pilot Study. Footnotes: (a) The Calibration sub-study was a random sample of the cohort, except for an overweighting of Black subjects so they formed 40% of the total. (b) The two pilot sub-studies were convenience samples of study subjects living in Texas (Black subjects) or Washington State designed to test our bio-sample acquisition strategies. (c) The BioMRS sub-study was a local sample of AHS-2 subjects who had responded to a request to complete a questionnaire containing psychosocial and religiosity questions. They lived within 50 miles of Loma Linda, Riverside, or downtown Los Angeles and were at least 50 years of age, Figure S2. Subclass loadings from first principal component (PC1). Principal component analysis was used to identify components explaining variation in metabolites comparing vegans and non-vegetarians. Metabolite loadings > 0.5 were extracted and averaged across each represented subclass.

Author Contributions

The authors’ contributions were as follows: F.L.M., G.E.F., M.J.O. and P.D.-H. designed research; G.E.F., F.L.M. and A.M. analyzed data, and P.D.C. and J.W.L. provided input regarding the analysis; F.L.M. and G.E.F. wrote the paper; F.L.M. had primary responsibility for final content. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Loma Linda University Health through a pilot grant given to PDH in support of AHS-2 research (Basic Science and Translational Research Pilot Funds). Additionally, the research was supported by the Ardmore Institute of Health (“The Adventist Health Study: The next generation in transformational health knowledge”), and NIH National Institute on Minority Health and Health Disparities K01MD015194.

Institutional Review Board Statement

The current study was conducted in accordance with the Declaration of Helsinki, and the protocol was reviewed by the Ethics Committee of Loma Linda University Health but determined to be exempt as outlined in federal regulations for protection of human subjects, 45 CFR Part 46.101(b)(4) XXX (5190502), as the study utilized biospecimens collected in a previously approved study, protocol #48134 on 5/28/03.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. All subjects gave their informed consent for inclusion before they participated in the original study approved by the Loma Linda University institutional review board, protocol #48134 in which plasma samples were obtained for the current study.

Data Availability Statement

Data supporting the results of this study will be available upon request, once reviewed and approved by the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rizzo, N.S.; Sabate, J.; Jaceldo-Siegl, K.; Fraser, G.E. Vegetarian dietary patterns are associated with a lower risk of metabolic syndrome: The adventist health study 2. Diabetes Care 2011, 34, 1225–1227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Tonstad, S.; Stewart, K.; Oda, K.; Batech, M.; Herring, R.P.; Fraser, G.E. Vegetarian diets and incidence of diabetes in the Adventist Health Study-2. Nutr. Metab. Cardiovasc. Dis. 2013, 23, 292–299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Orlich, M.J.; Singh, P.N.; Sabate, J.; Jaceldo-Siegl, K.; Fan, J.; Knutsen, S.; Beeson, W.L.; Fraser, G.E. Vegetarian dietary patterns and mortality in Adventist Health Study 2. JAMA Intern. Med. 2013, 173, 1230–1238. [Google Scholar] [CrossRef] [PubMed]
  4. Orlich, M.J.; Singh, P.N.; Sabate, J.; Fan, J.; Sveen, L.; Bennett, H.; Knutsen, S.F.; Beeson, W.L.; Jaceldo-Siegl, K.; Butler, T.L.; et al. Vegetarian dietary patterns and the risk of colorectal cancers. JAMA Intern. Med. 2015, 175, 767–776. [Google Scholar] [CrossRef] [PubMed]
  5. Tantamango-Bartley, Y.; Jaceldo-Siegl, K.; Fan, J.; Fraser, G. Vegetarian diets and the incidence of cancer in a low-risk population. Cancer Epidemiol. Biomarkers Prev. 2013, 22, 286–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Tantamango-Bartley, Y.; Knutsen, S.F.; Knutsen, R.; Jacobsen, B.K.; Fan, J.; Beeson, W.L.; Sabate, J.; Hadley, D.; Jaceldo-Siegl, K.; Penniecook, J.; et al. Are strict vegetarians protected against prostate cancer? Am. J. Clin. Nutr. 2016, 103, 153–160. [Google Scholar] [CrossRef] [Green Version]
  7. Bouchard-Mercier, A.; Rudkowska, I.; Lemieux, S.; Couture, P.; Vohl, M.C. The metabolic signature associated with the Western dietary pattern: A cross-sectional study. Nutr. J. 2013, 12, 158. [Google Scholar] [CrossRef] [Green Version]
  8. Chandler, P.D.; Balasubramanian, R.; Paynter, N.; Giulianini, F.; Fung, T.; Tinker, L.F.; Snetselaar, L.; Liu, S.M.; Eaton, C.; Tobias, D.K.; et al. Metabolic signatures associated with Western and Prudent dietary patterns in women. Am. J. Clin. Nutr. 2020, 112, 268–283. [Google Scholar] [CrossRef]
  9. Playdon, M.C.; Moore, S.C.; Derkach, A.; Reedy, J.; Subar, A.F.; Sampson, J.N.; Albanes, D.; Gu, F.; Kontto, J.; Lassale, C.; et al. Identifying biomarkers of dietary patterns by using metabolomics. Am. J. Clin. Nutr. 2017, 105, 450–465. [Google Scholar] [CrossRef] [Green Version]
  10. Tong, T.Y.N.; Koulman, A.; Griffin, J.L.; Wareham, N.J.; Forouhi, N.G.; Imamura, F. A Combination of Metabolites Predicts Adherence to the Mediterranean Diet Pattern and Its Associations with Insulin Sensitivity and Lipid Homeostasis in the General Population: The Fenland Study, United Kingdom. J. Nutr. 2020, 150, 568–578. [Google Scholar] [CrossRef] [Green Version]
  11. Vazquez-Fresno, R.; Llorach, R.; Urpi-Sarda, M.; Lupianez-Barbero, A.; Estruch, R.; Corella, D.; Fito, M.; Aros, F.; Ruiz-Canela, M.; Salas-Salvado, J.; et al. Metabolomic pattern analysis after mediterranean diet intervention in a nondiabetic population: A 1- and 3-year follow-up in the PREDIMED study. J. Proteome Res. 2015, 14, 531–540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Garcia-Perez, I.; Posma, J.M.; Gibson, R.; Chambers, E.S.; Hansen, T.H.; Vestergaard, H.; Hansen, T.; Beckmann, M.; Pedersen, O.; Elliott, P.; et al. Objective assessment of dietary patterns by use of metabolic phenotyping: A randomised, controlled, crossover trial. Lancet Diabetes Endocrinol. 2017, 5, 184–195. [Google Scholar] [CrossRef] [Green Version]
  13. Schmidt, J.A.; Rinaldi, S.; Ferrari, P.; Carayol, M.; Achaintre, D.; Scalbert, A.; Cross, A.J.; Gunter, M.J.; Fensom, G.K.; Appleby, P.N.; et al. Metabolic profiles of male meat eaters, fish eaters, vegetarians, and vegans from the EPIC-Oxford cohort. Am. J. Clin. Nutr. 2015, 102, 1518–1526. [Google Scholar] [CrossRef] [Green Version]
  14. Wu, G.D.; Compher, C.; Chen, E.Z.; Smith, S.A.; Shah, R.D.; Bittinger, K.; Chehoud, C.; Albenberg, L.G.; Nessel, L.; Gilroy, E.; et al. Comparative metabolomics in vegans and omnivores reveal constraints on diet-dependent gut microbiota metabolite production. Gut 2016, 65, 63–72. [Google Scholar] [CrossRef] [PubMed]
  15. Miles, F.L.; Lloren, J.I.C.; Haddad, E.; Jaceldo-Siegl, K.; Knutsen, S.; Sabate, J.; Fraser, G.E. Plasma, Urine, and Adipose Tissue Biomarkers of Dietary Intake Differ Between Vegetarian and Non-Vegetarian Diet Groups in the Adventist Health Study-2. J. Nutr. 2019, 149, 667–675. [Google Scholar] [CrossRef]
  16. Miles, F.L.; Mashchak, A.; Filippov, V.; Orlich, M.J.; Duerksen-Hughes, P.; Chen, X.; Wang, C.; Siegmund, K.; Fraser, G.E. DNA Methylation Profiles of Vegans and Non-Vegetarians in the Adventist Health Study-2 Cohort. Nutrients 2020, 12, 3697. [Google Scholar] [CrossRef]
  17. Jaceldo-Siegl, K.; Knutsen, S.F.; Sabate, J.; Beeson, W.L.; Chan, J.; Herring, R.P.; Butler, T.L.; Haddad, E.; Bennett, H.; Montgomery, S.; et al. Validation of nutrient intake using an FFQ and repeated 24 h recalls in black and white subjects of the Adventist Health Study-2 (AHS-2). Public Health Nutr. 2010, 13, 812–819. [Google Scholar] [CrossRef] [PubMed]
  18. Schakel, S.F. Maintaining a nutrient database in a changing marketplace: Keeping pace with changing food products—A research perspective. J. Food Compos. Anal. 2001, 14, 315–322. [Google Scholar] [CrossRef]
  19. Fraser, G.E.; Stram, D.O. Regression calibration when foods (measured with error) are the variables of interest: Markedly non-Gaussian data with many zeroes. Am. J. Epidemiol. 2012, 175, 325–331. [Google Scholar] [CrossRef] [Green Version]
  20. Fraser, G.E.; Yan, R. A multivariate method for measurement error correction using pairs of concentration biomarkers. Ann. Epidemiol. 2007, 17, 64–73. [Google Scholar] [CrossRef] [PubMed]
  21. Martins, M.C.T.; Jaceldo-Siegl, K.; Orlich, M.; Fan, J.; Mashchak, A.; Fraser, G.E. A New Approach to Assess Lifetime Dietary Patterns Finds Lower Consumption of Animal Foods with Aging in a Longitudinal Analysis of a Health-Oriented Adventist Population. Nutrients 2017, 9, 1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lee, J.W.; Morton, K.R.; Walters, J.; Bellinger, D.L.; Butler, T.L.; Wilson, C.; Walsh, E.; Ellison, C.G.; McKenzie, M.M.; Fraser, G.E. Cohort profile: The biopsychosocial religion and health study (BRHS). Int. J. Epidemiol. 2009, 38, 1470–1478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chan, J.; Knutsen, S.F.; Sabate, J.; Haddad, E.; Yan, R.; Fraser, G.E. Feasibility of running clinics to collect biological specimens in a nationwide cohort study--Adventist Health Study-2. Ann. Epidemiol. 2007, 17, 454–457. [Google Scholar] [CrossRef]
  24. Orlich, M.J.; Jaceldo-Siegl, K.; Sabate, J.; Fan, J.; Singh, P.N.; Fraser, G.E. Patterns of food consumption among vegetarians and non-vegetarians. Br. J. Nutr. 2014, 112, 1644–1653. [Google Scholar] [CrossRef] [Green Version]
  25. Mesnage, R.; Teixeira, M.; Mandrioli, D.; Falcioni, L.; Ducarmon, Q.R.; Zwittink, R.D.; Mazzacuva, F.; Caldwell, A.; Halket, J.; Amiel, C.; et al. Use of Shotgun Metagenomics and Metabolomics to Evaluate the Impact of Glyphosate or Roundup MON 52276 on the Gut Microbiota and Serum Metabolome of Sprague-Dawley Rats. Environ. Health Perspect. 2021, 129, 17005. [Google Scholar] [CrossRef]
  26. Storey, J.D.; Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. USA 2003, 100, 9440–9445. [Google Scholar] [CrossRef] [Green Version]
  27. Freedman, D.; Lane, D. A Nonstochastic Interpretation of Reported Significance Levels. J. Bus. Econ. Stat. 1983, 1, 292–298. [Google Scholar] [CrossRef]
  28. Hastie, T.; Friedman, J. Random Forests. In The Elements of Statistical Learning, 2nd ed.; Springer Publishers: New York, NY, USA, 2009. [Google Scholar]
  29. Setchell, K.D.; Clerici, C. Equol: History, chemistry, and formation. J. Nutr. 2010, 140, 1355S–1362S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Yang, Z.; Kulkarni, K.; Zhu, W.; Hu, M. Bioavailability and pharmacokinetics of genistein: Mechanistic studies on its ADME. Anticancer Agents Med. Chem. 2012, 12, 1264–1280. [Google Scholar] [CrossRef]
  31. Velenosi, T.J.; Hennop, A.; Feere, D.A.; Tieu, A.; Kucey, A.S.; Kyriacou, P.; McCuaig, L.E.; Nevison, S.E.; Kerr, M.A.; Urquhart, B.L. Untargeted plasma and tissue metabolomics in rats with chronic kidney disease given AST-120. Sci. Rep. 2016, 6, 22526. [Google Scholar] [CrossRef] [Green Version]
  32. Ma, X.; Laaksonen, O.; Heinonen, J.; Sainio, T.; Kallio, H.; Yang, B. Sensory profile of ethyl beta-d-glucopyranoside and its contribution to quality of sea buckthorn (Hippophae rhamnoides L.). Food Chem. 2017, 233, 263–272. [Google Scholar] [CrossRef]
  33. Xie, X.; Yang, C.; Cui, Q.; Ma, W.; Liu, J.; Yao, Q.; Zhang, Z.; Xiao, L.; Wang, N. Stachydrine Mediates Rapid Vascular Relaxation: Activation of Endothelial Nitric Oxide Synthase Involving AMP-Activated Protein Kinase and Akt Phosphorylation in Vascular Endothelial Cells. J. Agric. Food Chem. 2019, 67, 9805–9811. [Google Scholar] [CrossRef]
  34. Burstein, S.; McQuain, C.; Salmonsen, R.; Seicol, B. N-Amino acid linoleoyl conjugates: Anti-inflammatory activities. Bioorg. Med. Chem. Lett. 2012, 22, 872–875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. White, P.J.; Lapworth, A.L.; An, J.; Wang, L.; McGarrah, R.W.; Stevens, R.D.; Ilkayeva, O.; George, T.; Muehlbauer, M.J.; Bain, J.R.; et al. Branched-chain amino acid restriction in Zucker-fatty rats improves muscle insulin sensitivity by enhancing efficiency of fatty acid oxidation and acyl-glycine export. Mol. Metab. 2016, 5, 538–551. [Google Scholar] [CrossRef] [PubMed]
  36. Kaulmann, A.; Bohn, T. Carotenoids, inflammation, and oxidative stress-implications of cellular signaling pathways and relation to chronic disease prevention. Nutr. Res. 2014, 34, 907–929. [Google Scholar] [CrossRef] [PubMed]
  37. Saleh, H.A.; Yousef, M.H.; Abdelnaser, A. The Anti-Inflammatory Properties of Phytochemicals and Their Effects on Epigenetic Mechanisms Involved in TLR4/NF-kappaB-Mediated Inflammation. Front. Immunol. 2021, 12, 606069. [Google Scholar] [CrossRef] [PubMed]
  38. Wu, J.C.; Lai, C.S.; Tsai, M.L.; Ho, C.T.; Wang, Y.J.; Pan, M.H. Chemopreventive effect of natural dietary compounds on xenobiotic-induced toxicity. J. Food Drug Anal. 2017, 25, 176–186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Moller, N. Ketone Body, 3-Hydroxybutyrate: Minor Metabolite—Major Medical Manifestations. J. Clin. Endocrinol. Metab. 2020, 105, 2884–2892. [Google Scholar] [CrossRef]
  40. Puchalska, P.; Crawford, P.A. Multi-dimensional Roles of Ketone Bodies in Fuel Metabolism, Signaling, and Therapeutics. Cell Metab. 2017, 25, 262–284. [Google Scholar] [CrossRef] [Green Version]
  41. Ding, Y.; Svingen, G.F.; Pedersen, E.R.; Gregory, J.F.; Ueland, P.M.; Tell, G.S.; Nygard, O.K. Plasma Glycine and Risk of Acute Myocardial Infarction in Patients With Suspected Stable Angina Pectoris. J. Am. Heart Assoc. 2015, 5, e002621. [Google Scholar] [CrossRef] [Green Version]
  42. Floegel, A.; Stefan, N.; Yu, Z.; Muhlenbruch, K.; Drogan, D.; Joost, H.G.; Fritsche, A.; Haring, H.U.; Hrabe de Angelis, M.; Peters, A.; et al. Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach. Diabetes 2013, 62, 639–648. [Google Scholar] [CrossRef] [Green Version]
  43. Gall, W.E.; Beebe, K.; Lawton, K.A.; Adam, K.P.; Mitchell, M.W.; Nakhle, P.J.; Ryals, J.A.; Milburn, M.V.; Nannipieri, M.; Camastra, S.; et al. alpha-hydroxybutyrate is an early biomarker of insulin resistance and glucose intolerance in a nondiabetic population. PLoS ONE 2010, 5, e10883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Huffman, K.M.; Shah, S.H.; Stevens, R.D.; Bain, J.R.; Muehlbauer, M.; Slentz, C.A.; Tanner, C.J.; Kuchibhatla, M.; Houmard, J.A.; Newgard, C.B.; et al. Relationships Between Circulating Metabolic Intermediates and Insulin Action in Overweight to Obese, Inactive Men and Women. Diabetes Care 2009, 32, 1678–1683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. De Aguiar Vallim, T.Q.; Tarling, E.J.; Edwards, P.A. Pleiotropic roles of bile acids in metabolism. Cell Metab. 2013, 17, 657–669. [Google Scholar] [CrossRef] [Green Version]
  46. Ridlon, J.M.; Kang, D.J.; Hylemon, P.B. Bile salt biotransformations by human intestinal bacteria. J. Lipid Res. 2006, 47, 241–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Trefflich, I.; Marschall, H.U.; Giuseppe, R.D.; Stahlman, M.; Michalsen, A.; Lampen, A.; Abraham, K.; Weikert, C. Associations between Dietary Patterns and Bile Acids-Results from a Cross-Sectional Study in Vegans and Omnivores. Nutrients 2019, 12, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Xu, J.; Yang, S.; Cai, S.; Dong, J.; Li, X.; Chen, Z. Identification of biochemical changes in lactovegetarian urine using 1H NMR spectroscopy and pattern recognition. Anal. Bioanal. Chem. 2010, 396, 1451–1463. [Google Scholar] [CrossRef]
  49. Vessally, E.; Soleimani-Amiri, S.; Hosseinian, A.; Edjlali, L.; Bekhradnia, A. New protocols to access imidazoles and their ring fused analogues: Synthesis from N-propargylamines. Rsc Adv. 2017, 7, 7079–7091. [Google Scholar] [CrossRef] [Green Version]
  50. McKeran, R.O.; Halliday, D.; Purkiss, P.; Royston, P. 3-Methylhistidine excretion as an index of myofibrillar protein catabolism in neuromuscular disease. J. Neurol. Neurosurg. Psychiatry 1979, 42, 536–541. [Google Scholar] [CrossRef] [Green Version]
  51. Myint, T.; Fraser, G.E.; Lindsted, K.D.; Knutsen, S.F.; Hubbard, R.W.; Bennett, H.W. Urinary 1-methylhistidine is a marker of meat consumption in Black and in White California Seventh-day Adventists. Am. J. Epidemiol. 2000, 152, 752–755. [Google Scholar] [CrossRef]
  52. Sullards, M.C.; Lynch, D.V.; Merrill, A.H., Jr.; Adams, J. Structure determination of soybean and wheat glucosylceramides by tandem mass spectrometry. J. Mass. Spectrom. 2000, 35, 347–353. [Google Scholar] [CrossRef]
  53. Raj, D.S.; Ouwendyk, M.; Francoeur, R.; Pierratos, A. Plasma amino acid profile on nocturnal hemodialysis. Blood Purif. 2000, 18, 97–102. [Google Scholar] [CrossRef] [PubMed]
  54. Tuma, P.; Samcova, E.; Balinova, P. Determination of 3-methylhistidine and 1-methylhistidine in untreated urine samples by capillary electrophoresis. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2005, 821, 53–59. [Google Scholar] [CrossRef]
  55. Batch, B.C.; Shah, S.H.; Newgard, C.B.; Turer, C.B.; Haynes, C.; Bain, J.R.; Muehlbauer, M.; Patel, M.J.; Stevens, R.D.; Appel, L.J.; et al. Branched chain amino acids are novel biomarkers for discrimination of metabolic wellness. Metabolism 2013, 62, 961–969. [Google Scholar] [CrossRef] [Green Version]
  56. Weng, L.; Quinlivan, E.; Gong, Y.; Beitelshees, A.L.; Shahin, M.H.; Turner, S.T.; Chapman, A.B.; Gums, J.G.; Johnson, J.A.; Frye, R.F.; et al. Association of branched and aromatic amino acids levels with metabolic syndrome and impaired fasting glucose in hypertensive patients. Metab. Syndr. Relat. Disord. 2015, 13, 195–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yamaguchi, N.; Mahbub, M.H.; Takahashi, H.; Hase, R.; Ishimaru, Y.; Sunagawa, H.; Amano, H.; Kobayashi-Miura, M.; Kanda, H.; Fujita, Y.; et al. Plasma free amino acid profiles evaluate risk of metabolic syndrome, diabetes, dyslipidemia, and hypertension in a large Asian population. Environ. Health Prev. Med. 2017, 22, 35. [Google Scholar] [CrossRef] [Green Version]
  58. Guasch-Ferre, M.; Hruby, A.; Toledo, E.; Clish, C.B.; Martinez-Gonzalez, M.A.; Salas-Salvado, J.; Hu, F.B. Metabolomics in Prediabetes and Diabetes: A Systematic Review and Meta-analysis. Diabetes Care 2016, 39, 833–846. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Wang, T.J.; Larson, M.G.; Vasan, R.S.; Cheng, S.; Rhee, E.P.; McCabe, E.; Lewis, G.D.; Fox, C.S.; Jacques, P.F.; Fernandez, C.; et al. Metabolite profiles and the risk of developing diabetes. Nat. Med. 2011, 17, 448–453. [Google Scholar] [CrossRef] [PubMed]
  60. Wang, W.; Jiang, H.; Zhang, Z.; Duan, W.; Han, T.; Sun, C. Interaction between dietary branched-chain amino acids and genetic risk score on the risk of type 2 diabetes in Chinese. Genes Nutr. 2021, 16, 4. [Google Scholar] [CrossRef] [PubMed]
  61. Lachkar, F.; Ferre, P.; Foufelle, F.; Papaioannou, A. Dihydroceramides: Their emerging physiological roles and functions in cancer and metabolic diseases. Am. J. Physiol. Endocrinol. Metab. 2021, 320, E122–E130. [Google Scholar] [CrossRef]
  62. Meikle, P.J.; Summers, S.A. Sphingolipids and phospholipids in insulin resistance and related metabolic disorders. Nat. Rev. Endocrinol. 2017, 13, 79–91. [Google Scholar] [CrossRef]
  63. Sokolowska, E.; Blachnio-Zabielska, A. The Role of Ceramides in Insulin Resistance. Front. Endocrinol 2019, 10, 577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Haus, J.M.; Kashyap, S.R.; Kasumov, T.; Zhang, R.; Kelly, K.R.; Defronzo, R.A.; Kirwan, J.P. Plasma ceramides are elevated in obese subjects with type 2 diabetes and correlate with the severity of insulin resistance. Diabetes 2009, 58, 337–343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Filippov, V.; Song, M.A.; Zhang, K.; Vinters, H.V.; Tung, S.; Kirsch, W.M.; Yang, J.; Duerksen-Hughes, P.J. Increased ceramide in brains with Alzheimer’s and other neurodegenerative diseases. J. Alzheimers Dis. 2012, 29, 537–547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Grimm, M.O.W.; Michaelson, D.M.; Hartmann, T. Omega-3 fatty acids, lipids, and apoE lipidation in Alzheimer’s disease: A rationale for multi-nutrient dementia prevention. J. Lipid Res. 2017, 58, 2083–2101. [Google Scholar] [CrossRef] [Green Version]
  67. Kurz, J.; Parnham, M.J.; Geisslinger, G.; Schiffmann, S. Ceramides as Novel Disease Biomarkers. Trends Mol. Med. 2019, 25, 20–32. [Google Scholar] [CrossRef]
  68. Nakayama, H.; Nagafuku, M.; Suzuki, A.; Iwabuchi, K.; Inokuchi, J.I. The regulatory roles of glycosphingolipid-enriched lipid rafts in immune systems. FEBS Lett. 2018, 592, 3921–3942. [Google Scholar] [CrossRef] [Green Version]
  69. Aronov, P.A.; Luo, F.J.; Plummer, N.S.; Quan, Z.; Holmes, S.; Hostetter, T.H.; Meyer, T.W. Colonic contribution to uremic solutes. J. Am. Soc. Nephrol. 2011, 22, 1769–1776. [Google Scholar] [CrossRef] [Green Version]
  70. Ottosson, F.; Brunkwall, L.; Smith, E.; Orho-Melander, M.; Nilsson, P.M.; Fernandez, C.; Melander, O. The gut microbiota-related metabolite phenylacetylglutamine associates with increased risk of incident coronary artery disease. J. Hypertens. 2020, 38, 2427–2434. [Google Scholar] [CrossRef]
  71. Ussher, J.R.; Lopaschuk, G.D.; Arduini, A. Gut microbiota metabolism of L-carnitine and cardiovascular risk. Atherosclerosis 2013, 231, 456–461. [Google Scholar] [CrossRef] [PubMed]
  72. Makrecka-Kuka, M.; Sevostjanovs, E.; Vilks, K.; Volska, K.; Antone, U.; Kuka, J.; Makarova, E.; Pugovics, O.; Dambrova, M.; Liepinsh, E. Plasma acylcarnitine concentrations reflect the acylcarnitine profile in cardiac tissues. Sci. Rep. 2017, 7, 17528. [Google Scholar] [CrossRef] [Green Version]
  73. DiNicolantonio, J.J.; O’Keefe, J.H. Good Fats versus Bad Fats: A Comparison of Fatty Acids in the Promotion of Insulin Resistance, Inflammation, and Obesity. Mo. Med. 2017, 114, 303–307. [Google Scholar] [PubMed]
  74. Miller, M.; Stone, N.J.; Ballantyne, C.; Bittner, V.; Criqui, M.H.; Ginsberg, H.N.; Goldberg, A.C.; Howard, W.J.; Jacobson, M.S.; Kris-Etherton, P.M.; et al. Triglycerides and cardiovascular disease: A scientific statement from the American Heart Association. Circulation 2011, 123, 2292–2333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Kaneda, T. Iso- and anteiso-fatty acids in bacteria: Biosynthesis, function, and taxonomic significance. Microbiol. Rev. 1991, 55, 288–302. [Google Scholar] [CrossRef] [PubMed]
  76. Mika, A.; Stepnowski, P.; Kaska, L.; Proczko, M.; Wisniewski, P.; Sledzinski, M.; Sledzinski, T. A comprehensive study of serum odd- and branched-chain fatty acids in patients with excess weight. Obesity 2016, 24, 1669–1676. [Google Scholar] [CrossRef] [Green Version]
  77. Ran-Ressler, R.R.; Khailova, L.; Arganbright, K.M.; Adkins-Rieck, C.K.; Jouni, Z.E.; Koren, O.; Ley, R.E.; Brenna, J.T.; Dvorak, B. Branched chain fatty acids reduce the incidence of necrotizing enterocolitis and alter gastrointestinal microbial ecology in a neonatal rat model. PLoS ONE 2011, 6, e29032. [Google Scholar] [CrossRef] [Green Version]
  78. Su, X.; Magkos, F.; Zhou, D.; Eagon, J.C.; Fabbrini, E.; Okunade, A.L.; Klein, S. Adipose tissue monomethyl branched-chain fatty acids and insulin sensitivity: Effects of obesity and weight loss. Obesity 2015, 23, 329–334. [Google Scholar] [CrossRef] [Green Version]
  79. Jaceldo-Siegl, K.; Haddad, E.; Knutsen, S.; Fan, J.; Lloren, J.; Bellinger, D.; Fraser, G.E. Lower C-reactive protein and IL-6 associated with vegetarian diets are mediated by BMI. Nutr. Metab. Cardiovas. 2018, 28, 787–794. [Google Scholar] [CrossRef] [PubMed]
  80. Orlich, M.J.; Fraser, G.E. Vegetarian diets in the Adventist Health Study 2: A review of initial published findings. Am. J. Clin. Nutr. 2014, 100 (Suppl. S1), 353S–358S. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Random forest variable importance. Mean decrease accuracy from random forest analysis classifying vegans and non-vegetarians represents average decrease in accuracy in model prediction after permutation of each indicated variable.
Figure 1. Random forest variable importance. Mean decrease accuracy from random forest analysis classifying vegans and non-vegetarians represents average decrease in accuracy in model prediction after permutation of each indicated variable.
Nutrients 14 00709 g001
Figure 2. PCA plot of scores for vegan and non-vegetarian participants for principal components 1 and 3. Individual scores were obtained from covariance matrix of 930 log-transformed metabolites and plotted to analyze variation explained by diet group.
Figure 2. PCA plot of scores for vegan and non-vegetarian participants for principal components 1 and 3. Individual scores were obtained from covariance matrix of 930 log-transformed metabolites and plotted to analyze variation explained by diet group.
Nutrients 14 00709 g002
Table 1. Demographic, dietary, and lifestyle characteristics of study population 1,2.
Table 1. Demographic, dietary, and lifestyle characteristics of study population 1,2.
Non-VegetarianVeganp-Value
Participants, n4647
Demographic
Sex 1
Male23 (50)24 (51)
Female23 (50)23 (49)
Age (years)60.8 (11.1)66.5 (9.9)0.011
BMI (kg/m2) 331.3 (7.0)24.7 (3.4)<0.001
Ethnicity 1
Caucasian23 (50)23 (49)
Black23 (50) 24 (51)
Dietary
Total kca32520.8 (902.8)1776.8 (744.7)<0.001
Fiber (g/d)19.6 (6.5)40.5 (8.8)<0.001
Saturated fat (g/d)21.8 (7.9)9.0 (2.7)<0.001
Total meat (g/d)104.9 (41.0)0< 0.001
Red meat (g/d)54.9 (22.8)0< 0.001
Processed meat (g/d)6.4 (6.6)0<0.001
Poultry (g/d)31.2 (21.6)0< 0.001
Fish (g/d)20.7 (17.2)0< 0.001
Dairy (kcal/d)252.8 (151.7)17.0 (14.7) 3<0.001
Fruit (kcal/d)108.3 (64.9)320.9 (159.3)<0.001
Vegetables (kcal/d)67.6 (36.3)134.5 (115.4)<0.001
Soy (kcal/d)60.06 (42.3)156.3 (99.2)<0.001
Legumes (kcal/d)46.1 (56.4)73.2 (56.0)<0.001
Whole grains (kcal/d)172.4 (130.9)326.6 (189.8)<0.022
Lifestyle
Exercise (min/week) 463.6 (100.8)132.3 (110.6)0.003
Smoking (ever vs. never)12 (26.1)2 (4.3)0.008
Alcohol drinking (ever vs. never)24 (52.2)12 (25.5)0.015
Coffee (any vs. never/rarely)18 (39.1)0<0.001
Aspirin/NSAIDS (<1/wk vs. ≥1/wk)11 (23.9)3 (6.4)0.038
1 Values presented as n (%) or mean (SD). 2 Dietary variables adjusted for kcal/day. 3 Based on follow up phone calls in a sample, these small quantities of dairy are most likely artefactual coming from recipes of certain mixed foods where vegans can choose non-dairy options not reflected in these recipes. 4 Missing values: body mass index (BMI), n = 2; kcal, n = 1; Exercise, n = 2. NSAIDS = non-steroidal anti-inflammatory drugs.
Table 2. Top 40 (of 164) metabolites positively associated with a vegan relative to non-vegetarian dietary pattern at FDR < 0.05 1,2.
Table 2. Top 40 (of 164) metabolites positively associated with a vegan relative to non-vegetarian dietary pattern at FDR < 0.05 1,2.
Metabolite Fold ChangeFDR SubclassMajor Class
4-ethylphenyl sulfate6.585.2 × 10−4Benzoate MetabolismXenobiotics
s-methylmethionine4.612.0 × 10−4Methionine, Cysteine, SAM, and Taurine MetabolismAmino Acid
branched chain 14:0 dicarboxylic acid3.542.2 × 10−4Fatty Acid, DicarboxylateLipid
4-acetylphenyl sulfate3.202.6 × 10−4Benzoate MetabolismXenobiotics
glycohyocholate2.866.3 × 10−4Secondary Bile Acid MetabolismLipid
ethyl beta-glucopyranoside2.791.5 × 10−4Food Component/PlantXenobiotics
methyl glucopyranoside (alpha + beta)2.642.9 × 10−4Food Component/PlantXenobiotics
beta-cryptoxanthin2.521.0 × 10−3Vitamin A MetabolismCofactors and Vitamins
n-methylproline2.289.0 × 10−5Urea cycle; Arginine and Proline MetabolismAmino Acid
stachydrine2.271.6 × 10−4Food Component/PlantXenobiotics
4-allylcatechol sulfate2.216.7 × 10−5Benzoate MetabolismXenobiotics
indolepropionate2.039.8 × 10−5Tryptophan MetabolismAmino Acid
2-methylserine1.951.3 × 10−4Glycine, Serine, and Threonine MetabolismAmino Acid
1-linoleoyl-2-linolenoyl-gpc (18:2/18:3)1.934.5 × 10−4Phosphatidylcholine (PC)Lipid
s-methylcysteine sulfoxide1.931.1 × 10−4Methionine, Cysteine, SAM, and Taurine MetabolismAmino Acid
daidzein sulfate (2)1.919.8 × 10−3Food Component/PlantXenobiotics
gentisate1.907.8 × 10−5Tyrosine MetabolismAmino Acid
4-allylphenol sulfate1.892.4 × 10−3Food Component/PlantXenobiotics
dodecanedioate (c12:1-dc)1.881.8 × 10−4Fatty Acid, DicarboxylateLipid
glycochenodeoxycholate 3-sulfate1.872.3 × 10−3Primary Bile Acid MetabolismLipid
3-hydroxybutyrate (bhba)1.861.7 × 10−4Ketone BodiesLipid
octadecanedioate (c18:1-dc)1.833.5 × 10−4Fatty Acid, DicarboxylateLipid
cinnamoylglycine1.821.7 × 10−3Food Component/PlantXenobiotics
n-delta-acetylornithine1.791.6 × 10−4Urea cycle; Arginine and Proline MetabolismAmino Acid
3-hydroxydodecanedioate1.763.9 × 10−4Fatty Acid, DicarboxylateLipid
catechol sulfate1.751.2 × 10−4Benzoate MetabolismXenobiotics
octadecadienedioate (c18:2-dc)1.751.2 × 10−4Fatty Acid, DicarboxylateLipid
tryptophan betaine1.741.8 × 10−3Tryptophan MetabolismAmino Acid
n-linoleoylglycine1.706.2 × 10−5Fatty Acid Metabolism (Acyl Glycine)Lipid
4-methoxyphenol sulfate1.691.4 × 10−3Tyrosine MetabolismAmino Acid
octadecenedioylcarnitine (c18:1-dc)1.691.4 × 10−4Fatty Acid Metabolism (Acyl Carnitine, Dicarboxylate)Lipid
4-ethylcatechol sulfate1.683.5 × 10−3Benzoate MetabolismXenobiotics
12,13-dihome1.669.5 × 10−5Fatty Acid, DihydroxyLipid
carotene diol (2)1.658.0 × 10−5Vitamin A MetabolismCofactors and Vitamins
chiro-inositol1.651.2 × 10−3Inositol MetabolismLipid
2-acetamidophenol sulfate1.652.1 × 10−2Food Component/PlantXenobiotics
3-hydroxysebacate1.632.1 × 10−4Fatty Acid, MonohydroxyLipid
acetoacetate1.631.6 × 10−4Ketone BodiesLipid
Cis-4-decenoate (10:1n6)1.631.2 × 10−4Medium Chain Fatty AcidLipid
pentose acid1.627.7 × 10−5Partially Characterized MoleculesPartially Characterized Molecules
1 Fold change represents ratio of geometric means of vegans relative to non-vegetarians obtained from linear regression model with SmartSVA. 2 Adapted Storey et al. [26] permutation approach used to adjust for false discovery.
Table 3. Top 40 (of 422) metabolites inversely associated with a vegan relative to non-vegetarian dietary pattern at FDR < 3.9 × 10−5 1,2.
Table 3. Top 40 (of 422) metabolites inversely associated with a vegan relative to non-vegetarian dietary pattern at FDR < 3.9 × 10−5 1,2.
Metabolite Fold ChangeSubclassMajor Class
5-acetylamino-6-amino-3-methyluracil0.04Xanthine MetabolismXenobiotics
theobromine0.04Xanthine MetabolismXenobiotics
3-carboxy-4-methyl-5-propyl-2-furanpropanoate (cmpf)0.05Fatty Acid, DicarboxylateLipid
paraxanthine0.07Xanthine MetabolismXenobiotics
theophylline0.07Xanthine MetabolismXenobiotics
3-methylhistidine0.07Histidine MetabolismAmino Acid
1-methyl-5-imidazoleacetate0.09Histidine MetabolismAmino Acid
hydroxy-cmpf0.09Fatty Acid, DicarboxylateLipid
1,7-dimethylurate0.09Xanthine MetabolismXenobiotics
4-acetaminophen sulfate0.10Drug—Analgesics, AnestheticsXenobiotics
caffeine0.11Xanthine MetabolismXenobiotics
piperine0.12Food Component/PlantXenobiotics
7-methylxanthine0.17Xanthine MetabolismXenobiotics
2-hydroxyacetaminophen sulfate0.18Drug—Analgesics, AnestheticsXenobiotics
3-methylxanthine0.19Xanthine MetabolismXenobiotics
4-acetamidophenol0.19Drug—Analgesics, AnestheticsXenobiotics
sulfate of piperine metabolite c16h19no3 (2)0.19Food Component/PlantXenobiotics
1-methylurate0.20Xanthine MetabolismXenobiotics
5-acetylamino-6-formylamino-3-methyluracil0.20Xanthine MetabolismXenobiotics
3-bromo-5-chloro-2,6-dihydroxybenzoic acid0.22ChemicalXenobiotics
sulfate of piperine metabolite c16h19no3 (3)0.23Food Component/PlantXenobiotics
n,n,n-trimethyl-5-aminovalerate0.24Lysine MetabolismAmino Acid
glucuronide of piperine metabolite c17h21no3 (4)0.25Food Component/PlantXenobiotics
sulfate of piperine metabolite c18h21no3 (1)0.26Food Component/PlantXenobiotics
heptenedioate (c7:1-dc)0.27Fatty Acid, DicarboxylateLipid
glucuronide of piperine metabolite c17h21no3 (3)0.27Food Component/PlantXenobiotics
glucuronide of piperine metabolite c17h21no3 (5)0.28Food Component/PlantXenobiotics
1-methylxanthine0.29Xanthine MetabolismXenobiotics
(14 or 15)-methyl palmitate (a17:0 or i17:0)0.29Fatty Acid, BranchedLipid
sulfate of piperine metabolite c18h21no3 (3)0.29Food Component/PlantXenobiotics
ibuprofen0.29Drug—Analgesics, AnestheticsXenobiotics
(12 or 13)-methyl myristate (a15:0 or i15:0)0.30Fatty Acid, BranchedLipid
1,3-dimethylurate0.32Xanthine MetabolismXenobiotics
3-methyl catechol sulfate (1)0.32Benzoate MetabolismXenobiotics
1-margaroylglycerol (17:0)0.32MonoacylglycerolLipid
perfluorooctanesulfonate (pfos)0.33ChemicalXenobiotics
saccharin0.35Food Component/PlantXenobiotics
tridecenedioate (c13:1-dc)0.36Fatty Acid, DicarboxylateLipid
3,5-dichloro-2,6-dihydroxybenzoic acid0.38ChemicalXenobiotics
sphingomyelin (d18:1/25:0, d19:0/24:1, d20:1/23:0, d19:1/24:0)0.39SphingomyelinsLipid
1 Fold change represents ratio of geometric means of vegans relative to non-vegetarians obtained from linear regression model with SmartSVA. 2 Adapted Storey et al. [26] permutation approach used to adjust for false discovery.
Table 4. Metabolite subclasses associated with diet group (vegan vs. non-vegetarian) at FDR < 0.05 1.
Table 4. Metabolite subclasses associated with diet group (vegan vs. non-vegetarian) at FDR < 0.05 1.
SubclassFold Change (95% CI)FDR 2Significant
Metabolites (n)
#↓ 3#↑ 4n Total Metabolites 5
Ketone Bodies1.75 (2.16, 1.43)1.0 × 10−32022
Vitamin A Metabolism1.40 (1.57, 1.25)2.0 × 10−36246
Inositol Metabolism1.35 (1.57, 1.16)6.0 × 10−32022
Fatty Acid Metabolism (Acyl Glycine)1.21 (1.35, 1.08)2.6 × 10−25147
Lactosylceramides (LCER)1.17 (1.28, 1.07)2.8 × 10−22023
Benzoate Metabolism1.16 (1.27, 1.06)4.8 × 10−2167924
Amino Sugar Metabolism0.94 (1.00, 0.88)4.7 × 10−22205
Urea cycle; Arginine and Proline Metabolism0.94 (0.99, 0.90)2.3 × 10−21410421
Glutamate Metabolism0.94 (0.99, 0.89)3.4 × 10−285312
Tyrosine Metabolism0.92 (0.98, 0.86)1.1 × 10−2119222
Pyrimidine Metabolism, Uracil containing0.92 (0.97, 0.88)1.1 × 10−355012
Purine Metabolism, Adenine containing0.92 (0.97, 0.88)1.4 × 10−33306
Fatty Acid, Dicarboxylate0.91 (0.98, 0.85)1.5 × 10−222101234
Long-Chain Polyunsaturated Fatty Acid (n3 and n6)0.90 (0.97, 0.83)1.0 × 10−2129317
Phospholipid Metabolism0.90 (0.96, 0.85)1.2 × 10−34407
Sphingolipid Synthesis0.89 (0.98, 0.80)2.1 × 10−22203
Alanine and Aspartate Metabolism0.89 (0.94, 0.85)1.1 × 10−46519
Partially Characterized Molecules0.88 (0.98, 0.79)2.8 × 10−286213
Creatine Metabolism0.86 (0.92, 0.79)1.7 × 10−43213
Diacylglycerol0.84 (0.97, 0.72)2.5 × 10−266014
Pyrimidine Metabolism, Orotate containing0.84 (0.93, 0.76)9.9 × 10−43304
Purine Metabolism, (Hypo)Xanthine/Inosine containing0.84 (0.90, 0.77)<3.9 × 10−55507
Fructose, Mannose, and Galactose Metabolism0.83 (0.98, 0.70)3.4 × 10−22204
Secondary Bile Acid Metabolism0.83 (0.95, 0.72)8.1 × 10−3119221
Purine Metabolism, Guanine containing0.83 (0.90, 0.76)<3.9 × 10−53303
Sphingomyelins0.83 (0.89, 0.77)<3.9 × 10−52218429
Food Component/Plant0.83 (0.88, 0.78)<3.9 × 10−536221451
Fatty Acid Metabolism (also BCAA Metabolism)0.82 (0.93, 0.72)4.8 × 10−34405
Tryptophan Metabolism0.82 (0.88, 0.76)<3.9 × 10−51714320
Pantothenate and CoA Metabolism0.81 (0.97, 0.68)2.6 × 10−21102
Guanidino and Acetamido Metabolism0.80 (0.94, 0.69)1.0 × 10−22202
Lysine Metabolism0.80 (0.87, 0.74)<3.9 × 10−51111018
Plasmalogen0.80 (0.85, 0.75)<3.9 × 10−576111
Glycerolipid Metabolism0.79 (0.88, 0.71)4.3 × 10−53303
Leucine, Isoleucine, and Valine Metabolism0.79 (0.84, 0.75)<3.9 × 10−52423132
Dihydrosphingomyelins0.78 (0.89, 0.69)3.0 × 10−45415
Chemical0.78 (0.83, 0.73)<3.9 × 10−51310320
Monoacylglycerol0.77 (0.84, 0.69)<3.9 × 10−51111017
Lysoplasmalogen0.75 (0.82, 0.68)<3.9 × 10−54404
Phenylalanine Metabolism0.75 (0.80, 0.70)<3.9 × 10−56607
Long-Chain Monounsaturated Fatty Acid0.74 (0.80, 0.68)<3.9 × 10−54407
Long-Chain Saturated Fatty Acid0.73 (0.78, 0.68)<3.9 × 10−57708
Fatty Acid Metabolism (Acyl Carnitine, Long-Chain Saturated)0.72 (0.80, 0.65)<3.9 × 10−58808
Ceramides0.64 (0.72, 0.57)<3.9 × 10−51010011
Dihydroceramides0.62 (0.74, 0.52)<3.9 × 10−52202
Acetylated Peptides0.59 (0.74, 0.48)<3.9 × 10−54404
Histidine Metabolism0.59 (0.65, 0.55)<3.9 × 10−51110115
Fatty Acid, Branched0.33 (0.38, 0.29)<3.9 × 10−53303
Drug—Analgesics, Anesthetics0.22 (0.29, 0.16)<3.9 × 10−55505
Xanthine Metabolism0.14 (0.20, 0.10)<3.9 × 10−51313013
1 Linear regression analysis with SmartSVA based on composite t-statistics generated by dividing standard deviation of averaged numerators representing log transformed metabolites. 2 Adapted Storey et al. [26] permutation approach used to adjust for false discovery. 3 Number of metabolites within subclass that were inversely associated with a vegan diet among those significantly differential in the linear regression analysis. 4 Number of metabolites within subclass that were positively associated with a vegan diet among those significantly differential in the linear regression analysis. 5 Total number of metabolites measured.
Table 5. Adjusted linear regression predicting dietary and lifestyle characteristics from top principal components derived from principal component analysis 1.
Table 5. Adjusted linear regression predicting dietary and lifestyle characteristics from top principal components derived from principal component analysis 1.
βSET Valuep-ValueCorrelation
Coefficient
p-Value
(Correlation)
Vegan −0.57.9 × 10−7
PC1−7.331.96−3.73.3 × 10−4
PC38.521.884.61.7 × 10−5
PC4−5.231.44−3.65.0 × 10−4
BMI 0.612.1 × 10−10
PC10.990.156.54.7 × 10−9
PC2−0.280.13−2.14.0 × 10−2
PC3−0.350.15−2.42.0 × 10−2
Red meat 0.523.5 × 10−7
PC11.890.493.93.7 × 10−3
PC3−2.140.47−4.61.7 × 10−5
PC41.350.363.73.5 × 10−4
Total meat 0.514.0 × 10−7
PC11.610.423.82.7 × 10−4
PC3−1.860.4−4.61.4 × 10−5
PC41.170.313.83.2 × 10−4
Processed meat
PC11.510.374.18.9 × 10−50.442.04 × 10−5
PC3−0.790.39−2.14.0 × 10−2
PC40.750.282.79.6 × 10−3
Poultry
PC11.990.573.497.6 × 10−40.451.10 x 10−5
PC3−2.250.55−4.19.2 × 10−5
PC41.450.423.468.5 × 10−4
Fish 0.522.6 × 10−7
PC11.910.613.12.0 × 10−3
PC2−1.150.56−2.14.0 × 10−2
PC3−2.270.58−3.92.0 × 10−4
PC41.640.443.83.0 × 10−4
Fiber −0.417.2 × 10−5
PC1−6.672.18−3.12.9 × 10−3
PC39.132.044.52.2 × 10−5
PC4−3.321.65−2.05.0 × 10−2
Soy −0.220.04
PC4−1.70.72−2.42.0 × 10−2
PC71.210.532.32.5 × 10−2
Vegetables −0.230.03
PC1−2.71.18−2.32.5 × 10−2
PC8−2.40.65−3.74.3 × 10−4
Fruit
PC1−4.690.94−5.03.4 × 10−6−0.473.4 × 10−6
Dairy 0.639.6 × 10−11
PC12.270.564.19.7 × 10−5
PC2−1.030.53−2.05.4 × 10−2
PC3−2.060.56−3.73.8 × 10−4
PC41.260.423.03.9 × 10−3
Saturated fat
PC17.441.634.61.7 × 10−50.553.1 × 10−8
PC3−6.71.64−4.19.3 × 10−5
PC43.241.292.51.4 × 10−2
Whole grains
PC33.140.893.56.6 × 10−40.366.6 × 10−4
1 Individual principal components (PCs) regressed on indicated dietary variables in linear regression model adjusted for age, sex, race, and BMI.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Miles, F.L.; Orlich, M.J.; Mashchak, A.; Chandler, P.D.; Lampe, J.W.; Duerksen-Hughes, P.; Fraser, G.E. The Biology of Veganism: Plasma Metabolomics Analysis Reveals Distinct Profiles of Vegans and Non-Vegetarians in the Adventist Health Study-2 Cohort. Nutrients 2022, 14, 709. https://doi.org/10.3390/nu14030709

AMA Style

Miles FL, Orlich MJ, Mashchak A, Chandler PD, Lampe JW, Duerksen-Hughes P, Fraser GE. The Biology of Veganism: Plasma Metabolomics Analysis Reveals Distinct Profiles of Vegans and Non-Vegetarians in the Adventist Health Study-2 Cohort. Nutrients. 2022; 14(3):709. https://doi.org/10.3390/nu14030709

Chicago/Turabian Style

Miles, Fayth L., Michael J. Orlich, Andrew Mashchak, Paulette D. Chandler, Johanna W. Lampe, Penelope Duerksen-Hughes, and Gary E. Fraser. 2022. "The Biology of Veganism: Plasma Metabolomics Analysis Reveals Distinct Profiles of Vegans and Non-Vegetarians in the Adventist Health Study-2 Cohort" Nutrients 14, no. 3: 709. https://doi.org/10.3390/nu14030709

APA Style

Miles, F. L., Orlich, M. J., Mashchak, A., Chandler, P. D., Lampe, J. W., Duerksen-Hughes, P., & Fraser, G. E. (2022). The Biology of Veganism: Plasma Metabolomics Analysis Reveals Distinct Profiles of Vegans and Non-Vegetarians in the Adventist Health Study-2 Cohort. Nutrients, 14(3), 709. https://doi.org/10.3390/nu14030709

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop