Introduction: Obesity, characterized by excess adipose tissue, is a major public health problem worldwide. Brown adipose tissue (BAT) and beige adipose tissue participate in thermogenesis through uncoupling protein 1 (UCP1). Polyphenols including those from Calafate (a native polyphenol-rich Patagonian berry), are considered as potential anti-obesity compounds due to their pro-thermogenic characteristics. However, polyphenols are mainly metabolized by the gut microbiota (GM) that may influence their bioactivity and bioavailability. The aim of this study was to determine the impact of dietary administration with a Calafate polyphenol-rich extract on thermogenic activity of BAT and beige adipose tissue and GM composition. Methods: Eight-week-old C57BL6 mice (n = 30) were divided into 4 groups to receive for 24 weeks a control diet (C), a high-fat diet alone (HF), or high-fat diet supplemented with Calafate extract (HFC) or the same high-fat diet supplemented with Calafate extract but treated with antibiotics (HFCAB) from week 19–20. Administration with Calafate extract (50 mg/kg per day) was carried out for 3 weeks from week 21–23 in the HFC and HFCAB groups. After euthanasia, gene expression of thermogenic markers was analyzed in BAT and inguinal white adipose tissue (iWAT). Transmission electron microscopy was performed to assess mitochondrial morphology and cristae density in BAT. GM diversity and composition were characterized by deep sequencing with the MiSeq Illumina platform. Results: Calafate extract administration had no effect on weight gain in mice fed a high-fat diet. However, it prevented alterations in mitochondrial cristae induced by HFD and increased Dio2 expression in BAT and iWAT. The intervention also influenced the GM composition, preventing changes in specific bacterial taxa induced by the high-fat diet. However, the antibiotic treatment prevented in part these effects, suggesting the implications of GM. Conclusion: These results suggest that the acute administration of a Calafate extract modulates the expression of thermogenic markers, prevents alterations in mitochondrial cristae and intestinal microbiota in preclinical models. The study highlights the complex interaction between polyphenols, thermogenesis, and the GM, providing valuable insights into their potential roles in the treatment of obesity-related metabolic diseases.

Obesity, an abnormal increase in white adipose tissue (WAT) due to an energy imbalance, is a major health problem worldwide [1]. Alongside WAT, brown adipose tissue (BAT) contributes to the regulation of body temperature through thermogenesis [2]. The mitochondria of BAT express uncoupling protein 1 (UCP1), which dissipates energy in the form of heat [3]. Factors such as thyroid hormones, catecholamines, exposure to cold, and dietary phytochemicals can enhance the thermogenic activity of BAT [4]. It has also been reported that some stimuli may also promote the appearance of clusters of brown adipocytes in certain regions of the WAT [5].

Recently, polyphenols have been considered as possible anti-obesity compounds due to their thermogenic characteristics [6]. These secondary plant metabolites are poorly absorbed in the small intestine and pass into the colon where they are deglycosylated by the gut microbiota (GM) and their aglycone moiety subsequently metabolized by different bacterial taxa into a large number of metabolites that can be absorbed [7]. Although some polyphenol glycosides can be hydrolyzed by enzymes present in enterocytes, such as the lactase-phlorizin hydrolase, most of them are hydrolyzed by the GM [8].

The GM forms a complex bacterial ecosystem of almost 100 billion microorganisms. A greater richness of species is observed in the GM of healthy individuals, and the loss of this diversity is known as dysbiosis [9]. GM diversity is influenced by several factors, of which diet is one of the most important [10]. Recent data show that the metabolism of polyphenols by GM influences their bioavailability and bioactivity [11].

Calafate (Berberis microphylla) is a native Patagonian berry rich in anthocyanins, hydroxamic acid derivatives, and flavonols, which has been shown to exert a wide range of health benefits including anti-obesogenic and pro-thermogenic effects in animal models [12, 13]. The aim of this study was to determine the impact of acute treatment with a Calafate polyphenol-rich extract on the thermogenic activity of brown and beige adipose tissue and on the GM composition of mice fed a high-fat diet. The role of GM in this effect was also determined by using a group of animals treated with antibiotics.

Calafate Extract

Calafate extract derived from freeze-dried Calafate (INIA, Chile) was obtained using a low-toxicity solvent (ethanol:water). The crude extract was filtered and purified by HPLC using an Amberlite XAD7HP resin column. Characterization revealed a total polyphenol concentration of 4,257 mg GAE/100 g, total anthocyanins of 4,028 mg D-3-G/100 g, and a total antioxidant capacity of 59,009 mmol TE/100 g.

Animal Housing

Eight-week-old C57BL/6J male mice (n = 30; approximate weight 20 g) from the colony of the Nutrition Department of the University of Chile, were housed in cages of 3–4 animals, at controlled relative humidity, room temperature (21–23°C), with 12-h light-dark cycles, and were randomly assigned into 4 groups: (a) control group (C) fed a standard diet (D12328 diet, 11 kcal% fat, 16% protein, 73% carbohydrate with corn starch Surwit diet) (n = 5); (b) obese group (HF) fed a high-fat diet (D12330 diet, 58 kcal% fat, 16% protein, 26% carbohydrate with corn starch Surwit diet) (n = 10), (c) obese group fed the high-fat diet (D12330 diet) that received 50 mg/kg body weight per day of total polyphenols of Calafate extract in water for 3 weeks from week 21–23 (n = 10) (HFC), (d) obese group fed the high-fat diet (D12330 diet) that received a cocktail of broad-spectrum antibiotics [14] (ampicillin [0.5 mg/mL], neomycin [0.5 mg/mL], gentamicin [0.5 mg/mL], metronidazole [0.5 mg/mL], vancomycin [0.25 mg/mL]) in sterile drinking water ad libitum for 2 weeks (week 19–20) prior to Calafate extract administration, and then treated daily with the Calafate extract for 3 weeks (week 21–23) (n = 5) (HFCAB). The control and the HF diets were administered for 4 months. The dose of Calafate extract administered to the animals in the HFC and HFCAB groups was previously tested in our laboratory and the results published [15]. Antibiotic administration in the HFCAB group was performed to test whether GM was necessary to obtain the thermogenic effects of the Calafate extract. At the end of trial, animals were euthanized, their body weights recorded, and samples of inguinal (iWAT), epididymal (eWAT), and interscapular (BAT) adipose tissue collected, as well as the cecal contents, and stored at −80°C until subsequent analysis

