Gut Microbiota and Relevant Metabolites Analysis in Alcohol Dependent Mice.

Alcohol abuse is a major public health crisis. Relative evidences supported that the gut microbiota (GM) played an important role in central nervous system (CNS) function, and the composition of them had changed after alcohol drinking. We sought to explore the changes of GM in alcohol dependence. In our study, the GM of mice with alcohol administration was detected through analyzed 16S rRNA gene sequencing and the fecal metabolites were analyzed by LC-MS. The microbial diversity was significantly higher in the alcohol administration group, the abundance of phylum Firmicutes and its class Clostridiales were elevated, meanwhile the abundance of Lachnospiraceae, Alistipes, and Odoribacter showed significant differences among the three groups. Based on LC-MS results, bile acid, secondary bile acid, serotonin and taurine level had varying degrees of changes in alcohol model. From paraffin sections, tissue damage was observed in liver and colon. These findings provide direct evidence that alcohol intake affects the composition of GM, enable a better understanding of the function of GM in the microbiota-gut-brain (MGB) axis, and give a new thought for alcohol addiction treatment.

Introduction

Alcohol abuse lead to a series of healthy problems, like alcoholic liver disease, cardiovascular diseases and depression, some works reported that gut microbiota (GM) played an important role in treating these diseases (Bull-Otterson et al., 2013; Vassallo et al., 2015). Alcohol mainly metabolizes in liver (Cederbaum, 2012) and alcohol abuse could change the intestinal permeability which would lead bacteria displacement to mesenteric lymph nodes and liver exacerbates alcoholic liver disease (Leclercq et al., 2014; Wang et al., 2016). Again, alcohol could easily penetrate the blood brain barrier, which would affect the function of brain. Recent studies found that most alcohol addicted patients caused depression symptoms and cognitive dysfunction (Corrigan and Hutchinson, 2012; Retson et al., 2015; Hendricks et al., 2017). Long-term alcohol consumption leads to oxidative damage of brain (Wang et al., 2013), however, corticosterone, regulating the hypothalamic-pituitary-adrenal (HPA), could relieve these symptoms (Tilg and Mathurin, 2016). It was reported that alcohol abuse changed neurotransmitter systems including GABAergic and glutamatergic systems (Mukherjee et al., 2008), and many studies focused on cholinergic, dopaminergic, sero-tonergic, noradrenergic, corticotrophin releasing hormone, opioid, and neuropeptide Y (Leggio et al., 2012; Cui et al., 2013; Bell et al., 2016). Alcohol dependent patients produced severe alcohol withdrawal symptoms, such as heightened responses to sensory stimuli, tremors, hallucinations, increased levels of anxiety, generalized convulsions with pain and depression (Dina et al., 2006; Conner et al., 2009; Jordaan and Emsley, 2014). Relieving withdrawal symptoms was considered a key step to treat alcohol addiction. Since medicine do not bring good curative effect for patients (Berrettini, 2016), an urgent need of new thought to solve this problem was imperative to seek.

