Spatial and Temporal Distribution of Bacterioplankton Molecular Ecological Networks in the Yuan River under Different Human Activity Intensity.

Bacterioplankton communities play a crucial role in freshwater ecosystem functioning, but it is unknown how co-occurrence networks within these communities respond to human activity disturbances. This represents an important knowledge gap because changes in microbial networks could have implications for their functionality and vulnerability to future disturbances. Here, we compare the spatiotemporal and biogeographical patterns of bacterioplankton molecular ecological networks using high-throughput sequencing of Illumina HiSeq and multivariate statistical analyses from a subtropical river during wet and dry seasons. Results demonstrated that the lower reaches (high human activity intensity) network had less of an average degree (10.568/18.363), especially during the dry season, when compared with the upper reaches (low human activity intensity) network (10.685/37.552) during the wet and dry seasons, respectively. The latter formed more complexity networks with more modularity (0.622/0.556) than the lower reaches (high human activity intensity) network (0.505/0.41) during the wet and dry seasons, respectively. Bacterioplankton molecular ecological network under high human activity intensity became significantly less robust, which is mainly caused by altering of the environmental conditions and keystone species. Human activity altered the composition of modules but preserved their ecological roles in the network and environmental factors (dissolved organic carbon, temperature, arsenic, oxidation-reduction potential and Chao1 index) were the best parameters for explaining the variations in bacterioplankton molecular ecological network structure and modules.