design: Total DNA was extracted using the HiPure Soil DNA Kits (Magen, Guangzhou, China), then the bacterial V3-V4 hypervariable region of 16S rRNA was amplified by PCR. Finally, the amplicons were pooled in equimolar and paried-end sequenced (PE250) on an Illumina platform.