Simultaneous Intake of Euglena Gracilis and Vegetables Synergistically Exerts an Anti-Inflammatory Effect and Attenuates Visceral Fat Accumulation by Affecting Gut Microbiota in Mice

We determined whether the benefits provided by the consumption of Euglena gracilis (Euglena), which is a unicellular photosynthesizing green alga and rich in insoluble dietary fiber paramylon, can be enhanced by the co-consumption of vegetables that are rich in soluble dietary fiber. Nine-week-old male C57BL/6J mice were divided into four groups: group 1 received normal diet, whereas groups 2, 3 and 4 received normal diet containing 0.3% paramylon, 1.0% Euglena, or 1.0% Euglena plus 0.3% vegetables (barley leaf, kale and ashitaba), respectively. Mice were fed ad libitum until 18 weeks of age. Euglena intake significantly decreased serum markers of inflammation and co-consumption of vegetables enhanced this reduction. Notably, we observed an increase in the fraction of beneficial bacteria producing short-chain fatty acids, a reduction in harmful bacteria that cause inflammation and an increase in short-chain fatty acid production. Visceral fat accumulation was also reduced. Subsequent analyses showed that co-consumption of Euglena with vegetables reduced adipocyte area, suppressed the expression of genes related to fatty acid synthesis and increased the expression of genes related to adipocyte growth and lipolysis. Therefore, co-consumption of Euglena with vegetables enhanced its anti-inflammatory effect and the inhibitory effect on visceral fat accumulation likely by modulating the composition of gut microbiota.


Introduction
Euglena gracilis (Euglena) is a unicellular photosynthesizing green alga. It is rich in vitamins and minerals as well as in many other nutrients and it has recently attracted attention as a new health food [1]. Paramylon is an insoluble dietary fiber unique to Euglena with a triple-helical polymer structure composed of straight-chain β-1,3-glucans. Paramylon is known to attenuate manifestations of obesity and diabetes [2,3]. Paramylon also exhibits immunostimulatory activity by affecting cytokine production [4] and suppresses symptoms of atopic dermatitis and other chronic inflammatory disorders [5]. For these reasons, Euglena and paramylon have been the focus of many studies, although almost none of those revealed these beneficial effects in an efficient manner. It is therefore very important that effective methods of Euglena consumption are determined.
First, we focused on the properties of paramylon as a dietary fiber. Dietary fibers are broadly classified into soluble and insoluble dietary fibers [6]. Pectin, gum and storage polysaccharides are soluble dietary fibers, whereas cellulose, hemicellulose and lignin belong to insoluble dietary fibers [7]. Soluble dietary fiber interferes with the absorption of glucose and lipids from the intestinal tract, which reduces the concentration of sugars and lipids in blood [8]. Soluble dietary fiber promotes the production of short-chain fatty acids from fermentation by intestinal bacteria, which enhances the defensive function of the intestines because short-chain fatty acids are used as an energy source by colonic mucosal bacteria [9]. Due to its tendency to absorb water, increase stool volume and improve bowel movement, insoluble dietary fiber also prevents colon cancer by helping the body to get rid of harmful substances [10]. Paramylon and almost all other dietary fibers found in Euglena belong to insoluble dietary fibers. To obtain an optimal combination of dietary fibers, we decided to use vegetables with high soluble dietary fiber content. The soluble dietary fibers contained in vegetables such as barley leaf, kale and ashitaba have been reported to delay the speed of absorption of excess nutrients from the intestines, promote the production of short-chain fatty acids by intestinal fermentation and improve intestinal environment [11]. In addition, the intake of Euglena and vegetables is able to be expected not only to strengthen dietary fiber but also beneficial effects by micronutrients and phytochemicals. This study investigated whether combining Euglena with vegetables was an effective method of maintain health by enhancing the benefits of Euglena. The investigation was carried out in mice and measured body composition, serological and liver parameters. We also analyzed gut microbiota, performed histological examination of white adipose tissue and measured expression levels of genes related to lipid metabolism. In addition, in order to verify the utility of intake of Euglena and vegetables as an effective way to maintain health, normal mice were raised on the test diet that based on a normal diet.

