SRX3886242 T4_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125094 SAMN08848775 T4_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886243 T4_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125093 SAMN08848776 T4_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886244 T4_7 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125095 SAMN08848773 T4_7 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886245 T4_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125096 SAMN08848774 T4_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886246 T4_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125097 SAMN08848779 T4_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886247 H_7 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125098 SAMN08848689 H_7 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886248 T4_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125099 SAMN08848777 T4_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886249 T4_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125100 SAMN08848778 T4_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886250 T4_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125101 SAMN08848781 T4_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886251 H_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125102 SAMN08848690 H_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886252 T6_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125104 SAMN08848815 T6_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886253 T6_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125103 SAMN08848816 T6_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886254 T6_6 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125105 SAMN08848813 T6_6 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886255 T6_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125106 SAMN08848814 T6_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886256 T6_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125107 SAMN08848819 T6_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886257 T6_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125108 SAMN08848820 T6_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886258 T6_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125109 SAMN08848817 T6_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886259 T6_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125130 SAMN08848818 T6_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886260 T6_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125110 SAMN08848821 T6_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886261 T6_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125111 SAMN08848822 T6_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886262 T4_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125112 SAMN08848780 T4_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886263 H_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125113 SAMN08848691 H_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886264 T4_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125114 SAMN08848782 T4_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886265 H_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125116 SAMN08848692 H_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886266 H_1 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125115 SAMN08848683 H_1 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886267 H_2 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125117 SAMN08848684 H_2 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886268 T5_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125118 SAMN08848806 T5_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886269 T5_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125119 SAMN08848805 T5_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886270 H_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125120 SAMN08848687 H_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886271 T5_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125121 SAMN08848807 T5_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886272 T6_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125122 SAMN08848810 T6_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886273 T6_2 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125124 SAMN08848809 T6_2 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886274 T6_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125123 SAMN08848812 T6_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886275 T6_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125125 SAMN08848811 T6_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886276 H_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125126 SAMN08848701 H_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886277 H_20 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125127 SAMN08848702 H_20 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886278 H_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125128 SAMN08848695 H_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886279 H_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125129 SAMN08848696 H_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886280 H_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125138 SAMN08848693 H_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886281 H_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125140 SAMN08848694 H_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886282 H_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125131 SAMN08848699 H_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886283 H_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125132 SAMN08848700 H_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886284 H_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125133 SAMN08848697 H_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886285 H_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125134 SAMN08848698 H_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886286 T5_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125135 SAMN08848801 T5_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886287 T5_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125136 SAMN08848802 T5_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886288 T5_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125137 SAMN08848797 T5_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886289 T5_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125139 SAMN08848798 T5_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886290 T5_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125142 SAMN08848799 T5_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886291 T5_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125141 SAMN08848800 T5_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886292 T5_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125166 SAMN08848793 T5_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886293 T5_6 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125143 SAMN08848794 T5_6 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886294 T5_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125144 SAMN08848795 T5_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886295 T5_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125145 SAMN08848796 T5_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886296 T1_7 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125146 SAMN08848710 T1_7 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886297 T1_6 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125168 SAMN08848709 T1_6 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886298 T1_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125147 SAMN08848708 T1_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886299 T1_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125148 SAMN08848707 T1_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886300 T1_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125149 SAMN08848706 T1_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886301 T1_2 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125173 SAMN08848705 T1_2 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886302 T1_1 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125150 SAMN08848704 T1_1 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886303 H_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125176 SAMN08848703 H_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886304 T1_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125151 SAMN08848712 T1_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886305 T1_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125152 SAMN08848711 T1_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886306 T5_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125154 SAMN08848792 T5_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886307 T5_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125153 SAMN08848791 T5_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886308 T4_20 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125155 SAMN08848786 T4_20 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886309 T4_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125180 SAMN08848785 T4_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886310 T4_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125156 SAMN08848784 T4_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886311 T4_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125157 SAMN08848783 T4_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886312 T5_2 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125158 SAMN08848790 T5_2 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886313 T5_1 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125159 SAMN08848789 T5_1 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886314 T4_22 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125160 SAMN08848788 T4_22 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886315 T4_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125161 SAMN08848787 T4_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886316 T1_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125162 SAMN08848713 T1_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886317 T1_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125163 SAMN08848714 T1_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886318 T1_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125164 SAMN08848715 T1_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886319 T1_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125165 SAMN08848716 T1_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886320 T1_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125235 SAMN08848717 T1_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886321 T1_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125167 SAMN08848718 T1_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886322 T1_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125169 SAMN08848719 T1_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886323 T1_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125170 SAMN08848720 T1_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886324 T1_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125171 SAMN08848721 T1_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886325 T1_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125197 SAMN08848722 T1_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886326 T2_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125172 SAMN08848728 T2_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886327 T2_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125174 SAMN08848727 T2_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886328 T2_7 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125204 SAMN08848730 T2_7 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886329 T2_6 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125175 SAMN08848729 T2_6 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886330 T1_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125177 SAMN08848724 T1_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886331 T1_20 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125178 SAMN08848723 T1_20 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886332 T2_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125179 SAMN08848726 T2_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886333 T1_22 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125233 SAMN08848725 T1_22 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886334 T2_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125181 SAMN08848732 T2_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886335 T2_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125216 SAMN08848731 T2_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886336 T2_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125182 SAMN08848739 T2_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886337 T2_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125222 SAMN08848740 T2_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886338 T2_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125183 SAMN08848737 T2_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886339 T2_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125184 SAMN08848738 T2_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886340 T2_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125185 SAMN08848735 T2_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886341 T2_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125186 SAMN08848736 T2_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886342 T2_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125187 SAMN08848733 T2_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886343 T2_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125189 SAMN08848734 T2_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886344 T2_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125188 SAMN08848741 T2_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886345 T2_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125190 SAMN08848742 T2_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886346 T3_10 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125236 SAMN08848754 T3_10 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886347 T3_1 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125191 SAMN08848746 T3_1 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886348 T2_22 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125192 SAMN08848745 T2_22 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886349 T2_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125193 SAMN08848744 T2_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886350 T2_20 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125194 SAMN08848743 T2_20 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886351 T3_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125237 SAMN08848750 T3_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886352 T3_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125195 SAMN08848749 T3_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886353 T3_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125196 SAMN08848748 T3_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886354 T3_2 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125198 SAMN08848747 T3_2 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886355 T3_8 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125199 SAMN08848752 T3_8 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886356 T3_7 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125200 SAMN08848751 T3_7 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886357 T3_20 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125201 SAMN08848764 T3_20 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886358 T3_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125202 SAMN08848763 T3_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886359 T3_22 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125203 SAMN08848766 T3_22 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886360 T3_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125205 SAMN08848765 T3_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886361 T4_2 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125206 SAMN08848768 T4_2 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886362 T4_1 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125207 SAMN08848767 T4_1 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886363 T4_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125208 SAMN08848770 T4_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886364 T4_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125210 SAMN08848769 T4_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886365 T4_6 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125209 SAMN08848772 T4_6 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886366 T4_5 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125211 SAMN08848771 T4_5 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886367 T6_21 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125212 SAMN08848826 T6_21 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886368 T6_19 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125213 SAMN08848825 T6_19 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886369 T6_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125214 SAMN08848824 T6_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886370 T6_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125215 SAMN08848823 T6_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886371 T5_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125217 SAMN08848804 T5_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886372 T3_17 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125218 SAMN08848761 T3_17 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886373 T3_18 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125219 SAMN08848762 T3_18 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886374 T5_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125220 SAMN08848803 T5_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886375 T3_13 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125221 SAMN08848757 T3_13 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886376 T3_14 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125223 SAMN08848758 T3_14 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886377 T3_15 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125224 SAMN08848759 T3_15 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886378 T3_16 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125225 SAMN08848760 T3_16 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886379 T3_9 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125226 SAMN08848753 T3_9 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886380 H_3 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125227 SAMN08848685 H_3 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886381 T3_11 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125228 SAMN08848755 T3_11 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886382 T3_12 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125229 SAMN08848756 T3_12 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886383 H_4 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125230 SAMN08848686 H_4 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886384 T6_1 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125232 SAMN08848808 T6_1 AMPLICON GENOMIC PCR Illumina HiSeq 2500 SRX3886385 H_6 SRP137370 PRJNA445640 Extraction of bacterial DNA was performed from fecal samples by QIAamp Fast DNA Stool Mini Kit as previous reported (Wu et al., 2016). The freshly extracted DNA was purified by 0.5% melting point agarose gel followed by phenolchloroform-butanol extraction. The V3-V4 region of the bacterial 16S rRNA gene from each DNA sample was amplified using the bacterial universal primers (forward primer: 5 -ACTCCTACGGGRSGCAGCAG-3 ; reverse primer: 5 -GGACTACVVGGGTATCTAATC-3 ).The PCR amplification products were purified through AxyPrep DNA Gel Extraction Kit (Axygen, USA), eluted in 30 l water, and aliquoted into three PCR tubes. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm, and final DNA contents were quantified with a Qubit? dsDNA HS Assay Kit (Invitrogen, USE). Finally, we used of bacterial DNA amplicons from each fecal sample for 2 * 250 bp paired-end sequencing based on the Illumina Hiseq 2500. SRS3125231 SAMN08848688 H_6 AMPLICON GENOMIC PCR Illumina HiSeq 2500