Thermogenesis and Browning Markers

The AllPrep® DNA/RNA/Protein Mini Kit (Qiagen, Hilden, Germany) was used to obtain RNA, DNA, and proteins from samples of adipose tissue. RNA was treated with DNase (DNA-free, Amboina). 500 ng of total RNA was used for reverse transcription using a commercial kit (High-Capacity cDNA Reverse Transcription Kit, ThermoFisher Scientific, Waltham, MA, USA). The resulting cDNA was used for qPCR using TaqMan® probes (Applied Biosystems) in a Stratagene Mx3000P System (Agilent Technologies, CA, USA). Specific transcripts of thermogenesis-related markers (Ucp1, Dio2, Pgc1-α, Pparα, Pparγ, Prdm16, and Sirt1) were evaluated in inguinal WAT and interscapular BAT. The expression levels of genes were normalized using Gapdh and Ppia housekeeping genes. The reported gene expression changes between groups were calculated by the 2 (∆∆Ct) method [16].

Transmission Electron Microscopy

BAT was fixed in glutaraldehyde (2.5%) and dissected into bundles of fibers. Samples were then washed with water and stained with 1% uranyl acetate for 2 h. Finally, 80-nm sections were cut, mounted, and examined on a Tecnai G2 T12 transmission electron microscope (Philips-FEI). The mitochondrial size was calculated as the mean area of each organelle [17, 18]. Mitochondria cristae parameters were evaluated manually by trained personnel, as previously described [17, 18]. Morphometric analyses were formed by using ImageJ software [17, 18].

Evaluation of GM

Cecal samples were collected for GM analysis. Bacterial DNA was extracted using a commercial kit (QIAamp DNA Stool Mini Kit, Qiagen, Hilden, Germany), adding a bead-beater stage to optimize bacterial breakdown [19]. DNA concentrations (ng/μL) were measured on a NanoQuant Infinite 200Pro kit (Tecan, Männedorf, Switzerland), and DNA integrity was assessed by electrophoresis on 1% agarose gel. Samples were sent to the University of Illinois (Roy J. Carver Biotechnology Center - Urbana-Champaign, Champaign, IL, USA) for PCR amplification of the 16 s rDNA gene (v3–v4 region) and amplicon sequencing using the MiSeq Illumina platform (Illumina, San Diego, CA, USA). Results were expressed as relative abundances of each bacterial taxa (phyla, families, and genera) and the intra- and interindividual diversity indices were also determined.

Statistical Analysis

Statistical analysis was performed using SPSS Packaged 20.0 software (SPSS Inc., Chicago, IL, USA) and GraphPad Prism 6.0 software (GraphPad Software, Inc., San Diego, CA, USA) for graphing results. Normal distribution of the data was determined by using the Shapiro-Wilk test. Data were described as mean and standard deviations (for variables with parametric distribution) or as median and interquartile range (for these with nonparametric distribution). Differences between groups were determined by one-way ANOVA test followed by Tukey’s post hoc test in case of parametric data, and Kruskal-Wallis test followed by Dunn’s post hoc in case of nonparametric data. Relative bacterial abundances were expressed as proportions, i.e., percentages of the total sequences. Additionally, principal coordinate analysis with weighted UniFrac [20] was performed to analyze the clustering of samples based on their composition. A value of p < 0.05 was considered statistically significant.

The general characteristics of the animals in the four groups are shown in Table 1. Weight gain throughout the experiment was higher in the HF and HFC groups than in the HFCAB and C groups. The eWAT weight was significantly higher in the groups fed the high-fat diet compared to the C group. BAT weight was significantly higher in the HFCAB group than in the C group (Table 1).

Table 1.

General characteristics of adult male C57BL/6J mice

Groups/treatmentsOne-way ANOVA (p value)
CHFHFCHFCAB
Initial weight, g 24.5±1.6 25.7±1.9 23.3±2.8 24.6±1.5 0.541 
Final weight, g 28.2±1.8 (a) 36.9±1.4 (b) 32.8±1.4 (b) 34.9±3.4 (b) 0.003 
Weight gain, g 2.5±0.9 (a) 6.5±2.9 (b) 6.9±2.4 (b) 4.9±1.9 (a) 0.001 
Weight eWAT, g 0.5±0.1 (a) 1.9±0.2 (b) 1.2±0.3 (b) 1.6±0.5 (b) 0.003 
Weight iWAT, g 0.03±0.01 (a) 0.07±0.02 (a, b) 0.07±0.02 (b) 0.07±0.02 (a, b) 0.018 
Weight BAT, g 0.07±0.09 (a) 0.11±0.01 (a, b) 0.08±0.02 (a, b) 0.12±0.02 (b) 0.010 
Groups/treatmentsOne-way ANOVA (p value)
CHFHFCHFCAB
Initial weight, g 24.5±1.6 25.7±1.9 23.3±2.8 24.6±1.5 0.541 
Final weight, g 28.2±1.8 (a) 36.9±1.4 (b) 32.8±1.4 (b) 34.9±3.4 (b) 0.003 
Weight gain, g 2.5±0.9 (a) 6.5±2.9 (b) 6.9±2.4 (b) 4.9±1.9 (a) 0.001 
Weight eWAT, g 0.5±0.1 (a) 1.9±0.2 (b) 1.2±0.3 (b) 1.6±0.5 (b) 0.003 
Weight iWAT, g 0.03±0.01 (a) 0.07±0.02 (a, b) 0.07±0.02 (b) 0.07±0.02 (a, b) 0.018 
Weight BAT, g 0.07±0.09 (a) 0.11±0.01 (a, b) 0.08±0.02 (a, b) 0.12±0.02 (b) 0.010 

Values represent mean ± SD. Data analyzed with one-way ANOVA. Post hoc Tukey’s test, different letters mean differences of at least p < 0.05.

C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics.

As shown in Figures 1a, b, no differences were found in the mitochondrial size between groups. However, the relative distribution of mitochondrial size (Fig. 1e) showed a greater proportion of larger mitochondria (+1.3 µm2) in the HFC group than in the C, HF, and HFCAB groups. Compared to C group, HF and HFCAB groups had a lower number of mitochondrial cristae (Fig. 1c) and total length of the cristae (Fig. 1d). Animals from the HFC group displayed similar mitochondrial cristae number and total cristae length than in the C group (Fig. 1c, d).

