Next Article in Journal
A Multi-Year Examination of Gardening Experience and Fruit and Vegetable Intake During College
Next Article in Special Issue
New and Preliminary Evidence on Altered Oral and Gut Microbiota in Individuals with Autism Spectrum Disorder (ASD): Implications for ASD Diagnosis and Subtyping Based on Microbial Biomarkers
Previous Article in Journal
The Potential for Renal Injury Elicited by Physical Work in the Heat
Previous Article in Special Issue
The Interplay between Maternal and Post-Weaning High-Fat Diet and Gut Microbiota in the Developmental Programming of Hypertension
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Regular Kefir Consumption on Gut Microbiota in Patients with Metabolic Syndrome: A Parallel-Group, Randomized, Controlled Study

by
Ezgi BELLIKCI-KOYU
1,2,
Banu Pınar SARER-YUREKLI
3,
Yakut AKYON
4,
Fadime AYDIN-KOSE
5,
Cem KARAGOZLU
6,
Ahmet Gokhan OZGEN
3,
Annika BRINKMANN
7,
Andreas NITSCHE
7,
Koray ERGUNAY
4,
Engin YILMAZ
8 and
Zehra BUYUKTUNCER
1,*
1
Department of Nutrition and Dietetics, Faculty of Health Sciences, Hacettepe University, Ankara 06230, Turkey
2
Department of Nutrition and Dietetics, Faculty of Health Sciences, Izmir Katip Celebi University, Izmir 35620, Turkey
3
Department of Endocrinology, Faculty of Medicine, Ege University, Izmir 35040, Turkey
4
Department of Medical Microbiology, Faculty of Medicine, Hacettepe University, Ankara 06230, Turkey
5
Department of Biochemistry, Faculty of Pharmacy, Ege University, Izmir 35040, Turkey
6
Department of Dairy Technology, Faculty of Agriculture, Ege University, Izmir 35040, Turkey
7
Robert Koch Institute; Center for Biological Threats and Special Pathogens 1 (ZBS-1), Berlin 13353, Germany
8
Department of Medical Biology, Acıbadem Mehmet Ali Aydınlar University, Istanbul 34752, Turkey
*
Author to whom correspondence should be addressed.
Nutrients 2019, 11(9), 2089; https://doi.org/10.3390/nu11092089
Submission received: 13 August 2019 / Revised: 21 August 2019 / Accepted: 22 August 2019 / Published: 4 September 2019
(This article belongs to the Special Issue Nutrition, Microbiota and Noncommunicable Diseases)

Abstract

:
Several health-promoting effects of kefir have been suggested, however, there is limited evidence for its potential effect on gut microbiota in metabolic syndrome This study aimed to investigate the effects of regular kefir consumption on gut microbiota composition, and their relation with the components of metabolic syndrome. In a parallel-group, randomized, controlled clinical trial setting, patients with metabolic syndrome were randomized to receive 180 mL/day kefir (n = 12) or unfermented milk (n = 10) for 12 weeks. Anthropometrical measurements, blood samples, blood pressure measurements, and fecal samples were taken at the beginning and end of the study. Fasting insulin, HOMA-IR, TNF-α, IFN-γ, and systolic and diastolic blood pressure showed a significant decrease by the intervention of kefir (p ≤ 0.05, for each). However, no significant difference was obtained between the kefir and unfermented milk groups (p > 0.05 for each). Gut microbiota analysis showed that regular kefir consumption resulted in a significant increase only in the relative abundance of Actinobacteria (p = 0.023). No significant change in the relative abundance of Bacteroidetes, Proteobacteria or Verrucomicrobia by kefir consumption was obtained. Furthermore, the changes in the relative abundance of sub-phylum bacterial populations did not differ significantly between the groups (p > 0.05, for each). Kefir supplementation had favorable effects on some of the metabolic syndrome parameters, however, further investigation is needed to understand its effect on gut microbiota composition.

1. Introduction

Metabolic syndrome (MetS) is a pathologic condition that includes abdominal obesity, insulin resistance, dyslipidemia, and arterial hypertension [1]. Each component of MetS is known as a risk factor for the development of type 2 diabetes and cardiovascular diseases. It was found that the risk of type 2 diabetes was five times, the risk of cardiovascular disease was two times, and the risk of death was one-half times higher in individuals with MetS compared to those without the syndrome. Due to its high prevalence and related health problems, the MetS is currently considered as a significant public health problem [1,2].
MetS has a multi-factorial etiology comprising complex interactions between genetic predispositions and environmental factors including diet, physical activity, and other lifestyle factors [3,4]. Since Turnbaugh et al. showed the linked between gut microbiota and obesity, there has been growing evidence that suggests a causal relationship between gut microbiota and the components of MetS [5]. Primary, the low-grade chronic inflammation state in MetS has been explained by the metabolic endotoxemia that was a result of gut dysbiosis [6,7]. Most of the animal and human studies have reported that obesity and insulin resistance are associated with an altered ratio of Firmicutes and Bacteroidetes [8,9]. In addition to the effects on immune function, the gut microbiota also exerts its role through the influence on host energy metabolism and gut barrier integrity [10,11]. Therefore, the gut microbiota has been suggested as a potential target to modify the risk factors that contribute to conditions of MetS.
The modification of diet using prebiotics and probiotics has been suggested as a useful strategy to improve metabolic health via the modulation of gut microbiota. Although the effects of probiotic and prebiotic supplementation on metabolic health have been examined in previous studies, the results are inconsistent due to the choice of probiotic strain, formulation of the probiotic, outcome of interest, and duration of the intervention [12,13,14]. Furthermore, ingestion of probiotics through traditional fermented foods has not been widely examined in terms of their efficiency on MetS components. Kefir is a fermented milk product, traditionally produced with kefir grains that have a specific combination of bacteria and yeasts [15,16]. Microbial composition of kefir varies depending upon the type of kefir grains, the type and composition of milk, culture medium, fermentation period and temperature, and also storage conditions [17]. Lactobacillus, Lactococcus, Streptococcus Leuconostoc, and acetic acid bacteria are the most common bacteria; and Saccharomyces, Kluyveromyces, and Candida species are mostly found yeasts in kefir [18]. Animal studies have suggested that kefir has anticarcinogenic, antimicrobial, anti-inflammatory activities, and thus may ameliorate MetS components [19,20,21,22,23]. However, there is still limited clinical evidence for its potential effects on MetS patients. To our knowledge, especially, the effects of kefir on MetS components via the modulation of gut microbiota have not been examined widely in clinical settings. To address the research gap, this study aimed to investigate the effects of daily kefir consumption on gut microbiota composition and their relation with the components of metabolic syndrome in adults with MetS.

2. Materials and Methods

2.1. Subjects

Subjects with MetS, aged 18–65 years, were recruited from the outpatient clinic of the Department of Endocrinology and Metabolism at the Ege University, Izmir, Turkey. MetS was diagnosed using the IDF-2005 guidelines [24]. The eligibility of a subject was confirmed following a physical examination by the research endocrinologist and a nutritional assessment by the research dietitian in the screening period. Adults were excluded if they (1) were using antibiotics in the past 1 month or during the intervention period, (2) were using dietary supplements (probiotic, prebiotic, or symbiotic) during the past three months or during the intervention period, (3) were pregnant or lactating, (4) had severe liver, kidney, heart, or immune deficiency, (5) had chronic gastrointestinal system diseases, type 1 diabetes or cancer, (6) had allergy to the dairy products or lactose intolerance, (7) were currently taking prescribed drugs that can modulate lipid profile or glycaemic control, and (8) did not comply with the consumption of test drinks.
The compliance was assessed by interviewing the participants and reviewing the record of their consumption in each visit. Non-compliance was defined as consuming < 80% of the scheduled serving during the study period.
The study protocol was approved by the Ethics Committee of Clinical Research at Ege University Faculty of Medicine (15-2.1/14) and registered at clinicaltrials.gov (NCT03966846). All procedures were performed in accordance with the Declaration of Helsinki. Written informed consent was obtained from all participants.

2.2. Study Design

A parallel-group, randomized, controlled clinical trial was performed. A total of 40 eligible participants were randomized, and 20 participants in each group were allocated to intervention. Five participants in the kefir group and four participants in the unfermented milk group did not receive the allocated intervention due to medical conditions and not providing the fecal samples, and also three participants in the kefir group and six participants in the unfermented milk group discontinued the intervention due to taking antibiotics and declined consent. Therefore, the study was completed with 22 participants and an allocation ratio of 55%. The recruitment and follow up of participants were conducted between March 2015–July 2017. Participants were randomized into two groups (kefir group and unfermented milk group as control) by the research physicians using a stratified block randomization method. The random allocation sequence was provided by the Department of Biostatistics, Hacettepe University. Participants visited the research center 5 times in total. The first visit included the screening of individuals in terms of inclusion and exclusion criteria. The second visit (Week 0) included recording general characteristics, medical history, and lifestyle behaviors of participants, assessing the nutritional status of participants using 24-h dietary recall and anthropometrical measurements, collecting the initial blood and fecal samples, measuring the blood pressure, and also proving information about the consumption and storage of test drinks. The third (Week 4) and forth (Week 8) visits included the assessment of the compliance in terms of consumption of test drinks and dietary intake. The fifth visit (Week 12) included the assessment of the nutritional status of participants using 24-h dietary recall and anthropometrical measurements, collection of the final blood and fecal samples, and the measurement of blood pressure (Figure 1).
The primary outcome of the study was the change in the relative abundance of microorganisms in gut microbiota by regular kefir consumption. The potential correlations between the changes in dietary intake, anthropometrical measurements, biochemical parameters, or blood pressure and the change in microbiota composition were all secondary outcomes.

2.3. Intervention

During a 12-week intervention period, kefir group (n = 12) received kefir (180 mL/day) while control group (n = 10) received unfermented milk (180 mL/day) regularly. Participants were asked to maintain their habitual diet and physical activity. Additional products that contain probiotics were not allowed during the intervention period. No dietary supplement use was recorded before or during the study.

2.4. Test Drinks