Test Diets
In this study, the following test diets were prepared to examine in mice whether the effect of Euglena can be enhanced by its combination with vegetables: normal diet + paramylon, normal diet + Euglena (spray-dried powder of Euglena gracilis) and normal diet + Euglena + vegetables (mixture of spray-dried powder of Euglena gracilis, barley leaf (Hordeum vulgare), kale (Brassica oleracea var. acephala) and ashitaba (Angelica keiskei)). Diet-enriching samples were provided by Takeda Consumer Healthcare Company Limited and euglena Co., Ltd. (Tokyo, Japan). Paramylon is an indigestible dietary fiber, which is a major physiologically active substance of Euglena [1]. Vegetables that were combined with Euglena included barley leaf, kale and ashitaba, which all contain large amounts of water-soluble dietary fiber [11]. Paramylon, Euglena, or Euglena + vegetables samples were mixed with normal diet (98121701M, Research Diet Inc., New Brunswick, NJ, USA) ( Table 1). Four test diets were prepared as follows: control diet-normal diet; paramylon diet-normal diet containing 0.3% paramylon (as Euglena contains 30% paramylon [1][2][3][4][5], 0.3% paramylon content functionally corresponded to 1% of whole Euglena content); Euglena diet-normal diet containing 1.0% Euglena; and Euglena + vegetables diet-normal diet containing 1.0% Euglena and 0.3% vegetables. The composition of each test diet is shown in Table 1. The energy content per 100 g of each test diet was 388 or 389 kcal. Energy was calculated using the modified Atwater method (4 kcal of protein, 9 kcal of lipid, 4 kcal of carbohydrate per g) [12,13].

Animals
All animal procedures were performed in accordance with the Animal Experiment Guidelines of the Tohoku University and the animal protocol was approved by the Animal Use Committee at the Tohoku University (2016AgA-009) [14,15]. Male C57BL/6J mice (9 weeks old) were obtained from Clea Japan (Tokyo, Japan). The mice had access to their respective diet and distilled water ad libitum in a temperature-and humidity-controlled room with a 12/12-h light/dark cycle. After acclimatization to the normal diet for 1 week, 40 mice were randomly divided into four groups. Each group consisted of 10 mice and they were divided into two cages. Mice were fed control diet, paramylon diet, Euglena diet, or Euglena + vegetables diet (CO, paramylon, Euglena, or Euglena + V groups, respectively) for 8 weeks. At 18 weeks of age, after a 12-h fast, the mice were weighed and blood samples were collected after decapitation. To obtain serum, the blood was centrifuged (900× g, 5 • C, 15 min). The brain, cecum, heart, lung, liver, spleen, pancreas, kidney and white adipose tissues were removed and weighed. Serum and organs were stored at −80 • C until use.

Serum and Liver Biochemical Analyses
Biochemical analyses of serum and liver samples were performed as described previously [16,17]. Serum and liver triacylglycerol (TG), total cholesterol (TC) and serum phospholipid (PL), glucose, alanine aminotransferase (ALT) and aspartate transaminase (AST) levels were measured using commercial enzyme kits (Wako Pure Chemical, Osaka, Japan). Serum insulin and leptin were determined using an ELISA kit (Morinaga Institute of Biological Science, Yokohama, Japan). Liver phospholipids were measured using the procedure of Rouser et al. [18]. Serum IL-1β and IL-6 were determined using an ELISA kit (BD Biosciences, Franklin Lakes, NJ, USA). Serum and liver thiobarbituric acid reactive substances (TBARS) were measured as described previously [19]. A microplate reader (Infinite F200, Tecan Japan Co., Ltd., Kanagawa, Japan) was used for absorbance measurements.

