design: Total bacterial diversity in each sample was estimated by sequencing the V3 region of the 16S rRNA gene. The sequences were analyzed following the protocol described inhttps://github.com/GenomicaMicrob