Two dairy products (kefir and unfermented milk) were tested in parallel groups. Kefir was prepared using the culture of DC1500I (Danisco, Olsztzyn, Poland) containing Lactococcus lactis subsp. lactis, Lactococcus lactis subsp. cremoris, Lactococcus lactis subsp. diacetylactis, Leuconostoc mesenteroides subsp. cremoris, Lactobacillus kefir, Kluyveromyces marxianus, and Saccharomyces unisporus at Ege University Faculty of Agriculture, Department of Dairy Technology. Kefir was derived from the full-fat (3.5%) homogenized and pasteurized (at 85 °C) milk that was used as a control drink at the same time. The beverages were distributed and stored at 4 °C. The test drinks were received to participants twice a week, and they consumed the test drinks between 1 and 4 days of post-production.

2.5. Dietary Assessment

Dietary intake was assessed using 24-h dietary recall method by research dietician in each visit. A photographic atlas of food portion sizes was used to clarify the amounts of food items consumed. Dietary energy, macro- and micronutrient intakes were analyzed using BeBIS software (Ebispro for Windows, Stuttgart, Germany; Turkish Version BeBIS, Nutrition Information System, Version 8).

2.6. Anthropometrical Measurements

Body weight and composition (fat mass and fat-free mass) were measured by Tanita BC418 (USA), and height was measured by a calibrated stadiometer (Nan Tartı, TR). Body Mass Index (BMI) was calculated by dividing body weight (in kilograms) by the square of height (in meter). The waist circumference was measured at the midpoint between the lower ribs and the iliac crest, and hip circumference was measured horizontal at the largest circumference of hip. Waist-to-hip ratio (WHR) was calculated.

2.7. Biochemical Analysis and Blood Pressure

Venous blood samples were drawn after a 10-h overnight fasting excluding only the water at the Visit 1, Visit 2 (Week 0), and Visit 5 (Week 12). Serum glucose, insulin, HbA1c, total cholesterol, high-density lipoprotein (HDL-C), low-density lipoprotein (LDL-C), triglycerides, homocysteine, high-sensitivity C-reactive protein (hs-CRP), alanine aminotransferase (ALT), aspartate aminotransferase (AST), and gamma-glutamyl transferase (GGT) were analyzed at Ege University, Hospital of Medical School, Laboratory of Clinical Biochemistry. All biomarkers were analyzed using routine methods by Roche/Cobas analyzer series. Serum concentrations of tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), interleukin-10 (IL-10) and interferon-gamma (IFN-γ) were determined by enzyme-linked immunosorbent assay (ELISA) using standard kits, and the analyses were conducted as described by the manufacturer (DIAsource ImmunoAssays, Louvain-la-Neuve, Belgium). Insulin resistance was assessed using Homeostatic Model Assessment (HOMA-IR) model calculated with the equation of “the fasting insulin level (μU/L) × fasting plasma glucose (mg/dL)/405”. Systolic (SBP) and diastolic blood pressure (DBP) were measured at the brachial artery of right upper arm after 15 min rest. Both blood pressures were measured twice at 5-min intervals and recorded on average.

2.8. Specimen Processing, 16S rRNA Amplification and Sequencing

Fecal samples from individuals enrolled in the study were collected in sterile containers and kept frozen at −80 °C. A sterile spatula was used to obtain 4–5 pieces of frozen chunks from the surface and internal portions of the specimen. They were combined for a 150–200 mg total weight for each and mixed by vortexing. Following a bead-beated step described by Tomas et al. and Wu et al. [25,26], DNA was extracted using Qiagen Stool Mini Kit (Qiagen, Hilden, Germany) as directed by the manufacturer. DNA amount of 50 ng/µL was prepared for each specimen, using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Hennigsdorf, Germany).
The 16S rRNA sequences were amplified using previously-described primers targeting the V3-V4 region, frequently used to study bacterial diversity [27], with Illumina adapter overhang sequences added, as directed by the manufacturer. Attachment of sequencing adapters to PCR products, amplification and library preparation were performed using the Nextera XT Index and Nextera DNA Library Prep kits (Illumina, San Diego, CA, United States), as suggested by the manufacturer. Product clean-up, library quantification, and optimization were carried out using the Agencourt AMPure XP reagent (Beckman Coulter Biosciences, Krefeld, Germany) standard protocol and a Qubit 2.0 Fluorometer (Thermo Fisher Scientific, Rochester, NY, USA). The sequencing runs were performed in an Illumina MiSeq sequencer (Illumina Inc., New York, USA).

2.9. Data Handling, Phylogenetic and Statistical Analyses

The raw sequencing data were de-multiplexed and extracted in fastq format. Sequence data handling and taxonomic assignment were carried out using Geneious v11.1 (Biomatters Ltd., Auckland, New Zealand), MALT V0.3.8 and MEGAN v6.11 [28]. Trimming for read quality and length and adaptor sequence removal were performed using Trimmomatic v0.35 [29,30]. Trimmed reads were mapped to the NCBI-NT RefSeq 16S database via MALT V0.3.8, with hits down to 95% identity. For the operational taxonomic unit (OTU) identification and taxonomic binning, LCA-assignment algorithm (with 95% minimum identity) and 16s percent identity filter (species assignment at 99% identity) were employed. Relative bacterial abundance on the genus and species levels were calculated using the reads numbers of the corresponding OTUs.
Various alpha and beta diversity metrics were calculated for bacterial diversity and composition analyses. For this purpose, raw data were imported into QIIME2 [31], filtered and controlled for quality and chimeric sequences using DADA2, q2-demux, and dblur scripts [32,33]. The trimmed reads were subsequently mapped to the GreenGenes [34] and SILVA [35] databases for OTU identification and taxonomic binning. Faith phylogenetic diversity (PD), Pielou, Shannon, and Jaccard indices and Bray–Curtis and UniFrac distances were computed and evaluated using Kruskal–Wallis, Spearman, or permutational multivariate analysis of variance (PERMANOVA) tests as appropriate.
Data were analyzed using the Statistical Package for the Social Sciences (SPSS), version 22.0. (SPSS Inc., Chicago, IL, USA). Data normality was tested by Shapiro–Wilk test prior to further analyses. Kruskal–Wallis and Mann–Whitney U tests were employed for comparisons among groups where appropriate. Spearman’s rank-order correlation was used to analyze the correlation analysis between microbial taxa and biochemical and blood pressure measurements. A value of p < 0.05 was considered as significant.

3. Results

Baseline characteristics of the participants were summarized in Table 1. There were no differences in terms of age, dietary intake, anthropometrical measurements, and biochemical parameters except serum insulin levels between groups.
The changes in dietary intake, anthropometrical measurements, biochemical parameters, and blood pressure during the intervention period were given in Table 2. Intakes of energy and macronutrients did not change significantly during the intervention period in either kefir or unfermented milk groups (p > 0.05, for each). In terms of anthropometrical measurements, body weight and fat mass showed slight reductions after the 12-weeks intervention of kefir compared to unfermented milk, however, the changes in any of anthropometrical measurement from baseline to after the intervention did not differ significantly between the groups (p > 0.05, for each). Among the biochemical biomarkers, almost all parameters of lipid profile and glycaemic status showed amelioration by the intervention kefir group, however, only the difference in fasting insulin and thereby HOMA-IR from baseline to after intervention was significant (p = 0.050). Furthermore, TNF-α and IFN-γ showed a significant decrease after the intervention of kefir (p = 0.015 and p = 0.013, respectively), whereas IL-6 showed a slightly larger decrease in the unfermented milk group (p = 0.047). Both systolic blood pressure and diastolic blood pressure decreased significantly after the intervention in the kefir group (respectively p = 0.041 and p = 0.019), while only systolic blood pressure showed a modest decrease in the unfermented milk group (p = 0.047).
In regard to the analysis of gut microbiota composition, the mean number of total reads per sample was 90627 (Standart Deviation (SD): 44912; range: 31198–183068) at baseline and 118025 (SD:38831; range 39248–171915) after the intervention in kefir group, while it was recorded as 138775 (SD: 29961, range: 100922–195456) at the baseline and 95058 (SD: 23740, range: 55853–124889) after the intervention in the unfermented milk group. The gut microbiome of the participants was composed of phyla Bacteroidetes (51%), Firmicutes (30%), Proteobacteria (11%), Verrucomicrobia (0.02%) and Actinobacteria (0.003%) at the baseline in kefir group; the relative abundance of these phyla were detected as 39%, 39%, 10%, 0.03%, and 0.04% respectively after the intervention of kefir. Only the increase in the relative abundance of Actinobacteria was found to be statistically significant (p = 0.023). In the unfermented milk group, the phyla Bacteroidetes (66%), Firmicutes (27%), Verrucomicrobia (0.01%), Actinobacteria (0.01%), and Proteobacteria (0.03%) were detected at the baseline, however, the relative abundance of these phyla were changed to 33%, 56%, 0.03%, 0.04%, and 0.02% respectively after the intervention (Figure 2). When the changes in the relative abundance of each phyla distribution from after the intervention to baseline were compared between the groups, no significant difference was obtained (p > 0.05, for each). In the kefir group, the median of Firmicutes/Bacteroidetes ratio was 0.62 (range:0.06–10.01) at the baseline, and 1.77 (range:0.14–46.16) after the intervention (p = 0.388). This ratio was 0.30 (range:0.03–25.27) and 2.22 (range:0.34–12.04) respectively at baseline and end of the intervention in the unfermented milk group (p = 0.333). No significant difference was obtained between the changes in Firmicutes to Bacteroidetes ratio of groups (p > 0.05).
The changes in the relative abundance of Bacteriodetes and Firmucutes by the consumption of test drinks were given in Figure 3. The phyla Bacteroidetes was composed of five dominant genera; Bacteroides, Odoribacteraceae, Porphyromonadaceae, Prevotellaceae, and Alistipes in the gut microbiome of participants. Bacteroides (54% in kefir group and 44% in unfermented milk group), Prevotellacea (26% in kefir group and 40% in unfermented milk group) and Alistipes (11% in kefir group and 7% in unfermented milk group) were respectively most abundant families among Bacteroidetes at baseline. The relative abundance of Bacteroides changed to 59% in kefir group and 50% in unfermented milk group; Prevotellacea changed to 25% in kefir group and 29% in unfermented milk group; Alistipes changed to 9% in kefir group and 10% in unfermented milk group after the intervention. Despite the modest changes in the relative abundance of some genera, no significant difference was obtained when the changes in the relative abundances were compared between groups (p > 0.05, for each). Among Firmicutes, Clostridia, Erysipelotrichaceae, Veillonellaceae, and Lactobacillales were obtained in the fecal samples of participants. Although an increase in the relative abundance of Clostridia (from 73% to 85%) and Lactobacillales (2% to 5%), and also a decrease in the relative abundance of Veillonellaceae (from 9% to 6%) were obtained from baseline to after the intervention in kefir group, none of these changes were statistically significant (p > 0.05, for each). In the unfermented milk group, the relative abundance of Clostridia was increased from 75% to 89%, whereas Lactobacillales (5% to 2%) and Veillonellaceae (from 5% to 4%) were decreased from baseline to after the intervention. Similar to the phyla Bacteriodetes, no significant difference was obtained in the changes of the relative abundance of Firmicutes at genus level between the kefir and unfermented milk groups (p > 0.05, for each) (Figure 3). Among the phyla Actinobacteria, the relative abundance of Bifidobacterium was increased from 31% to 39% by the intervention of kefir, and from 23% to 32% by the intervention of unfermented milk. However, these changes were not found significant when compared between the groups (p > 0.05). Furthermore, Bifidobacterium species were detectable in only 50% of participants’ the fecal samples at the baseline, whereas they could be detected in 91.7% after the intervention of kefir (data not shown). On the contrary, Verrucomicrobia was obtained less frequently from baseline (75% of participants) to after the intervention (58.3% of participants) in the kefir group.
The correlation between the change in gut microbiota composition and the change in physiological characteristics, including anthropometrical measurements, biochemical markers, or blood pressure were conducted to examine the potential associations. Correlations between the changes in anthropometrical measurements and fecal microbiota composition at phylum and subphylum level were summarized in Table 3. The body weight and BMI were positively correlated with the relative abundance of Firmicutes and Proteobacteria. However, they were negatively correlated with the relative abundance of Clostridia (p < 0.05, for each). The body fat mass was negatively correlated with the relative abundance of Bacteroidetes (p < 0.01), and positively correlated with the relative abundance of Porphyromonadaceae, Firmicutes, and Actinobacteria (p < 0.05, for each). The waist circumference was negatively correlated with the relative abundance of Clostridia and positively correlated with the relative abundance Veillonellaceae (p < 0.05, for each).
Table 4 summarizes the correlation between the change in gut microbiota composition and biochemical markers. The changes in the relative abundance of Bacteroides was negatively correlated with both total and LDL cholesterol while the changes in the relative abundance of Veillonellaceae was negatively correlated with only LDL cholesterol (p < 0.05, for each). The change in the relative abundance of Bacteroidetes was negatively, and the change in the relative abundance of Odoribacteraceae and Alistipes groups was positively correlated with the change in serum glucose. The change in the relative abundance of Verrucomicrobia was positively correlated with changes in serum homocysteine and insulin (p < 0.05, for each). When we analyzed the correlation between the change in gut microbiota composition and blood pressure in the phylum level, the change in the relative abundance of Actinobacteria and Proteobacteria was positive, whereas Bacteroidetes was negatively correlated with the change in blood pressure. In the sub–phylum level, the change in the relative abundance Lactobacillales was positively correlated with the change in systolic and diastolic blood pressure (p < 0.05).
We further assessed several alpha and beta diversity metrics to assess bacterial biodiversity in specimens collected in week 0 and 12 from kefir and unfermented milk groups. No statistically significant differences were observed in OTU counts (Figure 4). Alpha diversity indices, indicating species richness and evenness with/without phylogenetic relations; namely, Shannon, Jaccard, and Faith PD indices, were similar between the study groups (Figure 5) (Jaccard plots not provided). Biodiversity between study cohorts, assessed by Bray–Curtis and weighted/unweighted UniFrac distances revealed no significant variation among study groups or in different time points. The PCoA plots of the unweighted UniFrac distances were given in Figure 5. No differences in OTU counts, alpha or beta diversity measures were observed when specimens from week 0 and 12 were assessed, regardless of the study group.
No side effect was reported by the participants during or after the intervention period that included the consumption of kefir.

