key: cord-0070019-1v0eg49x
authors: Zhao, Xinyue; Zhou, Ruilin; Li, Hanyu; Fan, Yue; Sun, Yueshen; Hu, Xiaomin; Zhang, Shuyang
title: The Effects of Moderate Alcohol Consumption on Circulating Metabolites and Gut Microbiota in Patients With Coronary Artery Disease
date: 2021-11-02
journal: Front Cardiovasc Med
DOI: 10.3389/fcvm.2021.767692
sha: 8d434800a98110748f0268ca959ff2b86aae9881
doc_id: 70019
cord_uid: 1v0eg49x

Background: Epidemiological studies confirmed that moderate alcohol consumption was associated with a reduced risk of adverse cardiovascular events. It is increasingly recognized that the composition of gut microbiota and metabolites is involved in modulating the cardiovascular health of the host. However, the association of moderate alcohol consumption with serum metabolites and gut microbiome and its impact on coronary artery disease (CAD) is not fully investigated. Method: Serum untargeted metabolomics analysis and fecal 16S rRNA sequencing were performed on 72 male patients with CAD having various alcohol consumption (36 non-drinkers, 18 moderate drinkers, and 18 heavy drinkers) and 17 matched healthy controls. MetaboAnalyst and PICRUSt2 were utilized to analyze the possible involved metabolic pathways. Multi-omics analysis was achieved by Spearman correlation to reveal the interactions of alcohol consumption with gut microbiome and serum metabolites in patients with CAD. Results: We noted distinct differences between patients with CAD, with varying levels of alcohol consumption and healthy controls in aspects of serum metabolome and the gut microbiome. Moderate alcohol consumption significantly changed the lipidomic profiles, including reductions of sphingolipids and glycerophospholipids in moderate drinkers with CAD when compared with non and heavy drinkers with CAD. Moreover, we also found the reduction of microbial-derived metabolites in moderate drinkers with CAD, such as 2-phenylacetamide and mevalonic acid. To be noted, the gut microbiota of moderate drinkers with CAD tended to resemble that of healthy controls. Compared with non-drinkers, the relative abundance of genus Paraprevotella, Lysinibacillus was significantly elevated in moderate drinkers with CAD, while the genus Bifidobacterium, Megasphaera, and Streptococcus were significantly reduced in moderate drinkers with CAD. Multi-omics analysis revealed that specific metabolites and microbes associated with moderate alcohol consumption were correlated with the severity of CAD. Conclusions: Our study revealed that the impact of moderate alcohol consumption on serum metabolites and gut microbiota in patients with CAD seemed to be separated from that of heavy and non-alcohol consumption. Moderate drinking tended to have more positive effects on metabolic profiles and commensal flora, which may explain its beneficial effects on cardiovascular health. Overall, our study provides a novel insight into the effects of moderate alcohol consumption in patients with CAD.

Coronary artery disease is one of the most common causes of death in the general population (1) , and alcohol abuse is a widely recognized risk factor in adverse cardiovascular events. Alcohol affects the cardiovascular system of human beings in several ways. It was reported that alcohol abuse can be associated with various adverse cardiovascular events, such as stroke, acute heart failure, and sudden cardiac death (2) . However, clinical studies demonstrated that moderate drinking presented an inverse correlation with the incidence of coronary artery disease, perhaps due to its effects on lipoprotein levels (3) . Epidemiological studies illustrated that moderate alcohol intake (one to two cups a day, or 100 grams per week) was associated with a decreased incidence of CAD (4, 5) . What is more, the association between the severity of cardiovascular diseases and different alcohol consumption presented a U-shaped curve. In other words, light to moderate alcohol intake had a reduced risk of cardiovascular events, while the relationship reversed in chronic alcohol overconsumption (4) (5) (6) (7) (8) . In addition, recent studies have reported possible mechanisms for the protective effects of moderate alcohol consumption on patients with CAD. Some research reported that moderate alcohol consumption may elevate the level of high-density lipoprotein cholesterol (HDL-C), apolipoprotein A, and adiponectin levels (9, 10) , which were considered as protective factors in cardiovascular disease. It was also reported that moderate alcohol consumption may affect the modulation of insulin resistance and chronic inflammation. By suppressing inflammation and promoting vasodilation, adverse cardiovascular events can be reduced (11) . Furthermore, by moderate alcohol consumption, the hypercoagulable state of patients with CAD may be reduced by the anticoagulation effect of alcohol (12) . Overall, it was believed that moderate alcohol consumption was associated with a reduced risk of cardiovascular events by various possible mechanisms.

Meanwhile, it is increasingly recognized that the composition of gut microbiota and serum metabolites is involved in modulating the cardiovascular health of the host. In recent years, more and more studies have confirmed that alterations in the gut microbiome were associated with the severity of CAD (13) (14) (15) . Gut microbiota participates in modulating the cardiovascular health of the host through its metabolites, such as short-chain fatty acids (SCFA), trimethylamine Noxide (TMAO), and phenylacetylglutamine (PAGln) (16) (17) (18) . Alterations in serum metabolome, such as C34:2 hydroxyphosphatidylcholine, N-acetylneuraminic acid, sphingomyelin (SM), were reported to be significantly associated with CAD and were considered as potential biomarkers (19) (20) (21) (22) . Moreover, the metabolic features of alcohol drinkers were investigated and showed profound modifications of lipidomic profiles and amino acids in humans (23) (24) (25) . Alcohol use was associated with lower levels of phosphatidylcholine acyl-alkyls, hydroxysphingomyelin, glutamine, and citrate, and higher concentration of fatty acids, phosphatidylcholine diacyls, tyrosine, and alanine (23) . Many metabolic profiles showed U-shape associations with alcohol consumption, such as total triglycerides and phenylalanine (24) . Thus, there may exist the possibility that moderate alcohol consumption may have a protective effect on cardiovascular health by modulating the gut microbiome and serum metabolome.