Analysis of Gut Microbiota
DNA of intestinal bacteria was extracted from mouse stool samples and the gene encoding 16S rRNA was amplified by PCR. The composition of intestinal bacteria in feces was examined by meta 16S analysis as described previously [12,20]. Mouse feces were collected every day during the last week of the study period, DNA was extracted from all recovered feces and pooled for each group to make a stool sample. DNA extraction from mouse feces was carried out using a QIAmp DNA Stool Mini Kit (Qiagen, Hilden, Germany). PCR-based 16S rRNA gene amplicon sequencing was performed using 10 ng DNA, 10 µM of each barcoded forward and reverse primer, 2× Gflex PCR Buffer (Mg 2+ , dNTP plus) (Takara Bio, Otsu, Japan) and Tks Gflex DNA Polymerase (Takara Bio) in a total volume of 25 µL. To target 16S rRNA variable regions 3 and 4 (V3 and V4), a forward primer 341F (5-TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG CCT ACG GGN GGC WGC AG-3) and reverse primer 806R (5-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GGG ACT ACH VGG GTW TCT AAT-3) were used at 94 • C for 1 min followed by 28 cycles of 98 • C for 10 s, 50 • C for 15 s and 68 • C for 15 s. PCR products were tagged using a Nextera XT Index Kit (Illumina, San Diego, CA, USA) to distinguish sample IDs using 10 ng DNA, Nextera XT Index Primer1, Nextera XT Index Primer2, 2× Gflex PCR Buffer (Mg 2+ , dNTP plus) and Tks Gflex DNA Polymerase (Takara Bio, Otsu, Japan) in a total volume of 25 µL at 94 • C for 1 min followed by 8 cycles of 98 • C for 10 s, 60 • C for 15 s and 68 • C for 15 s. The final PCR products were mixed to give 2 ng of DNA in each sample. This final mixture was sent to Takara Bio for metagenomic analysis.

Short-Chain Fatty Acids (SCFAs) and Gamma-Aminobutyric Acid (GABA) Analyses
SCFAs were derivatized and analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) (AB SCIEX, Framingham, USA). Acetic, propionic and butyric acid were quantified and the results were expressed in µmol/g of wet weight cecum content or µmol/L of serum. The details of the derivatization procedure and LC-MS/MS analysis of SCFAs were as described by Zeng et al. [21]. Serum GABA was determined using an ELISA kit (ImmuSmol SAS, Pessac, France) as described previously [16].

Histological Analysis
For histological analysis, epididymal adipose tissue was fixed in 10% formalin and embedded in paraffin [22,23]. Vertical sections (5 µm) were cut, mounted on a glass slide, stained with hematoxylin and eosin (HE) and observed using a microscope (BZ-9000; Keyence, Osaka, Japan). The mean area of adipocytes was calculated in each group [22,23].

mRNA Expression Analysis
For real-time quantitative reverse transcriptase PCR (qRT-PCR), total RNA was isolated from epididymal adipose tissue using an RNeasy Mini Kit (Qiagen, Valencia, CA, USA) [24], eluted with 40 µL of RNase-free water and stored at −80 • C until use. To quantify gene expression, mRNA levels for fatty acid synthase (Fans), glucose-6-phosphate dehydrogenase X-linked (G6pdx), malic enzyme (Me), peroxisome proliferator activated receptor gamma (Pparg), sterol regulatory element binding factor 1 (Srebf1) and -actin (Actb) in epididymal adipose tissue were determined with a Thermal Cycler Dice Real Time System ® (Takara Bio, Otsu, Japan). This system allows real-time quantitative detection of PCR products by measuring the increase in fluorescence caused by binding of SYBR green to double-stranded DNA [25]. In brief, cDNA was made using a Prime Script ® RT Master Mix (Perfect Real Time; Takara Bio, Otsu, Japan) from total RNA in epididymal adipose tissue. The cDNA was subjected to PCR amplification using SYBR ® Premix Ex TaqTM (Perfect Real Time; Takara Bio) and gene-specific primers for Fas, G6pdx, Me, Pparg, Srebf1 and Actb (Table 2). PCR amplification was performed with activation at 95 • C for 10 s, followed by 40 cycles at 95 • C for 5 s (denaturation) and 60 • C for 31 s (extension) and dissociation at 95 • C for 15 s, 60 • C for 30 s and 95 • C for 15 s for each gene. Melting curve analysis was performed following each reaction to confirm the presence of only a single reaction product. The threshold cycle (Ct) represents the PCR cycle at which an increase in reporter fluorescence above a baseline signal can first be detected. The levels of mRNAs of genes of interest were normalized to that of Actb in test samples. Table 2. Primer pairs used for the real time qRT-PCR analysis.

Statistical Analysis
Results are expressed as the mean ± standard error of the mean (SE). Data were analyzed by the one-way analysis of variance followed by the Tukey-Kramer post hoc test, if appropriate. All analyses were performed with a significance level of α = 0.05 (p < 0.05) using BellCurve for Excel (Social Survey Research Information Co., Ltd., Tokyo, Japan).