4. Discussion

In this parallel–group randomized controlled study, regular kefir consumption during 12 weeks provided some improvements in anthropometrical measurements, lipid profile, glycaemic status, and inflammation in participants with MetS. In particular, insulin and HOMA–IR levels were significantly decreased, and also pro–inflammatory cytokines (TNF–α and IFN–γ) and blood pressure were ameliorated by kefir consumption. However, the magnitude of the improvements stayed insignificant when compared to unfermented milk. The effects of kefir on metabolic status were previously investigated in both animal models and human studies [37,38,39]. Some animal models suggested that kefir might have a potential to benefit the management of MetS by reducing body weight, fasting blood glucose, insulin, total, and LDL cholesterol, triacylglycerol, and pro–inflammatory cytokines, including IL–1β and IL–6 [38]. However, the evidence from human studies has been controversial. For instance, Ostadrahimi et al. reported that consumption of 600 mL/d kefir containing Lactobacillus casei, Lactobacillus acidophilus, and Bifidobacteria species had beneficial effects on fasting blood glucose and HbA1c compared to the control drink in patients with type 2 diabetes [37]. On the other hand, St–Onge et al. showed that 500 mL/day of kefir consumption for four weeks had no effect on lipid profile [40]. Furthermore, Fathi et al. showed that two servings of kefir in a day during eight weeks led a similar improvement both in lipid profile and weight management compared with milk [41,42]. The variation in response to the kefir consumption could be mainly explained by the variation of kefir composition, and the characteristics of study samples in different studies. Many different bacteria and yeast might be used for kefir production, and this might lead to distinct effects on metabolism and gut microbiota. Kefir used in this study contained Lactococcus lactis subsp. lactis, Lactococcus lactis subsp. cremoris, Lactococcus lactis subsp. diacetylactis, Leuconostoc mesenteroides subsp. cremoris, Lactobacillus kefir, Kluyveromyces marxianus, and Saccharomyces unisporus, and differed from the kefir samples used in other studies [37,38]. Furthermore, the initial metabolic profile of the participants was suggested an essential factor for the efficacy of probiotic interventions. Fuentes et al. showed that probiotics are more effective in patients with high baseline total cholesterol levels (251–300 mg/dL) compared to the patients with low baseline total cholesterol levels (200–250 mg/dL) [43]. Similarly, Nikbakht found out that probiotic supplementation was only effective in patients with baseline fasting blood glucose level above 126 mg/dL [44]. In our study, the median of total cholesterol levels was 243.50 mg/dL and 220.00 mg/dL, and the baseline glucose levels were 105.00 mg/dL and 101.50 mg/dL for kefir group and unfermented milk group, respectively. This may partly explain the lack of efficacy of kefir on metabolic status in our study.
Alterations in gut microbiota diversity, composition, and function were suggested to play a significant role in the development of MetS [12]. Ameliorating the intestinal dysbiosis with prebiotics and probiotics have gained considerable attention in recent years for the management of MetS [45]. However, the studies have yielded inconsistent results regarding the influence of probiotics on fecal microbial diversity and composition [46,47,48,49,50,51,52]. Furthermore, the effects of kefir as a probiotic on gut microbiota have been examined very limited and mainly with animal studies. Kim et al. revealed that three-week oral administration of kefir provided a decrease in the number of Firmicutes, Proteobacteria, and Enterobacteriaceae, and an increase in the number of Bacteroidetes, Lactobacillus, Lactococcus, and total yeast compared to milk group in mice [53]. However, in their follow–up study, no significant difference apart from the increase in Lactobacillus/Lactococcus populations was observed in the kefir group compared to control [39]. Similarly, an increase in Lactobacillus and Bifidobacterium populations and a reduction in Clostridium populations by consumption of kefir have been reported in mice previously [54,55,56]. The present study is one of the first reports showing the impact of kefir on human microbiota composition in patients with MetS. In this study, regular kefir and milk consumption for 12 weeks resulted in some alterations in the gut microbiota composition. For instance, Lactobacillus and Bifidobacterium spp. were increased by kefir consumption. However, apart from the increase in the relative abundance of Actinobacteria, no significant change by kefir consumption was recorded. Furthermore, the changes in the relative abundance of bacterial populations did not differ significantly between the groups. In some studies, following the probiotic supplementation, increases in the supplemented genera without an additional impact on the main microbial groups were observed [57,58]. Previously, Yılmaz et al. showed that 400 mL/day kefir consumption for four weeks in patients with inflammatory bowel diseases resulted in the significant increase of Lactobacillus bacterial load in feces [59]. Our study also showed an increase from 2% to 5% in the relative abundance of Lactobacillales by kefir consumption, albeit lacking statistical significance. It has suggested that the change of the microbiota composition may be related to several factors such as age, gender, initial microbiota composition, dietary intake, lifestyle factors, menopausal status, and medical therapy of the individuals [60,61]. Moreover, the microbiota composition of the product (kefir) that was tested [17,18] and the consumption pattern, including the period and frequency of consumption and amount of the product should be considered as the factors that have potential to influence the magnitude of the changes in gut microbiota [62]. Therefore, evaluating the effect of probiotics or fermented foods such as kefir on an individual basis may be set as a goal for future studies.
In this study, the correlations between changes in microbiota and anthropometric measurements or biochemical status were demonstrated. Our results pointed out a negative correlation between body fat mass and abundance of Bacteroidetes, whereas a positive correlation with the abundance of Firmicutes and Actinobacteria were observed. Although the data regarding the abundance of Bacteriodetes and Firmicutes phyla in obese and lean individuals is inconsistent, an overall analysis of results indicates an increase in Firmicutes with obesity [63]. Our results supported the previous reports, which revealed increased Firmicutes and decreased the abundance of Bacteroidetes are associated with obesity [64]. Turnbaugh et al. revealed a higher proportion of Actinobacteria in obese individuals compared to lean individuals [65]. In this study, we observed a positive correlation between body fat mass and abundance of Actinobacteria, which is in line with Turnbaugh et al.’s work. Members of the phylum of Proteobacteria are gram–negative bacteria and include several common human pathogens. An association between the increased relative abundance of Proteobacteria and increased risk of cardio–metabolic disorders was suggested previously [66]. In parallel with these findings, our study showed a positive correlation between an increased relative abundance of Proteobacteria and both body weight and diastolic blood pressure.
In terms of glycaemic status, Larsen et al. showed a lower abundance of Firmicutes and Clostridia, and a higher abundance of Bacteroidetes and Betaproteobacteria in diabetic patients compared to the non–diabetics [67]. Accordingly, we observed a lower abundance of Firmicutes and a higher abundance of Bacteroidetes in both groups at the beginning of the study. However, only the change in the relative abundance of Bacteroidetes by the dietary intervention was negatively correlated with the change in fasting blood glucose. This was parallel to results of the study conducted by Egshatyan et al., which found that microbiota of glucose–intolerant subjects were represented by Firmicutes phylum and to a lesser degree by Bacteroidetes phylum [68]. In the subphylum level, studies mostly indicate higher levels of Bacteroides and Prevotella and lower levels of butyrate producing–bacteria in type 2 diabetic patients [68,69]. In this study, we have observed a positive correlation between the relative abundance of Odoribacteraceae and Alistipes, and fasting plasma glucose. The correlation between Alistipes and blood glucose was also observed in a previous study [70].
Many researchers demonstrated a link between dysbiosis of gut microbiota and blood pressure. Yang et al. reported an increase in Firmicutes/Bacteroidetes ratio in hypertensive rats and humans. They also recorded a lower abundance of Actinobacteria as well as acetate– and butyrate–producing bacteria [71]. Yan et al. indicated higher levels of Proteobacteria but lower levels of Actinobacteria in hypertensive subjects [72]. In this study, the change in the relative abundance of Bacteroidetes was negatively correlated with the change in systolic and diastolic blood pressure as reported previously. Surprisingly, we observed a strong positive correlation with systolic and a weak positive correlation with diastolic blood pressure, respectively, with Actinobacteria abundance. The phylum Actinobacteria includes Bifidobacterium genera, possessing probiotic features [73]. Studies that report reductions of Actinobacteria in hypertensive patients explained this relationship mostly with Bifidobacterium levels [71,74]. In our study, when we analyzed the association between Bifidobacterium and blood pressure, no significant correlation was noted. This may be due to the lack of significant changes in Bifidobacterium abundance after the intervention. The species other than Bifidobacterium within the Actinobacteria phylum might be further investigated in terms of their contributions to hypertension.
Apart from fermentation, unfermented dairy products may also affect the gut microbiota [75,76]. In our study, regular milk consumption that was used as control also led to some changes in microbiota composition compared to the baseline. Firmucutes and Verrucomicrobiota were increased with milk consumption. However, Bacteroidetes group was decreased compared to the baseline. In accordance with our results, Ntemiri et al. also found out that whole milk consumption was associated with an increase in taxons belonging to Firmicutes and a higher Firmicutes/Bacteroidetes ratio [76]. In another randomized cross–over study, consumption of probiotic yogurt and milk acidified with D-(+)-glucono-δ-lactone showed some distinct effects on microbiota composition. In both groups, the abundance of Bilophila wadsworthia was reduced. However, only the abundance of Bifidobacterium species was increased with acidified milk intake, and it was suggested that gluconic acid in milk might possess prebiotic activity [75]. In a like manner, exopolysaccharides such as kefiran derived from kefir also suggested as bioactive compounds due to their potential prebiotic effects and relation to alteration of intestinal microbiota [77]. These results support a strong interaction between diet and microbiota even without probiotic intervention.
Studies showed that not only the composition of microbiota but also its functionality plays a role in the metabolic status [78,79]. This study focused only on the composition of gut microbiota; any consideration of its functionality was not taken into account. Using the metabolites of microbiota, such as postbiotics, as markers of the efficiency might have provided a better understanding. This should be noted as the main limitation of the study. The small number of participants in each arm could also be considered as the other limitation of the study.

