Abstract
The emergence and development of next-generation sequencing technologies (NGS) has made the analysis of the water microbiome in drinking water distribution systems (DWDSs) more accessible and opened new perspectives in microbial ecology studies. The current study focused on the characterization of the water microbiome employing a gene- and genome-centric metagenomic approach to five waterworks in Finland with different raw water sources, treatment methods, and disinfectant. The microbial communities exhibit a distribution pattern of a few dominant taxa and a large representation of low-abundance bacterial species. Changes in the community structure may correspond to the presence or absence and type of disinfectant residual which indicates that these conditions exert selective pressure on the microbial community. The Archaea domain represented a small fraction (up to 2.5%) and seemed to be effectively controlled by the disinfection of water. Their role particularly in non-disinfected DWDS may be more important than previously considered. In general, non-disinfected DWDSs harbor higher microbial richness and maintaining disinfectant residual is significantly important for ensuring low microbial numbers and diversity. Metagenomic binning recovered 139 (138 bacterial and 1 archaeal) metagenome-assembled genomes (MAGs) that had a >50% completeness and <10% contamination consisting of 20 class representatives in 12 phyla. The presence and occurrence of nitrite-oxidizing bacteria (NOB)-like microorganisms have significant implications for nitrogen biotransformation in drinking water systems. The metabolic and functional complexity of the microbiome is evident in DWDSs ecosystems. A comparative analysis found a set of differentially abundant taxonomic groups and functional traits in the active community. The broader set of transcribed genes may indicate an active and diverse community regardless of the treatment methods applied to water. The results indicate a highly dynamic and diverse microbial community and confirm that every DWDS is unique, and the community reflects the selection pressures exerted at the community structure, but also at the levels of functional properties and metabolic potential.