The genome DNA was used as template for PCR amplification with the barcoded primers and Tks Gflex DNA Polymerase (Takara). For bacterial diversity analysis, V3-V4 variable regions of 16S rRNA genes was amplified with universal primers 343F and 798R. The library sequencing and data processing were conducted by OE biotech Co., Ltd. (Shanghai, China).