To address the question above, we recruited a total of 89 patients with CAD enrolled at Peking Union Medical College Hospital (PUMCH) and divided them into four groups (detailed in the "Materials and Methods" section): healthy controls (HC) (n = 17), non-drinkers with CAD (ND-CAD) (n = 36), moderate drinkers with CAD (MD-CAD) (n = 18), and heavy drinkers with CAD (HD-CAD) (n = 18). Detailed clinical data of all 89 subjects were collected. Moreover, 16S rRNA sequencing and liquid chromatography-mass spectrometry (LC-MS) were applied. The aim of this study was to reveal the association between moderate drinking and gut microbiome, serum metabolome in patients with CAD, as well as key alterations related to clinical benefits.

The participants in the study were recruited consecutively at the department of cardiology in Peking Union Medical College Hospital from 2016 to 2018. The inclusion criteria were as follows: (1) male patients (to exclude the influence of hormone); (2) patients who exhibited ≥ 50% stenosis in at least one main coronary artery in coronary angiography; exclusion criteria included the complication of gastrointestinal diseases, malignant tumor, autoimmune disorders, infectious diseases, renal dysfunction (creatine > 3 mg/dl), a history of gastrointestinal surgery within a year, and had antibiotics over 3 days in the last 3 months. In addition, the study focused on current drinkers and current non-drinkers, excluding the abstinent drinkers, who may interfere with the results. The subjects enrolled in the study had detailed questionnaires about their demographic features, living habits, and clinical information. Peripheral venous blood and stool samples were collected the next morning after admission. The preparation and the storage of blood and stool samples were described in our previous study (26) . The freshly collected samples from each participant were immediately transported to our laboratory and stored in a −80 • C refrigerator.

The level of alcohol consumption was assessed based on the information collected in the questionnaire, including drinking history, frequency, amount, and the most often consumed type of alcoholic beverage. The alcohol intake (g/d) was obtained by calculating drinking frequency, amount, ethanol density (0.8 g/L), and the alcohol content of each beverage (%v/v): 50.0% for liquor, 12.9% for wine, and 5.3% for beer (27) . Based on the cut-off value for different drinking categories in other research (28) (29) (30) , 72 patients with CAD were divided into three groups: (1) non-drinkers with CAD (ND-CAD): patients with CAD who denied a drinking history; (2) moderate drinkers with CAD (MD-CAD): patients with CAD who had an ongoing regular light-tomoderate drinking habit (the average alcohol intake was 0-40 g pure alcohol per day); (3) heavy drinkers with CAD (HD-CAD): patients with CAD who had alcoholism habits or chronic alcohol over-intake (more than 40 g pure alcohol per day).

In addition, 17 healthy volunteers who met the following criteria were enrolled as the healthy control (HC): (1) male, (2) did not suffer from CAD, (3) and did not meet any of the exclusion criteria above.

The statistical analysis of the characteristics in the study population was described in detail in our previous work (31) . The study complied with the principles of the Declaration of Helsinki. All the participants in the study provided written informed consent.