With the interest in GM in recently years, many evidences indicated that GM were closely related to the immune system, nervous system, obesity and diabetes, and damaging GM would lead to severe pathologies, metabolic disease, cancer and irritable bowel syndrome (Ijssennagger et al., 2015; Nakatsu et al., 2015; Kim et al., 2017; Li et al., 2017; Postler and Ghosh, 2017; Sekera et al., 2017; Silverman et al., 2017; Zmora et al., 2017). Generally, the initial development and maturation of the neonatal microbiome is largely determined by maternal–offspring exchanges of microbiota (Mueller et al., 2015). It was reported that high-fat mother diet could change the GM, which would lead to mental retardation of the fetus, resulting in social barriers (Buffington et al., 2016). Some researches also indicated that GM influenced host social behaviors like stress, cognition and anxiety (Parashar and Udayabanu, 2016). The higher diversity of GM was more conducive to adapt the changes in the external environment (Rosshart et al., 2017). GM and host influenced each other through a variety of ways, including vagus nerve, microbiota-hormonal signaling and short-chain fatty acids (SCFAs) produced by GM consuming fiber (Yang et al., 2013; Ijssennagger et al., 2015; Koh et al., 2016). SCFAs like propionate and butyrate could influence intestinal gluconeogenesis (IGN), which have benefit on glucose and energy homeostasis (De Vadder et al., 2014). Diet was an important factor on the change of GM, and high dietary fiber diet had proven beneficial to the brain through butyric acid that playing an important role in microbiota-gut-brain (MGB) axis, a bidirectional neurohumoral communication system (Bourassa et al., 2016). The changes in the GM affected both stress reactivity and stress-related behaviors (Luna and Foster, 2015; Tarr et al., 2015), and mental disorders could be treated by the MGB axis (Foster and McVey Neufeld, 2013) that connected to the brain in three ways: the metabolic substrate such as SCFAs, bile acid (Devlin and Fischbach, 2015) and succinic acid (Watanabe et al., 2012) that directly affects the brain through the peripheral circulation, the vagus nerve (Sampson and Mazmanian, 2015) and the immune system (Filiano et al., 2015; Kelly et al., 2015). Serotonin, an important neurotransmitter, is related to depression, and ninety percent of serotonin biosynthesizes from colonic enterochromaffin cells (ECs) (Bellono et al., 2017). Relatively studies have reported that GM is important modulators of serotonin (Yano et al., 2015).

Many studies showed that GM composition in alcohol abusers was different from healthy people (Dubinkina et al., 2017), and GM played an important role in alcohol-dependence (Leclercq et al., 2014; Gorky and Schwaber, 2016). Admittedly, GM and depression were inseparable. Moreover, depression had been known as one of the reasons leading to addiction (Skosnik and Cortes-Briones, 2016). In this study, we established a mouse model of alcohol addiction through different feeding methods, open field exploration and light-dark transition test demonstrated the reliability of the model. High-throughput sequencing of 16S rRNA and metabolomics analysis were performed to analyze the function of GM in alcohol dependence mice. The results of these experiments demonstrate that alcohol administration takes profoundly influence on the composition of GM.

Materials and Methods

Animals

Thirty female BALB/c mice (6 weeks old, 18–23 g) were purchased from Jiesijie (Shanghai, China). All mice were housed at ambient temperature (21°C) in a room maintained on a reversed 12L:12D cycle (lights on at 9:00 AM, lights off at 9:00 PM). Purchased mice were taken 1 week to adapt the environment before starting the experiments. All animal studies have been approved by China Ethics Committee and performed in accordance with the ethical standards.

Chemicals and Alcohol

Alcohol was purchased from Maotai (Guizhou, China), executive standard GB/T26760. ELISA kits were purchased from Fankewei (Shanghai, China), chemicals of LC-MS were purchased from Thermo Fisher Scientific (San Jose, CA, United States).

Chronic Alcohol Consumption

Thirty female mice were randomly equally divided into three groups, two of them are experimental groups, another group is viewed as control group (CT). One of the experimental group named active drinking group (ADG) was treated in two-bottle drinking mode. One bottle with increasing alcohol concentration (3%, 6%, 10%, v/v) was given for each 2 days to train mice adapting alcohol, and another bottle filled with water. The placement of two bottles was exchanged everyday to avoid side preferences (Hwa et al., 2011). Another experimental group named forced drinking group (FDG) was treated in one-bottle drinking mode with increasing alcohol concentration (3%, 6%, 10%, v/v) for each 2 days to train mice adapting alcohol. After 6 days, the concentration of alcohol solution was access to 20% to feed mice for 7 weeks. Control group took no treatment. All three groups were given free food intake. At the same time of the following day, food and bottles were weighed to calculate intake.

Alcohol Withdrawal Assessments

After 8 weeks feeding, open field exploration and light-dark transition test were taken to evaluate the anxiety levels by measuring the general locomotor activity with CT mice. ADG and FDG mice were extracted alcohol solution for 24 h, and then tasting anxiety levels in the withdrawn mice by open field exploration and light-dark transition test.

