title: data
design: Paired-end reads were assigned to samples based on unique barcodes and merged using FLASH. Then FqTrim (v. 0.94) was used to clean raw reads while Vsearch (v2.3.4) was utilized to filter chimeric sequences. Next, DADA2 was used for dereplication to obtain a feature table and a sequence followed by the calculation of alpha and beta diversities utilizing QIIME2. Subsequently, R version 3.5.2 was used to construct the corresponding phylogenetic tree, and the Basic Local Alignment Search Tool (BLAST) was utilized for sequence alignment and species annotation, whereas sequence alignment was performed using SILVA and NT-16S.