Bilberry Anthocyanins Ameliorate NAFLD by Improving Dyslipidemia and Gut Microbiome Dysbiosis

Non-alcoholic fatty liver disease (NAFLD) is a manifestation of metabolic syndrome closely linked to dyslipidemia and gut microbiome dysbiosis. Bilberry anthocyanins (BA) have been reported to have preventive effects against metabolic syndrome. This study aimed to investigate the protective effects and mechanisms of BA in a Western diet (WD)-induced mouse model. The results revealed that supplementation with BA attenuated the serum levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT), low-density lipoprotein cholesterol (LDL-c), fat content in liver, 2-thiobarbituric acid reactive substances (TBARS) and α-smooth muscle actin (α-SMA) caused by WD. Furthermore, gut microbiota characterized by 16S rRNA sequencing revealed that BA reduced remarkably the ratio of Firmicutes/Bacteroidetes (F/B) and modified gut microbiome. In particular, BA increased the relative abundance of g_Akkermansia and g_Parabacteroides. Taken together, our data demonstrated that BA might ameliorate WD-induced NAFLD by attenuating dyslipidemia and gut microbiome dysbiosis.


Introduction
Non-alcoholic fatty liver disease (NAFLD) is a condition in which excess fat is stored in the liver without significant alcohol consumption. It embraces a pathological spectrum from simple steatosis to steatohepatitis; 10-20% have steatohepatitis (non-alcoholic steatohepatitis: NASH), which can progress to cirrhosis and hepatocellular carcinoma [1,2]. NAFLD prevalence is reported from 21% to 24.7% in USA and North America, 24% in Europe, and from 23% to 26% in Japan [3]. NAFLD is closely associated with type 2 diabetes and obesity, and both of these conditions drive progressive disease toward the more advanced stages. The mechanisms that govern hepatic lipid accumulation and the predisposition to inflammation and fibrosis are still not fully understood but reflect a complex interplay between metabolic target tissues including adipose, immune and inflammatory cells [4]. and proteinase inhibitor cocktail were purchased from Nacalai Tesque, Inc. (Kyoto, Japan). AIN-93G mineral mix and AIN-93G vitamin mix were purchased from Oriental Yeast Co., Ltd. (Tokyo, Japan). Corn starch was purchased from Sanwa Starch Co., Ltd. (Nara, Japan). Edible Acid Casein 30-60 Mesh was purchased from Meggle GmbH & Co. KG (Wasserburg am Inn., Germany). Sucrose was purchased from Hayashi Pure Chemical Industry., Ltd. (Osaka, Japan). Paraformaldehyde, hexane and 2-ethyl butyric acid were purchased from Wako Pure Chemical Co. (Osaka, Japan). Oil red O reagent was purchased from Muto Pure Chemicals Co., Ltd. (Tokyo, Japan). Antibodies against NF-E2-related factor 2 (Nrf2), Kelch-like ECH-associated protein 1 (Keap1), superoxide dismutase-2 (SOD2), and α-smooth muscle actin (α-SMA) were purchased from Santa Cruz Biotechnology (Santa Cruz, CA, USA). Antibodies against β-actin were purchased from Cell Signaling Technology (Beverly, MA, USA).

Animal Experiment Design
The animal experimental protocol was drafted according to the Guidelines of the Animal Care and Use Committee of Kagoshima University (Permission NO. A12005). Male C57BL/6N mice (5 weeks age) from Japan SLC Inc. (Shizuoka, Japan) were housed separately in cages with wood shavings bedding under controlled light (12 h light/day) and temperature (23.5 • C), and free access to water and feed. After acclimatization for 1 week, the mice were randomly divided into 4 groups (n = 5): normal diet (ND) group, ND + 2% BA group (NDBA), Western diet (WD) group, WD + 2% BA group (WDBA) (Appendix A Table A1). ND contained 3% lard and 3% soybean oil, WD contained 30% lard, 3% soybean oil and 1.5% cholesterol. Normal water was provided to ND and NDBA groups, and 4% fructose water was provided to WD and WDBA groups [14]. The supplemented dose of BA was determined according to our preliminary study without toxicity (data not shown). Mice were sacrificed at 18 weeks age after overnight fasting.