5. Conclusions

In conclusion, to our knowledge, this was the first report exploring the effect of kefir on microbiota composition in patients with metabolic syndrome. This study indicated that kefir consumption could provide some potential improvements, especially in glycaemic status, inflammation–related indicators, and blood pressure, however, none of these improvements might stay significant when compared the changes led by unfermented milk consumption. Regarding to microbiota composition, the relative abundance of Actinobacteria phylum were increased in the kefir group compared to the baseline, even though a similar change by unfermented milk was also reported. Furthermore, this study underlined the potential alterations in gut microbiota composition that can be correlated with some indicators of the metabolic status led by both kefir and milk consumption, even if the magnitude of the efficiency remained limited. Further studies, especially randomized controlled trials, are needed to clarify the efficiency of kefir on gut microbiota and its link to metabolic status.

Author Contributions

Conceptualization, E.B.-K. and Z.B.; Data curation, E.B.-K. and Z.B.; Formal analysis, E.B.-K., K.E. and Z.B.; Funding acquisition, Y.A., E.Y. and Z.B.; Investigation, E.B.-K., B.P.S.-Y., C.K. and A.G.O.; Methodology, E.B.-K., B.P.S.-Y., Y.A., F.A.-K., C.K., A.B., A.N., K.E. and E.Y.; Writing—original draft, E.B.-K., K.E. and Z.B.; Writing—review & editing, C.K., K.E. and Z.B.

Funding

This research was funded by the Turkish Council of Higher Education.

Acknowledgments