Fig. 1.

Effect of Calafate on mitochondrial size and mitochondrial cristae density in mice with high-fat diet. a Representative image of brown adipose tissue (BAT). b Mitochondrial size from BAT. c Mitochondrial cristae number from BAT. d Mitochondrial cristae length from BAT. e Histogram plot of relative distributions of mitochondrial size from BAT. n = 3 mice per group. Data were analyzed using the Kruskal-Wallis test, followed by Dunn’s posttest. Differences from the control group are indicated by * (p < 0.05). C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics.

Fig. 1.

Effect of Calafate on mitochondrial size and mitochondrial cristae density in mice with high-fat diet. a Representative image of brown adipose tissue (BAT). b Mitochondrial size from BAT. c Mitochondrial cristae number from BAT. d Mitochondrial cristae length from BAT. e Histogram plot of relative distributions of mitochondrial size from BAT. n = 3 mice per group. Data were analyzed using the Kruskal-Wallis test, followed by Dunn’s posttest. Differences from the control group are indicated by * (p < 0.05). C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics.

Close modal

Changes in gene expression of the biomarkers of thermogenesis in BAT and iWAT are described in Table 2. In both tissues, a higher expression of Dio2 transcript was observed in the HFC group compared to the C and HFCAB groups, while no differences were observed with Ucp1, Pgc1-α, Pparϒ, Pparα, Prdm16, and Sirt1. Therefore, antibiotic administration prevented the Calafate-induced increase of Dio2 in BAT and iWAT.

Table 2.

Genetic markers of thermogenesis in brown and inguinal adipose tissue of adult male C57BL/6J mice

Groups/treatmentsOne-way ANOVA (p value)
CHFHFCHFCAB
BAT gene expression 
 Ucp1/Ppia (AU) 1.000±0.21 1.218±0.44 1.835±1.14 1.659±0.73 0.2141 
 Pgc1-α/Ppia (AU) 1.000±0.19 2.005±0.44 2.036±0.85 1.627±1.12 0.1341 
 Pparγ/Ppia (AU) 1.000±0.30 0.997±0.20 0.940±0.20 0.916±0.22 0.8971 
 Pparα/Ppia (AU) 1.000±0.15 1.100±0.09 0.966±0.22 1.066±0.18 0.4239 
 Prdm16/Ppia (AU) 1.000±0.61 1.375±0.69 1.584±0.87 1.696±0.75 0.4667 
 Sirt1/Ppia (AU) 1.000±0.44 0.7291±0.13 0.928±0.13 0.877±0.13 0.1640 
 Dio2/Ppia (AU) 1.000(a)±0.47 2.311(a,b)±1.23 3.314(b)±1.05 0.870(a)±0.62 0.0000 
iWAT gene expression 
 Ucp1/Gapdh (AU) 1.000±0.20 0.851±0.55 0.650±0.32 0.538±0.44 0.3489 
 Pgc1-α/Gapdh (AU) 1.000±0.19 2.005±0.44 2.036±0.85 1.627±1.27 0.1341 
 Pparγ/Gapdh (AU) 1.000±0.46 0.868±0.31 0.671±0.29 0.428±0.21 0.0514 
 Pparα/Gapdh (AU) 1.000±0.32 0.786±0.38 0.595±0.24 0.589±0.36 0.1973 
 Prdm16/Gapdh (AU) 1.000±0.18 1.174±0.53 0.870±0.43 0.764±0.48 0.4045 
 Sirt1/Gapdh (AU) 1.000±0.29 1.258±0.38 1.504±0.59 1.114±0.55 0.3343 
 Dio2/Gapdh (AU) 1.000(a)±0.30 2.837(a)±1.74 6.465(b)±5.09 1.548(a)±0.66 0.0247 
Groups/treatmentsOne-way ANOVA (p value)
CHFHFCHFCAB
BAT gene expression 
 Ucp1/Ppia (AU) 1.000±0.21 1.218±0.44 1.835±1.14 1.659±0.73 0.2141 
 Pgc1-α/Ppia (AU) 1.000±0.19 2.005±0.44 2.036±0.85 1.627±1.12 0.1341 
 Pparγ/Ppia (AU) 1.000±0.30 0.997±0.20 0.940±0.20 0.916±0.22 0.8971 
 Pparα/Ppia (AU) 1.000±0.15 1.100±0.09 0.966±0.22 1.066±0.18 0.4239 
 Prdm16/Ppia (AU) 1.000±0.61 1.375±0.69 1.584±0.87 1.696±0.75 0.4667 
 Sirt1/Ppia (AU) 1.000±0.44 0.7291±0.13 0.928±0.13 0.877±0.13 0.1640 
 Dio2/Ppia (AU) 1.000(a)±0.47 2.311(a,b)±1.23 3.314(b)±1.05 0.870(a)±0.62 0.0000 
iWAT gene expression 
 Ucp1/Gapdh (AU) 1.000±0.20 0.851±0.55 0.650±0.32 0.538±0.44 0.3489 
 Pgc1-α/Gapdh (AU) 1.000±0.19 2.005±0.44 2.036±0.85 1.627±1.27 0.1341 
 Pparγ/Gapdh (AU) 1.000±0.46 0.868±0.31 0.671±0.29 0.428±0.21 0.0514 
 Pparα/Gapdh (AU) 1.000±0.32 0.786±0.38 0.595±0.24 0.589±0.36 0.1973 
 Prdm16/Gapdh (AU) 1.000±0.18 1.174±0.53 0.870±0.43 0.764±0.48 0.4045 
 Sirt1/Gapdh (AU) 1.000±0.29 1.258±0.38 1.504±0.59 1.114±0.55 0.3343 
 Dio2/Gapdh (AU) 1.000(a)±0.30 2.837(a)±1.74 6.465(b)±5.09 1.548(a)±0.66 0.0247 

Values represent mean ± SD Data analyzed with one-way ANOVA *p < 0.05. Post hoc Tukey’s test, different letters mean differences of at least p < 0.05.

C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics; BAT, brown adipose tissue; UA, iWAT, inguinal white adipose tissue; AU, arbitrary units.