Measurement of Serum Biochemical Indexes
Blood was obtained from mice orbital veins and collected into a tube with coagulant (Separable microtubes, FUCHIGAMI, Kyoto, Japan) for 30 min at room temperature. Serum were separated from blood after centrifugation at 4000 rpm for 5 min and stored at −80 • C until use. The serum levels of alanine transaminase (ALT), aspartate aminotransferase (AST), total cholesterol (T-Cho), total triacylglycerol (TG), high density lipoprotein cholesterol (HDL-c), and glucose were measured with an automated analyzer for clinical chemistry (SPOTCHEM EZ, Arkray, Kyoto, Japan). The level of low-density lipoproteins (LDLs) was calculated using the Friedewald equation (LDL = T-Cho -HDL-c -TG/5) [15]. The serum concentrations of insulin and monocyte chemotactic protein-1 (MCP-1) were measured with an enzyme-linked immunosorbent assay (ELISA) kit (Invitrogen, San Diego, CA, USA) according to the manufacturer's instructions. The index of homeostatic model assessment for insulin resistance (HOMA-IR) was calculated with the function of fasting glucose × fasting insulin/405 [16].

Measurement of 2-Thiobarbituric Acid Reactive Substances (TBARS)
The liver TBARS concentration was measured according to the previous method modified slightly [17]. The reaction mixture contained 30 µL of liver homogenate (liver:PBS = 1:10), 30 µL of 8.1% sodium dodecyl sulfate (SDS), 225 µL of 20% acetic acid, 225 µL of 0.8% thiobarbituric acid. After incubation for 60 min at 100 • C, 150 µL of water was added and soon moved on to ice to stop the reaction. After cooling, 750 µL of the butanol-pyridine mixture (15:1) was added to the reaction mixture, and then centrifuged at 8000 rpm for 10 min. The absorbance of supernatant was read at 532 nm. The concentration of the TBARS was expressed in nmol/mg liver proteins. The range of 10~100 µM of 1,1,3,3-tetraethoxypropane was used as the standards.

Histomorphological Analysis
Mouse liver was collected at 18 weeks and fixed in 4% paraformaldehyde dissolved in phosphate buffered (PBS) at 4 • C overnight. After fixation, the liver was washed with PBS and preserved at Nutrients 2020, 12, 3252 4 of 16 4 • C in 20% sucrose dissolved with PBS overnight. After dehydration, liver was mounted with O.C.T. Compound (Sakura Finetek Japan Co., Ltd., Tokyo, Japan) under dry ice -isopentane, and preserved at −80 • C. The mounted liver was sliced at 5 µm by microtome (Yamato, Saitama, Japan). The sliced liver section was stained by oil red O for 10 min, washed by 60% isopropanol, and then observed under a fluorescence microscope (Keyence, Tokyo, Japan).

Extraction of Hepatic and Fecal Lipid
One hundred mg of liver or feces from each mouse were homogenized in hexane (200 mg/mL) with a Speed-Mill PLUS homogenizer (Analytik Jena, Jena, Germany). The supernatant was obtained by centrifugation at 8000 rpm for 10 min at room temperature. Three hundred µL of supernatant were transferred to new tube and evaporated by centrifugal concentration equipment, the residue was then weighed.

Measurement of Short-Chain Fatty Acid (SCFA) Concentration in Caecum
One hundred mg of caecum contents from each were homogenized in MQ water (100 mg/mL) by vortex for 30 s, and then centrifuged at 8000 rpm for 10 min at room temperature to get supernatant. A mixture of 100 µL of the caecum supernatant and 200 µL of 0.25 mM 2-ethylbutyric acid as internal standard were labelled with 2-nitrophenyl hydrazide (2-NPH) using a short-and long-chain fatty acid analysis kit (YMC Co., Ltd., Kyoto, Japan) according to the manufacturer's manual. 2-NPH labelled caecum supernatant was filtered through a 0.45 µm syringe filter (Millex Syringe Driven Filter Unit, MERCK, Darmstadt, Germany), and then injected into a high-performance liquid chromatography (HPLC) system (JASCO Corporation, Tokyo, Japan) with a YMC-Pack FA column (250 × 4.5 mm; YMC Co., Ltd. Tokyo, Japan). Mobile phase was acetonitrile-methanol-water (30:16:54, v/v/v), which adjusted to pH 4.5 by 0.01 M of hydrochloric acid. The injection volume was 20 µL and flow rate was kept at 1 mL/min and column oven temperature was set as 50 • C.