Anxiety-Like Behavior Measurements

Open Field Exploration Test

Spontaneous exploration of sports were analyzed using automated activity chambers with camera (25 × 25 × 30 cm, W × L × H) (Zheng et al., 2016). Before the test, every mouse was brought to the center of automated activity chambers for 30 min everyday continued 1 week to adapt environment. Total distance of exercise and the time of stay in the center of the open field box were recorded within 5 min. A decrease in time spent in the center and total distance of exercise were considered as anxiety-like behaviors. The chamber was cleaned with 75% ethanol solution and dried before the next mouse test.

Light-Dark Transition Test

The light-dark apparatus was made up of an automated activity monitor with a dark box and insert to create an equally spaced light and dark compartment (20 × 20 × 25 cm, W × L × H). The entire apparatus was positioned in a sound-attenuating chamber. The light side was illuminated to a degree of 60 Lx, and 5 Lx in the dark side. Each mouse was placed into the light side and allowed to freely explore the chamber for 5 min. The light-dark apparatus was cleaned with 75% ethanol solution and dried before the next mouse was tested. A photo beam-based tracking system was used to track the movement and calculate the time spent in each area. Anxiety-like effects were detected by increasing time spent in the dark compartment (Alongkronrusmee et al., 2016).

Fecal Sample Collection and DNA Extraction

Fecal sample was collected in metabolism cages. After the behavioral test, mice were placed in the metabolism cages sterilized by absolute alcohol before collection. Fecal samples were stored immediately at −80°C and extracted using OMEGA soil DNA kit (OMEGA, United States). DNA quantity and quality was assessed using a NanoDrop 2000 (Thermo Fisher Scientific, United States).

High-Through Sequencing

The bacterial communities in the fecal samples were investigated by Illumina MiSeq high-throughput sequencing (US-Kyrgyzstan Biotechnology Company, Shanghai, China). The V3 and V4 regions of the 16S rDNA gene were selected for PCR. The primers were barcoded 338F (50-ACTCCTACGGGAGGCAGCA-30) and 806R (50-GGACTACHVGGGTWTCTA AT-30; H, W, and V were degenerate bases; H represented A, T or C; V represented G, A or C; W represented A or T), where the barcode was an eight-base sequence unique to each sample. The 20 μL PCR reaction mixture was composed of 4 μL of 5× FastPfu buffer, 2 μL of 2.5 μM dNTPs, 5 μM each of forward and reverse primers, 0.4 μL FastPfu Polymerase, 10 ng Template DNA, and ddH2O making up to 20 μL. The following cycling parameters were used: maintain at 95°C for 3 min, 25 cycles (95°C for 30 s, 55°C for 30 s, and 72°C for 45 s), and a final extension at 72°C for 10 min (Ning et al., 2017).

Serotonin Measurements

Serotonin levels were detected in sera by ELISA method according to the manufacturer instructions (Fankewei, Shanghai, China).

Liver, Colon Paraffin Sections

The livers and colons were soaked in 10% formalin for 24 h (Thunnissen et al., 2018), dehydrated with increasing concentrations of ethanol, embedded in paraffin and cut into 5-μm sections (Mueller-Ortiz et al., 2014). The liver and colon sections were then stained with hematoxylin and eosin (H&E), analyzed and photographed using a Leica DM2500 microscope (Leica Microsystems, United States).

UHPLC–MS Analysis of Fecal Samples