Regarding the GM, Shannon index as a marker of alpha diversity is shown in Figure 2. The HFCAB group displayed a lower alpha diversity than the other groups, confirming the impact of antibiotic administration on GM composition. The principal component analysis (Fig. 3) showed a significant difference in β-diversity. Both the HF and the HFC groups clustered together and differed from the C and HFCAB groups, showing that the antibiotic treatment and high-fat diet administration affect the GM composition.

Fig. 2.

Alpha diversity using Shannon index of the GM of adult male C57BL/6J mice. Differences in alpha diversity between groups were analyzed by One-way ANOVA. Post hoc Tukey’s test, different letters mean differences of at least p < 0.05. C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics.

Fig. 2.

Alpha diversity using Shannon index of the GM of adult male C57BL/6J mice. Differences in alpha diversity between groups were analyzed by One-way ANOVA. Post hoc Tukey’s test, different letters mean differences of at least p < 0.05. C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics.

Close modal
Fig. 3.

Comparison of the intestinal microbiota of adult male C57BL/6J mice. Principal coordinate analysis of beta diversity values based on weighted UniFrac distances. Each color represents a population from a specific geographic area. Ellipses were drawn using a 95% confidence interval for each group.

Fig. 3.

Comparison of the intestinal microbiota of adult male C57BL/6J mice. Principal coordinate analysis of beta diversity values based on weighted UniFrac distances. Each color represents a population from a specific geographic area. Ellipses were drawn using a 95% confidence interval for each group.

Close modal

Changes in the composition of bacterial phyla, families, and genera are described in Table 3. Results show that the high-fat diet, alone or with Calafate extract, had no effect on the microbiota at the phylum level. Animals from the HF group had lower abundances of the Burkholderiaceae, Bacteroidales F082 and Coriobacteriales families, and a higher abundance of Lachnospiraceae, compared to C group. At the genus level, animals fed the high-fat diet had a lower abundance of Parasutterella, Muribaculum,PrevotellaceaeUCG-003, and Rikenellaceae RC9 gut group, and a higher abundance of Tyzzerella 3, Desulfovibrio, and Lachnospiraceae NK4A136 group, versus the C group. Administration with Calafate extract prevented (fully or in part) the increase of Tyzzerella 3, Desulfovibrio, and Lachnospiraceae NK4A136 group, and the decrease of Prevotellaceae UCG-003 and Rikenellaceae RC9 gut group, induced by the high-fat diet. Finally, HFC group also showed a higher abundance of Ruminococcaceae UCG-010 genus versus the HF group.

Table 3.

GM composition of adult male C57BL/6J mice by 16S rRNA V3–V4 gene sequencing using the MiSeq Illumina system

