title: S2
design: Total genomic DNA for each isolated strain was extracted according to QIAamp DNA Mini Kit Cat# 51304, the quality and concentration of DNA were assessed in NanoDrop 2000C, Spectrophotometer, Thermoscientific, Germany. The variable regions 34 of the 16S rRNA gene were amplified using specific forward (5 TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG 3) and reverse primers (5 GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC 3) containing the Illumina adapter overhang nucleotide sequences. The 16S V3-V4 amplicons were quantified using TapeSation (Agilent, Santa Clara, CA, USA) using Ampure XP bits to purify them. The amplicons were roughly 550 bp in size. The Nextera XT V2 Kit was then used to perform an index PCR in order to affix dual indices. The fluorometer Qubit 3.0 (Invitrogen, Waltham, MA, USA) was used to quantify the concentration of double-stranded DNA, and the amplicon size was roughly 630 bp. Equimolar concentrations of the final amplicon library were combined. The Illumina NovaSeq platform (Illumina, San Diego, CA, USA) was used for the sequencing to produce paired-end reads that were 250 bases long in each direction.