The extraction procedure was performed as previously described (Cesbron et al., 2017). Briefly, fecal samples were soaked in a solvent mixture of methanol/acetonitrile for 24 h, then ultrasound for 20 min and filtration. UHPLC-MS analysis was performed using a Q Exactive HF-X instrument (Thermo Fisher, Carlsbad, CA, United States) in combination with an UHPLC system (Thermo Fisher, Carlsbad, CA, United States). Fecal extracts were separated in Thermo Dionex Ultimate 3000 (Thermo Fisher, Carlsbad, CA, United States) with Amide (2.1 × 100 mm, 1.7 μm) column in positive and negative ionization mode respectively. In Amide column a solvent mixture of ammonium acetate/acetonitrile in pH 9.0 (A: 10 millimolar ammonium acetate/0.1% formic acid in 95% acetonitrile; B: 10 millimolar ammonium acetate/0.1% formic acid in 50% acetonitrile). Column temperature, injection volume at partial loop condition and flow rate was 40°C, 2 μL, and 0.3 mL/min, respectively. The source parameters of Q Exactive HF-X were as follows: spray voltage of 3500/3200 V (+/−); capillary temperature of 320°C; vaporiser temperature of 350°C; sheath gas flow rate of 40 Arb; auxiliary gas flow rate of 10 Arb. Data analysis and identification were performed by the software Compound Discoverer 3.0 (Thermo Fisher, Carlsbad, CA, United States). The software integrates retention time alignment, pick detection, group unknown compounds, identification and statistics node into the workflow.

Statistical Analysis

According to the different level of similarity, all the sequences clustered into operational taxonomic units (OTUs) which is a unity mark set artificially for a taxonomic unit (strain, genera, species, grouping, etc.), and the statistical analysis of biological information is usually carried out at the similar level of 97% of OTU. Similarity or difference in composition of sample communities was analyzed using principal component analysis (PCA). Statistical analysis between multiple samples was performed using One-way ANOVA. Significant differences between the two microbial communities were analyzed using Student’s t-test. Linear discriminant analysis effect size (LEfSe) was used for identifying differences in population abundance, and assessing the magnitude of the effect of each species abundance on the differences (Segata et al., 2011). Differences were considered significant if the p-value was less than 0.05.

Results

Histological Examination

The liver is the main site for alcohol metabolism and the colon is a habitat for the GM. To investigate functional and pathologic outcomes in recipient mice after alcohol administration, colon and liver sections were examined after 8 weeks feeding (Figure ). Staining of ADG showed the slight damage of colon mucosa, staining of FDG showed the shedding of colon villus, partial necrosis and colon chronic inflammatory cell infiltration. Moreover, staining of ADG and FDG showed slight congestion and loosening of liver cytoplasm.

(A–C) Representative H&E staining of colon sections. (D–F) Representative H&E staining of liver sections. (A,D) Representative H&E staining of colon and liver sections from ADG. (B,E) Representative H&E staining of colon and liver sections from FDG. (C,F) Representative H&E staining of colon and liver sections from CT.

Acquisition and Maintenance of Alcohol

Mice of ADG and FDG had less feed intake than CT group. After 8 weeks feeding, food (Figure ) and alcohol (Figure ) intake of ADG had increased and tended to stability. The results suggested that mice of ADG adapted to alcohol and took the initiative to drink, mice of FDG gradually adapted to the lifestyle of only drinking alcohol. Weight of FDG mice showed significant decrease compared with CT group (Figure ).

(A) Food intake (g) over 24 h of the three groups. (B) Alcohol solution intake (ml) over 24 h for ADG and FDG, at the first week the concentrations of alcohol were from 3, 6, 10% to 20%, and then the concentration kept 20% for 7 weeks. (C) Body weight of the three groups after 8 weeks feeding. ∗P < 0.05 and ∗∗P < 0.01.

Mice Withdrawal Performed Anxiety and Depression

We used open field exploration test and light-dark transition test to investigate the differences of the physiological state in the three groups (Figures ). After 8 weeks feeding, the mice of two groups were withdrawn from alcohol for 24 h. Both of them were performed open field exploration test and light-dark transition test to confirm whether the mice of ADG or FDG were addicted in alcohol. Mice of ADG (p-value between the two groups was 0.0004, as determined by a two-tailed t-test) and FDG (p-value between the two groups was 0.0086, as determined by a two-tailed t-test) spent more time in the dark side, the distances traveled and time in the center of the open field exploration box significantly reduced (p-value between the two groups was less than 0.0001, as determined by a two-tailed t-test) compared to the CT group. The result indicated that the alcohol groups performed anxiety and depression compared with no alcohol group and alcohol-withdrawn mice shown more anxiety and depression than before withdrawal. Compared with the CT group (Figure ), the serotonin levels of the experimental groups were increased, and FDG had a significant differences (p-value between the two groups was 0.0175, as determined by a two-tailed t-test).

