Statistical analysis
Microbiota alpha diversity Shannon and Chao1 were calculated using the R program package ‘vegan’ (version 2.5.6). β-diversity metrics were obtained with rda and PERMANOVA with the adonis function. Principal Components analysis (PCA) was performed using the package vegan. Different analysis was performed to identify taxa with differentiating abundance in the different groups (