Hepatic Protein Extraction and Western Blot Assay
The total proteins of liver were extracted as described in our previous study [18]. Liver proteins (50 µg) were run on sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) gel before transferring to poly vinylidene di-fluoride (PVDF) membrane (GE Healthcare, Buckinghamshire, UK). The membrane was incubated with specific primary antibody and then with horseradish peroxidase (HRP)-conjugated secondary antibody. The detection was performed with a LumiVision PRO system (TAITEC Co., Saitama, Japan).

Gut Microbiota Analysis by 16S rRNA Gene Sequencing
The feces DNA were extracted by FastDNA SPIN kit for Feces (MP Bio Japan K. K., Tokyo). The composition of gut bacterial communities was analyzed by sequencing 16S rRNA genes as described in our previous study [19]. The sequences were grouped in operational taxonomic units (OTUs) with 97% similarity by QIIME 2.0. Relative abundance is calculated as the percentage of a certain bacteria OTUs to the total number of OTUs.

Statistical Analysis
Significant differences between groups were determined by Duncan's multiple range tests (IBM SPSS Statistics 24, IBM Japan, Ltd., Tokyo, Japan). A probability of p < 0.05 was considered significant. Pearson correlation was analyzed by IBM SPSS Statistics 24. The gplots package was used to draw a heatmap in RStudio (version 1.2.5001).

Body Weight and Indexes of Liver Injury
The final body weight, liver weight, epididymis fat weight, the ratio of liver weight/body weight, liver fat weight as well as total feces fat of mice fed with WD at 18 weeks was significantly higher than the mice fed with ND (Table 1). Supplementation with BA decreased significantly WD-enhanced all of these items except epididymis fat weight ( Table 1). The calories from food intake had no significant difference in all groups ( Figure A1). The serum levels of liver injury indexes including AST ( Figure 1A) and ALT ( Figure 1B) were significantly decreased by supplementation with BA in both the ND and WD groups. The serum level of MCP-1 ( Figure 1C), an inflammatory cytokine, was significantly increased in WD, and significantly reduced by supplementation with BA. The serum level of T-Cho ( Figure 1D) was also significantly increased in WD, and significantly reduced by supplementation with BA. LDL-c ( Figure 1F) was significantly decreased by supplementation with BA in both the ND and WD groups. However, the serum levels of HDL-c ( Figure 1E), TG ( Figure 1G) or glucose ( Figure 1H) were not significantly different in all groups. Additionally, the serum levels of WD-enhanced insulin ( Figure 1I) and insulin resistance ( Figure 1J) showed reduced trend by supplementation with BA.
To assess fat accumulation and excretion, we investigated the liver lipid droplets stained by Oil red O, and simultaneously measured lipid contents in both liver and feces. Oil red O stain on the liver section showed that WD caused accumulation of many lipid droplets, which were reduced by supplementation with BA ( Figure 2). Quantitative results revealed that fat content in the liver of WD was 3.6-fold higher than that of ND, and BA significantly decreased the fat contents, compared with ND and WD, respectively (Table 1).