We thank all study participants for their cooperation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alberti, K.G.; Eckel, R.H.; Grundy, S.M.; Zimmet, P.Z.; Cleeman, J.I.; Donato, K.A.; Fruchart, J.C.; James, W.P.; Loria, C.M.; Smith, S.C., Jr.; et al. Harmonizing the metabolic syndrome: a joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation 2009, 120, 1640–1645. [Google Scholar] [CrossRef] [PubMed]
  2. Mottillo, S.; Filion, K.B.; Genest, J.; Joseph, L.; Pilote, L.; Poirier, P.; Rinfret, S.; Schiffrin, E.L.; Eisenberg, M.J. The metabolic syndrome and cardiovascular risk a systematic review and meta–analysis. J. Am. Coll. Cardiol. 2010, 56, 1113–1132. [Google Scholar] [CrossRef] [PubMed]
  3. Elder, S.J.; Lichtenstein, A.H.; Pittas, A.G.; Roberts, S.B.; Fuss, P.J.; Greenberg, A.S.; McCrory, M.A.; Bouchard, T.J., Jr.; Saltzman, E.; Neale, M.C. Genetic and environmental influences on factors associated with cardiovascular disease and the metabolic syndrome. J. Lipid. Res. 2009, 50, 1917–1926. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Stancakova, A.; Laakso, M. Genetics of metabolic syndrome. Rev. Endocr. Metab. Disord. 2014, 15, 243–252. [Google Scholar] [CrossRef] [PubMed]
  5. Turnbaugh, P.J.; Ley, R.E.; Mahowald, M.A.; Magrini, V.; Mardis, E.R.; Gordon, J.I. An obesity–associated gut microbiome with increased capacity for energy harvest. Nature 2006, 444, 1027–1031. [Google Scholar] [CrossRef] [PubMed]
  6. Burcelin, R.; Luche, E.; Serino, M.; Amar, J. The gut microbiota ecology: a new opportunity for the treatment of metabolic diseases? Front Biosci 2009, 14, 5107–5117. [Google Scholar] [CrossRef]
  7. Cani, P.D.; Amar, J.; Iglesias, M.A.; Poggi, M.; Knauf, C.; Bastelica, D.; Neyrinck, A.M.; Fava, F.; Tuohy, K.M.; Chabo, C.; et al. Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes 2007, 56, 1761–1772. [Google Scholar] [CrossRef]
  8. Ley, R.E.; Backhed, F.; Turnbaugh, P.; Lozupone, C.A.; Knight, R.D.; Gordon, J.I. Obesity alters gut microbial ecology. Proc. Natl. Acad. Sci. USA 2005, 102, 11070–11075. [Google Scholar] [CrossRef] [Green Version]
  9. Backhed, F.; Ley, R.E.; Sonnenburg, J.L.; Peterson, D.A.; Gordon, J.I. Host–bacterial mutualism in the human intestine. Science 2005, 307, 1915–1920. [Google Scholar] [CrossRef]
  10. Shen, J.; Obin, M.S.; Zhao, L. The gut microbiota, obesity and insulin resistance. Mol. Asp. Med. 2013, 34, 39–58. [Google Scholar] [CrossRef]
  11. Esteve, E.; Ricart, W.; Fernandez-Real, J.M. Gut microbiota interactions with obesity, insulin resistance and type 2 diabetes: did gut microbiote co–evolve with insulin resistance? Curr. Opin. Clin. Nutr. Metab. Care 2011, 14, 483–490. [Google Scholar] [CrossRef]
  12. Festi, D.; Schiumerini, R.; Eusebi, L.H.; Marasco, G.; Taddia, M.; Colecchia, A. Gut microbiota and metabolic syndrome. World J. Gastroenterol. 2014, 20, 16079–16094. [Google Scholar] [CrossRef] [PubMed]
  13. Le Barz, M.; Anhe, F.F.; Varin, T.V.; Desjardins, Y.; Levy, E.; Roy, D.; Urdaci, M.C.; Marette, A. Probiotics as Complementary Treatment for Metabolic Disorders. Diabetes Metab. J. 2015, 39, 291–303. [Google Scholar] [CrossRef] [PubMed]
  14. Cani, P.D.; Van Hul, M. Novel opportunities for next–generation probiotics targeting metabolic syndrome. Curr. Opin. Biotechnol. 2015, 32, 21–27. [Google Scholar] [CrossRef] [PubMed]
  15. Ahmed, Z.; Wang, Y.; Ahmad, A.; Khan, S.T.; Nisa, M.; Ahmad, H.; Afreen, A. Kefir and health: A contemporary perspective. Crit. Rev. Food Sci. Nutr. 2013, 53, 422–434. [Google Scholar] [CrossRef] [PubMed]
  16. Koyu, E.B.; Buyuktuncer Demirel, Z. A functional food: Kefir. J. Nutr. Diet. 2018, 46, 166–175. [Google Scholar] [CrossRef]
  17. Rosa, D.D.; Dias, M.M.S.; Grzeskowiak, L.M.; Reis, S.A.; Conceicao, L.L.; Peluzio, M. Milk kefir: nutritional, microbiological and health benefits. Nutr. Res. Rev. 2017, 30, 82–96. [Google Scholar] [CrossRef]
  18. Bourrie, B.C.; Willing, B.P.; Cotter, P.D. The microbiota and health promoting characteristics of the fermented beverage kefir. Front. Microbiol. 2016, 7, 647. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, J.R.; Chen, M.J.; Lin, C.W. Antimutagenic and antioxidant properties of milk–kefir and soymilk–kefir. J. Agric. Food Chem. 2005, 53, 2467–2474. [Google Scholar] [CrossRef] [PubMed]
  20. Lee, M.Y.; Ahn, K.S.; Kwon, O.K.; Kim, M.J.; Kim, M.K.; Lee, I.Y.; Oh, S.R.; Lee, H.K. Anti–inflammatory and anti–allergic effects of kefir in a mouse asthma model. Immunobiology 2007, 212, 647–654. [Google Scholar] [CrossRef]
  21. Turan, I.; Dedeli, O.; Bor, S.; Ilter, T. Effects of a kefir supplement on symptoms, colonic transit, and bowel satisfaction score in patients with chronic constipation: a pilot study. Turk. J. Gastroenterol. 2014, 25, 650–656. [Google Scholar] [CrossRef] [PubMed]
  22. Friques, A.G.F.; Arpini, C.M.; Kalil, I.C.; Gava, A.L.; Leal, M.A.; Porto, M.L.; Nogueira, B.V.; Dias, A.T.; Andrade, T.U.; Pereira, T.M.C.; et al. Chronic administration of the probiotic kefir improves the endothelial function in spontaneously hypertensive rats. J. Transl. Med. 2015, 13, 390. [Google Scholar] [CrossRef] [PubMed]
  23. Hadisaputro, S.; Djokomoeljanto, R.R.; Judiono; Soesatyo, M.H. The effects of oral plain kefir supplementation on proinflammatory cytokine properties of the hyperglycemia Wistar rats induced by streptozotocin. Acta Med. Indones. 2012, 44, 100–104. [Google Scholar] [PubMed]
  24. Zimmet, P.; KG, M.M.A.; Serrano Rios, M. A new international diabetes federation worldwide definition of the metabolic syndrome: the rationale and the results. Rev. Esp. Cardiol. 2005, 58, 1371–1376. [Google Scholar] [CrossRef]
  25. Thomas, V.; Clark, J.; Dore, J. Fecal microbiota analysis: an overview of sample collection methods and sequencing strategies. Future Microbiol. 2015, 10, 1485–1504. [Google Scholar] [CrossRef]
  26. Wu, W.K.; Chen, C.C.; Panyod, S.; Chen, R.A.; Wu, M.S.; Sheen, L.Y.; Chang, S.C. Optimization of fecal sample processing for microbiome study—The journey from bathroom to bench. J. Formos. Med. Assoc. 2019, 118, 545–555. [Google Scholar] [CrossRef]
  27. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glockner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next–generation sequencing–based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef]
  28. Huson, D.H.; Beier, S.; Flade, I.; Górska, A.; El-Hadidi, M.; Mitra, S.; Ruscheweyh, H.J.; Tappu, R. MEGAN Community Edition– Interactive exploration and analysis of large–scale microbiome sequencing data. PLoS Comput. Biol. 2016, 12, e1004957. [Google Scholar] [CrossRef]
  29. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  30. Langmead, B.; Salzberg, S.L. Fast gapped–read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef]
  31. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. QIIME 2: Reproducible, interactive, scalable, and extensible microbiome data science. PeerJ Prepr. 2018, 6, e27295v27292. [Google Scholar] [CrossRef] [PubMed]
  32. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.; Holmes, S.P. DADA2: High–resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [PubMed]
  33. Amir, A.; McDonald, D.; Navas-Molina, J.A.; Kopylova, E.; Morton, J.T.; Zech Xu, Z.; Kightley, E.P.; Thompson, L.R.; Hyde, E.R.; Gonzalez, A.; et al. Deblur rapidly resolves single–nucleotide community sequence patterns. MSystems 2017, 2, e00191-16. [Google Scholar] [CrossRef] [PubMed]
  34. McDonald, D.; Price, M.N.; Goodrich, J.; Nawrocki, E.P.; DeSantis, T.Z.; Probst, A.; Andersen, G.L.; Knight, R.; Hugenholtz, P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, 610–618. [Google Scholar] [CrossRef] [PubMed]
  35. Yilmaz, P.; Parfrey, L.W.; Yarza, P.; Gerken, J.; Pruesse, E.; Quast, C.; Schweer, T.; Peplies, J.; Ludwig, W.; Glockner, F.O. The SILVA and “All–species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014, 42, D643–D648. [Google Scholar] [CrossRef] [PubMed]
  36. Vazquez-Baeza, Y.; Pirrung, M.; Gonzalez, A.; Knight, R. EMPeror: a tool for visualizing high–throughput microbial community data. GigaScience 2013, 2, 16. [Google Scholar] [CrossRef]
  37. Ostadrahimi, A.; Taghizadeh, A.; Mobasseri, M.; Farrin, N.; Payahoo, L.; Beyramalipoor Gheshlaghi, Z.; Vahedjabbari, M. Effect of probiotic fermented milk (kefir) on glycemic control and lipid profile in type 2 diabetic patients: a randomized double–blind placebo–controlled clinical trial. Iran. J. Public Health 2015, 44, 228–237. [Google Scholar]
  38. Rosa, D.D.; Grzeskowiak, L.M.; Ferreira, C.L.; Fonseca, A.C.; Reis, S.A.; Dias, M.M.; Siqueira, N.P.; Silva, L.L.; Neves, C.A.; Oliveira, L.L.; et al. Kefir reduces insulin resistance and inflammatory cytokine expression in an animal model of metabolic syndrome. Food Funct. 2016, 7, 3390–3401. [Google Scholar] [CrossRef]
  39. Kim, D.H.; Kim, H.; Jeong, D.; Kang, I.B.; Chon, J.W.; Kim, H.S.; Song, K.Y.; Seo, K.H. Kefir alleviates obesity and hepatic steatosis in high–fat diet–fed mice by modulation of gut microbiota and mycobiota: targeted and untargeted community analysis with correlation of biomarkers. J. Nutr. Biochem. 2017, 44, 35–43. [Google Scholar] [CrossRef]
  40. St-Onge, M.P.; Farnworth, E.R.; Savard, T.; Chabot, D.; Mafu, A.; Jones, P.J. Kefir consumption does not alter plasma lipid levels or cholesterol fractional synthesis rates relative to milk in hyperlipidemic men: A randomized controlled trial [ISRCTN10820810]. BMC Complement. Altern. Med. 2002, 2, 1. [Google Scholar] [CrossRef]
  41. Fathi, Y.; Ghodrati, N.; Zibaeenezhad, M.J.; Faghih, S. Kefir drink causes a significant yet similar improvement in serum lipid profile, compared with low–fat milk, in a dairy–rich diet in overweight or obese premenopausal women: A randomized controlled trial. J. Clin. Lipidol. 2017, 11, 136–146. [Google Scholar] [CrossRef] [PubMed]
  42. Fathi, Y.; Faghih, S.; Zibaeenezhad, M.J.; Tabatabaei, S.H. Kefir drink leads to a similar weight loss, compared with milk, in a dairy–rich non–energy–restricted diet in overweight or obese premenopausal women: A randomized controlled trial. Eur. J. Nutr. 2016, 55, 295–304. [Google Scholar] [CrossRef] [PubMed]
  43. Fuentes, M.C.; Lajo, T.; Carrion, J.M.; Cune, J. Cholesterol–lowering efficacy of Lactobacillus plantarum CECT 7527, 7528 and 7529 in hypercholesterolaemic adults. Br. J. Nutr. 2013, 109, 1866–1872. [Google Scholar] [CrossRef] [PubMed]
  44. Nikbakht, E.; Khalesi, S.; Singh, I.; Williams, L.T.; West, N.P.; Colson, N. Effect of probiotics and synbiotics on blood glucose: a systematic review and meta–analysis of controlled trials. Eur. J. Nutr. 2018, 57, 95–106. [Google Scholar] [CrossRef] [PubMed]
  45. Mallappa, R.H.; Rokana, N.; Duary, R.K.; Panwar, H.; Batish, V.K.; Grover, S. Management of metabolic syndrome through probiotic and prebiotic interventions. Indian J. Endocrinol. Metab. 2012, 16, 20–27. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, H.; Sun, J.; Liu, X.; Hong, C.; Zhu, Y.; Liu, A.; Li, S.; Guo, H.; Ren, F. Lactobacillus paracasei subsp. paracasei LC01 positively modulates intestinal microflora in healthy young adults. J. Microbiol. 2013, 51, 777–782. [Google Scholar] [CrossRef] [PubMed]
  47. Toscano, M.; De Grandi, R.; Stronati, L.; De Vecchi, E.; Drago, L. Effect of Lactobacillus rhamnosus HN001 and Bifidobacterium longum BB536 on the healthy gut microbiota composition at phyla and species level: A preliminary study. World J. Gastroenterol. 2017, 23, 2696–2704. [Google Scholar] [CrossRef] [PubMed]
  48. Larsen, N.; Vogensen, F.K.; Gobel, R.J.; Michaelsen, K.F.; Forssten, S.D.; Lahtinen, S.J.; Jakobsen, M. Effect of Lactobacillus salivarius Ls–33 on fecal microbiota in obese adolescents. Clin. Nutr. 2013, 32, 935–940. [Google Scholar] [CrossRef] [PubMed]
  49. Stadlbauer, V.; Leber, B.; Lemesch, S.; Trajanoski, S.; Bashir, M.; Horvath, A.; Tawdrous, M.; Stojakovic, T.; Fauler, G.; Fickert, P.; et al. Lactobacillus casei Shirota Supplementation Does Not Restore Gut Microbiota Composition and Gut Barrier in Metabolic Syndrome: A Randomized Pilot Study. PLoS ONE 2015, 10, e0141399. [Google Scholar] [CrossRef]
  50. Nova, E.; Pérez de Heredia, F.; Gómez-Martínez, S.; Marcos, A. The Role of Probiotics on the Microbiota: Effect on Obesity. Nutr. Clin. Pract. 2016, 31, 387–400. [Google Scholar] [CrossRef]
  51. Costabile, A.; Buttarazzi, I.; Kolida, S.; Quercia, S.; Baldini, J.; Swann, J.R.; Brigidi, P.; Gibson, G.R. An in vivo assessment of the cholesterol–lowering efficacy of Lactobacillus plantarum ECGC 13110402 in normal to mildly hypercholesterolaemic adults. PLoS ONE 2017, 12, e0187964. [Google Scholar] [CrossRef]
  52. Kristensen, N.B.; Bryrup, T.; Allin, K.H.; Nielsen, T.; Hansen, T.H.; Pedersen, O. Alterations in fecal microbiota composition by probiotic supplementation in healthy adults: a systematic review of randomized controlled trials. Genome Med. 2016, 8, 52. [Google Scholar] [CrossRef] [PubMed]
  53. Kim, D.-H.; Chon, J.-W.; Kim, H.; Seo, K.-H. Modulation of intestinal microbiota in mice by kefir administration. Food Sci. Biotechnol. 2015, 24, 1397–1403. [Google Scholar] [CrossRef]
  54. Liu, J.-R.; Wang, S.-Y.; Chen, M.-J.; Yueh, P.-Y.; Lin, C.-W. The anti–allergenic properties of milk kefir and soymilk kefir and their beneficial effects on the intestinal microflora. J. Sci. Food Agric. 2006, 86, 2527–2533. [Google Scholar] [CrossRef]
  55. Carasi, P.; Racedo, S.M.; Jacquot, C.; Romanin, D.E.; Serradell, M.A.; Urdaci, M.C. Impact of kefir derived Lactobacillus kefiri on the mucosal immune response and gut microbiota. J. Immunol. Res. 2015, 2015, 361604. [Google Scholar] [CrossRef] [PubMed]
  56. Hamet, M.F.; Medrano, M.; Perez, P.F.; Abraham, A.G. Oral administration of kefiran exerts a bifidogenic effect on BALB/c mice intestinal microbiota. Benef. Microbes 2016, 7, 237–246. [Google Scholar] [CrossRef] [PubMed]
  57. Lahtinen, S.J.; Forssten, S.; Aakko, J.; Granlund, L.; Rautonen, N.; Salminen, S.; Viitanen, M.; Ouwehand, A.C. Probiotic cheese containing Lactobacillus rhamnosus HN001 and Lactobacillus acidophilus NCFM(R) modifies subpopulations of fecal lactobacilli and Clostridium difficile in the elderly. Age 2012, 34, 133–143. [Google Scholar] [CrossRef] [PubMed]
  58. Bogovic Matijasic, B.; Obermajer, T.; Lipoglavsek, L.; Sernel, T.; Locatelli, I.; Kos, M.; Smid, A.; Rogelj, I. Effects of synbiotic fermented milk containing Lactobacillus acidophilus La–5 and Bifidobacterium animalis ssp. lactis BB–12 on the fecal microbiota of adults with irritable bowel syndrome: A randomized double–blind, placebo–controlled trial. J. Dairy Sci. 2016, 99, 5008–5021. [Google Scholar] [CrossRef]
  59. Yilmaz, I.; Dolar, M.E.; Ozpinar, H. Effect of administering kefir on the changes in fecal microbiota and symptoms of inflammatory bowel disease: A randomized controlled trial. Turk. J. Gastroenterol. 2019. [Google Scholar] [CrossRef]
  60. Uyeno, Y.; Sekiguchi, Y.; Kamagata, Y. Impact of consumption of probiotic lactobacilli–containing yogurt on microbial composition in human feces. Int. J. Food Microbiol. 2008, 122, 16–22. [Google Scholar] [CrossRef]
  61. Santos-Marcos, J.A.; Rangel-Zuniga, O.A.; Jimenez-Lucena, R.; Quintana-Navarro, G.M.; Garcia-Carpintero, S.; Malagon, M.M.; Landa, B.B.; Tena-Sempere, M.; Perez-Martinez, P.; Lopez-Miranda, J.; et al. Influence of gender and menopausal status on gut microbiota. Maturitas 2018, 116, 43–53. [Google Scholar] [CrossRef] [PubMed]
  62. Gerritsen, J.; Smidt, H.; Rijkers, G.T.; de Vos, W.M. Intestinal microbiota in human health and disease: The impact of probiotics. Genes Nutr. 2011, 6, 209–240. [Google Scholar] [CrossRef] [PubMed]
  63. Chakraborti, C.K. New–found link between microbiota and obesity. World J. Gastrointest. Pathophysiol. 2015, 6, 110–119. [Google Scholar] [CrossRef] [PubMed]
  64. Ley, R.E.; Turnbaugh, P.J.; Klein, S.; Gordon, J.I. Human gut microbes associated with obesity. Nature 2006, 444, 1022–1023. [Google Scholar] [CrossRef] [PubMed]
  65. Turnbaugh, P.J.; Hamady, M.; Yatsunenko, T.; Cantarel, B.L.; Duncan, A.; Ley, R.E.; Sogin, M.L.; Jones, W.J.; Roe, B.A.; Affourtit, J.P.; et al. A core gut microbiome in obese and lean twins. Nature 2008, 457, 480. [Google Scholar] [CrossRef]
  66. Rizzatti, G.; Lopetuso, L.; Gibiino, G.; Binda, C.; Gasbarrini, A.J.B.r.i. Proteobacteria: A common factor in human diseases. Biomed. Res. Int. 2017, 2017. [Google Scholar] [CrossRef] [PubMed]
  67. Larsen, N.; Vogensen, F.K.; van den Berg, F.W.; Nielsen, D.S.; Andreasen, A.S.; Pedersen, B.K.; Al-Soud, W.A.; Sorensen, S.J.; Hansen, L.H.; Jakobsen, M. Gut microbiota in human adults with type 2 diabetes differs from non–diabetic adults. PLoS ONE 2010, 5, e9085. [Google Scholar] [CrossRef]
  68. Egshatyan, L.; Kashtanova, D.; Popenko, A.; Tkacheva, O.; Tyakht, A.; Alexeev, D.; Karamnova, N.; Kostryukova, E.; Babenko, V.; Vakhitova, M.; et al. Gut microbiota and diet in patients with different glucose tolerance. Endocr. Connect. 2016, 5, 1–9. [Google Scholar] [CrossRef]
  69. Qin, J.; Li, Y.; Cai, Z.; Li, S.; Zhu, J.; Zhang, F.; Liang, S.; Zhang, W.; Guan, Y.; Shen, D.; et al. A metagenome–wide association study of gut microbiota in type 2 diabetes. Nature 2012, 490, 55. [Google Scholar] [CrossRef]
  70. Huang, G.; Xu, J.; Lefever, D.E.; Glenn, T.C.; Nagy, T.; Guo, T.L. Genistein prevention of hyperglycemia and improvement of glucose tolerance in adult non–obese diabetic mice are associated with alterations of gut microbiome and immune homeostasis. Toxicol. Appl. Pharmacol. 2017, 332, 138–148. [Google Scholar] [CrossRef]
  71. Yang, T.; Santisteban, M.M.; Rodriguez, V.; Li, E.; Ahmari, N.; Carvajal, J.M.; Zadeh, M.; Gong, M.; Qi, Y.; Zubcevic, J.; et al. Gut dysbiosis is linked to hypertension. Hypertension 2015, 65, 1331–1340. [Google Scholar] [CrossRef] [PubMed]
  72. Yan, Q.; Gu, Y.; Li, X.; Yang, W.; Jia, L.; Chen, C.; Han, X.; Huang, Y.; Zhao, L.; Li, P.; et al. Alterations of the Gut Microbiome in Hypertension. Front. Cell. Infect. Microbiol. 2017, 7, 381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Ventura, M.; Canchaya, C.; Tauch, A.; Chandra, G.; Fitzgerald, G.F.; Chater, K.F.; van Sinderen, D. Genomics of Actinobacteria: tracing the evolutionary history of an ancient phylum. Microbiol. Mol. Biol. Rev. 2007, 71, 495–548. [Google Scholar] [CrossRef] [PubMed]
  74. Santisteban, M.M.; Qi, Y.; Zubcevic, J.; Kim, S.; Yang, T.; Shenoy, V.; Cole-Jeffrey, C.T.; Lobaton, G.O.; Stewart, D.C.; Rubiano, A.; et al. Hypertension–Linked Pathophysiological Alterations in the Gut. Circulation Res. 2017, 120, 312–323. [Google Scholar] [CrossRef] [PubMed]
  75. Burton, K.J.; Rosikiewicz, M.; Pimentel, G.; Butikofer, U.; von Ah, U.; Voirol, M.J.; Croxatto, A.; Aeby, S.; Drai, J.; McTernan, P.G.; et al. Probiotic yogurt and acidified milk similarly reduce postprandial inflammation and both alter the gut microbiota of healthy, young men. Br. J. Nutr. 2017, 117, 1312–1322. [Google Scholar] [CrossRef] [Green Version]
  76. Ntemiri, A.; Ribiere, C.; Stanton, C.; Ross, R.P.; O’Connor, E.; O’Toole, P.W. Retention of microbiota diversity by lactose–free milk in a mouse model of elderly gut microbiota. J. Agric. Food Chem. 2019. [Google Scholar] [CrossRef] [PubMed]
  77. Lim, J.; Kale, M.; Kim, D.H.; Kim, H.S.; Chon, J.W.; Seo, K.H.; Lee, H.G.; Yokoyama, W.; Kim, H. Antiobesity Effect of Exopolysaccharides Isolated from Kefir Grains. J. Agric. Food Chem. 2017, 65, 10011–10019. [Google Scholar] [CrossRef]
  78. Hemarajata, P.; Versalovic, J. Effects of probiotics on gut microbiota: mechanisms of intestinal immunomodulation and neuromodulation. Therap. Adv. Gastroenterol. 2013, 6, 39–51. [Google Scholar] [CrossRef]
  79. Sanchez, B.; Delgado, S.; Blanco-Miguez, A.; Lourenco, A.; Gueimonde, M.; Margolles, A. Probiotics, gut microbiota, and their influence on host health and disease. Mol. Nutr. Food Res. 2017, 61. [Google Scholar] [CrossRef]
