Estimation of microbial beta diversity.
We calculated weighted UniFrac distances to obtain the pairwise beta diversity, which we further evaluated by the adonis2 function using the R package vegan to test the significance with 999 permutations (