Activation of Nrf2-Antioxidant Pathway and Amelioration of Non-Alcoholic Fatty Liver Disease (NAFLD)-Associated Markers
To evaluate the oxidative status in the liver, we measured the levels of some antioxidant factors and enzymes in the liver. The typical Western blots of three mice randomly selected from each group were shown in Figure 3A. The level of Nrf2, a key transcription factor that regulates the expression of antioxidant proteins, decreased to 0.7-fold in WD compared with ND, and recovered to 1.3-fold in WDBA ( Figure 3B). Interestingly, the level of ubiquitinated Nrf2 (Ub-Nrf2) was increased to 3.6 in the NDBA group ( Figure 3C). Additionally, Keap-1 ( Figure 3D), a negative factor of Nrf2 by ubiquitination, showed the opposite result with Nrf2. The level of SOD2 ( Figure 3E), a Nrf2 downstream target, also showed similar result with Nrf2. Next, we measured the level of liver TBARS, which is a secondary product of lipid peroxidation and has been widely adopted as a sensitive marker for lipid peroxidation in animal tissues [17]. As shown in Figure 3G, TBARS increased 6.4 times in the WD group compared with the ND group, and decreased significantly by supplementation with BA. To evaluate the oxidative status in the liver, we measured the levels of some antioxidant factors and enzymes in the liver. The typical Western blots of three mice randomly selected from each group were shown in Figure 3A. The level of Nrf2, a key transcription factor that regulates the expression of antioxidant proteins, decreased to 0.7-fold in WD compared with ND, and recovered to 1.3-fold in WDBA ( Figure 3B). Interestingly, the level of ubiquitinated Nrf2 (Ub-Nrf2) was increased to 3.6 in the NDBA group ( Figure 3C). Additionally, Keap-1 ( Figure 3D), a negative factor of Nrf2 by ubiquitination, showed the opposite result with Nrf2. The level of SOD2 ( Figure 3E), a Nrf2 downstream target, also showed similar result with Nrf2. Next, we measured the level of liver TBARS, which is a secondary product of lipid peroxidation and has been widely adopted as a sensitive marker for lipid peroxidation in animal tissues [17]. As shown in Figure 3G, TBARS increased 6.4 times in the WD group compared with the ND group, and decreased significantly by supplementation with BA. Finally, we assayed the level of α-SMA, an indicator of NAFLD-associated protein. The level was 1.3-fold higher in the WD group than that in the ND group, and decreased significantly to 1.0fold by supplementation with BA ( Figure 3F).   Finally, we assayed the level of α-SMA, an indicator of NAFLD-associated protein. The level was 1.3-fold higher in the WD group than that in the ND group, and decreased significantly to 1.0-fold by supplementation with BA ( Figure 3F).

The Concentration of Caecal Organic Acids
To elucidate the effect of BA on gut environment, we measured caecal weight and short-chain fatty acids (SCFAs) concentrations. The caecal weight was significantly increased by supplementation with BA in both the ND and WD groups ( Table 2). The concentration of lactic acid was significantly increased in both the ND and WD groups, and the concentration of butyric acid was significantly decreased in the NDBA group.

Modulation of Mouse Gut Microbiome
To further understand the effect of BA on gut environment, the composition and relative abundance of gut microbiota were determined, using high throughput 16S rRNA gene sequencing. As shown in Figure 4A,B, Chao1 (species richness) and Shannon (species evenness) decreased in the WD group compared with the ND group. Supplementation with BA decreased them in both the ND and WD groups.
Furthermore, we used principal coordinate analysis (PCoA) plots (β-diversity: between-habitat diversity) based on Jaccard to investigate the similarities in gut microbial community structure among the different groups ( Figure 4C). The percent of dataset variability explained by each principal coordinate is shown in the axis's titles (PC1: 22.6%, PC2: 12.3%). The PCoA plot indicated that the structure of gut microbiota was divided 3 groups: (1) ND (6 weeks) and ND (18 weeks), (2) NDB and WDB (both at 18 weeks); (3) WD (18 weeks). These data suggested that the BA may regulate gut microbiome. Therefore, we further investigated the changes in individual microbial species at the phylum level ( Table 3). The ratio of p_Firmicutes/p_Bacteroidetes (F/B) increased in WD, and significantly decreased by supplementation with BA. Moreover, BA significantly increased the abundance of p_Verrucomicrobia and decreased abundance of p_Deferribactere. Furthermore, BA significantly increased the abundance of g_Bacteroides;s_acidifaciens ( Figure

3.5.A Correlation Between Gut Microbiota and NAFLD-Associated Factors
To further explore the relationships between gut microbiome and NAFLD biomarkers, we generated a heatmap of Spearman correlations between gut microbial population and biomarkers ( Figure 6). Both AST and ALT, the liver injury markers, were positively correlated with f_Ruminococcaceae, and f_Deferribacteraceae, and negatively correlated with

A Correlation Between Gut Microbiota and NAFLD-Associated Factors
To further explore the relationships between gut microbiome and NAFLD biomarkers, we generated a heatmap of Spearman correlations between gut microbial population and biomarkers ( Figure 6). Both AST and ALT, the liver injury markers, were positively correlated with f_Ruminococcaceae, and f_Deferribacteraceae, and negatively correlated with f_Porphyromonadaceae and f_Verrucomicrobiaceae. T-Cho and LDL-c were positively correlated with f_ Desulfovibrionaceae. Liver fat was positively correlated with f_S24-7, f_Desulfovibrionaceae and f_Ruminococcaceae. TBARS was positively correlated with f_Erysipelotrichaceae, and negatively correlated with f_Peptococcaceae. MCP-1 was positively correlated with f_S24-7 and f_Enterococcaceae. Serum glucose level was positively correlated with f_Desulufovibrionaceae, and serum insulin level was positively correlated with f_Enterococcaceae. Fecal fat was positively correlated with f_S24-7, and negatively correlated with f_Peptococcaceae, f_Paraprevotellaceae; f_Rikenellaceae.