(A) The serotonin levels measured after 8 weeks feeding. (B) The mice of three groups spend time in dark side. (C) The mice of three groups traveled total distance in open field exploration. (D) The mice of three groups spend time in central of open field exploration. ∗P < 0.05, ∗∗P < 0.01 and ∗∗∗P < 0.001.

Richness and Diversity of the Microbial Community

In total, approximately 1,505,392 sequence reads of 16S rRNA genes were obtained after feeding and withdrawn alcohol, and each full length was 437 bp on average. A 97% similarity cut-off was used to delineate OTUs in the downstream analyses. After subsampling, a total of 617 OTUs were acquired.

Rarefaction curves indicated that the bacterial community was well represented because the OTU level had no changes as randomly selected number of sequencing analyzed increased. The curves characterize species abundance and species uniformity, in the horizontal direction, abundance of species is reflected by the width of the curve. The higher the abundance of species, the greater the range of the curve on the horizontal axis, and the shape of the curve(smoothness) reflects the homogeneity of the species in the sample, the smoother the curve, as well as the more uniform the species distribution. Community richness is the number of bacterial species assigned by OTUs detected in the samples. Richness estimates were obtained from the observed number of species by extrapolation using estimators such as the ACE and Chao1 indices (Table ). ACE and Chao1 were estimated to be 472.9335/477.2634 in ADG, 444.8102/449.9084 in FGD and 449.4612/456.777 in CT group. Respectively, there were little differences between ADG and CT group. However, the observed richness (Sobs) was different from 409.9 in ADG and 433.9 in FDG (p-value for the Sobs indices between the two groups was 0.03368, as determined by a two-tailed t-test). The Shannon and Simpson indices were 4.441124/0.027667 in ADG, 4.37675/0.026892 in FDG and 4.411108/0.027046 in the CT group, the higher Shannon and lower Simpson index with the same change in the three groups indicated higher community diversity. The results suggested that there were no significant change in community diversity in the three groups, but FDG had an increasing tendency compared with the other two groups.

Altered Microbiota Composition in the Alcohol Group

The analysis of species composition demonstrated that bacteria from 30 samples had different richness and diversity. Hundred and twenty one genus in 30 samples were detected, and each relative abundance of genus were shown in Figure . Based on similarity in spaces abundance, cluster analysis by community Heatmap was shown in Figure , and the results suggested that FDG differed from the other two groups, ADG and CT group had a high similarity in species abundance. Interestingly, the acetate-producing Bifidobacterium was one of the increased species in this study. The species composition of the three groups could be separated clearly by PCA (Principal Component Analysis) (Figure ). Circos (Sui et al., 2016) is a visual loop diagram that describes the correspondence between a sample and a species (Figures ), which suggested that Bacteroidetes, Firmicutes, Proteobacteria and Deferribacteres had a high abundance in the three groups and each group had different abundance in the four dominant species, and also indicated changes in Clostridiales and Bacteroidales, the two orders with the highest abundance. Ratio of Firmicutes to Bacteroidetes was increased in the experimental group, abundance of Clostridiales was also observed higher in FDG than that in CT group.

Structures of bacterial community in all samples at the genus level. The abundance is presented in terms of the percentage of the total effective bacterial sequences in the sample.

Community Heatmap based on the abundance of the microbiota in the samples. The columns correspond to the samples. Each row corresponds to the microbiota in genus level. The hierarchical clustering was performed using the Euclidean metric and complete linkage.

(A) Principal component analysis (PCA) of the samples on OTU level. The fecal microbiotas of the three groups could be divided into clusters according to community composition. (B,C) The abundance was presented by the percentage on the total effective bacterial sequences in the samples.

