SRX18319702 TaC21 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE250 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810159 TaC21 TaC21 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319703 TaC22 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE251 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810160 TaC22 TaC22 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319704 EaW2 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE260 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810161 EaW2 EaW2 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319705 EaW3 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE261 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810162 EaW3 EaW3 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319706 TaC23 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE252 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810163 TaC23 TaC23 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319707 TaW21 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE253 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810164 TaW21 TaW21 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319708 TaW22 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE254 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810165 TaW22 TaW22 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319709 TaW213 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE255 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810166 TaW213 TaW213 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319710 EaC1 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE256 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810167 EaC1 EaC1 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319711 EaC2 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE257 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810168 EaC2 EaC2 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319712 EaC3 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE258 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810169 EaC3 EaC3 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000 SRX18319713 EaW1 SRP408743 bp0 DNA extraction, amplification, and sequencing were conducted by PersonalBio Biotechnology Co., Ltd. (Shanghai, China). In brief, complete DNA samples were extracted using the E.Z.N.ATM Mag-Bind Soil DNA Kit (M5635, OMEGA, USA), following the manufacturers protocol. The quantity and quality of the extracted DNA were assessed using a fluorescence spectrophotometer (Quantifluor-ST fluorometer, Promega, E6090; Quant-iT PicoGreen dsDNA Assay Kit, Invitrogen, P7589) at 260 nm and 280 nm, respectively, and was also detected via 1.2% agarose gel electrophoresis. PCR amplification of the bacterial 16S rRNA gene V3-V4 hypervariable region was performed using a forward primer (338F:5-ACTCCTACGGGAGGCAGCA-3) and reverse primer (806R:5-GGACTACHVGGGTWTCTAAT-3). The PCR reaction (25 L) was prepared as follows: template DNA 1 L, amplicon PCR forward primer (10 M) 1 L, amplicon PCR reverse primer (10 M) 1 L, dNTP (2.5 mM) 2 L, Fast Pfu DNA Polymerase (0.25 L), 2 buffer (5 L), and ddH2O (14.75 L). PCR was performed using the following program: denaturation at 98 C for 5 min followed by 25 cycles consisting of denaturation at 98 C for 30 s, annealing at 52 C for 30 s, and extension at 72 C for 45 s, with a final extension of 5 min at 72 C. PE259 paired-end sequencing was performed according to the Illumina MiSeq (Illumina San Diego, CA, USA) instrument manual after the DNA libraries were mixed. SRS15810170 EaW1 EaW1 AMPLICON METAGENOMIC PCR Illumina NovaSeq 6000