design: DNA for each sample was extracted with FastDNA SPIN Kit (MP Biomedicals), then V3-V4(338F-806R) and of 16S rDNA was amplified by two-step PCR.    Finally, the amplicons were sequenced on Illumina MiSeq PE300 platform in pair-end 300 bp mode