Growth Parameters
Body weights, food intake and tissue weights of experimental animals are shown in Table 3. There were no significant differences in body weights, food intake, or energy consumption between the groups. In addition, there were no significant differences in brain, heart, lung, liver, pancreas, spleen, kidney, or mesenteric adipose tissue weights between the groups. However, the weights of perinephric, epididymal and total white adipose tissue significantly decreased in Euglena and Euglena + V groups compared to those of CO group that fed the control diet, with Euglena + V group having the lowest white adipose tissue weight (Euglena + V vs. CO: p < 0.001, 0.001 and 0.001 for perinephric, epididymal and total white adipose tissues, respectively; Euglena vs. CO: p < 0.001, 0.005 and 0.005, respectively). Although white adipose tissue weights showed nominally lower values following the intake of paramylon, the effect of diet did not reach statistical significance in those cases. These results show that visceral fat accumulation was suppressed by ingesting Euglena and that simultaneous intake of Euglena and vegetables had a particularly strong effect.

Biochemical Parameters in Serum and Liver
Biochemical parameters in serum and liver are shown in Figure 1 and Table 4. Serum IL-6, an index of inflammation, was significantly decreased in Euglena +V group compared to that in CO group (Euglena + V vs. CO, p < 0.05). Serum IL-1β, another index of inflammation, was also significantly lower in Euglena and Euglena +V groups compared to that in CO group, with Euglena + V group having the lowest IL-1β level (Euglena + V vs. CO, p < 0.05; Euglena vs. CO, p < 0.05, respectively). Liver TBARS, an index of lipid peroxidation, was significantly lower in Euglena +V group compared to that in CO group (Euglena + V vs. CO, p < 0.05). Although these parameters showed nominally lower values following the intake of Euglena or paramylon, the effect of diet did not reach statistical significance in those cases. Therefore, we concluded that inflammation was suppressed by the simultaneous intake of Euglena and vegetables. CO, control diet group; Paramylon, the group that received the normal diet containing 0.3% paramylon; Euglena, the group that received the normal diet containing 1.0% Euglena; Euglena + V, the group that received the normal diet containing 1.0% Euglena and 0.3% vegetables.

Gut Microbiota
The above results showed that simultaneous intake of Euglena and vegetables had a greater impact than consumption of Euglena only. Thus, the beneficial effects appeared to be enhanced by vegetables. To clarify the factors involved in the enhancing effects of vegetables, changes in gut microbiota in Euglena and Euglena + V groups were examined. The 16S rRNA region of DNA extracted from feces was amplified and sequence analysis was performed using the next generation sequencer. Homology searches and lineage classification analyzes were performed on the 16S rRNA database using the obtained sequences and the genus levels of the microorganisms in the samples were determined. Detailed analysis of gut microbiota identified nine taxonomic units whose proportion was at least twice higher or lower in Euglena + V group compared with that in Euglena group (Table 5). Erysipelotrichaceae; Other, YS2;f__;g__ and Lactobacillus were twofold more abundant in Euglena + V group compared with their levels in Euglena group. Furthermore, Staphylococcus, Jeotgalicoccus, Sporosarcina, Streptococcus, Sutterella and Clostridiaceae; Other were one-half or even less frequent in Euglena + V group compared with their levels in the Euglena group. It was not possible to determine the family and the genus of YS2;f__;g__. These differences in gut microbiota composition caused by the simultaneous intake of Euglena and vegetables likely were one of the factors underlying different effects of diets on physiological functions.

