whole DNA was extracted with QIAamp Fast DNA Stool Mini Kit; 16S rRNA gene V3V4 region was PCR-amplified to generate amplicon metagenomic libraries; sequencing was performed on an Illumina NovaSeq 6000 platform with paired-end 250 bp reads