Figure 1. Timeline of the study.
Figure 1. Timeline of the study.
Nutrients 11 02089 g001
Figure 2. Gut microbiota composition before and after the intervention in each group.
Figure 2. Gut microbiota composition before and after the intervention in each group.
Nutrients 11 02089 g002
Figure 3. Bacterial changes in the relative abundance of Bacteriodetes (a) and Firmucutes (b) by the consumption of test drinks.
Figure 3. Bacterial changes in the relative abundance of Bacteriodetes (a) and Firmucutes (b) by the consumption of test drinks.
Nutrients 11 02089 g003
Figure 4. Box–and–whisker plots of the alpha diversity metrics calculated for the study groups. The box represents the median and the interquartile range, whereas whiskers indicate the 90th and 10th percentiles (A: Observed OTUs, B: Shannon index, C: Faith phylogenetic diversity index).
Figure 4. Box–and–whisker plots of the alpha diversity metrics calculated for the study groups. The box represents the median and the interquartile range, whereas whiskers indicate the 90th and 10th percentiles (A: Observed OTUs, B: Shannon index, C: Faith phylogenetic diversity index).
Nutrients 11 02089 g004
Figure 5. Principal coordinate analysis (PCoA) plot of the unweighted (A) and weighted (B) UniFrac distance matrices in the study groups. The plots were generated using EMPeror [36]. Axis titles indicate the percentage variations. The colors indicate sampling time (red: week 0, blue: week 12).
Figure 5. Principal coordinate analysis (PCoA) plot of the unweighted (A) and weighted (B) UniFrac distance matrices in the study groups. The plots were generated using EMPeror [36]. Axis titles indicate the percentage variations. The colors indicate sampling time (red: week 0, blue: week 12).
Nutrients 11 02089 g005
Table 1. Baseline characteristics of kefir and unfermented milk groups.
Table 1. Baseline characteristics of kefir and unfermented milk groups.
CharacteristicsKefir GroupUnfermented Milk Group
BaselineBaselinep
Sex (Female/Male)10/26/40.348
Age (year)52.00 (47.50–60.50)53.00 (45.00–60.00)0.821
Dietary intake
Energy (kcal/day)1694.16 (1590.92–1936.72)1655.35 (1423.52–2026.12)0.821
Carbohydrate (g)182.87 (166.79–205.74)155.79 (141.61–224.90)0.283
Protein (g)73.12 (59.29–83.91)65.23 (47.55–75.36)0.254
Fat (g)73.89 (68.16–97.67)85.21 (68.06–105.53)0.418
Fibre (g)26.11 (18.42–36.90)23.28 (17.11–26.23)0.418
Anthropometrical measurements
Weight (kg)84.05 (69.23–88.78)87.65 (75.60–100.60)0.180
Body mass index (kg/m2)30.67 (26.94–34.66)32.38 (29.18–34.59)0.381
Body fat mass (%)37.05 (31.33–44.05)37.45 (27.05–41.45)0.582
Waist circumference (cm)100.50 (90.75–110.00)106.75 (102.25–119.00)0.228
Hip circumference (cm)111.50 (106.00–116.50)112.00 (106.00–119.25)0.771
Waist–to–hip ratio0.92 (0.86–0.99)0.97 (0.92–1.00)0.203
Lipid profile
Total cholesterol (mg/dL)243.50 (217.25–265.25)220.00 (199.75–249.00)0.228
HDL cholesterol (mg/dL)45.00 (39.00–55.75)42.50 (34.50–56.25)0.456
LDL cholesterol (mg/dL)154.50 (135.75–177.00)141.00 (114.50–177.50)0.283
Triglycerides (mg/dL)185.00 (114.50–216.75)164.50 (126.25–220.75)0.923
Homocysteine (µmoL/L)10.01 (8.64–12.40)13.10 (10.73–15.25)0.050
Glycaemic status
Glucose (mg/dL)105.00 (93.75–109.75)101.50 (97.00–107.25)>0.99
Insulin (mU/L)15.94 (11.75–17.64)19.04 (18.09–25.49)0.011 *
HbA1c (%)5.60 (5.25–5.88)5.65 (5.20–6.03)0.872
HOMA–IR4.18 (2.86–4.59)4.52 (4.29–6.65)0.180
Inflammation–related indicators
hs–CRP (mg/dL)0.22 (0.69–0.80)0.27 (0.21–0.41)0.722
TNF–α (pg/mL)12.01 (0.76–43.05)8.51 (0.49–25.85)0.418
IL–6 (pg/mL)15.82 (11.52–29.75)19.73 (13.85–28.71)0.418
IL–10 (pg/mL)4.38 (1.13–32.90)1.45 (1.13–9.34)0.456
IFN–γ (IU/mL)1.23 (0.12–2.19)0.56 (0.02–3.04)>0.99
ALT (U/L)18.50 (16.50–24.00)25.00 (20.75–31.25)0.140
AST (U/L)19.00 (18.00–20.00)19.00 (18.00–20.25)0.923
GGT (U/L)15.00 (10.75–23.00)19.00 (16.00–40.50)0.169
Blood pressure
Systolic blood pressure (mmHg)134.50 (115.25–140.50)132.50 (123.75–144.00)0.722
Diastolic blood pressure (mmHg)85.00 (77.50–92.00)89.00 (81.00–92.00)0.497
Data are given as median (25th percentile–75th percentile). Mann–Whitney U test was used to compare differences between groups. Fisher’s exact test was used to compare gender between groups. HDL: High Density Lipoprotein, LDL: Low Density Lipoprotein, HbA1c: Glycosylated Hemoglobin, HOMA-IR: Homeostasis Model of Assessment Insulin Resistance, hs-CRP: High-sensitivity C-reactive Protein, TNF: Tumor Necrosis Factor, IL: Interleukin, IFN: Interferon, ALT: Alanine Aminotransferase, AST: Aspartate Aminotransferase, GGT: Gamma-glutamyl Transferase. * p < 0.05.
Table 2. Dietary intake, anthropometrical measurements, biochemical parameters, and blood pressure in kefir and unfermented milk groups.
Table 2. Dietary intake, anthropometrical measurements, biochemical parameters, and blood pressure in kefir and unfermented milk groups.
CharacteristicsKefir Group Unfermented Milk Group
BaselineWeek 12p1BaselineWeek 12p2p3
Dietary intake
Energy (kcal/day)1694.16 (1590.92–1936.72)1995.73 (1567.23–2351.80)0.3471655.35 (1423.52–2026.12)1979.09 (1606.15–2123.62)0.5750.821
Carbohydrate (%)44.00 (38.00–45.75)42.50 (37.75–48.75)0.96441.00 (32.00–45.25)46.00 (38.50–52.50)0.2140.283
Protein (%)16.50 (15.25–19.00)13.50 (12.00–18.50)0.06615.50 (12.75–17.25)15.00 (12.00–16.50)0.7170.254
Fat (%)39.50 (37.00–44.75)41.00 (38.00–46.00)0.66643.50 (41.75–48.00)40.00 (32.75–44.25)0.1670.228
Fibre (g)26.11 (18.42–36.90)26.81 (21.58–32.65)0.81423.28 (17.11–26.23)25.17 (16.67–33.61)0.6460.722
Anthropometrical measurements
Weight (kg)84.05 (69.23–88.78)83.50 (66.90–88.75)0.69587.65 (75.60–100.60)88.55 (74.33–96.65)0.2070.418
Body mass index (kg/m2)30.67 (26.94–34.66)30.58 (26.24–34.31)0.75432.38 (29.18–34.59)31.90 (29.05–33.71)0.2410.418
Body fat mass (%)37.05 (31.33–44.05)35.85 (30.58–44.23)0.24837.45 (27.05–41.45)38.30 (29.63–43.98)0.2410.069
Waist circumference (cm)100.50 (90.75–110.00)102.25 (90.00–109.00)0.407106.75 (102.25–119.00)106.75 (100.50–118.50)0.9520.722
Hip circumference (cm)111.50 (106.00–116.50)110.00 (106.25–118.63)0.813112.00 (106.00–119.25)111.75 (105.13–116.25)0.4830.228
Waist–to–hip ratio0.92 (0.86–0.99)0.92 (0.86–0.95)0.9290.97 (0.92–1.00)0.99 (0.92–1.03)0.3860.497
Lipid profile
Total cholesterol (mg/dL)243.50 (217.25–265.25)222.00 (201.25–275.00)0.209220.00 (199.75–249.00)226.50 (198.75–240.25)0.9530.539
HDL cholesterol (mg/dL)45.00 (39.00–55.75)46.00 (41.00–63.00)0.27142.50 (34.50–56.25)43.50 (36.00–58.00)0.4120.346
LDL cholesterol (mg/dL)154.50 (135.75–177.00)144.00 (116.50–188.75)0.170141.00 (114.50–177.50)147.50 (115.75–167.50)0.9590.314
Triglycerides (mg/dL)185.00 (114.50–216.75)152.50 (116.50–191.25)0.530164.50 (126.25–220.75)161.50 (117.00–236.75)0.8781.000
Homocysteine (µmol/L)10.01 (8.64–12.40)9.31 (7.45–12.70)0.18213.10 (10.73–15.25)12.00 (11.05–14.35)0.2130.710
Glycaemic status
Glucose (mg/dL)105.00 (93.75–109.75)100.50 (96.50–103.00)0.157101.50 (97.00–107.25)98.50 (97.50–116.25)0.9180.159
Insulin (mU/L)15.94 (11.75–17.64)13.64 (7.33–16.36)0.050 *19.04 (18.09–25.49)22.08 (15.05–28.54)0.3860.123
HbA1c (%)5.60 (5.25–5.88)5.65 (5.50–5.98)0.1575.65 (5.20–6.03)5.70 (5.10–5.90)0.9180.123
HOMA–IR4.18 (2.86–4.59)3.42 (1.93–4.22)0.050 *4.52 (4.29–6.65)5.52 (3.38–8.49)0.4450.159
Inflammation–related indicators
hs–CRP (mg/dL)0.22 (0.69–0.80)0.16 (0.10–0.46)0.5330.27 (0.21–0.41)0.24 (0.13–0.50)0.9170.733
TNF–α (pg/mL)12.01 (0.76–43.05)1.13 (0.49–8.33)0.015 *8.51 (0.49–25.85)4.12 (0.49–13.03)0.4010.123
IL–6 (pg/mL)15.82 (11.52–29.75)13.47 (5.65–21.39)0.09919.73 (13.85–28.71)10.03 (6.16–16.45)0.047 *0.872
IL–10 (pg/mL)4.38 (1.13–32.90)1.91 (1.13–14.77)0.3861.45 (1.13–9.34)1.13 (1.13–15.95)0.7350.539
IFN–γ (IU/mL)1.23 (0.12–2.19)0.38 (0.04–0.85)0.013 *0.56 (0.02–3.04)0.49 (0.18–1.19)0.0860.628
ALT (U/L)18.50 (16.50–24.00)22.00 (19.50–24.00)0.28825.00 (20.75–31.25)24.50 (18.75–29.00)0.2150.180
AST (U/L)19.00 (18.00–20.00)19.00 (17.00–22.50)0.88719.00 (18.00–20.25)17.50 (17.00–19.50)0.1360.203
GGT (U/L)15.00 (10.75–23.00)14.50 (12.00–23.75)0.37119.00 (16.00–40.50)19.00 (14.25–29.25)0.0650.169
Blood pressure
Systolic blood pressure (mmHg)134.50 (115.25–140.50)118.00 (103.25–137.75)0.041 *132.50 (123.75–144.00)118.00 (105.75–137.00)0.047 *0.974
Diastolic blood pressure (mmHg)85.00 (77.50–92.00)78.50 (69.00–80.00)0.019 *89.00 (81.00–92.00)78.50 (66.75–89.50)0.0591.000
Data are given as median (25th percentile–75th percentile). HDL: High Density Lipoprotein, LDL: Low Density Lipoprotein, HbA1c: Glycosylated Hemoglobin, HOMA-IR: Homeostasis Model of Assessment Insulin Resistance, hs-CRP: High-sensitivity C-reactive Protein, TNF: Tumor Necrosis Factor, IL: Interleukin, IFN-γ: Interferon-γ, ALT: Alanine Aminotransferase, AST: Aspartate Aminotransferase, GGT: Gamma-glutamyl Transferase. p1 gives the differences between baseline and after intervention in the kefir group. p2 gives the differences between baseline and after intervention in the unfermented group. p3 gives the comparison of the changes from baseline to week 12th between groups. p1 and p2 were analysed by the Wilcoxon test while p3 was analyzed by the Mann–Whitney U test. * p ≤ 0.05.
Table 3. Correlations between changes in anthropometrical measurements and fecal microbiota composition at phylum and subphylum level .
Table 3. Correlations between changes in anthropometrical measurements and fecal microbiota composition at phylum and subphylum level .
BacteroidetesBacteroidesOdoribacteraceaePorphyromonadaceaePrevotellaceaeAlistipesFirmicutesClostridiaErysipelotrichaceaeVeillonellaceaeLactobacillalesVerrucomicrobiaActinobacteriaBifidobacteriumProteobacteria
Anthropometric measurements
Weight (kg)−0.3840.0700.1760.160−0.1310.1330.433 *−0.432 *−0.3380.2950.3520.1320.4000.0300.456 *
Body mass index (kg/m2)−0.3830.1030.1890.172−0.1680.1390.434 *−0.455 *−0.3410.2910.3470.0950.3820.0470.461 *
Body fat mass (%)−0.563 **0.2230.2850.468 *−0.3880.4020.599 **0.017−0.0700.2100.0380.4200.536 *−0.0180.371
Waist circumference (cm)−0.1510.267−0.321−0.009−0.217−0.0470.242−0.505 *−0.2870.432 *0.1350.2470.1280.0440.332
Hip circumference (cm)0.0820.0290.3260.164−0.149−0.0780.070−0.096−0.2640.2720.215−0.1150.127−0.3670.081
Waist–to–hip ratio−0.2390.112−0.526 *−0.086−0.0380.0140.168−0.2510.0460.066−0.0430.3760.0900.2230.208
Values indicate Spearman correlation coefficients (n = 22). * Correlation is significant at the 0.05 level (2–tailed). ** Correlation is significant at the 0.01 level (2–tailed).
Table 4. Correlations between changes in biochemical markers and fecal microbiota composition at phylum and sub–phylum level .
Table 4. Correlations between changes in biochemical markers and fecal microbiota composition at phylum and sub–phylum level .
BacteroidetesBacteroidesOdoribacteraceaePorphyromonadaceaePrevotellaceaeAlistipesFirmicutesClostridiaErysipelotrichaceaeVeillonellaceaeLactobacillalesVerrucomicrobiaActinobacteriaBifidobacteriumProteobacteria
Lipid profile
Total cholesterol (mg/dL)−0.151−0.444 *−0.025−0.0260.334−0.007−0.1010.3890.205−0.4080.0090.0590.190−0.164−0.096
HDL cholesterol (mg/dL)−0.003−0.4180.130−0.0660.414−0.004−0.248−0.2240.1610.1540.2590.175−0.016−0.0120.130
LDL cholesterol (mg/dL)−0.016−0.535 *−0.021−0.0230.3050.075−0.2980.3220.330−0.431 *0.0050.1520.042−0.034−0.094
Triglycerides (mg/dL)−0.080−0.007−0.118−0.0790.197−0.2630.2510.334−0.211−0.242−0.082−0.1760.199−0.292−0.210
Homocysteine (µmoL/L)−0.1620.077−0.1860.141−0.135−0.023−0.057−0.0440.2320.2240.1490.477 *0.1100.1270.359
Glycaemic status
Glucose (mg/dL)−0.590 **0.0920.423 *0.370−0.3740.629 **0.3870.3170.056−0.037−0.0330.3280.3650.0900.187
Insulin (mU/L)−0.3310.0690.3440.244−0.2650.3610.357−0.098−0.1590.3120.0100.428 *0.284−0.2180.123
HbA1c (%)0.1890.062−0.168−0.2660.004−0.1230.0210.126−0.3720.211−0.282−0.003−0.2800.068−0.277
HOMA–IR−0.3560.0300.4060.207−0.2550.458 *0.346−0.119−0.1480.2690.0670.3950.294−0.2090.123
Inflammation Related Indicators
hs–CRP (mg/dL)−0.1450.110−0.2940.238−0.255−0.0320.0960.1450.3580.078−0.0150.3800.321−0.2890.078
TNF–α (pg/mL)−0.2970.087−0.305−0.034−0.2240.0510.0100.3170.0930.023−0.0620.3920.281−0.3230.322
IL–6 (pg/mL)−0.0160.3810.0010.013−0.287−0.1420.313−0.320−0.3330.399−0.1700.039−0.0430.0910.274
IL–10 (pg/mL)0.086−0.130−0.414−0.3140.165−0.216−0.054−0.289−0.0400.0150.2730.151−0.0060.008−0.054
IFN–γ (IU/mL)0.066−0.241−0.089−0.0740.1110.291−0.182−0.0240.074−0.1620.0840.093−0.0290.086−0.076
AST−0.3430.2920.0160.058−0.159−0.1780.283−0.160−0.031−0.0060.203−0.3280.190−0.1870.169
ALT−0.2510.235−0.0890.067−0.081−0.423 *0.3910.138−0.1200.2530.0980.1580.252−0.4020.175
GGT0.255−0.2640.112−0.0990.493 *−0.385−0.0930.025−0.285−0.089−0.093−0.457 *−0.016−0.169−0.208
Blood Pressure
Systolic blood pressure (mmHg)−0.531 *−0.148−0.0800.1370.0920.0760.243−0.1720.144−0.1980.536 *0.2050.710 **−0.1680.379
Diastolic blood pressure (mmHg)−0.491 *0.1690.1830.271−0.3480.1710.244−0.2050.239−0.1280.561 **0.0960.452 *0.0130.469 *
Values indicate Spearman correlation coefficients (n = 22). * Correlation is significant at the 0.05 level (2–tailed). ** Correlation is significant at the 0.01 level (2–tailed). HDL: High Density Lipoprotein, LDL: Low Density Lipoprotein, HbA1c: Glycosylated Hemoglobin, HOMA-IR: Homeostasis Model of Assessment Insulin Resistance, hs-CRP: High-sensitivity C-reactive Protein, TNF: Tumor Necrosis Factor, IL: Interleukin, IFN: Interferon, ALT: Alanine Aminotransferase, AST: Aspartate Aminotransferase, GGT: Gamma-glutamyl Transferase.

