Eukaryotes are ubiquitous in natural environments such as soil and freshwater. Little is known of their presence in drinking water distribution systems (DWDSs) or of the environmental conditions that affect their activity and survival. Eukaryotes were characterized by Illumina high-throughput sequencing targeting 18S rRNA gene (DNA) that estimates the total community and the 18S rRNA gene transcript (RNA) that is more representative of the active part of the community. DWDS cold water (N = 124), hot water (N = 40), and biofilm (N = 16) samples were collected from four cities in Finland. The sampled DWDSs were from two waterworks A–B with non-disinfected, recharged groundwater as source water and from three waterworks utilizing chlorinated water (two DWDSs of surface waterworks C–D and one of ground waterworks E). In each DWDS, samples were collected from three locations during four seasons of 1 year. A beta-diversity analysis revealed that the main driver shaping the eukaryotic communities was the DWDS (A–E) (R = 0.73, P < 0.001, ANOSIM). The kingdoms Chloroplastida (green plants and algae), Metazoa (animals: rotifers, nematodes), Fungi (e.g., Cryptomycota), Alveolata (ciliates, dinoflagellates), and Stramenopiles (algae Ochrophyta) were well represented and active—judging based on the rRNA gene transcripts—depending on the surrounding conditions. The unchlorinated cold water of systems (A–B) contained a higher estimated total number of taxa (Chao1, average 380–480) than chlorinated cold water in systems C–E (Chao1 ≤ 210). Within each DWDS, unique eukaryotic communities were identified at different locations as was the case also for cold water, hot water, and biofilms. A season did not have a consistent impact on the eukaryotic community among DWDSs. This study comprehensively characterized the eukaryotic community members within the DWDS of well-maintained ground and surface waterworks providing good quality water. The study gives an indication that each DWDS houses a unique eukaryotic community, mainly dependent on the raw water source and water treatment processes in place at the corresponding waterworks. In particular, disinfection as well as hot water temperature seemed to represent a strong selection pressure that controlled the number of active eukaryotic species.

Bacteriophage control of harmful or pathogenic bacteria has aroused growing interest, largely due to the rise of antibiotic resistance. The objective of this study was to test phages as potential agents for the biocontrol of an opportunistic pathogen Pseudomonas aeruginosa in water. Two P. aeruginosa bacteriophages (vB_PaeM_V523 and vB_PaeM_V524) were isolated from wastewater and characterized physically and functionally. Genomic and morphological characterization showed that both were myoviruses within the Pbunavirus genus. Both had a similar latent period (50–55 min) and burst size (124–134 PFU/infected cell), whereas there was variation in the host range. In addition to these environmental phages, a commercial Pseudomonas phage, JG003 (DSM 19870), was also used in the biocontrol experiments. The biocontrol potential of the three phages in water was tested separately and together as a cocktail against two P. aeruginosa strains; PAO1 and the environmental strain 17V1507. With PAO1, all phages initially reduced the numbers of the bacterial host, with phage V523 being the most efficient (>2.4 log₁₀ reduction). For the environmental P. aeruginosa strain (17V1507), only the phage JG003 caused a reduction (1.2 log₁₀) compared to the control. The cocktail of three phages showed a slightly higher decrease in the level of the hosts compared to the use of individual phages. Although no synergistic effect was observed in the host reduction with the use of the phage cocktail, the cocktail-treated hosts did not appear to acquire resistance as rapidly as hosts treated with a single phage. The results of this study provide a significant step in the development of bacteriophage preparations for the control of pathogens and harmful microbes in water environments.