Relative abundance of bacteria in groups/treatmentsOne-way ANOVA (p value)
CHFHFCHFCAB
Bacterial phyla
Proteobacteria 0.141±0.05 (a) 0.157±0.08 (a) 0.116±0.05 (a) 0.001±0.01 (b) 0.003 
Deferribacteres 0.042±0.04 (a) 0.020±0.01 (a) 0.021±0.01 (a) 0.000±0.01 (b) 0.005 
Spirochetes 0.001±0.002 (a, b) 0.005±0.006 (a) 0.003±0.002 (a) 0.000±0.000 (b) 0.007 
Patescibacteria 0.010±0.01 (a, b) 0.031±0.02 (a) 0.023±0.02 (a) 0.000±0.001 (b) 0.005 
Verrumicrobia 0.065±0.06 (a) 0.010±0.02 (a) 0.013±0.02 (a) 0.410±0.08 (b) 0.001 
Bacterial families 
Burkholderiaceae 0.042±0.02 (a) 0.001±0.002 (b) 0.003±0.003 (b) 0.000±0.00 (b) 0.000 
Bacteroidales F082 0.014±0.01 (a) 0.002±0.002 (b) 0.004±0.005 (b) 0.000±0.00 (b) 0.003 
Coriobacteriales 0.007±0.01 (a) 0.000±0.00 (b) 0.000±0.00 (b) 0.000±0.00 (a, b) 0.030 
Lachnospiraceae 0.143±0.09 (a) 0.410±0.13 (b) 0.438±0.12 (b) 0.288±0.06 (a, b) 0.003 
Rikenellaceae 0.017 ± 0.012 (a) 0.011 ± 0.005 (a) 0.009 ± 0.004 (a) 0.003 ± 0.004 (b) 0.020 
Deferribacteraceae 0.042 ± 0.04 (a) 0.020 ± 0.01 (a, b) 0.021 ± 0.01 (a, b) 0.000 ± 0.001 (b) 0.037 
Peptococcaceae 0.000 ± 0.000 (a) 0.000 ± 0.00 (a) 0.000 ± 0.001 (a) 0.004 ± 0.005 (b) 0.015 
Akkermansiaceae 0.065 ± 0.06 (a) 0.010 ± 0.082 (a) 0.013 ± 0.02 (a) 0.410 ± 0.08 (b) 0.000 
Prevotellaceae 0.020 ± 0.02 (a) 0.003 ± 0.005 (a, b) 0.010 ± 0.009 (a) 0.000 ± 0.00 (b) 0.039 
Bacterial genera 
Parasutterella 0.042±0.02 (a) 0.001±0.002 (b) 0.002±0.002 (b) 0.000±0.00 (b) 0.000 
Muribaculum 0.002±0.002 (a) 0.000±0.00 (b) 0.000±0.00 (b) 0.001±0.01 (a, b) 0.001 
Mucispirillum 0.042±0.04 (a) 0.020±0.01 (a, b) 0.021±0.01 (a, b) 0.000±0.001 (b) 0.037 
Alistipes 0.010±0.008 (a) 0.009±0.005 (a) 0.006±0.002 (a, b) 0.000±0.001 (b) 0.003 
Ruminococcaceae UCG−004 0.000 ± 0.00 (a, b) 0.000 ± 0.00 (a) 0.003 ± 0.001 (a, b) 0.001 ± 0.002 (b) 0.018 
Erysipelatoclostridium 0.000 ± 0.00 (a) 0.000 ± 0.00 (a) 0.001 ± 0.003 (a) 0.005 ± 0.003 (b) 0.002 
Akkermansia 0.065 ± 0.06 (a) 0.010 ± 0.02 (a) 0.0013 ± 0.03 (a) 0.410 ± 0.08 (b) 0.000 
Tyzzerella3 0.000 ± 0.00 (a) 0.001 ± 0.001 (b) 0.000 ± 0.00 (a) 0.000 ± 0.00 (a) 0.002 
Prevotellaceae UCG−003 0.017± 0.02 (a) 0.001 ± 0.002 (b) 0.005 ± 0.006 (a, b) 0.000 ± 0.00 (a, b) 0.040 
Desulfovibrio 0.026 ± 0.01 (a, c) 0.066 ± 0.02 (b) 0.044 ± 0.01 (a) 0.000 ± 0.00 (c) 0.000 
Ruminococcaceae UCG−010 0.000 ± 0.00 (a, b) 0.000 ± 0.00 (b) 0.003 ± 0.00 (a) 0.000 ± 0.00 (a, b) 0.018 
Lachnospiraceae NK4A136 group 0.004 ± 0.004 (a) 0.202 ± 0.133 (b) 0.121 ± 0.07 (a, b) 0.179 ± 0.14 (b) 0.011 
Rikenellaceae RC9 gut group 0.004 ± 0.003 (a) 0.001 ± 0.001 (b) 0.002 ± 0.002 (a, b) 0.000 ± 0.000 (b) 0.013 
Bacterial species 
Mucispirillum schaedleri 0.037±0.04 (a) 0.018±0.01 (a, b) 0.017±0.01 (a, b) 0.000±0.001 (b) 0.048 
Akkermansia muciniphila 0.064±0.062 (a) 0.010±0.02 (a) 0.013±0.02 (a) 0.346±0.07 (b) 0.000 
Relative abundance of bacteria in groups/treatmentsOne-way ANOVA (p value)
CHFHFCHFCAB
Bacterial phyla
Proteobacteria 0.141±0.05 (a) 0.157±0.08 (a) 0.116±0.05 (a) 0.001±0.01 (b) 0.003 
Deferribacteres 0.042±0.04 (a) 0.020±0.01 (a) 0.021±0.01 (a) 0.000±0.01 (b) 0.005 
Spirochetes 0.001±0.002 (a, b) 0.005±0.006 (a) 0.003±0.002 (a) 0.000±0.000 (b) 0.007 
Patescibacteria 0.010±0.01 (a, b) 0.031±0.02 (a) 0.023±0.02 (a) 0.000±0.001 (b) 0.005 
Verrumicrobia 0.065±0.06 (a) 0.010±0.02 (a) 0.013±0.02 (a) 0.410±0.08 (b) 0.001 
Bacterial families 
Burkholderiaceae 0.042±0.02 (a) 0.001±0.002 (b) 0.003±0.003 (b) 0.000±0.00 (b) 0.000 
Bacteroidales F082 0.014±0.01 (a) 0.002±0.002 (b) 0.004±0.005 (b) 0.000±0.00 (b) 0.003 
Coriobacteriales 0.007±0.01 (a) 0.000±0.00 (b) 0.000±0.00 (b) 0.000±0.00 (a, b) 0.030 
Lachnospiraceae 0.143±0.09 (a) 0.410±0.13 (b) 0.438±0.12 (b) 0.288±0.06 (a, b) 0.003 
Rikenellaceae 0.017 ± 0.012 (a) 0.011 ± 0.005 (a) 0.009 ± 0.004 (a) 0.003 ± 0.004 (b) 0.020 
Deferribacteraceae 0.042 ± 0.04 (a) 0.020 ± 0.01 (a, b) 0.021 ± 0.01 (a, b) 0.000 ± 0.001 (b) 0.037 
Peptococcaceae 0.000 ± 0.000 (a) 0.000 ± 0.00 (a) 0.000 ± 0.001 (a) 0.004 ± 0.005 (b) 0.015 
Akkermansiaceae 0.065 ± 0.06 (a) 0.010 ± 0.082 (a) 0.013 ± 0.02 (a) 0.410 ± 0.08 (b) 0.000 
Prevotellaceae 0.020 ± 0.02 (a) 0.003 ± 0.005 (a, b) 0.010 ± 0.009 (a) 0.000 ± 0.00 (b) 0.039 
Bacterial genera 
Parasutterella 0.042±0.02 (a) 0.001±0.002 (b) 0.002±0.002 (b) 0.000±0.00 (b) 0.000 
Muribaculum 0.002±0.002 (a) 0.000±0.00 (b) 0.000±0.00 (b) 0.001±0.01 (a, b) 0.001 
Mucispirillum 0.042±0.04 (a) 0.020±0.01 (a, b) 0.021±0.01 (a, b) 0.000±0.001 (b) 0.037 
Alistipes 0.010±0.008 (a) 0.009±0.005 (a) 0.006±0.002 (a, b) 0.000±0.001 (b) 0.003 
Ruminococcaceae UCG−004 0.000 ± 0.00 (a, b) 0.000 ± 0.00 (a) 0.003 ± 0.001 (a, b) 0.001 ± 0.002 (b) 0.018 
Erysipelatoclostridium 0.000 ± 0.00 (a) 0.000 ± 0.00 (a) 0.001 ± 0.003 (a) 0.005 ± 0.003 (b) 0.002 
Akkermansia 0.065 ± 0.06 (a) 0.010 ± 0.02 (a) 0.0013 ± 0.03 (a) 0.410 ± 0.08 (b) 0.000 
Tyzzerella3 0.000 ± 0.00 (a) 0.001 ± 0.001 (b) 0.000 ± 0.00 (a) 0.000 ± 0.00 (a) 0.002 
Prevotellaceae UCG−003 0.017± 0.02 (a) 0.001 ± 0.002 (b) 0.005 ± 0.006 (a, b) 0.000 ± 0.00 (a, b) 0.040 
Desulfovibrio 0.026 ± 0.01 (a, c) 0.066 ± 0.02 (b) 0.044 ± 0.01 (a) 0.000 ± 0.00 (c) 0.000 
Ruminococcaceae UCG−010 0.000 ± 0.00 (a, b) 0.000 ± 0.00 (b) 0.003 ± 0.00 (a) 0.000 ± 0.00 (a, b) 0.018 
Lachnospiraceae NK4A136 group 0.004 ± 0.004 (a) 0.202 ± 0.133 (b) 0.121 ± 0.07 (a, b) 0.179 ± 0.14 (b) 0.011 
Rikenellaceae RC9 gut group 0.004 ± 0.003 (a) 0.001 ± 0.001 (b) 0.002 ± 0.002 (a, b) 0.000 ± 0.000 (b) 0.013 
Bacterial species 
Mucispirillum schaedleri 0.037±0.04 (a) 0.018±0.01 (a, b) 0.017±0.01 (a, b) 0.000±0.001 (b) 0.048 
Akkermansia muciniphila 0.064±0.062 (a) 0.010±0.02 (a) 0.013±0.02 (a) 0.346±0.07 (b) 0.000 