Discussion
In this study, bilberry anthocyanins showed preventive effects on WD-induced NAFLD by modulating serum biochemical markers, liver oxidative stress, caecal SCFA, and the gut microbiome structure and functions. In a previous study with the same bilberry anthocyanins extract (Mirtoselect), BA showed attenuation effects against the development of NASH by reducing hepatic lipid accumulation, inflammation and fibrosis [13]. Our data partially confirmed the effects and mechanisms of BA on NAFLD prevention in the previous report, and further suggested a novel mechanism that BA might prevent WD-induced NAFLD by modulating gut microbiome structure and functions.
BA has been reported to attenuate hepatic lipid accumulation although the serum markers of lipid metabolism were not investigated [13]. In the present study, we observed that BA attenuated significantly the serum levels of T-Cho ( Figure 1D) and LDL-c ( Figure 1F), and liver fat droplets ( Figure 2) and fat amounts (Table 1). These data confirmed that bilberry anthocyanins extract could attenuate hepatic lipid accumulation as reported previously [13]. However, the previous study was designed to investigate the treatment effects of Mirtoselect on NASH [13] because all of the female E3L mice received WD for 4 weeks, and only the positive female E3L mice who matched the serum levels of cholesterol and triglycerides after 4 weeks WD intake were selected to further feed with

