design: The V4 region of the 16S rRNA gene was amplified using primers 515F and 806R. Resulting PCR products were subjected to pair-end 250 bp sequencing using Illumina NovaSeq 6000 platform (Illumina, San Diego, CA, USA) at Novogene (Beijing, PR China) at read depth of 30,000 tags per sample.