Compared with the CT group, Odoribacter and Ruminococcaceae_UCG-014 abundance had a significant change in ADG (Figure ) (p-value for the species between the two groups was 0.04535 and 0.0004385, respectively, as determined by a two-tailed t-test). The abundance of Lachnospiraceae, Alistipes, and Odoribacter had a significant difference between FDG and the CT group (p-value for the species between the two groups was 0.00004845, 0.0005043, and 0.003309, respectively, as determined by a two-tailed t-test, Figure ). LEfSe was used to determine the taxa that best characterized each population. LEfSe scores measure the consistency of differences in the relative abundance between taxa in the groups analyzed with a higher score, thus indicating higher consistency. LDA showed distinct taxa in the microbiome of the CT group versus the alcohol administration group (Figures ).

Excessive alcohol consumption shaped the gut microbiota composition. (A,B) Species abundance calculated by Student’s t-test bar plot on genus level.

Linear discriminant analysis (LDA) gave distinct community composition on genus level. (A) LEfSe cladogram indicated differentially abundant taxa. (B) LDA scores were calculated by LEfSe of taxa differentially abundant.

Analysis of LC-MS

The result of LC-MS showed that the principal component analysis (PCA) nicely clustered FDG and the pooled QCs in alkaline condition both positive and negative ion mode. There was clearly a difference between FDG and the CT group, and ADG showed no significant difference with CT group (Figures ). As shown in Figure , the content of taurine, butyric acid, bile acid and secondary bile acids was increased in ADG. Higher concentrations of serotonin than the CT group were found in FDG. Moreover, taurine was found increased in ADG while decreased in FDG.

Principal component analysis (PCA) of fecal samples. The fecal metabolites of the three groups could be divided into clusters according to community composition. (A,B) Principal component analysis (PCA) of fecal samples.

Analysis of taurine, serotonin, bile acid and secondary bile acids in the fecal matter of the three groups. (A–C) Representative LC-MS spectrum of the butanol derivative of taurine, serotonin, bile acid and secondary bile acids in the fecal sample. (D–F) Relative abundance of taurine, serotonin, bile acid and secondary bile acids in the fecal matter of the three groups tested through LC-MS.

Discussion

Alcohol abuse caused body disease and large social burden. Generally, alcohol addicted patients have alcohol withdrawal symptoms such as anxiety, depression and other brain diseases (Becker, 2012). Currently, increasing evidences have shown that GM played an important role in brain disease, and chronic alcohol drinking could change GM composition in human body (Dubinkina et al., 2017). Herein, we developed and validated two robust mouse models with alcohol withdrawal syndrome in different levels to explore the mechanism of GM in alcohol addicts.

Notably, the data from ADG and FDG demonstrated that the diversity of GM had increased, while some commensal gut taxa had reduced. Interestingly, the results indicated that the Firmicutes (Watanabe et al., 2012), the largest phylum of bacteria, was increased. However, Bacteroidetes had decreased in ADG and FDG (Ning et al., 2017). Firmicutes is considered to be the anti-inflammatory bacteria (Natividad et al., 2015). The ratio of Firmicutes to Bacteroidetes increased in the alcohol groups, which might bring the protection to intestinal lesions. Obviously, the increase of Firmicutes abundance and relevant metabolites in the external environment could make the pH value decreased, resulting in a lower tolerance for an acidified gut environment spaces would outcompeted. Based on the data, lower feed intake in the alcohol group might be a reason to the weight decrease, furthermore efficient use of substrates was a key factor for explaining the advantage of bacteria in the competition.