SCFAs and GABA Levels in Feces and Serum
Given that gut microbiota substantially changed due to simultaneous intake of Euglena and vegetables, the levels of SCFAs and GABA, whose production is promoted by intestinal bacteria, were measured (Figures 2 and 3). Acetic, propionic and butyric acids levels in cecum content were significantly increased in Euglena + V group compared with the levels in CO group (Euglena + V vs. CO, p < 0.01, 0.01 and 0.05, respectively). Although SCFAs levels showed nominally higher values following the intake of Euglena or paramylon, the effect of diet did not reach statistical significance in those cases.
Acetic, propionic and butyric acids levels in serum significantly increased in the Euglena + V group compared with those in CO group (Euglena + V vs. CO; p < 0.05, 0.05 and 0.01, respectively). Although SCFAs levels showed nominally higher values following the intake of Euglena or paramylon, the effect of diet did not reach statistical significance in those cases. GABA level in serum was significantly higher in Euglena + V group than in CO group (Euglena + V vs. CO; p < 0.01). These results showed that the production of SCFAs and GABA was promoted by the simultaneous intake of Euglena and vegetables.  GABA level are presented as the mean ± standard error of the mean, n = 10. Different superscript letters indicate significantly different means at p < 0.05 (refer to "Result" in detail). CO, control diet group; Paramylon, the group that received the normal diet containing 0.3% paramylon; Euglena, the group that received the normal diet containing 1.0% Euglena; Euglena + V, the group that received the normal diet containing 1.0% Euglena and 0.3% vegetables.

Histological Analysis of Epididymal Adipose Tissue
The large changes in adipose tissue weight suggested that the groups might differ in adipocyte size. Histological observations of epididymal adipose tissue, using HE staining (Figure 4a), showed that in Euglena and Euglena + V groups, adipocytes were smaller than in CO group. Quantification of adipocyte area showed a significantly lower value in Euglena + V group compared to that in CO group (p < 0.05) (Figure 4b). Although adipocyte area showed nominally smaller values following the intake of Euglena or paramylon, the effect of diet did not reach statistical significance in those cases. These results show that the expansion of adipocytes area was suppressed by Euglena consumption, with a particularly strong effect of the simultaneous intake of Euglena and vegetables.

mRNA Levels of Lipid Metabolism-Related Genes in Epididymal Adipose Tissue
The apparent suppression of visceral fat accumulation in the Euglena + V group prompted us to examine mRNA levels for lipid metabolism-related genes in epididymal adipose tissue that play a central role in visceral fat accumulation (Table 6). mRNA levels of G6pdx (fatty acid synthesis pathway enzyme) was lower in Euglena + V groups compared to that in CO group, which suggests that simultaneous intake of Euglena and vegetables suppressed the fatty acid synthesis system in the adipose tissue (Euglena + V vs. CO; p < 0.01). mRNA level of Me (fatty acid synthesis pathway enzyme) and Srebf1 (transcription factor for these genes) were lower in Euglena + V group compared to those in paramylon group (Euglena + V vs. paramylon, p < 0.05 and 0.05, respectively). Euglena + V group had the lowest mRNA levels for these genes, showing that simultaneous intake of Euglena and vegetables suppressed the fatty acid synthesis system in the adipose tissue. mRNA levels for Pparg, which promotes the cell division and Hsl, which promotes lipolysis, increased in Euglena + V group compared to that in CO group (p < 0.05 in both cases), which suggests that simultaneous intake of Euglena and vegetables promoted cell division and lipolysis in the adipose tissue. There were no significant differences in Fasn mRNA level between the groups. (b) Adipocyte size ratio values are presented as the mean ± standard error of the mean, n = 10. Different superscript letters indicate significantly different means at p < 0.05 (refer to "Result" in detail). CO, control diet group; Paramylon, the group that received the normal diet containing 0.3% paramylon; Euglena, the group that received the normal diet containing 1.0% Euglena; Euglena + V, the group that received the normal diet containing 1.0% Euglena and 0.3% vegetables. Table 6. mRNA expression level in white adipose tissue (ratio).