Share and Cite

MDPI and ACS Style

BELLIKCI-KOYU, E.; SARER-YUREKLI, B.P.; AKYON, Y.; AYDIN-KOSE, F.; KARAGOZLU, C.; OZGEN, A.G.; BRINKMANN, A.; NITSCHE, A.; ERGUNAY, K.; YILMAZ, E.; et al. Effects of Regular Kefir Consumption on Gut Microbiota in Patients with Metabolic Syndrome: A Parallel-Group, Randomized, Controlled Study. Nutrients 2019, 11, 2089. https://doi.org/10.3390/nu11092089

AMA Style

BELLIKCI-KOYU E, SARER-YUREKLI BP, AKYON Y, AYDIN-KOSE F, KARAGOZLU C, OZGEN AG, BRINKMANN A, NITSCHE A, ERGUNAY K, YILMAZ E, et al. Effects of Regular Kefir Consumption on Gut Microbiota in Patients with Metabolic Syndrome: A Parallel-Group, Randomized, Controlled Study. Nutrients. 2019; 11(9):2089. https://doi.org/10.3390/nu11092089

Chicago/Turabian Style

BELLIKCI-KOYU, Ezgi, Banu Pınar SARER-YUREKLI, Yakut AKYON, Fadime AYDIN-KOSE, Cem KARAGOZLU, Ahmet Gokhan OZGEN, Annika BRINKMANN, Andreas NITSCHE, Koray ERGUNAY, Engin YILMAZ, and et al. 2019. "Effects of Regular Kefir Consumption on Gut Microbiota in Patients with Metabolic Syndrome: A Parallel-Group, Randomized, Controlled Study" Nutrients 11, no. 9: 2089. https://doi.org/10.3390/nu11092089

APA Style

BELLIKCI-KOYU, E., SARER-YUREKLI, B. P., AKYON, Y., AYDIN-KOSE, F., KARAGOZLU, C., OZGEN, A. G., BRINKMANN, A., NITSCHE, A., ERGUNAY, K., YILMAZ, E., & BUYUKTUNCER, Z. (2019). Effects of Regular Kefir Consumption on Gut Microbiota in Patients with Metabolic Syndrome: A Parallel-Group, Randomized, Controlled Study. Nutrients, 11(9), 2089. https://doi.org/10.3390/nu11092089

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

Article Metrics

Back to TopTop