Serum metabolome analysis was conducted on a Waters ACQUITY ultra-high-performance liquid chromatography system (Milford, MA), coupled with a Waters Q-TOF Micromass system (Manchester, UK) in positive and negative ionization modes. Polar ionic and lipid modes were performed based on the properties of metabolites. Sample preparation and LC-MS experiment procedure were detailed and described previously (26) . The peak-ion intensity matrix was then filtered by removing peaks with a zero value in more than 80% of samples. The threshold of the coefficient variation value of quality control samples was set at 30%. Wilcoxon rank-sum test was adopted to identify the metabolites with significant differences between groups. Partial least squares discriminant analysis (PLS-DA) was conducted by SIMCA software (v14.1, Umetrics, Sweden). Variable importance in the projection (VIP)-value > 1 and P < 0.05 was adopted for select important peaks. Online databases Human Metabolome Database (https://hmdb.ca), LipidMaps (https://www.lipidmaps.org), and PubChem (https:// pubchem.ncbi.nlm.nih.gov/) were used for classifying peaks according to the molecular mass data (m/z). MetaboAnalyst (https://www.metaboanalyst.ca) was used to conduct pathway enrichment analysis. Pathways were regarded as potential targets with a threshold of impact-value > 0.10 when utilizing MetaboAnalyst (32) .

Microbial DNAs were extracted from the stool samples using the bead-beating method (33) . Then, the amplification of the V3-V4 region of 16S rRNA genes was performed using PCR (34) . The sequencing library was established as described previously (35) , and purified products were sequenced with the Illumina Miseq system (Illumina Inc., USA). The downstream amplicon analysis was conducted by EasyAmplicon v1.0 (36) . Dereplication was performed by the -derep_fullength command of VSEARCH (v2.15) (37). Operational taxonomic units (OTU) were clustered via the -cluster_otus command of USEARCH (v10.0) at the cutoff of 97% (38) . The feature table was created with vsearch-usearch_global. Taxonomic annotation was generated according to Greengenes database using usearch-otutab (39) .

The sequences of all samples were downsized to the sample with the least sequences to calculate the diversity indices. Alpha diversity was evaluated with Shannon's index and Chao1 index. Beta diversity was performed by principal coordinate analysis (PCoA) and constrained PCoA (CPCoA) with Bray-Curtis distances. The composition of each group was represented as a boxplot plot at the phylum level and as a Chord diagram at the genus level with R package ggplot2. As for the comparison of differences, edgeR was applied to identify the differences between groups and the Benjamini-Hochberg method to control the FDR (40) . A threshold of P < 0.05 with FDR < 0.2 was considered statistically significant. Bugbase and Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) were utilized to perform functional predication of the gut microbiota (41, 42) . Furthermore, PICRUSt2 was utilized to predict the metagenomic pathways based on the MetaCyc database (43) . Pathways that were significantly different between the moderate drinkers and non-drinkers were identified by Welch's t-test, and Storey FDR was utilized for multiple pathways. STAMP software (v2.1.3) was applied for statistical analysis and visualization of the identified pathways.

For metabolomics analysis, Wilcoxon rank-sum test was conducted, and the threshold of VIP > 1 and P < 0.05 was used to identify differential metabolites. To be detailed, the metabolites that were significantly higher or lower in MD-CAD than those in both ND-CAD and HD-CAD groups were defined as key metabolites associated with moderate alcohol consumption. In terms of the microbiota, using edgeR with a threshold of P < 0.05 and FDR < 0.2, there were a total of 21 differential OTUs between MD-CAD and ND-CAD. Among the 21 OTUs, key OTUs were defined as those whose average relative abundances represented a U-shaped or convex-shaped relationship with different levels of alcohol consumption (i.e., the key OTUs present either the highest or the lowest average relative abundance in MD-CAD when comparing with ND-CAD and HD-CAD).

Spearman correlation analysis was conducted between key bacterial taxa, metagenomic pathway, serum metabolites, and clinical parameters with SPSS (v. 24.0), and visualized as a heatmap using the R package pheatmap. A Sankey plot was drawn by R package networkD3 to illustrate the interrelationship among the differential OTUs, metabolic features, and clinical parameters.

Continuous normally distributed data among three groups were analyzed by one-way ANOVA. Kruskal-Wallis H-test was employed to analyze continuous data with non-normal distribution among three groups, and Mann-Whitney U test was applied for that between two groups. Categorical variables were analyzed by χ2 test or Fisher's exact test. Data were analyzed using SPSS (v.24.0). Figures were made by utilizing R 4.0.4.

The clinical characteristics of the study cohort were summarized in Table 1 . In general, among different alcohol consumption groups, there was no significant difference in age, body mass index (BMI), and systolic blood pressure. In terms of alcohol consumption, there were significant differences between the three disease groups. The amount of alcohol intake was significantly higher in HD-CAD when compared with MD-CAD (p < 0.001), and liquor was the most common alcohol intake type. Drinking history, drinking years, and drinking frequency showed no significant difference. The three disease groups showed no significant difference in hypertension (HTN), type 2 diabetes mellitus (T2DM), hyperlipidemia (HLP), and fatty liver disease (FLD). In terms of CAD severity, although there was no significant difference in Gensini score. Gensini score is an assessment of the CAD severity using coronary angiography and positively associated with disease severity (44) . The average Gensini score of MD-CAD [27.50 (16.00, 42.75) ] was the lowest when compared with ND-CAD [39.50 (19.25, 66.75) ] and HD-CAD [34.75 (24.75, 56 .00)]. Also, the average level of CK-MB [CK-MB is a cardiac biomarker that represents cardiovascular injury (45) ] in MD-CAD was the lowest when compared with non-drinkers and heavy drinkers, so as the level of hs-CRP [The level of hs-CRP represents the inflammation status in the human body (46) , and chronic inflammation plays a central role in the development of cardiovascular disease (47) ]. The other laboratory tests, such as total cholesterol (TC), triacylglycerol (TG), low-density lipoprotein cholesterol (LDL-C), HDL-C, alanine aminotransferase (ALT), and aspartate aminotransferase (AST), showed no significant difference among the three disease groups. Thus, moderate alcohol consumption was associated with the improvement of CAD severity, that is, MD-CAD presented relatively better status than other patients with CAD.

Since it was reported that alcohol may affect serum metabolic profiles, serum metabolome in ND-CAD, MD-CAD, and HD-CAD was explored through the untargeted LC-MS method. After QC and removal of peaks with low abundance, a total of 10,416 metabolites in polar ionic mode and 3,131 in lipid mode were obtained. PLS-DA modeling was utilized to reveal the alteration of serum metabolites. The score scatter plots of four modes were shown in Figures 1A-D , and the parameters of each mode were provided in Supplementary Table 1. It was obvious that the metabolic feature of three diseased groups (ND-CAD, MD-CAD, and HD-CAD) showed significant differences in score scatter plot, especially in polar ionic positive mode ( Figure 1B ) and lipid positive mode ( Figure 1D ).

A total of 203 serum metabolites were identified as key serum metabolites associated with moderate alcohol consumption with a threshold of VIP > 1 and Wilcoxon rank-sum P < 0.05 (Supplementary Table 2 ). All the 203 key metabolites in MD-CAD were significantly different from those in ND-CAD or HD-CAD, either the most abundant or the least in the three disease groups (Supplementary Table 2 ). In other words, the relative abundance of these 203 metabolites showed a "U" shaped or convex curve with alcohol consumption. Since MD-CAD tended to have better cardiovascular health according to the previous studies, we speculated that the most abundant metabolites (i.e., convex-shaped metabolites) in MD-CAD may have protective effects in cardiovascular health, while the most deprived metabolites (i.e., "U" shaped metabolites) may be harmful. Among these 203 metabolites, we observed unignorable quantities of metabolites were lipids and lipidslike molecules, such as sphingolipids, glycerophospholipids, prenol lipids, and fatty acyls (Figure 2A) , which were consistent with other existing reports (23) (24) (25) . It was reported that a high plasma level of SM was an independent risk factor in cardiovascular diseases and predicted a poor prognosis of CAD (48, 49) . Consistently, we observed sphingolipids, such as SM (d18:0/22:1) (LPN77) and SM (d18:1/20:0) (LPN288), were depleted in MD-CAD when compared with the other two diseased groups. Moreover, we observed all annotated glycerophospholipid showed significant depletion in MD-CAD when compared with other two CAD groups (ND-CAD and HD-CAD) (Figure 2A) (51), the depletion of glycerophospholipid in MD-CAD may be associated with cardiovascular health. What is more, PC can be degraded by the gut microbiome and thus converted to TMAO, a well-known risk factor in the major adverse cardiovascular events (17, 52, 53 ). Thus, a higher level of PC was often associated with higher cardiovascular risk (17, 52, 53) . Consistently, PC was observed to be depleted in MD-CAD when compared with the other two diseased groups.

To further elucidate the cause of serum metabolome change, metabolic pathway enrichment was performed by MetaboAnalyst. The results showed that these differential metabolites mainly clustered in 11 pathways (Figure 2B,  Supplementary Table 3) . Among the 11 pathways, terpenoid backbone synthesis (P = 0.28955, impact-value = 0.0.11429), and lipoic acid metabolism (P = 0.14053, impact-value = 0.16667) were regarded as potential target pathways. Terpenoid had significant effects on the prevention and treatment of cardiovascular disease (54), while lipoic acid was reported to be beneficial in cardiovascular disease by reducing atherosclerosis (55, 56) .

Since some of the differential metabolites may be microbiome derived, such as PC, we further elucidate the taxonomic features of the gut microbiome by 16S rRNA sequencing with sufficient depths (Supplementary Figure 1A) . A total of 626 OTUs were obtained. In terms of alpha-diversity, no significant difference was observed in HC, ND-CAD, MD-CAD, and HD-CAD (Supplementary Figures 1B,C) . As for beta-diversity, although no significant difference was observed (Supplementary Figure 1D) among three diseased groups in PCoA analysis, the taxonomic features of MD-CAD were almost overlapping with HC while the other two disease groups separated clearly from HC in CPCoA analysis (Figure 3A) . Firmicutes and Bacteroides were the two most dominant phyla in all groups (Figure 3B, Supplementary Table 4) . We noted that MD-CAD and HC presented similar proportions of Actinobacteria and Proteobacteria, while ND-CAD and HD-CAD presented a much higher proportion of the two phyla ( Figure 3B ). It could be inferred that the taxonomic features of the gut microbiome in MD-CAD were more similar with healthy people when compared with the other two disease groups. What is more, Bacteroides, Prevotella, and Faecalibacterium were the top three abundant genera ( Figure 3C ) in all four groups. Venn plot was utilized to further illustrate the taxonomic features with a relative abundance >0.1% in HC, ND-CAD, MD-CAD, and HD-CAD (Figure 3D) . A total of 82 OTUs coexisted in four groups, and 94 were specifically shared by three CAD groups. Twelve OTUs were unique in MD-CAD, while those in ND-CAD and HD-CAD were 15 and 16. We further predicted the phenotypes of the gut microbiome based on their 16S rRNA sequence utilizing BugBase. Notably, HC and MD-CAD had a similar abundance of "biofilm-forming" bacteria, while ND-CAD and HD-CAD had a significantly higher abundance of "biofilmforming" bacteria ( Figure 3E) . Thus, it was clear that MD-CAD presented more similar taxonomic and functional features with HC than ND-CAD and HD-CAD. This taxonomic similarity may (F) An overview of the relative abundance of the differential OTUs associated with moderate alcohol consumption. The differential OTUs were filtered with P < 0.05 as well as FDR < 0.2. ns, not significant. *P < 0.05 and FDR < 0.2, **P < 0.01 and FDR < 0.2, ***P < 0.001 and FDR < 0.2, analyzed by edgeR.

partially explain the healthier cardiovascular status discovered in MD-CAD.

To further elucidate the alterations of the gut microbiome in MD-CAD, edgeR was utilized with a threshold of P < 0.05 and FDR < 0.02. Only five OTUs showed significant differences between HC and MD-CAD ( Figure 4A) . However, there were 21 significant differentials OTUs between MD-CAD and ND-CAD (Figure 4B) , while there were 27 significant differentials OTUs between MD-CAD and HD-CAD ( Figure 4C) . It was obvious that, regarding MD-CAD, more differential OTUs were found when comparing with ND-CAD and HD-CAD than comparing with HC. This further confirmed the similarity of gut microbiome taxonomic features between HC and MD-CAD. We suspected that it was moderate alcohol consumption that tended to alter the gut microbiome of patients with CAD to a more similar feature with HC. We further looked into the differential microbiome between ND-CAD and MD-CAD. At the phylum level, Actinobacteria, Firmicutes, and Bacteroidetes were found significantly different (Figure 4D) . At the genus level, Paraprevotella and Lysinibacillus were significantly elevated in MD-CAD while Bifidobacterium, Megasphaera, and Streptococcus were significantly depleted in MD-CAD when compared with ND-CAD (Figure 4E) . At the species level, we paid attention to 15 key OTUs associated with moderate alcohol consumption (key OTUs were defined in the "Materials and Methods" section) that showed significant differences between MD-CAD and ND-CAD (Figure 4F, Supplementary Table 5 ). Bacteroides ovatus (OTU 135, depleted in MD-CAD) was found to drive the production of IgA (57) , an autoimmune disease inflammatory mediator (58) , and was thus regarded as a driver of cardiovascular diseases. Prevotella stercorea (OTU 155, enriched in MD-CAD) was reported to be beneficial in overweight patients (59) . Bacteroides coprocola (OTU 14, depleted in MD-CAD) was reported to be enriched in patients with hypertension (60). Our result indicated that potentially beneficial bacteria were possibly promoted, and potentially pathogenic bacteria were possibly suppressed by moderate alcohol consumption.

In addition, we predicted possible functional pathways by utilizing PICRUSt2. We found 12 significantly affected pathways in MD-CAD when compared with ND-CAD (Supplementary Figure 2, Supplementary Table 6 ). Moreover, we found that the microbes depleted in MD-CAD (OTU 19, OTU 13, OTU 12, OTU 51, and OTU 135) were all positively correlated with pathways that were depleted in MD-CAD, especially the methionine-related pathways [superpathway of L-methionine biosynthesis (transsulfuration), L-methionine biosynthesis I, and superpathway of S-adenosyl-L-methionine biosynthesis]. In addition, we found most of these pathways were closely associated with metabolites degradation or biosynthesis, such as glucose and glucose-1-phosphate degradation, lactose, and galactose degradation I. Thus, we inferred that there may be exhibited functional changes of the gut microbiome that led to the alterations of serum metabolome in MD-CAD.

We analyzed the association between key OTUs and serum metabolites by utilizing the Spearman correlation ( Figure 5A) .

Several key OTUs were significantly associated with serum metabolites. We inferred that it may be the alterations in the gut microbiome that may affect the serum metabolites features. Furthermore, we correlated the key metabolites with clinical indexes and found several significant Spearman correlations ( Figure 5B) . We represented these interrelationships in a Sankey map, shown in Figure 6 . We noted that Bifidobacterium longum (OTU 12, depleted in MD-CAD), Lachnospiraceae gen. (OTU 127, depleted in MD-CAD), and Oscillibacter valericigenes (OTU 57, enriched in MD-CAD) were the three main species that impacted the serum metabolome. In terms of clinical indexes, Gensini score was positively correlated with 9R-hydroxy-octadecadienoic acid (PLP716, elevated in MD-CAD) and syringic acid (PLP7081, elevated in MD-CAD), and negatively correlated with β-citronellol (PLP3157, depleted in MD-CAD) and trans,cis-3,6-nonadien-1ol (PLP4864, depleted in MD-CAD). Moreover, cTnI negatively correlated with 12S-hydroperoxy-eicosapentaenoic acid (12-HPEPE, PLP2108, elevated in MD-CAD). Correlation analysis showed that the relative better cardiovascular health of MD-CAD may be associated with serum metabolome and gut microbiome alterations.

Alcohol presented complex effects on the cardiovascular system that vary with dose (61) . Observational and prospective studies consistently showed a lower risk of cardiovascular and all-cause mortality in people with low levels of alcohol consumption when compared with people with alcohol abuse and nondrinkers (62) . In our study, we found patients with CAD with moderate alcohol consumption presented a relatively healthier cardiovascular status, that is, a lower level of Gensini score, CK-MB, and hs-CRP. Thus, we looked further into the serum metabolome profiles and gut microbiome taxonomic features to reveal the possible associations between moderate alcohol consumption and healthier cardiovascular status.

Serum metabolome presented significant alterations in MD-CAD when compared with ND-CAD and HD-CAD. The representative metabolites were sphingolipids and glycerophospholipids, which were both significantly depleted in MD-CAD. Sphingolipids accumulate in the formation of atherosclerotic plaque both in humans and primates (63, 64) . Emerging studies uncovered that a high plasma level of sphingolipids was an independent risk factor in CAD (48, 49) . Moreover, studies on animal models confirmed that sphingolipids contributed to the development of obesity and insulin resistance (65) . The medication targeting sphingolipid metabolism may improve the condition and prognosis of metabolic disorders (66) . Consistently, in our study, the FIGURE 5 | Spearman correlations between differential serum metabolites and OTUs associated with moderate alcohol consumption or major clinical parameters for CAD. (A) Spearman correlations between differential serum metabolites and OTUs associated with moderate alcohol consumption. (B) Spearman correlations between differential serum metabolites and major clinical parameters associated with moderate alcohol consumption. The IDs of metabolites or OTUs are highlighted in red (enriched in MD-CAD) and blue (depleted in MD-CAD). *P < 0.05, **P < 0.01. ApoB, apolipoprotein B; TC, total cholesterol; LDL-C, low-density lipoprotein cholesterol; family history: having a family history of CAD; cTnI, cardiac tropinin I; HDL-C, high-density lipoprotein cholesterol; ApoA1, apolipoprotein A1; FFA, free fatty acid; Lpa, lipoprotein (a); HsCRP, high-sensitivity C-reactive protein; Cr, creatine; TG, triglyceride; Glu, fasting blood glucose; UA, uric acid; BUN, blood urea nitrogen; no. vessel, number of stenosed vessels; HTN, hypertension; T2DM, type 2 diabetes mellitus; HLP, hyperlipidemia; CK-MB, creatine kinase-MB; CK, creatine kinase; NYHA, grade of NYHA classification.

depletion of sphingolipids associated with moderate alcohol consumption tended to relieve cardiovascular stress. In addition, C19 sphingosine-1-phosphate (S1P, LPP138, depleted in MD-CAD) was reported to contribute to the pathogenesis of multiple cardiovascular disorders (67) by activating the proliferation and migration of vascular smooth muscle cells (68) . Hence, we suspected that moderate alcohol consumption may suppress the level of sphingolipids and thus present a cardiovascular protective effect. As for glycerophospholipids, including PC, PI, and PS, accumulating evidence showed that glycerophospholipids promoted the progression of cardiovascular diseases (69, 70) . Studies showed that PC was positively associated with the prevalence of CAD (71), while moderate alcohol consumption was found to reduce the level of PC (25) . Consistently, in our study, we found a decrease of PC (LPP1438, HMDB0008182) in MD-CAD when compared with HD-CAD and ND-CAD, suggesting that moderate alcohol consumption may reduce PC and lead to better cardiovascular health. In addition, it was reported that, in patients with hypertension, PI was found positively related to systolic blood pressure (SBP) (72) . What is more, researchers found that alcohol-treated rats exhibited a decrease in PI and PS than FIGURE 6 | Interrelationship between intestinal flora, serum metabolic features, and major clinical parameters associated with moderate alcohol consumption with CAD. A Sankey plot was utilized to examine the relationship among the differential OTUs and serum metabolites associated with moderate alcohol consumption, and major clinical parameters of CAD. Red connections indicate significant positive correlations, and blue connections indicate significant negative correlations (Spearman correlation analysis, P < 0.05). In the left column, light blue boxes indicate OTUs that are significantly depleted in MD-CAD when compared with ND-CAD and HD-CAD, and light orange boxes indicate OTUs that are significantly elevated in MD-CAD when compared with ND-CAD and HD-CAD. In the middle column, dark blue boxes indicate metabolic features that are significantly decreased in MD-CAD when compared with ND-CAD and HD-CAD, and dark orange boxes indicate metabolic features that are significantly increased in MD-CAD when compared with ND-CAD and HD-CAD. The right column exhibits some major clinical parameters of CAD. CK-MB, creatine kinase-MB; cTnI, cardiac tropinin I; Cr, creatine; ApoA1, apolipoprotein A1; Glu, fasting blood glucose; TG, triglyceride; FFA, free fatty acid; UA, uric acid; ApoB, apolipoprotein B; LDL-C, low-density lipoprotein cholesterol; TC, total cholesterol. control (25) . Consistently, in our study, PI (LPN 207, LPN208, both reduced in MD-CAD) was positively correlated with CK-MB, a biomarker for myocardial injury. Overall, the reduction of glycerophospholipids associated with moderate alcohol consumption tended to be beneficial for cardiovascular health.

In addition to sphingolipids and glycerophospholipids, several other metabolites also caught our attention. Firstly, it was worth mentioning that 2-phenylacetamide (PLP4198, reduced in MD-CAD) positively correlated with Bifidobacterium longum (OTU 12, depleted in MD-CAD) (P < 0.029, Rho = 0.302) and uric acid (UA) (P = 0.003, Rho = 0.409). A high concentration of UA was associated with cardiovascular diseases (73) . It was reported that 2-phenylacetamide (PLP4198) was the intermediate product in the bacteria fermentation of aromatic amino acid phenylalanine (Phe) into phenylacetic acid (PAA) (74) . The downstream product of PAA in the host is PAGln, which acts via adrenergic receptors and increases the risk of thrombosis and CAD (18) . Thus, we may infer that the reduction of 2-phenylacetamide (PLP4198) in MD-CAD may be associated with relevant bacteria, and thus improve the health of the cardiovascular system. Secondly, mevalonic acid (PLP5058, depleted in MD-CAD) was the key intermediate product in terpenoid backbone synthesis catalyzed by key enzyme 3-hydroxy-3-methyl glutaryl coenzyme A (HMG-CoA) reductase. Wang et al. observed a considerable elevation of mevalonic acid in hyperlipidemic rats and a decrease after lipid-lowering therapy with a folk medicine in China (75) . As for the key enzyme in terpenoid backbone synthesis, genomic and phylogenetic analysis revealed that multiple bacteria had the similar HMG-CoA reductase as in human beings (76). According to the results of our study, we speculated that moderate drinking might inhibit cholesterol synthesis by restraining the proliferation of bacteria with the ability to synthesize mevalonic acid. Thirdly, 1-methylhistamine (PLP6139, depleted in MD-CAD) was a biomarker of allergic response, which was the main catabolite of histamine. Depleted 1-methylhistamine concentrations were observed in moderate drinkers when compared with heavy drinkers (77) , and increased 1-methylhistamine was found in alcohol-preferring rats (78) . Besides, 1-methylhistamine was a vital cardiac substrate that activated monoamine oxidase and induced reactive oxygen species production in the heart during oxidative stress (79) . Our result was consistent with these findings and indicated that the reduction of 1-methylhistamine (PLP6139, depleted in MD-CAD) in moderate drinkers may play a role in the amelioration of allergic response and oxidative stress, owing to moderate alcohol consumption. In addition, 12-HPEPE (PLP2108, increased in MD-CAD) drew our special attention. An early study in vitro confirmed the potent platelet aggregation inhibitor 12-HPEPE could inhibit platelet aggregation and serotonin release (80) . The rise of 12-HPEPE observed in MD-CAD may suggest a positive effect on the progression of atherosclerotic plaque.

Since several metabolites were microbiome derived, we paid attention to gut microbiome alterations in MD-CAD. The taxonomic features of MD-CAD resembled that of HC, and several gut microbes in MD-CAD altered significantly when compared with HD-CAD and ND-CAD. Oscillibacter valericigenes (OTU 57, elevated in MD-CAD) caught our special attention. It was reported that Oscillibacter valericigenes was depleted in people with Crohn's disease (81) , and the decrease of the Oscillibacter genus may promote inflammation (82) . Furthermore, genus Oscillibacter was reported to be able to produce butyric acid (83) and valeric acid (84) . Butyric acid and valeric acid were both well-known beneficial substances for gut health. Thus, Oscillibacter valericigenes was considered a beneficial species. Our results indicated that Oscillibacter valericigenes (OTU 57, elevated in MD-CAD) may be promoted by alcohol and thus be beneficial to cardiovascular health. The relative abundance of Enterococcus villorum (OTU 19) was reduced in MD-CAD when compared with HD-CAD and ND-CAD. It was reported that Enterococcus faecalis was identified as a cause of hepatocyte death and liver injury (85) . Since Enterococcus villorum (OTU 19) belonged to the same genus as Enterococcus faecalis, we suspected that the depletion may be beneficial to the liver. Our result indicated that moderate alcohol consumption may possibly promote potentially beneficial bacteria and may inhibit potentially pathogenic bacteria.

In addition, the multi-omics correlation revealed that moderate alcohol consumption was beneficial for cardiovascular health possibly by affecting the gut microbiome and serum metabolome. In our study, by moderate alcohol consumption, potential beneficial gut microbes (such as Oscillibacter valericigenes OTU 57 and Paraprevotella clara OTU 95) presented the highest abundance in MD-CAD among three disease groups, while potential pathogenic gut microbes (such as Bacteroides ovatus OTU 135 and Bacteroides coprocola OTU 14) showed the lowest abundance in MD-CAD when compared with ND-CAD and HD-CAD. Furthermore, the elevation of Oscillibacter valericigenes OTU 57 and Paraprevotella clara OTU 95 were both negatively associated with sphingolipids and glycerophospholipids (such as LPP138, LPP863, and LPP911), which were known to promote the progression of cardiovascular diseases (69, 70) by activating the proliferation and migration of vascular smooth muscle cells (68) . Moreover, the depletion of potential sphingolipids and glycerophospholipids was positively associated with clinical indexes that indicate the severity of disease, such as Gensini score, CK-MB, FFA, and UA.

In a word, we came up with a new concept of "alcohol modulation through gut microbiome" in human cardiovascular health. It was found that alcohol consumption impacted many diseases causally, including ischemic heart disease, hypertensive heart disease, and stroke (86) . We speculated that a moderate amount of alcohol consumption may modulate the intestinal microecology and serum metabolites, which may play an important role in the amelioration of CAD. Alcohol consumption should be taken into considerations as an important modulator in the cardiovascular system. However, there is still a need for functional validation and exploration studies to verify the cause-and-effect relationship in the future.

Several limitations of the study need to be acknowledged. First, although the clinical features of our cohort were consistent with common patients with CAD, the sample size was relatively small. Second, the study population only included men; hence, it was uncertain whether the finding could apply to women with CAD as well. Third, untargeted metabolomics had limited accuracy in the annotation of serum metabolites. Future work needs to be done to establish the exact relationship between moderate drinking and these alterations in serum metabolites and gut microbiota.

In general, our study provided a novel insight into the effect of moderate alcohol consumption on cardiovascular health by affecting gut microbiota and serum metabolome. The impact on metabolites and microbiota in patients with CAD with moderate drinking seems to be separated from those in patients with CAD with heavy drinking or non-drinking. Drinking moderately may have more positive effects on the metabolic profiles and commensal flora of patients with CAD, which may explain how moderate drinking affects cardiovascular health.

The dataset supporting the results of this article has been deposited in the Sequence Read Archive under BioProject accession code SRP167862.

The studies involving human participants were reviewed and approved by the Ethics Review Board at the Peking Union Medical College Hospital. The patients/participants provided their written informed consent to participate in this study.

XZ, RZ, and XH conceived and designed the study and wrote the manuscript. RZ and HL contributed to the bioinformatics analysis and made the tables and figures. YF and YS conducted the literature search. HL contributed to the collections of samples and data acquisition. XH and SZ critically revised the manuscript. All authors contributed to the article and approved the submitted version.

Inflammation, atherosclerosis, and coronary artery disease

The clinical profile and pathophysiology of atrial fibrillation: relationships among clinical features, epidemiology, and mechanisms

Cardiovascular risks and benefits of moderate and heavy alcohol consumption

Association of alcohol consumption with selected cardiovascular disease outcomes: a systematic review and meta-analysis

Risk thresholds for alcohol consumption: combined analysis of individualparticipant data for 599 912 current drinkers in 83 prospective studies

Alcohol use and burden for 195 countries and territories, 1990-2016: a systematic analysis for the global burden of disease study 2016

Alcohol and coronary heart disease: a meta-analysis

Association between clinically recorded alcohol consumption and initial presentation of 12 cardiovascular diseases: population based cohort study using linked health records

Moderate alcohol intake and lower risk of coronary heart disease: meta-analysis of effects on lipids and haemostatic factors

The effect of moderate alcohol consumption on adiponectin oligomers and muscle oxidative capacity: a human intervention study

Moderate alcohol consumption, estrogen replacement therapy, and physical activity are associated with increased insulin sensitivity: is abdominal adiposity the mediator? Diabetes Care

Wine, alcohol, platelets, and the French paradox for coronary heart disease

Transplantation of an obesity-associated human gut microbiota to mice induces vascular dysfunction and glucose intolerance

Gut microbiome and its meta-omics perspectives: profound implications for cardiovascular diseases

Modulating the gut microbiome with dietary interventions to reduce cardiometabolic disease risk

Short-Chain fatty acid propionate protects from hypertensive cardiovascular damage

Gut microbial metabolite TMAO enhances platelet hyperreactivity and thrombosis risk

A cardiovascular disease-linked gut microbial metabolite acts via adrenergic receptors

Metabolic predictors of incident coronary heart disease in women

Functional metabolomics characterizes a key role for -acetylneuraminic acid in coronary artery diseases

Plasma ceramides, mediterranean diet, and incident cardiovascular disease in the PREDIMED trial (prevención con dieta mediterránea)

Sphingomyelin and progression of renal and coronary heart disease in individuals with type 1 diabetes

Changes in the human metabolome associated with alcohol use: a review

Metabolic profiling of alcohol consumption in 9778 young adults

Lipidomic signature of serum from the rats exposed to alcohol for one year

Alterations in the gut microbiome and metabolism with coronary artery disease severity

China Food Composition

The relation of alcohol intake to coronary heart disease and all-cause mortality in a beer-drinking population

Alcohol-related morbidity and mortality

Genderdependent impacts of body mass index and moderate alcohol consumption on serum uric acid-an index of oxidant stress status?

Multi-omics study reveals that statin therapy is associated with restoration of gut microbiota homeostasis and improvement in outcomes in patients with acute coronary syndrome

Pattern recognition approaches and computational systems tools for ultra performance liquid chromatographymass spectrometry-based comprehensive metabolomic profiling and pathways analysis of biological data sets

Molecular microbial diversity of an anaerobic digestor as determined by smallsubunit rDNA sequence analysis

Interactions between gut microbiota, host genetics and diet relevant to development of metabolic syndromes in mice

Accelerated dysbiosis of gut microbiota during aggravation of DSS-induced colitis by a butyrateproducing bacterium

A practical guide to amplicon and metagenomic analysis of microbiome data

VSEARCH: a versatile open source tool for metagenomics

UPARSE: highly accurate OTU sequences from microbial amplicon reads

An improved greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea

edgeR: a bioconductor package for differential expression analysis of digital gene expression data

BugBase predicts organism-level microbiome phenotypes. bioRxiv

PICRUSt2 for prediction of metagenome functions

The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases

Plasma nesfatin-1 and DDP-4 levels in patients with coronary artery disease: kozani study

Cardiac biomarkers, cardiac injury, and comorbidities associated with severe illness and mortality in coronavirus disease 2019 (COVID-19): a systematic review and metaanalysis

Anti-Inflammatory potentials of the n-hexane fraction of alstonia boonei stem bark in lipopolysaccharide-induced inflammation in wistar rats

Colchicine for prevention of atherothrombotic events in patients with coronary artery disease: review and practical approach for clinicians

Plasma sphingomyelin and subclinical atherosclerosis: findings from the multi-ethnic study of atherosclerosis

Plasma sphingomyelin level as a risk factor for coronary artery disease

Protective effect of isosteviol sodium against LPS-induced multiple organ injury by regulating of glycerophospholipid metabolism and reducing macrophage-driven inflammation

Salidroside simultaneously reduces de novo lipogenesis and cholesterol biosynthesis to attenuate atherosclerosis in mice

Gut microbiota metabolites and risk of major adverse cardiovascular disease events and death: a systematic review and meta-analysis of prospective studies

New insights into immunomodulation via overexpressing lipoic acid synthase as a therapeutic potential to reduce atherosclerosis

The beneficial effects of alpha lipoic acid supplementation on Lp-PLA2 mass and its distribution between HDL and apoBcontaining lipoproteins in type 2 diabetic patients: a randomized, double-blind, placebo-controlled trial

Fecal IgA levels are determined by strain-level differences in bacteroides ovatus and are modifiable by gut microbiota manipulation

IgA vasculitis: etiology, treatment, biomarkers and epigenetic changes

Impact of a moderately hypocaloric mediterranean diet on the gut microbiota composition of italian obese patients

Gut metagenomic and short chain fatty acids signature in hypertension: a cross-sectional study

Alcohol use disorders and the heart

Alcohol and cardiovascular health: the dose makes the poison. . . or the remedy

Intimal and medial lipids in human aortas

Arterial metabolism in primates

Untargeted metabolomics identifies a plasma sphingolipidrelated signature associated with lifestyle intervention in prepubertal children with obesity

Inhibition of de novo ceramide synthesis reverses diet-induced insulin resistance and enhances whole-body oxygen consumption

Sphingosine-1-phosphate receptors and the development of the vascular system

A new role of lipid receptors in vascular and cardiac morphogenesis

Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease

Intestinal microbial metabolism of phosphatidylcholine and cardiovascular risk

Trimethylamine N-oxide promotes apoE(-/-) mice atherosclerosis by inducing vascular endothelial cell pyroptosis via the SDHB/ROS pathway

Gut microbiota and host plasma metabolites in association with blood pressure in Chinese adults

Uric acid and cardiovascular disease: a clinical review

A gut bacterial pathway metabolizes aromatic amino acids into nine circulating metabolites

Serum metabonomics coupled with HPLC-LTQ/orbitrap MS and multivariate data analysis on the ameliorative effects of Bidens bipinnata L. in hyperlipidemic rats

Changes in the serum metabolite profile correlate with decreased brain gray matter volume in moderate-to-heavy drinking young adults

Increased brain histamine in an alcohol-preferring rat line and modulation of ethanol consumption by H(3) receptor mechanisms

Monoamine oxidase-dependent histamine catabolism accounts for post-ischemic cardiac redox imbalance and injury

Comparison of the in vitro effect of eicosapentaenoic acid (EPA)-derived lipoxygenase metabolites on human platelet function with those of arachidonic acid

Highlighting new phylogenetic specificities of Crohn's disease microbiota

A tryptophan-deficient diet induces gut microbiota dysbiosis and increases systemic inflammation in aged mice

Oscillibacter ruminantium sp. nov., isolated from the rumen of Korean native cattle

Oscillibacter valericigenes gen. nov., sp. nov., a valerate-producing anaerobic bacterium isolated from the alimentary canal of a Japanese corbicula clam

Bacteriophage targeting of gut bacterium attenuates alcoholic liver disease

The relation between different dimensions of alcohol consumption and burden of disease: an overview

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.

Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.The reviewer WL declared a shared affiliation with all authors to the handling editor at the time of the review.