Taxonomic relatedness and environmental pressure synergistically drive the primary succession of biofilm microbial communities in reclaimed wastewater distribution systems

Compared to drinking water, the higher bacterial abundance, diversity, and organic matter concentration in reclaimed wastewater suggest that it is more likely to form biofilms. Nevertheless, little is known regarding many important aspects of the biofilm ecology in reclaimed wastewater distribution systems (RWDS), such as the long-term microbial community succession and the underlying driving factors. In the present study, by sampling and analysing microbial compositions of pipe wall biofilms from six frequently used pipe materials under NaClOdisinfection (sodium hypochlorite-treated), NONdisinfection (without disinfection), and UVdisinfection (UV-treated) treatments over one year, it was found that the succession of microbial community structure followed a primary succession pattern. This primary succession pattern was reflected as increases in live cell number and α-diversity, along with metagenic succession in taxonomic composition. Proteobacteria, Nitrospirae, Bacteroidetes, Acidobacteria, Planctomycetes, Actinobacteria, and Verrucomicrobia comprised the dominant phyla in biofilm samples. Compared to biofilms in the NaClOdisinfection reactor, the bacterial communities of biofilms in NONdisinfection and UVdisinfection reactors were distributed more evenly among different bacterial phyla. Principal component analysis revealed a clear temporal pattern of microbial community structures in six kinds of pipe wall biofilms albeit a difference in microbial community structures among the three reactors. Adonis testing indicated that the microbial community composition variation caused by disinfection methods (R2 = 0.283, P < 0.01) was more pronounced than that from the time variable (R2 = 0.070, P < 0.01) and pipe material (R2 = 0.057, P < 0.01). Significantly positive correlation between average local abundance and occupancy was observed in biofilm communities of the three reactors, suggesting that the ‘core-satellite’ model could be applied to identify biofilm-preferential species under specific disinfection conditions in RWDS. The prevalence of family Sphingomonadaceae, known to show chlorine tolerance and powerful biofilm-forming ability in NaClOdisinfection reactors, evidenced the habitat filtering consequent to environment pressure. Correlation-based network analysis revealed that taxonomic relatedness such as similar niches, cooperation, taxa overdispersion, and competition all functioned toward driving the bacterial assembly succession in RWDS. Keywords: Reclaimed wastewater, Pipe wall biofilm, Primary succession, Core-satellite model, Network analysis