S1_Fig.pdf

The horizontal axle represents the 13,004 qualified OTUs. The singletons in sequenced data were removed and the remained singletons were caused by the removal of samples with <200 reads. The vertical axle is the number of sequences in each OTU. To fit graphing, the reads numbers were transformed to logarithm.

10.1371/journal.pone.0141842.s001 S2_Fig.pdf

The branch color means different sites, red branch: Na, yellow branch: Gb and Ba, blue branch: Is and Af, green branch: Hb and Vf. The sample names from different sites with same gender, age and place of residences were showed with same color.

10.1371/journal.pone.0141842.s002 S3_Fig.pdf 10.1371/journal.pone.0141842.s003 S1_Table.xlsx 10.1371/journal.pone.0141842.s004 S2_Table.xlsx

All P values were corrected by the Bonferroni procedure for multiple comparisons.

10.1371/journal.pone.0141842.s005 S3_Table.xlsx

The table matched to Fig 4.

10.1371/journal.pone.0141842.s006 S4_Table.xlsx 10.1371/journal.pone.0141842.s007 S5_Table.xlsx

All P values were corrected by the Bonferroni procedure for multiple comparisons.

10.1371/journal.pone.0141842.s008 S6_Table.xlsx 10.1371/journal.pone.0141842.s009