In addition, Ruminococcaceae, positively with anxiety and negatively related to memory, was reported decreased in patients of hepatic encephalopathy and inflammation (Bajaj et al., 2012), and the abundance of ADG species showed decreased. Alistipes and Odoribacter, members of the Bacteroidetes phylum, were both changed in ADG and FDG. Alistipes was known to be correlated to pain in Irritable Bowel Syndrome (IBS) patients (Saulnier et al., 2011). In FDG, the abundance of Alistipes had a significant decrease compared with that in CT group, however, Alistipes was found to be more abundant in depression mice model (Maria et al., 2012). Interestingly, the abundance of Odoribacter was increased in ADG and decreased in FDG. Odoribacter was involved in the inflammatory process, IBD, Crohn’s disease, ulcerative colitis and colon cancer (Zackular et al., 2013; Jiang et al., 2018). Noteworthy, the abundance of Bifidobacterium which was involved in weight gain, body fat, fasting glucose, and insulin resistance was increased in FDG.

A series of evidences indicated that brain function and social behavior were influenced by microbial metabolites such as butyric acid (Stilling et al., 2016). Butyric acid is the main product of Firmicutes by fermenting dietary fiber. Indeed, butyric acid was a mediator of host-microbe crosstalk through energy metabolism and immune functions. A study evidenced the butyric acid could influence immune milieu of brain through changing peripheral immune system function (Filiano et al., 2015). Again, butyric acid regulated serotonin and gut hormone level in the enteric nervous system, stimulated the vagus nerve and elicited endocrine signaling. However, SCFAs are mostly absorbed in the colon and less than 5% remaining in the fecal (Topping and Clifton, 2001). In this study, we found that there was no significant change about SCFAs.

Serotonin (5-hydroxytryptamine, 5-HT), 90% 5-HT synthesized in the gut by enterochromaffin cells (ECs), myenteric neurons and mucosal mast cells (Côté et al., 2007), was a brain neurotransmitter that could regulate gastrointestinal (GI) tract and other organ systems. In addition, it was reported that 14 different 5-HT receptor subtypes (Gershon and Tack, 2007) had observed in enterocytes (Hoffman et al., 2012), enteric neurons (Mawe and Hoffman, 2013) and immune cells (Baganz and Blakely, 2013). Moreover, including immune responses (Baganz and Blakely, 2013), cardiac function (Côté et al., 2003), platelet aggregation (Mercado et al., 2013), bone development (Chabbi-Achengli et al., 2012), enteric motor and secretory reflexes were influenced by gut-derived 5-HT. As the reported shown that depression patients had low 5-HT level and GM played a key role to regulate 5-HT. In our study, 5-HT level was increased in FDG, the results suggested that GM of FDG mice might play a positive role in regulating 5-HT.

Taurine (2-aminoethane-sulfonic acid), an organic osmolyte, regulated cell volume and maintained cellular integrity in the heart, muscle, retina, and throughout the CNS (Ripps and Shen, 2012). Again, Taurine was a substrate for the production of bile salts and a key in modulating intracellular free calcium concentration (Voss et al., 2004). Taurine may be present in fecal as a result of bacterial deconjugation of bile acids (Ridlon et al., 2006). Moreover, taurine was observed in each region of the brain including the pineal (Omura et al., 1997), pons medulla, hypothalamus (Junyent et al., 2011), striatum (Fordahl et al., 2010), and cerebellum (Taranukhin et al., 2010). It was reported that taurine could ameliorate certain forms of neuropathology, and lacked of taurine would slow down cell differentiation and migration in cerebellum, pyramidal cells and visual cortex (Maar et al., 1995). In our study, the results of LC-MS indicated that taurine level was increased in the ADG group and decreased in the FDG group, the result might explain that moderate drinking apparently took positive influence on brain while excessive drinking had negative influence on brain, respectively.

