SRX3908239 GD-C1-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144108 SAMN08875358 GD-C1-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908240 HN-W1-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144109 SAMN08875356 HN-W1-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908241 HN-C1-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144107 SAMN08875359 HN-C1-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908242 JX-C1-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144110 SAMN08875357 JX-C1-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908243 JX-W2-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144111 SAMN08875360 JX-W2-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908244 JX-W1-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144112 SAMN08875363 JX-W1-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908245 GD-W2-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144114 SAMN08875370 GD-W2-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908246 HN-W2-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144113 SAMN08875371 HN-W2-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908247 HN-C1-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144115 SAMN08875368 HN-C1-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908248 JX-W2-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144116 SAMN08875369 JX-W2-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908249 JX-C1-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144117 SAMN08875366 JX-C1-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908250 GD-C1-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144118 SAMN08875367 GD-C1-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908251 GD-W1-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144119 SAMN08875364 GD-W1-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908252 HN-W1-F SRP139050 PRJNA449026 The primer pair used to amplify the fungal ITS region was ITS3F (5'-GATGAAGAACGYAGYRAA-3') and ITS4R (5'-TCCTCCGCTTATTGATATGC-3'). The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144120 SAMN08875365 HN-W1-F AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908253 JX-W1-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144121 SAMN08875354 JX-W1-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908254 HN-W2-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144122 SAMN08875362 HN-W2-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908255 GD-W2-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144124 SAMN08875361 GD-W2-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500 SRX3908256 GD-W1-B SRP139050 PRJNA449026 The primer pair used to amplify the V3-V4 hypervariable regions of the bacterial 16SrRNA genes was 341F (5 -CCTACGGGNGGCWGCAG-3 ) and 785R (5 -GACTACHVGGGTATCTAATCC-3 ) . The PCR program was as follows: 2 min at 95 C; 25 cycles of 30 s at 95 C, 30 s at 55 C, and 45 s at 72 C; and a final extension step at 72 C for 10 min. SRS3144123 SAMN08875355 GD-W1-B AMPLICON METAGENOMIC PCR Illumina HiSeq 2500