Discussion
The objective of this study was to determine whether it is possible to enhance the benefits provided by consumption of Euglena, which is rich in insoluble dietary fiber, by consuming it with vegetables that are rich in soluble dietary fiber. Euglena intake caused a significant reduction in serum inflammatory markers in mice and this reduction was enhanced by co-ingestion of Euglena with vegetables. In the search for the mechanism of these diet-induced changes, we examined the effect of Euglena and vegetable co-ingestion on gut microbiota and found an increase in beneficial bacteria that produce short-chain fatty acids, a reduction in harmful bacteria that cause inflammation and an increase in the actual amount of short-chain fatty acids produced. We also observed a reduction in visceral fat in mice and this reduction was enhanced by co-ingestion of Euglena with vegetables. Histological examination of white adipose tissue and analysis of the expression of genes related to lipid metabolism revealed that co-ingestion of Euglena with vegetables caused a reduction in the area of adipocytes, suppressed the expression of genes related to fatty acid synthesis and increased the expression of genes related to adipocyte growth and lipolysis. Based on these findings, we propose that co-ingestion of Euglena with vegetables potentiates anti-inflammatory effect and enhances its suppressive effect on visceral fat accumulation. In addition, these results suggest that changes in gut microbiota are one of the mechanisms involved in these effects.
Although there were no significant differences in final body weight, mice that consumed paramylon tended to have a lower final body weight compared to the mice in CO group. Because the same result was also observed in mice that consumed Euglena only, this phenomenon presumably was due to paramylon presence. White adipose tissue weight also generally correlated with final body weight, with particularly large inter-group changes being observed around the kidneys and testicles. In a study of a type 2 diabetes model, visceral fat accumulation in rats that ingested paramylon or Euglena was suppressed and the effect was more prominent in rats that ingested Euglena compared to those that consumed paramylon [1,3]. This difference is presumably due to the presence in Euglena of various other beneficial substances in addition to paramylon. The same tendency was also observed in the present study, as with the effects of Euglena were further enhanced by the co-ingestion of Euglena with vegetables. Based on these findings, enhancement of the effect of Euglena by its co-ingestion with vegetables was likely due to the presence of various additional beneficial constituents. Euglena is rich in vitamins and minerals as well as in many other nutrients [1]. Although these ingredients have a beneficial effect, they had not reached the amount that each ingredient exerts its effect in this study. Therefore, it was thought that the interaction of many ingredients is important. We previously showed that the diet of a Japanese person in 1975 was more beneficial for health maintenance than the diet of a modern day Japanese person [13,[26][27][28]. The 1975 diet was characterized in one aspect by a diverse variety of food constituents and by the consumption of these constituents in small quantities. This led us to presume that consuming a variety of constituents in small amounts was very effective for health maintenance and therefore, co-consumption of Euglena with vegetables in the present study represented a positive enhancement of the diversity of food constituents.
In the present study, serum concentration of the inflammatory marker IL-6 and IL-1β [29] was reduced by Euglena intake and this reduction was enhanced by the co-consumption of Euglena with vegetables. This finding showed that vegetables enhanced the anti-inflammatory effect of Euglena [1,5]. Maintaining inflammatory markers at low levels prevents the onset of aging-associated diseases and is useful for health maintenance [29]. Thus, this finding suggested that the co-consumption of Euglena with vegetables is useful for achieving a long and healthy life. TBARS levels were reduced significantly by the co-consumption of Euglena with vegetables, indicating the possibility that this diet also reduced oxidative stress.
Adipocyte hypertrophy is the basis of lifestyle diseases such as diabetes mellitus, hyperlipidemia and arteriosclerosis [30]. Therefore, it is very important to maintain small adipocyte size. Histological examination of epididymal adipose tissue revealed the tendency of beneficial adipocyte-related changes in visceral fat weight. Although it is necessary to confirm our present results in animal models of obesity, it is possible that the co-consumption of Euglena with vegetables can prevent lifestyle diseases. Analysis of gene expression related to lipid metabolism in epididymal adipose tissue further confirmed the same trend: expression levels of G6pdx and Srebf1 mRNAs encoding lipid accumulation-promoting molecules [31] were lower in Euglena + V group, whereas expression levels of Pparg [32], which promotes adipocyte growth and induces adipocyte size reduction and of Hsl [33], which induces lipolysis and adipocyte size reduction, exhibited the opposite tendency. The above findings showed that simultaneous intake of Euglena with vegetables enhanced the effects of Euglena.
The effect of the co-consumption of Euglena with vegetables on gut microbiota was examined to elucidate part of the mechanism of the anti-inflammatory effect of Euglena and its suppressive effect on visceral fat accumulation. The co-consumption of Euglena and vegetables caused a substantial increase in the following taxonomic groups: Erysipelotrichaceae; Other, YS2;f__;g__ and Lactobacillus. The fraction of the bacterial subclass Erysipelotrichaceae; Other has been shown to increase following the consumption of β-glucans and inulin [34]. These bacteria participate in the production of short-chain fatty acids [35]. In the present study, the fraction of Erysipelotrichaceae; Other was increased by the consumption of vegetables rich in soluble dietary fiber, which also caused an actual increase in short-chain fatty acid concentrations. The beneficial effects of short-chain fatty acids include the anti-inflammatory action due to strengthened intestinal tract barrier function and anti-obesity effect due to increased energy consumption through GPR41 and GPR43 in the liver [36][37][38]. In this study, the co-consumption of Euglena with vegetables was associated with higher serum concentrations of short chain fatty acids and lower visceral fat accumulation. These findings suggest that the effects of vegetables observed in this study may be attributed to enhanced production of short-chain fatty acids. YS2;f__;g__ bacteria were shown to promote the biosynthesis of vitamins B and K [39,40]. It is possible that ingestion of vegetables led to greater uptake of vitamins with beneficial effects for the body due to the associated increase in YS2;f__;g__ amount. Lactobacillus produces lactic acid by sugar metabolism in the intestines and promotes the production of GABA and other useful substances [41,42]. Lactic acid suppresses the growth of harmful bacteria susceptible to acidic conditions and exerts an anti-inflammatory effect. We confirmed that the co-consumption of Euglena with vegetables resulted in elevated levels of GABA, a substance with a stress-reducing effect [16,43]. It is possible that the increase in Lactobacillus observed in this study was one of important effects of vegetable consumption. Furthermore, the co-consumption of Euglena with vegetables caused a substantial reduction in the genera Staphylococcus, Jeotgalicoccus, Sporosarcina, Streptococcus, Sutterella as well as in the representatives of the family Clostridiaceae; Other. Staphylococcus aureus causes a variety of diseases and belongs to the genus Staphylococcus [44]. In the present study, compared to the fraction of Staphylococcus in control group, the representatives of this genus were more abundant in mice that consumed paramylon but less abundant in mice that consumed Euglena and even less frequent after the co-consumption of Euglena with vegetables. The fractions of Jeotgalicoccus and Sporosarcina in the gut microbiota were shown to increase during hepatitis induced by artificial sweeteners [45]. An increase in the number of microorganisms of these genera may therefore promote inflammation. In the present study, differences in the fractions of these genera between experimental groups exhibited similar pattern to that of Staphylococcus: their levels were higher after paramylon ingestion, lower after Euglena ingestion and decreased even further after the co-consumption of Euglena with vegetables. Streptococcus causes streptococcal infection [46]. Therefore, an increase in the abundance of this genus is also presumed to promote inflammation. In the present study, changes in Streptococcus levels were similar to those in Staphylococcus: Streptococcus levels were (relative to control group) higher after paramylon consumption but lower in animals that received Euglena in the diet and minimal after the co-consumption of Euglena with vegetables. Based on the changes observed in these four genera, the co-consumption of Euglena with vegetables likely protected against the negative effects of paramylon. Sutterella fraction was shown to correlate positively with intestinal inflammation [47]. Furthermore, the family Clostridiaceae; Other includes harmful bacteria that cause food poisoning and other conditions [48]. In the present study, the abundance values of these taxa after feeding mice with paramylon, Euglena and Euglena + V vegetables were progressively lower than in control group. These findings suggest that the co-consumption of Euglena with vegetables enhanced the positive effects of paramylon. In particular, the abundance of beneficial bacteria was higher, whereas the fractions of harmful bacteria were lower in mice that consumed Euglena with vegetables, which may be one of the mechanisms of the positive effect of vegetables.
In this study, no significant effect of paramylon intake was observed compared to the parameters measured after control diet intake. The apparent lack of effect could be due to the fact that the dose of paramylon in this study was small compared to those used in previous studies. It is likely that increasing paramylon intake would make its effect more pronounced, because the tendencies of lower inflammatory marker concentration and suppressed visceral fat accumulation were observed following paramylon intake. Therefore, studies similar to ours in which paramylon effects are shown more efficiently, will be very important.

Conclusions
The above findings showed that the co-consumption of Euglena with vegetables enhanced its anti-inflammatory effect and suppressive effect on visceral fat accumulation and suggested that modulation of gut microbiota could be one of the mechanisms of these changes. It is therefore extremely important that more effective methods of consuming Euglena-based products are determined. In future, we would like to establish more effective methods of introducing Euglena and to examine potential positive effects of Euglena and vegetable co-consumption in animal models of obesity and verify these effects in humans.