Bile acids, produced in liver from cholesterol, modulated gut microbial composition through activation of innate immune genes, and also metabolized by the GM (Tq et al., 2013). Moreover, there was a study identified that Clostridium and Eubacterium, which belonged to the Firmicutes phylum, had capability to produce secondary bile acids (Kitahara et al., 2001). Deoxycholate, a secondary bile acid produced by microbial biotransformation of cholate, was reported facilitated lipid absorption and had endocrine, immunological, antibiotic effects (Islam et al., 2011). But deoxycholate exhibited negative effects on carcinogenic properties (Yoshimoto et al., 2013). Notably, deoxycholate promoted GI motility through regulated TGR5G protein-coupled receptors on ECs (Alemi et al., 2013). Particular Clostridium species were reported to possess high 7a-dehydroxylation activity to produce deoxycholate based on cholate (Atarashi et al., 2013), which was consistent with our results that GM of FGD was high abundance of Clostridia and increasing deoxycholate levels. Moreover, levels of bile acids and secondary bile acids were increased in ADG and FDG, which apparently indicated the changes of GM composition might contribute to body metabolism under alcohol damage. In summary, based on genomics and metabolomics, the gut microbial community and relevant metabolomics in mice model showed significant difference in alcohol addiction mice. Overall, our associative findings pave the way to understanding the changes of microbiota composition.

Author Contributions

QL and GW contributed to the conception and design of the work and specific experiments. GW, ZF, QZ, and JP performed the experiments. GW, JQ, DX, WZ, and XW contributed to analysis and interpretation of the data. GW, HZ, CD, QD, and LG wrote the manuscript.

Conflict of Interest Statement

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.

Table 1

Comparison of richness and diversity estimation of the 16S rRNA gene libraries for individuals at 97% similarity.

Sample ID	Reads	97%
		Sobs	Shannon	Simpson	Ace	Chao	Coverage
ADG1	43734	406	4.496882	0.022641	432.3101	444.0769	0.998621
ADG 2	51362	445	4.594691	0.020657	481.0918	508.3704	0.998278
ADG 3	46696	412	4.564669	0.019709	444.5324	450.2813	0.998448
ADG 4	47165	409	4.468058	0.023304	444.8268	440.0976	0.998509
ADG 5	49626	434	4.616726	0.020018	457.7531	460.0909	0.998802
ADG 6	48767	397	4.457343	0.025063	447.4724	446.4595	0.998176
ADG 7	42536	375	3.425182	0.135337	417.1205	413.8864	0.998265
ADG 8	42122	398	4.120826	0.036646	427.2269	427.6842	0.998608
ADG 9	49195	432	4.473285	0.027802	479.2604	483.6923	0.998153
ADG 10	53172	391	4.549837	0.021098	416.5079	424.4444	0.998933
FDG1	49139	424	4.280943	0.030396	455.2955	469	0.998385
FDG 2	55265	424	4.656459	0.019121	448.6656	449.3235	0.998756
FDG 3	54593	468	4.552514	0.024291	523.3777	528.8571	0.997838
FDG 4	51115	442	4.579632	0.023342	475.3675	483.7576	0.998438
FDG 5	49700	426	4.419177	0.026768	468.6892	470.6757	0.998262
FDG 6	51197	385	3.950514	0.049724	430.0143	423.8864	0.998245
FDG 7	52662	450	4.576815	0.021395	494.6033	492.0222	0.998125
FDG 8	51217	466	4.548253	0.023622	521.0144	543.0833	0.997885
FDG 9	55292	440	4.52235	0.02292	467.1975	470.0278	0.998769
FDG 10	52095	414	4.324587	0.035095	445.1103	442	0.998687
CT1	53334	415	4.592401	0.021985	459.2259	472.0333	0.998217
CT2	46953	386	4.283976	0.033218	414.2329	411.8	0.998759
CT3	55815	409	4.579986	0.02141	425.8144	426.0323	0.999081
CT4	61086	427	4.369427	0.040109	452.9092	462.3571	0.998896
CT5	47333	394	4.076004	0.031373	459.6481	460.7333	0.997572
CT6	42400	436	4.414414	0.033331	478.614	478.1429	0.998079
CT7	42974	401	4.417089	0.026979	435.1644	435.5882	0.998421
CT8	51645	427	4.667595	0.018428	457.6022	474	0.998665
CT9	64283	438	4.420129	0.025028	465.2132	476.6071	0.998988
CT10	42919	419	4.290059	0.037375	446.1881	470.4762	0.998704