Sample size and statistical analysis
Sample size calculation was estimated by performing a statistical power analysis [The normality of distribution was assessed by the Shapiro–Wilk test or Kolmogorov–Smirnov test. Unpaired Student T-test or Mann–Whitney U-test were used for normally or not normally distributed data, respectively. Differences between the groups at each time point were examined post hoc using an independent sample t-test or Chi-square test. Bray–Curtis, unweighted, weighted and UniFrac distances matrix, Bokulich longitudinal, Robust Aitchison Principal Component Analysis with DEICODE and Principal Coordinates Analyses (PCoA) were performed using QIIME 2 software [