Failures in the drinking water distribution system cause gastrointestinal outbreaks with multiple pathogens. A water distribution pipe breakage caused a community-wide waterborne outbreak in Vuorela, Finland, July 2012. We investigated this outbreak with advanced epidemiological and microbiological methods. A total of 473/2931 inhabitants (16%) responded to a web-based questionnaire. Water and patient samples were subjected to analysis of multiple microbial targets, molecular typing and microbial community analysis. Spatial analysis on the water distribution network was done and we applied a spatial logistic regression model. The course of the illness was mild. Drinking untreated tap water from the defined outbreak area was significantly associated with illness (RR 5.6, 95% CI 1.9–16.4) increasing in a dose response manner. The closer a person lived to the water distribution breakage point, the higher the risk of becoming ill. Sapovirus, enterovirus, single Campylobacter jejuni and EHEC O157:H7 findings as well as virulence genes for EPEC, EAEC and EHEC pathogroups were detected by molecular or culture methods from the faecal samples of the patients. EPEC, EAEC and EHEC virulence genes and faecal indicator bacteria were also detected in water samples. Microbial community sequencing of contaminated tap water revealed abundance of Arcobacter species. The polyphasic approach improved the understanding of the source of the infections, and aided to define the extent and magnitude of this outbreak.

Many factors, including microbiome structure and activity in the drinking water distribution system (DWDS), affect the colonization potential of opportunistic pathogens. The present study aims to describe the dynamics of active bacterial communities in DWDS and identify the factors that shape the community structures and activity in the selected DWDSs. Large-volume drinking water and hot water, biofilm, and water meter deposit samples were collected from five DWDSs. Total nucleic acids were extracted, and RNA was further purified and transcribed into its cDNA from a total of 181 water and biofilm samples originating from the DWDS of two surface water supplies (disinfected with UV and chlorine), two artificially recharged groundwater supplies (non-disinfected), and a groundwater supply (disinfected with UV and chlorine). In chlorinated DWDSs, concentrations of <0.02-0.97 mg/l free chlorine were measured. Bacterial communities in the RNA and DNA fractions were analysed using Illumina MiSeq sequencing with primer pair 341F-785R targeted to the 16S rRNA gene. The sequence libraries were analysed using QIIME pipeline, Program R, and MicrobiomeAnalyst. Not all bacterial cells were active based on their 16S rRNA content, and species richness was lower in the RNA fraction (Chao1 mean value 490) than in the DNA fraction (710). Species richness was higher in the two DWDSs distributing non-disinfected artificial groundwater (Chao1 mean values of 990 and 1 000) as compared to the two disinfected DWDSs using surface water (Chao1 mean values 190 and 460) and disinfected DWDS using ground water as source water (170). The difference in community structures between non-disinfected and disinfected water was clear in the beta-diversity analysis. Distance from the waterworks also affected the beta diversity of community structures, especially in disinfected distribution systems. The two most abundant bacteria in the active part of the community (RNA) and total bacterial community (DNA) belonged to the classes Alphaproteobacteria (RNA 28 %, DNA 44 %) and Gammaproteobacteria (RNA 32 %, DNA 30 %). The third most abundant and active bacteria class was Vampirovibrionia (RNA 15 %), whereas in the total community it was Paceibacteria (DNA 11 %). Class Nitrospiria was more abundant and active in both cold and hot water in DWDS that used chloramine disinfection compared to non-chlorinated or chlorine-using DWDSs. Thirty-eight operational taxonomic units (OTU) of Legionella, 30 of Mycobacterium, and 10 of Pseudomonas were detected among the sequences. The (RT)-qPCR confirmed the presence of opportunistic pathogens in the DWDSs studied as Legionella spp. was detected in 85 % (mean value 4.5 × 10⁴ gene copies/100 ml), Mycobacterium spp. in 95 % (mean value 8.3 × 10⁶ gene copies/100 ml), and Pseudomonas spp. in 78 % (mean value 1.6 × 10⁵ gene copies/100 ml) of the water and biofilm samples. Sampling point inside the system (distance from the waterworks and cold/hot system) affected the active bacterial community composition. Chloramine as a chlorination method resulted in a recognizable community composition, with high abundance of bacteria that benefit from the excess presence of nitrogen. The results presented here confirm that each DWDS is unique and that opportunistic pathogens are present even in conditions when water quality is considered excellent.

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.