design: 350 mg soil were extracted DNA by the procedure of MP BIO FASTDNA SPIN Kit (MP Bio, Santa Ana, CA). The V5V7 region of 16S rRNA gene was amplified with primers [799F (5-AACMGGATTAGATACCCKG-3) and 1193R (5-ACGTCATCCCCACCTTCC-3)] by PCR.  The PCR products were sent to Genesky Biotechnologies Inc., Shanghai, 201, 315 (China) to high throughput sequence by Illumina NovaSeq 6000 platform. We removed the primer sequences and adaptor with the cutadapt44 plug-in of QIIME245. We used DADA2 plug-in of QIIME2 to filter the data, denoise, merge and remove chimera to obtain high-quality sequencing46. The clean reads were processed into amplified sequence variant by the QIIME2.