Statistical analysis
Data are presented as means ± standard deviations unless otherwise stated. All variables were tested for normality using D’Agostino and Pearson omnibus. Most of the variables displayed a non-normal distribution and, thus, non-parametric tests were used for all analyses. We did not detect any sex interaction across the variables studied (Fig. S1 and Table S1), nor differences between the status of BMI (data not shown); therefore, all main data for men and women, as well as the status of BMI, were pooled together. Partial Spearman correlations were used to investigate the correlation of fecal microbiota composition with BAT volume, 