Discussion
In this study, bilberry anthocyanins showed preventive effects on WD-induced NAFLD by modulating serum biochemical markers, liver oxidative stress, caecal SCFA, and the gut microbiome structure and functions. In a previous study with the same bilberry anthocyanins extract (Mirtoselect), BA showed attenuation effects against the development of NASH by reducing hepatic lipid accumulation, inflammation and fibrosis [13]. Our data partially confirmed the effects and mechanisms of BA on NAFLD prevention in the previous report, and further suggested a novel mechanism that BA might prevent WD-induced NAFLD by modulating gut microbiome structure and functions.
BA has been reported to attenuate hepatic lipid accumulation although the serum markers of lipid metabolism were not investigated [13]. In the present study, we observed that BA attenuated significantly the serum levels of T-Cho ( Figure 1D) and LDL-c ( Figure 1F), and liver fat droplets ( Figure 2) and fat amounts (Table 1). These data confirmed that bilberry anthocyanins extract could attenuate hepatic lipid accumulation as reported previously [13]. However, the previous study was designed to investigate the treatment effects of Mirtoselect on NASH [13] because all of the female E3L mice received WD for 4 weeks, and only the positive female E3L mice who matched the serum levels of cholesterol and triglycerides after 4 weeks WD intake were selected to further feed with Mirtoselect [13]. In the present study, we supplied BA in diet from the beginning of WD intake to estimate the preventive effects on NAFLD induction by WD. Thus, the symptoms of NASH in previous study was more critical than that of NAFLD in the present study, especially in inflammation and fibrosis. Additionally, our data also agreed with other reports in which blueberry polyphenol extract (anthocyanins compositions are similar as BA with lower contents [20]) inhibited serum LDL-c level, fecal TG and cholesterol production, but did not affect serum TG in high-fat diet (HFD)-fed C57BL/6 mice [21].
Oxidative stress is associated with the progression of NAFLD [22]. It has been reported that activation of Nrf2 inhibited NAFLD and NASH through upregulating heme oxygenase-1 (HO-1) /NAD(P)H:quinone oxidoreductase 1 (NQO1)/glutathione S-transferase (GST), then promoted hepatocellular carcinoma by competing with Keap1 or through FGFR4-GSK3β signal pathway [23]. In this study, we investigated the protein levels related to Nrf2/Keap1 antioxidant pathway. The results revealed that BA enhanced the levels of Nrf2 (a positive transcriptional factor) and SOD2 (a downstream antioxidant enzyme), and reduced the level of Keap1 (a negative transcriptional factor) and TBARS (a lipid peroxidation marker) in the liver. Therefore, the antioxidant action of BA may be involved in the preventive effects of NAFLD progression. On the other hand, it has been reported that the function of Nrf2 in lipid metabolism differs with age. In young Nrf2 knockout mice (12 weeks age) fed HFD (35% fat + 0.15% cholesterol) for 1 month, the expressions of hydroxymethylglutaryl-CoA reductase (a rate-limiting enzyme of cholesterol synthesis), acetyl-coenzyme A carboxylase (an enzyme for fatty acid synthesis), and sterol regulatory element-binding protein (a transcriptional factor of fatty acid synthesis) gene in liver were increased, compared with that in wild-type fed HFD [24]. However, in aged Nrf2 knockout mice (6 months age) fed HFD (35% fat + 0.15% cholesterol) for 3 months, their expressions were decreased [25]. In present study, all of biochemical indicators were measured at 18 weeks age, which belongs to young mouse. Therefore, an increase in Nrf2 level in WDBA might have a positive effect on the decrease in T-Cho ( Figure 1D).
It has been reported that SCFAs have two major effects on host animals: one is to increase the calorie intake, another is to suppress the inflammatory reactions in the intestines [26]. In our study, lactic acid increased with BA administration (Table 2). Since lactic acid has a much lower pKa than acetic acid, propionic acid and butyric acid, thus, pH in the gut might be lower by BA administration. Lower pH is considered to be important factor in bacterial growth inhibition [27]. In fact, anthocyanins-rich food decreased butyric acid [28] and propionic acid in a rodent model [29]. Our data indicated that BA increased lactic acid, and decreased butyric acid. Therefore, BA administration may decrease harmful bacteria (e.g., o_Clostridiales) by selectively changing SCFAs composition.
The ratio of p_Firmicutes and p_Bacteroidetes (F/B) is an important marker representing the status of gut microbiome. The high ratio of F/B are often observed in obese and HFD-fed mice [30,31] with disrupted gut microbiome. In the present study, the ratio of F/B was significantly increased in the WD group, and reduced by BA administration (Table 3) These data demonstrated that BA attenuate WD-disrupted ratio of F/B in the gut microbiome. Moreover, BA increased the relative abundances of g_Bacteroides; s_acidifaciens, g_Parabacteroides ( Figure 5B). It has reported that g_Bacteroides;s_acidifaciens may have potential for treating metabolic diseases such as diabetes and obesity [32], and g_Parabacteroides; s_distasonis attenuated toll-like receptor 4 signaling and Akt activation to block inflammation and colon tumor formation in the mice fed with HFD-azoxymethane diet [33]. Thus, BA may potentially have an anti-inflammatory effect by regulating these gut bacteria. Finally, BA significantly increased the relative abundances of g_Akkermansia; s_muciniphila ( Figure 5C) g_Akkermansia, which is a Gram-negative, strict anaerobe and mucin-degrading bacterium [34]. Interestingly, treatment with g_Akkermansia;s_muciniphila could improve mucus layer thickness and metabolic disorders caused by HFD, including fat-mass gain, metabolic endotoxemia, adipose tissue inflammation, and insulin resistance [35]. These data suggested that BA may have protective effect on gut mucosa layer by enhancing the relative abundances of g_Akkermansia.
On the other hand, PCoA analysis of the gut microbiome showed that NDBA and WDBA groups were assigned in the same group. These data suggested that BA might also selectively regulate some gut bacteria in a diet-independent manner. A heatmap of Spearman correlation between the alterations in gut microbial population and the changes in host parameters revealed that both f_Porphyromonadaceae; g_Parabateroides, and f_Verrucomicrobiaceae; g_Akkermansia were negatively correlated with AST and ALT, serum markers of liver damage. It is known that f_Porphyromonadaceae; g_Parabateroides contributed to anti-inflammation [33] and f_Verrucomicrobiaceae; g_Akkermansia contributed to enhance mucus layer thickness [35]. Thus, it is possible that BA exert hepatic protective effects by increasing the relative abundance of these bacterial to protect the damage by intestinal toxins or inflammatory factors from enterohepatic circulation.

Conclusions
Bilberry anthocyanins ameliorated WD-induced NAFLD by attenuating gut microbiome dysbiosis and dyslipidemia. These data provide insights to fully understand the bioavailability and biological function of bilberry anthocyanins.

Conflicts of Interest:
The authors declare no conflict of interest.