The CAP ordinations demonstrated that, in terms of both com-munity structure and species composition, the microbial eukaryotecommunities showed a similar pattern along the gradient of stress,i.e., The station A and D were clearly separated by the first canonical axis (CAP 1), while these two stations may easily distinguished fromstations B and C by the second canonical axis (CAP 2) (Fig. 3a and c).The PERMDISP signified a difference among the four stations andbetween each pair of stations (P < 0.05), apart from pairs of stationsB/C (P > 0.05). In terms of the vectors for top species with correla-tion coefficients >0.2 pointed toward the samples at each station,those for the species abundance resolution mainly were dominantforms with both high abundance and occurrence (Fig. 3b; Table S2),while those for the species composition resolution were the taxawith only high occurrence frequency (Fig. 3d).Based on the modified Gower dissimilarity matrix, both theCAP ordination and the PERMDIDP analyses did not exhibit a clearthe spatial variation in response to the environmental conditionsamong the four stations (Fig. 4e). The vectors for top species withcorrelation coefficients >0.2 mainly were dominant forms with onlyhigh abundance (Fig. 3f; Table S2).