Values represent mean ± SD. Data analyzed with one-way ANOVA. Post hoc Tukey’s test, different letters mean differences of at least p < 0.05.

C, control; HF, high-fat diet; HFC, high-fat diet + Calafate; HFCAB, high-fat diet + Calafate + Antibiotics.

The administration of antibiotics affects the relative abundance of various bacterial groups in the HFCAB group. This group exhibited a lower abundance of the Proteobacteria and Deferribacteres phyla, and a higher abundance of Verrumicrobia, compared to C group. At the family level, it also had a lower abundance of Burkhohlderiaceae, Bacteroidales F082, Rikenellaceae, Deferribactetaceae, and Prevotellaceae while that of Peptococcaceae and Akkermansiaceae was higher, compared to C group. Moreover, in the HFCAB group, the genera Erysipelatoclostridium, Akkermansia, Lachnospiraceae NK4A136 group, and the species Akkermansia muciniphila were more abundant, while Parasutterella, Mucispirillum, Alistipes, and Rikenellaceae RC9 gut group were less abundant, compared to C group.

In this study, we were able to demonstrate that acute administration of Calafate increased the relative expression of Dio2, a thermogenesis marker, and prevented alterations in mitochondrial cristae and GM composition induced by the consumption of a high-fat diet. Contrarily to previous studies [12, 13] which observed a decreased weight gain after 16 weeks of administration of the same dose of Calafate extract, the current study showed no effect of the administration of Calafate extract for 3 weeks to mice fed a high-fat diet on this parameter. Probably, this lower impact can be attributed to the shorter duration of Calafate extract administration in the present study.

Interestingly, one of the most notable findings of our study is that mitochondria from the HF group presented lower mitochondrial cristae density compared to control. Mitochondrial cristae is a dynamic energetic compartment involved in oxidative phosphorylation and ATP production, whose shape changes according to physiological conditions, as shown in mononuclear cells [17] and skeletal muscle cells [21]. Decrease in mitochondrial cristae density has been associated with oxidative stress, an early risk factor for metabolic alterations in the obesity [22]. Thus, administration of Calafate may have prevented mitochondrial cristae alterations in BAT of obese mice, probably due to the antioxidant properties of its polyphenols.

Type II iodothyronine deiodinase is a protein encoded by the DIO2 gene that is involved in BAT activation and thermogenesis [23‒25]. We observed increased expression of the Dio2 transcript in the HFC group compared to the C group in both BAT and iWAT. This increase was not observed in the HFCAB group, although they also received Calafate. This could be caused by the administration of antibiotics and the consequent decrease in GM alpha diversity in the HFCAB group. Previous studies indicate that GM could stimulate Type II iodothyronine deiodinase activity and expression in BAT and beige adipose tissue through the stimulation of the TGR5 receptor by secondary bile acids formed from primary bile acids by specific bacterial taxa, resulting in increased thermogenesis and browning [26].

When analyzing GM composition, a lower abundance of Burkholderiaceae, Bacteroidales F082, and Parasutterella was observed in all HFD groups, confirming results from previous studies in mice fed a similar diet [27, 28]. Additionally, we observed that HFD also decreased the abundance of Coriobacteriales and Muribaculum and increased that of Lachnospiraceae, compared to C group, but these observations were not observed in the HFCAB group. These results confirmed those from previous studies [29].

In our study, Calafate administration prevented the HFD-induced increase in Tyzzerella 3, Lachnospiraceae NK4A136 group y Desulfovibrio genera. Previous studies have described that higher Tyzzerella abundances correlated with body and liver weight, perirenal fat, mesenteric fat, triglycerides, blood glucose, HOMA-IR, TNF-α, and IL-1β [30]. HFD has also been associated with an increase in bacteria from the Lachnospiraceae NK4A136 genus, this change being prevented by the 10 week-administration of a rosemary extract [31]. In contrast to our results, the relative abundance of bacteria from the Desulfovibrio genus has been shown to increase with the administration of polyphenols [26]. However, a diminution of Desulfovibrio might be considered as beneficial for the host’s health, since they are sulfate-reducing bacteria that are involved in production of H2S, a potentially harmful metabolite [32].

Our research shows that the administration of Calafate in the HFC group prevented the decrease in the relative abundance of the Prevotellaceae family and the genera Prevotellaceae UCG-003, Rikenellaceae RC9 gut group induced by HFD. Similarly, the administration of procyanidins (100 mg/kg body weight) for 12 weeks to HFD-fed mice has been shown to increase the abundance of Rikenellaceae RC9 [33]. Furthermore, Calafate extract administration increased the abundance of the Ruminococcaceae UCG-010 genus, compared to the HF group, therefore confirming previous studies in HFD-fed animals [34].

Despite the fact that the HFCAB group was treated with antibiotics, it presented a higher abundance of Akkermansiaceae and Peptococcaceae families, and A. muciniphila species, compared to the other groups. Previous investigations with animals treated with broad-spectrum antibiotics have reported similar expansion of A. muciniphila [35]. Increased abundance of Akkermansiaceae and Peptococcaceae has also been observed in germ-free animals after fecal transplant [36, 37], and it has been proposed that it could be due to an imbalance in the immune system of these animals [38].

Our study has some limitations. First, it is not possible to determine the mechanisms that best explain the results. Furthermore, the antibiotic combinations used in the HFCAB group cannot completely deplete GM, and we cannot reject the idea that these antibiotics may have an effect on thermogenesis. Finally, the animals were housed at room temperature throughout the experiment. Recent research has indicated that housing mice at room temperature (below their thermoneutrality) can lead to a significant increase in metabolic rate since mice must defend their core temperature by activating both shivering and non-shivering thermogenesis [39]. Despite these limitations, our results confirm the acute administration of Calafate extract induce the expression of thermogenic markers, prevents alterations in mitochondrial cristae and modulates GM in mice fed a high-fat diet.

This study was approved by the Institutional Committee for the Care and Use of Animals (CICUA) of the University of Chile, Santiago, Chile (protocol 17049-MED-UCH, accepted on June 27, 2017). All procedures were performed following the institutional policies and guidelines of the CICUA of the Universidad de Chile.

The authors have no conflicts of interest to declare.

The National Agency for Research and Development (ANID, Chile) funded this work (Grant FONDECYT#1171550 to D.F.G-D.).

Lissette Duarte: contributed to the design and the implementation of the research, the analysis of the results, and the writing of the manuscript; Vannesa Villanueva, Robert Barroux, Juan Francisco Orellana, Carlos and Poblete-Aro: contributed to the data curation of the research and the writing of the manuscript; Martin Gotteland: contributed to the analysis of the results and the writing of the manuscript; Mauricio Castro and Fabien Magne: contributed to the data curation, the analysis of the results, and the writing of the manuscript; Diego Garcia-Diaz: conceived the study and were in charge of direction and planning of project, directed the project and were involved in planning and supervision the work, and contributed to the analysis of the results and the writing of the manuscript.

All data generated in this study are included in this article. For more information, please contact the corresponding author.

1.
Hartroft
WS
.
The pathology of obesity
.
Bull N Y Acad Med
.
1960
;
36
(
5
):
313
22
.
2.
Hull
D
,
Segall
MM
.
The contribution of brown adipose tissue to heat production in the new-born rabbit
.
J Physiol
.
1965
;
181
(
3
):
449
57
.
3.
Cannon
B
,
Nedergaard
J
.
Brown adipose tissue: function and physiological significance
.
Physiol Rev
.
2004
;
84
(
1
):
277
359
.
4.
Concha
F
,
Prado
G
,
Quezada
J
,
Ramirez
A
,
Bravo
N
,
Flores
C
, et al
.
Nutritional and non-nutritional agents that stimulate white adipose tissue browning
.
Rev Endocr Metab Disord
.
2019
;
20
(
2
):
161
71
.
5.
Klingenspor
M
,
Herzig
S
,
Pfeifer
A
.
Brown fat develops a brite future
.
Obes Facts
.
2012
;
5
(
6
):
890
6
.
6.
Neyrinck
AM
,
Bindels
LB
,
Geurts
L
,
Van Hul
M
,
Cani
PD
,
Delzenne
NM
.
A polyphenolic extract from green tea leaves activates fat browning in high-fat-diet-induced obese mice
.
J Nutr Biochem
.
2017
;
49
:
15
21
.
7.
Manach
C
,
Scalbert
A
,
Morand
C
,
Remesy
C
,
Jimenez
L
.
Polyphenols: food sources and bioavailability
.
Am J Clin Nutr
.
2004
;
79
(
5
):
727
47
.
8.
Marin
L
,
Miguelez
EM
,
Villar
CJ
,
Lombo
F
.
Bioavailability of dietary polyphenols and gut microbiota metabolism: antimicrobial properties
.
BioMed Res Int
.
2015
;
2015
:
905215
.
9.
Le Chatelier
E
,
Nielsen
T
,
Qin
J
,
Prifti
E
,
Hildebrand
F
,
Falony
G
, et al
.
Richness of human gut microbiome correlates with metabolic markers
.
Nature
.
2013
;
500
(
7464
):
541
6
.
10.
Heiman
ML
,
Greenway
FL
.
A healthy gastrointestinal microbiome is dependent on dietary diversity
.
Mol Metab
.
2016
;
5
(
5
):
317
20
.
11.
Rechner
AR
,
Smith
MA
,
Kuhnle
G
,
Gibson
GR
,
Debnam
ES
,
Srai
SK
, et al
.
Colonic metabolism of dietary polyphenols: influence of structure on microbial fermentation products
.
Free Radic Biol Med
.
2004
;
36
(
2
):
212
25
.
12.
Duarte
LQ
,
Quezada
J
,
Ramirez
LA
,
Vasquez
K
,
Orellana
JF
,
Villanueva
V
, et al
.
The treatment with an extract from Calafate (Berberis microphylla) induces transcript and protein expression of molecules involved in thermogenesis and adipocyte browning in adipose tissue from obese mice
.
J Berry Res
.
2021
;
11
(
2
):
267
77
.
13.
Ramirez
LAQJ
,
Quezada
J
,
Duarte
L
,
Concha
F
,
Escobillana
L
,
Rincon-Cervera
MA
, et al
.
The administration of an extract from Berberis microphylla stimulates energy expenditure, thermogenesis and mitochondrial dynamics in mice brown adipose tissue
.
Food Biosci
.
2021
;
41
:
100988
.
14.
Abt
MC
,
Osborne
LC
,
Monticelli
LA
,
Doering
TA
,
Alenghat
T
,
Sonnenberg
GF
, et al
.
Commensal bacteria calibrate the activation threshold of innate antiviral immunity
.
Immunity
.
2012
;
37
(
1
):
158
70
.
15.
Soto-Covasich
JR-FM
,
Reyes-Farias
M
,
Torres
R
,
Vasquez
K
,
Duarte
L
,
Quezada
J
, et al
.
A polyphenol-rich Calafate (Berberis microphylla) extract rescues glucose tolerance in mice fed with cafeteria diet
.
J Funct Foods
.
2020
;
67
:
103856
.
16.
Livak
KJ
,
Schmittgen
TD
.
Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method
.
Methods
.
2001
;
25
(
4
):
402
8
.
17.
Castro-Sepulveda
M
,
Morio
B
,
Tunon-Suarez
M
,
Jannas-Vela
S
,
Diaz-Castro
F
,
Rieusset
J
, et al
.
The fasting-feeding metabolic transition regulates mitochondrial dynamics
.
FASEB J
.
2021
;
35
(
10
):
e21891
.
18.
Echeverria
F
,
Jimenez Patino
PA
,
Castro-Sepulveda
M
,
Bustamante
A
,
Garcia Concha
PA
,
Poblete-Aro
C
, et al
.
Microencapsulated pomegranate peel extract induces mitochondrial complex IV activity and prevents mitochondrial cristae alteration in brown adipose tissue in mice fed on a high-fat diet
.
Br J Nutr
.
2021
;
126
(
6
):
825
36
.
19.
Caporaso
JG
,
Kuczynski
J
,
Stombaugh
J
,
Bittinger
K
,
Bushman
FD
,
Costello
EK
, et al
.
QIIME allows analysis of high-throughput community sequencing data
.
Nat Methods
.
2010
;
7
(
5
):
335
6
.
20.
Lozupone
C
,
Knight
R
.
UniFrac: a new phylogenetic method for comparing microbial communities
.
Appl Environ Microbiol
.
2005
;
71
(
12
):
8228
35
.
21.
Castro-Sepulveda
M
,
Fernandez-Verdejo
R
,
Tunon-Suarez
M
,
Morales-Zuniga
J
,
Troncoso
M
,
Jannas-Vela
S
, et al
.
Low abundance of Mfn2 protein correlates with reduced mitochondria-SR juxtaposition and mitochondrial cristae density in human men skeletal muscle: examining organelle measurements from TEM images
.
FASEB J
.
2021
;
35
(
4
):
e21553
.
22.
Eramo
MJ
,
Lisnyak
V
,
Formosa
LE
,
Ryan
MT
.
The 'mitochondrial contact site and cristae organising system’ (MICOS) in health and human disease
.
J Biochem
.
2020
;
167
(
3
):
243
55
.
23.
de Jesus
LA
,
Carvalho
SD
,
Ribeiro
MO
,
Schneider
M
,
Kim
SW
,
Harney
JW
, et al
.
The type 2 iodothyronine deiodinase is essential for adaptive thermogenesis in brown adipose tissue
.
J Clin Invest
.
2001
;
108
(
9
):
1379
85
.
24.
Ohno
H
,
Shinoda
K
,
Spiegelman
BM
,
Kajimura
S
.
PPARγ agonists induce a white-to-brown fat conversion through stabilization of PRDM16 protein
.
Cell Metab
.
2012
;
15
(
3
):
395
404
.
25.
Seale
P
,
Kajimura
S
,
Yang
W
,
Chin
S
,
Rohas
LM
,
Uldry
M
, et al
.
Transcriptional control of brown fat determination by PRDM16
.
Cell Metab
.
2007
;
6
(
1
):
38
54
.
26.
Duarte
L
,
Gasaly
N
,
Poblete-Aro
C
,
Uribe
D
,
Echeverria
F
,
Gotteland
M
, et al
.
Polyphenols and their anti-obesity role mediated by the gut microbiota: a comprehensive review
.
Rev Endocr Metab Disord
.
2021
;
22
(
2
):
367
88
.
27.
Yin
J
,
Li
Y
,
Han
H
,
Chen
S
,
Gao
J
,
Liu
G
, et al
.
Melatonin reprogramming of gut microbiota improves lipid dysmetabolism in high-fat diet-fed mice
.
J Pineal Res
.
2018
;
65
(
4
):
e12524
.
28.
Singh
RP
,
Halaka
DA
,
Hayouka
Z
,
Tirosh
O
.
High-fat diet induced alteration of mice microbiota and the functional ability to utilize fructooligosaccharide for ethanol production
.
Front Cell Infect Microbiol
.
2020
;
10
:
376
.
29.
Rabot
S
,
Membrez
M
,
Blancher
F
,
Berger
B
,
Moine
D
,
Krause
L
, et al
.
High fat diet drives obesity regardless the composition of gut microbiota in mice
.
Sci Rep
.
2016
;
6
:
32484
.
30.
Li
L
,
Ma
L
,
Wen
Y
,
Xie
J
,
Yan
L
,
Ji
A
, et al
.
Crude polysaccharide extracted from moringa oleifera leaves prevents obesity in association with modulating gut microbiota in high-fat diet-fed mice
.
Front Nutr
.
2022
;
9
:
861588
.
31.
He
X
,
Zhang
M
,
Li
ST
,
Li
X
,
Huang
Q
,
Zhang
K
, et al
.
Alteration of gut microbiota in high-fat diet-induced obese mice using carnosic acid from rosemary
.
Food Sci Nutr
.
2022
;
10
(
7
):
2325
32
.
32.
Beaumont
M
,
Andriamihaja
M
,
Lan
A
,
Khodorova
N
,
Audebert
M
,
Blouin
J-M
, et al
.
Detrimental effects for colonocytes of an increased exposure to luminal hydrogen sulfide: The adaptive response
.
Free Radic Biol Med
.
2016
;
93
:
155
164
.
33.
Zheng
SJ
,
Huang
KL
,
Zhao
CH
,
Xu
WT
,
Sheng
Y
,
Luo
YB
, et al
.
Procyanidin attenuates weight gain and modifies the gut microbiota in high fat diet induced obese mice
.
J Funct Foods
.
2018
;
49
:
362
8
.
34.
Molinari
R
,
Merendino
N
,
Costantini
L
.
Polyphenols as modulators of pre-established gut microbiota dysbiosis: state-of-the-art
.
Biofactors
.
2022
;
48
(
2
):
255
73
.
35.
Leclercq
S
,
Le Roy
T
,
Furgiuele
S
,
Coste
V
,
Bindels
LB
,
Leyrolle
Q
, et al
.
Gut microbiota-induced changes in β-hydroxybutyrate metabolism are linked to altered sociability and depression in alcohol use disorder
.
Cell Rep
.
2020
;
33
(
2
):
108238
.
36.
Manca
C
,
Boubertakh
B
,
Leblanc
N
,
Deschenes
T
,
Lacroix
S
,
Martin
C
, et al
.
Germ-free mice exhibit profound gut microbiota-dependent alterations of intestinal endocannabinoidome signaling
.
J Lipid Res
.
2020
;
61
(
1
):
70
85
.
37.
Le Roy
T
,
Debedat
J
,
Marquet
F
,
Da-Cunha
C
,
Ichou
F
,
Guerre-Millo
M
, et al
.
Comparative evaluation of microbiota engraftment following fecal microbiota transfer in mice models: age, kinetic and microbial status matter
.
Front Microbiol
.
2018
;
9
:
3289
.
38.
Zhang
H
,
Sparks
JB
,
Karyala
SV
,
Settlage
R
,
Luo
XM
.
Host adaptive immunity alters gut microbiota
.
ISME J
.
2015
;
9
(
3
):
770
81
.
39.
Lac
M
,
Tavernier
G
,
Moro
C
.
Does housing temperature influence glucose regulation and muscle-fat crosstalk in mice
.
Biochimie
.
2023
;
210
:
35
9
.